bzoj4032: [HEOI2015]最短不公共子串(SAM+DP)

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

题目：传送门

---

题解：

　　 陈年老题良心%你赛膜爆嘎爷

　　 当初做题...一眼SAM...结果只会两种直接DP的情况...

　　 情况1: 直接设f[i][j] 表示的是a串的第i个位置和b串的第j个位置开始的最长公共前缀（灵感来源于SA）。然后就枚举开头直接瞎搞啊。。。

　　 情况2: 定义一个last[i][j] 表示b串的位置i后面的第一个字符j，贪心的思想直接做，匹配不到位置就直接记录答案了嘛

　　 情况3: 因为是在b串当中找子串，a中找子序列，那就考虑对b建SAM，然后定义da[] 表示的是a串的第i个位置结尾的串在b串中跑到第j个状态的最短长度（具体看代码）、

　　 情况4: 利用情况2预处理的last，贪心思想，仿照情况三把last当作自动机来跑

　　

---

代码：

 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<cmath>
 #include<algorithm>
 #define inf 1e9
 using namespace std;
 struct SAM
 {
     int son[];
 }ch[];int cnt,dep[],a[],fail[],da[],la,root,ss;、
 //da 表示的是a串的第i个位置结尾的串在b串中跑到第j个状态的最短长度
 void add(int k)
 {
     int x=a[k];
     int p=la,np=++cnt;dep[np]=k;
     while(p!= && ch[p].son[x]==)ch[p].son[x]=np,p=fail[p];
     if(p==)fail[np]=root;
     else
     {
         int q=ch[p].son[x];
         if(dep[p]+==dep[q])fail[np]=q;
         else
         {
             int nq=++cnt;dep[nq]=dep[p]+;
             ch[nq]=ch[q];fail[nq]=fail[q];fail[np]=fail[q]=nq;
             while(p && ch[p].son[x]==q)ch[p].son[x]=nq,p=fail[p];
         }
     }
     la=np;
 }
 int len1,len2,opt;
 int c[],last[][];//位置i后面的第一个字符j
 char sa[],sb[];
 int f[][];//分别以i,j位置开头的最长公共前缀
 int main()
 {
     //freopen("a.in","r",stdin);
     //freopen("a.out","w",stdout);
     //scanf("%d",&opt);
     scanf("%s",sa+);len1=strlen(sa+);
     scanf("%s",sb+);len2=strlen(sb+);
     //if(opt==1)
     //{
         memset(f,,sizeof(f));int ans=inf,len=,maxx=;
         for(int i=len1;i>=;i--)for(int j=len2;j>=;j--)if(sa[i]==sb[j])f[i][j]=f[i+][j+]+;
         for(int i=;i<=len1;i++)
         {
             maxx=;for(int j=;j<=len2;j++)maxx=max(maxx,f[i][j]);
             if(maxx!=len1-i+)ans=min(ans,maxx+);//加一定不相等
         }
         if(ans==inf)printf("-1\n");
         else printf("%d\n",ans);
     //}
     //else if(opt==2)
     //{
         len=,ans=inf;for(int i=;i<=;i++)c[i]=inf;sb[]='a';
         for(int i=len2;i>=;i--){for(int j=;j<=;j++)last[i][j]=c[j];c[sb[i]-'a'+]=i;}
         for(int i=;i<=len1;i++)
         {
             len=;
             for(int j=i;j<=len1;j++)
             {
                 len=last[len][sa[j]-'a'+];
                 if(len>len1){ans=min(ans,j-i+);break;}//跳出去证明没有了
             }
         }
         if(ans==inf)printf("-1\n");
         else printf("%d\n",ans);
     //}
 //    else if(opt==3)
 //    {
         cnt=;root=la=++cnt;memset(da,,sizeof(da));
         for(int i=;i<=len2;i++)a[i]=sb[i]-'a';for(int i=;i<=len2;i++)add(i);ss=;da[]=;ans=inf;
         for(int i=;i<=len1;i++)
             for(int j=;j<=cnt;j++)
             {
                 ss=ch[j].son[sa[i]-'a'];
                 if(ss==)ans=min(ans,da[j]+);else da[ss]=min(da[ss],da[j]+);
             }
         if(ans==inf)printf("-1\n");
         else printf("%d\n",ans);
 //    }
     //else if(opt==4)
     //{
         ans=inf;//for(int i=1;i<=50;i++)c[i]=inf;sb[0]='a';
         //for(int i=len2;i>=0;i--){for(int j=1;j<=26;j++)last[i][j]=c[j];c[sb[i]-'a'+1]=i;}
         memset(da,,sizeof(da));da[]=;ss=;
         for(int i=;i<=len1;i++)
             for(int j=len2;j>=;j--)
             {
                 ss=last[j][sa[i]-'a'+];
                 if(ss>len1)ans=min(ans,da[j]+);else da[ss]=min(da[ss],da[j]+);
             }
         if(ans==inf)printf("-1\n");
         else printf("%d\n",ans);
     //}
     return ;
 }