http://www.elijahqi.win/archives/3195
Description
在虐各种最长公共子串、子序列的题虐的不耐烦了之后,你决定反其道而行之。
一个串的“子串”指的是它的连续的一段,例如bcd是abcdef的子串,但bde不是。
一个串的“子序列”指的是它的可以不连续的一段,例如bde是abcdef的子串,但bdd不是。
下面,给两个小写字母串A,B,请你计算:
(1) A的一个最短的子串,它不是B的子串
(2) A的一个最短的子串,它不是B的子序列
(3) A的一个最短的子序列,它不是B的子串
(4) A的一个最短的子序列,它不是B的子序列
Input
有两行,每行一个小写字母组成的字符串,分别代表A和B。
Output
输出4行,每行一个整数,表示以上4个问题的答案的长度。如果没有符合要求的答案,输出-1.
Sample Input
aabbcc
abcabc
Sample Output
2
4
2
4
HINT
对于100%的数据,A和B的长度都不超过2000
我好菜..
orz题解http://jiruyi910387714.is-programmer.com/posts/90425.html
第一问最长公共后缀 n^2枚举所有a的后缀然后最长公共后缀+1即是答案
第二问 枚举在A中的起点然后 设转移数组nxt 贪心转移即可 尽量多匹配 使得 当前起点的A以后不可能满足答案
第三问 对B 建SAM 然后n^2递推dp 表示dp[i][j]当前在i位置 在SAM上j位置的最短匹配 然后如果不匹配了+1即可
第四问 同第三问但需要加入贪心转移即可 使得当前A结尾的序列于B开头的后缀序列尽可能匹配 如果不匹配了就一定是答案
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=2200;
const int inf=0x3f3f3f3f;
int last=1,cnt=1,root=1;
int fa[N<<1],ch[N<<1][26],len[N<<1],a[N],b[N],nxt[N][26],n,m,lst[N],dp[N][N],dp1[N<<1],ans;
char s1[N],s2[N];
inline void insert1(int x){
int p=last,np=++cnt;len[np]=len[p]+1;
for (;p&&!ch[p][x];p=fa[p]) ch[p][x]=np;
if (!p) fa[np]=root;else{
int q=ch[p][x];if (len[p]+1==len[q]) fa[np]=q;else{
int nq=++cnt;memcpy(ch[nq],ch[q],sizeof(ch[q]));fa[nq]=fa[q];fa[q]=fa[np]=nq;
len[nq]=len[p]+1;for (;p&&ch[p][x]==q;p=fa[p]) ch[p][x]=nq;
}
}last=np;
}
int main(){
freopen("bzoj4032.in","r",stdin);
scanf("%s",s1+1);scanf("%s",s2+1);
n=strlen(s1+1);m=strlen(s2+1);ans=inf;
for (int i=1;i<=n;++i) a[i]=s1[i]-'a';
for (int i=1;i<=m;++i) b[i]=s2[i]-'a';
for (int i=0;i<26;++i) lst[i]=m+1;
for (int i=m;~i;--i){
for (int j=0;j<26;++j) nxt[i][j]=lst[j];lst[b[i]]=i;
}
for (int i=n;i;--i)
for (int j=m;j;--j)
if (a[i]==b[j]) dp[i][j]=dp[i+1][j+1]+1;
for (int i=1;i<=n;++i){static int ans1;ans1=0;
for (int j=1;j<=m;++j) ans1=max(ans1,dp[i][j]);
if (ans1!=n-i+1) ans=min(ans,ans1+1);
}printf("%d\n",ans==inf?-1:ans);ans=inf;
for (int i=1,now;i<=n;++i){now=0;
for (int j=i;j<=n;++j){
int y=nxt[now][a[j]];
if (y>m){ans=min(j-i+1,ans);break;}now=nxt[now][a[j]];
}
}
printf("%d\n",ans==inf?-1:ans);
memset(dp1,0x3f,sizeof(dp1));dp1[1]=0;ans=inf;
for (int i=1;i<=m;++i) insert1(b[i]);
for (int i=1;i<=n;++i)
for (int j=1,y;j<=cnt;++j){
y=ch[j][a[i]];
if (!y) ans=min(ans,dp1[j]+1);else dp1[y]=min(dp1[y],dp1[j]+1);
}
printf("%d\n",ans==inf?-1:ans);
memset(dp1,0x3f,sizeof(dp1));dp1[0]=0;ans=inf;
for (int i=1;i<=n;++i){
for (int j=m;~j;--j){
int y=nxt[j][a[i]];
if (y>m) ans=min(ans,dp1[j]+1);else dp1[y]=min(dp1[j]+1,dp1[y]);
}
}printf("%d\n",ans==inf?-1:ans);
return 0;
}