这其实是道水题。。。

题目链接: (bzoj)https://www.lydsy.com/JudgeOnline/problem.php?id=4032

(luogu)https://www.luogu.org/problemnew/show/P4112

题解:

Task 1

\(O(n^2)\)做法无数(也有不用SAM的\(O(n^2)\)做法),我讲一下我的做法。

直接对B串建后缀自动机,用A串在上面跑,而且是像求两个串的最长公共子串那样跑,如果遇到失配了就拿当前长度\(+1\)更新答案。但是要注意失配后当前的长度不能设成当前节点的\(len\)(最大长度),而应该是\(len[fail[u]]+1\)(最小长度)。

时间复杂度\(O(n)\).

由于没有在网上看到\(O(n)\)做法,所以此做法正确性未知(能过,但是数据很水,一开始写成\(len\)了依然能过10个点中的9个)。

Task 2

设\(nxt[i][j]\)表示B序列中第\(i\)个位置之后最靠前的字符\(j\)的位置。枚举A子串的起始位置,贪心即可。

时间复杂度\(O(n^2)\).

听说有人把这个东西称作“序列自动机”,感觉也有道理啊,这玩意确实像个自动机。

Task 3

子串就用后缀自动机,子序列就用“序列自动机”美滋滋。

设\(dp[i][j]\)表示A串前\(i\)个位置匹配B串后缀自动机的节点\(j\),最少用多少长度。如果能转移就转移,不能转移就用\(dp[i][j]+1\)更新答案。注意要在每个\(i\)都更新(因为子序列可以不到头)。

然后第一维可以像01背包一样省掉,我觉得第二维要按先儿子后父亲的顺序循环。但是网上有人直接从\(1\)到\(siz\)循环也能过?

注意\(dp\)数组要开两倍。

Task 4

设\(dp[i][j]\)表示A串前\(i\)个位置匹配B串前\(j\)个位置的最小长度。

转移同Task 3. 如果用一维数组,从\(n\)到\(1\)循环。

代码

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std; const int N = 2000;
const int S = 26;
int len[(N<<1)+3];
int fa[(N<<1)+3];
int son[(N<<1)+3][S+3];
char a[N+3],b[N+3];
int pos[S+3];
int nxt[N+3][S+3];
int dp[(N<<1)+3];
int ord[(N<<1)+3];
int buc[N+3];
int n,m,siz,rtn,lstpos; void initSAM()
{
siz = rtn = lstpos = 1;
} void insertchar(char ch)
{
int p = lstpos,np; siz++; np = lstpos = siz; len[np] = len[p]+1;
for(; p && son[p][ch]==0; p=fa[p]) {son[p][ch] = np;}
if(!p) {fa[np] = rtn;}
else
{
int q = son[p][ch];
if(len[p]+1==len[q]) {fa[np] = q;}
else
{
siz++; int nq = siz; len[nq] = len[p]+1;
memcpy(son[nq],son[q],sizeof(son[q]));
fa[nq] = fa[q]; fa[np] = fa[q] = nq;
for(; p && son[p][ch]==q; p=fa[p]) {son[p][ch] = nq;}
}
}
} void update(int &x,int y) {x = min(x,y);} int main()
{
initSAM();
scanf("%s",a+1); n = strlen(a+1); for(int i=1; i<=n; i++) a[i] -= 96;
scanf("%s",b+1); m = strlen(b+1); for(int i=1; i<=m; i++) b[i] -= 96;
for(int i=1; i<=m; i++)
{
insertchar(b[i]);
}
for(int i=1; i<=m; i++)
{
for(int j=pos[b[i]]; j<i; j++)
{
nxt[j][b[i]] = i;
}
pos[b[i]] = i;
}
for(int i=1; i<=S; i++)
{
for(int j=0; j<=m; j++) {if(!nxt[j][i]) nxt[j][i] = m+1;}
}
for(int i=1; i<=siz; i++) buc[len[i]]++;
for(int i=1; i<=n; i++) buc[i] += buc[i-1];
for(int i=siz; i>=1; i--) ord[buc[len[i]]--] = i;
//Question 1
int u = rtn,cur = 0,ans1 = m+1;
for(int i=1; i<=n; i++)
{
while(u && son[u][a[i]]==0) {ans1 = min(ans1,cur+1); u = fa[u]; cur = len[fa[u]]+1;}
if(son[u][a[i]]!=0) {cur++; u = son[u][a[i]];}
else {ans1 = 1; u = rtn; cur = 0;}
}
if(ans1==m+1) printf("-1\n");
else printf("%d\n",ans1);
//Question 2
int ans2 = n+1;
for(int i=1; i<=n; i++)
{
int j = 0;
for(int k=i; k<=n; k++)
{
if(nxt[j][a[k]]<=m)
{
j = nxt[j][a[k]];
}
else
{
ans2 = min(ans2,k-i+1);
}
}
}
if(ans2==n+1) printf("-1\n");
else printf("%d\n",ans2);
//Question 3
int ans3 = n+1;
for(int i=1; i<=siz; i++) dp[i] = n+1; dp[rtn] = 0;
for(int i=1; i<=n; i++)
{
for(int j=siz; j>=1; j--)
{
int u = ord[j];
if(son[u][a[i]])
{
update(dp[son[u][a[i]]],dp[u]+1);
}
else
{
update(ans3,dp[u]+1);
}
}
}
if(ans3==n+1) printf("-1\n");
else printf("%d\n",ans3);
//Question 4
int ans4 = n+1;
for(int i=1; i<=n; i++) dp[i] = n+1; dp[0] = 0;
for(int i=1; i<=n; i++)
{
for(int j=m; j>=0; j--)
{
if(nxt[j][a[i]]<=m)
{
update(dp[nxt[j][a[i]]],dp[j]+1);
}
else
{
update(ans4,dp[j]+1);
}
}
}
if(ans4==n+1) printf("-1\n");
else printf("%d\n",ans4);
return 0;
}

BZOJ 4032 Luogu P4112 [HEOI2015]最短不公共子串 (DP、后缀自动机)的更多相关文章

  1. BZOJ.4032.[HEOI2015]最短不公共子串(DP 后缀自动机)

    题目链接 1.求A的最短子串,它不是B的子串. 子串是连续的,对B建SAM,枚举起点,在SAM上找到第一个无法匹配点即可.O(n)用SAM能做吗..开始想错了. 2.求A的最短子串,它不是B的子序列. ...

  2. 【BZOJ4032】[HEOI2015]最短不公共子串(后缀自动机,序列自动机)

    [BZOJ4032][HEOI2015]最短不公共子串(后缀自动机,序列自动机) 题面 BZOJ 洛谷 题解 数据范围很小,直接暴力构建后缀自动机和序列自动机,然后直接在两个自动机上进行\(bfs\) ...

  3. BZOJ4032: [HEOI2015]最短不公共子串(后缀自动机+序列自动机)

    题目描述 在虐各种最长公共子串.子序列的题虐的不耐烦了之后,你决定反其道而行之. 一个串的“子串”指的是它的连续的一段,例如bcd是abcdef的子串,但bde不是. 一个串的“子序列”指的是它的可以 ...

  4. BZOJ4032 [HEOI2015]最短不公共子串 【后缀自动机 + 序列自动机 + dp】

    题目链接 BZOJ4032 题解 首先膜\(hb\) 空手切神题 一问\(hash\),二问枚举 三问\(trie\)树,四问\(dp\) 南二巨佬神\(hb\) 空手吊打自动机 \(orz orz ...

  5. 洛谷 P4112 [HEOI2015]最短不公共子串 解题报告

    P4112 [HEOI2015]最短不公共子串 题目描述 在虐各种最长公共子串.子序列的题虐的不耐烦了之后,你决定反其道而行之. 一个串的"子串"指的是它的连续的一段,例如bcd是 ...

  6. luoguP4112 [HEOI2015]最短不公共子串 SAM,序列自动机,广搜BFS

    luoguP4112 [HEOI2015]最短不公共子串 链接 luogu loj 思路 子串可以用后缀自动机,子序列可以用序列自动机. 序列自动机是啥,就是能访问到所有子序列的自动机. 每个点记录下 ...

  7. BZOJ 4032: [HEOI2015]最短不公共子串 (dp*3 + SAM)

    转博客大法好 第4个子任务中,为什么只转移最近的一个位置,自己YY吧(多YY有益身体健康). #include <bits/stdc++.h> using namespace std; t ...

  8. BZOJ 4032: [HEOI2015]最短不公共子串

    4032: [HEOI2015]最短不公共子串 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 446  Solved: 224[Submit][Sta ...

  9. BZOJ 4032: [HEOI2015]最短不公共子串 后缀自动机 暴力

    4032: [HEOI2015]最短不公共子串 题目连接: http://www.lydsy.com/JudgeOnline/problem.php?id=4032 Description 在虐各种最 ...

随机推荐

  1. YTU 2705:用重载求距离

    2705: 用重载求距离. 时间限制: 1 Sec  内存限制: 128 MB 提交: 208  解决: 114 题目描述 使用函数重载的方法定义两个重名函数,分别求出整型数的两点间距离和浮点型数的两 ...

  2. E20171006-hm

    trace  vt. 跟踪,追踪; 追溯,探索; 探索; 查找;          vi. 沿着一小径或道路前进; 可以追溯的;            n. 痕迹; 痕迹,踪迹; 微量,极少量; [植 ...

  3. Akka源码分析-Remote-Actor创建

    在之前的博客中,我们分析过local模式下Actor的创建过程,最终还是调用了provider的actorOf的函数创建了Actor,在remote模式下provider就是RemoteActorRe ...

  4. 基于Spark Streaming预测股票走势的例子(一)

    最近学习Spark Streaming,不知道是不是我搜索的姿势不对,总找不到具体的.完整的例子,一怒之下就决定自己写一个出来.下面以预测股票走势为例,总结了用Spark Streaming开发的具体 ...

  5. 微信小程序获取地理位置

    小程序只支持获取当前位置的经纬度,并不能直接获取到地理名称,需要通过第三方来逆地址解析,这里我选择的是腾讯位置服务 在使用前需要去申请key,这里是地址:https://lbs.qq.com/cons ...

  6. Java初级进阶中高级工程师必备技能

    很多人学了javase以为自己学的已经很OK了,但是其实javase里边有很多的知识点是你不知道的,不管你找的是哪里的javase的视频,大多数是不会讲这些东西,而这些东西你平时业务又不会主动去接触, ...

  7. ACM_新七步诗(深搜)

    新七步诗 Time Limit: 2000/1000ms (Java/Others) Problem Description: 突然的一天,小锴做了一个梦,梦见自己来到了三国,而自己也成了梦寐以求的帅 ...

  8. EF在应用程序配置文件中找不到名为“XXX”的连接字符串。

    现象: 在配置EF的时候需要如题所述的问题,仔细检查了在EF实体模型对应程序集下的APP.Config文件中的ConnectionString配置项有了XXX项的数据库字符串的配置: <conn ...

  9. Spring Cloud (7) 服务容错保护-Hystrix服务降级

    在微服务架构中,根据业务来拆分成一个个的服务,服务与服务之间可以互相调用,在Spring Cloud可以用RestTemplate+Ribbon和Feign来调用.为了保证其高可用,单个服务通常会集群 ...

  10. ssh项目导入报the import javax.servlet cannot be resolved

    在做javaWeb项目时,我们经常会出现丢失包的情况,如下图所示的错误,我们应该怎么解决呢? 根据网上教程向工程中加入tomcat的servlet-api.jar和jsp-api.jar的包 此时项目 ...