[acm/icpc2016ChinaFinal][CodeforcesGym101194] Mr. Panda and Fantastic Beasts

地址：http://codeforces.com/gym/101194

题目：略

思路：

　　这题做法挺多的，可以sam也可以后缀数组，我用sam做的。

　　1.我自己yy的思路（瞎bb的）

　　　　把第一个串建立sam，然后让其他串在上面跑。

　　　　每走到一个位置p，当前长度为num，就代表这个endpos集合里的长度小于等于num的字符串是被包含了，同时parent树的所有祖先节点也要标记为被包含。

　　　　这一步的具体做法是：用一个mi数组表示到达该节p点时的长度大于mi[p]时，才算未被包含（到达长度指的是从root出发经过某条路径到p的长度），所以其他串在sam上走到每个节点p，长度为num时，记录下mi[p]=max(mi[p],num);

　　　　再跑完所有串后，来更新parent树上的节点。先拓扑排序，然后如果mi[p]>0，mi[fa[p]]=len[p]，否则mi[p]=len[fa[p]]。(mi数组初始时为0）

　　　　这样可以在O(n)更新parent树.

　　　　然后再扫一遍节点，如果mi[p]<len[p]，anslen=min(anslen,mi[p]+1)，这样就得到了anslen。（anslen为最小串的长度）

　　　　在扫一遍节点把所有可行答案求个字典序最小的就行了。（这一步远远没有O(anslen*|p|）,p是可行答案集合的大小，时间复杂度可以看下一个做法的时间）

　　2.网上看到的做法

　　　　把其他串建立广义sam，然后让第一个串在上面跑。

　　　　这样每当失配的时候，沿parent树走到符合条件的节点p，则此时最小可行字符串长度为len[fa[p]]+2，拿他去更新答案即可。

　　　　这做法时间复杂度更高，因为字符串比较次数比做法一多（它可能比较了长度大于anslen的字符串），且建立sam的时间复杂度也高。

贴个做法2的代码把，一的不想写了

 #include <bits/stdc++.h>

 using namespace std;

 struct SAM
 {
     static const int MAXN = <<;//大小为字符串长度两倍
     static const int LetterSize = ;

     int tot, last, ch[MAXN][LetterSize], fa[MAXN], len[MAXN];
     int sum[MAXN], tp[MAXN], cnt[MAXN]; //sum,tp用于拓扑排序，tp为排序后的数组

     void init( void)
     {
         last = tot = ;
         len[] = ;
         memset( ch[], , sizeof ch[]);
     }

     void add( int x)
     {
         int p = last, np = last = ++tot;
         len[np] = len[p] + , cnt[last] = ;
         memset( ch[np], , sizeof ch[np]);
         while( p && !ch[p][x]) ch[p][x] = np, p = fa[p];
         if( p == )
             fa[np] = ;
         else
         {
             int q = ch[p][x];
             if( len[q] == len[p] + )
                 fa[np] = q;
             else
             {
                 int nq = ++tot;
                 memcpy( ch[nq], ch[q], sizeof ch[q]);
                 len[nq] = len[p] + , fa[nq] = fa[q], fa[q] = fa[np] = nq;
                 while( p && ch[p][x] == q)  ch[p][x] = nq, p = fa[p];
             }
         }
     }

     void toposort( void)
     {
         for(int i = ; i <= len[last]; i++)   sum[i] = ;
         for(int i = ; i <= tot; i++)   sum[len[i]]++;
         for(int i = ; i <= len[last]; i++)   sum[i] += sum[i-];
         for(int i = ; i <= tot; i++)   tp[sum[len[i]]--] = i;
     }

     void go(char *ss)
     {
         int mi=0x3f3f3f3f,l;
         for(int i=,p=,num=,slen=strlen(ss);i<slen;i++)
         {
             int c=ss[i]-'a';
             if(ch[p][c])    p=ch[p][c],num++;
             else
             {
                 while(p&&!ch[p][c]) p=fa[p];
                 if(p)
                     num=len[p]+,p=ch[p][c];
                 else
                     p=,num=;
                 if(num+<mi)    mi=num+,l=i-num;
                 else if(num+==mi)
                 for(int j=l,k=i-num;k<=i;k++,j++)
                 if(ss[j]!=ss[k])
                 {
                     if(ss[j]>ss[k]) l=i-num;
                     break;
                 }
             }
         }
         if(mi==0x3f3f3f3f)
             printf("Impossible");
         else
             for(int i=;i<mi;i++)
                 printf("%c",ss[i+l]);
     }
 } sam;

 char sa[],sb[];
 int main(void)
 {
     int t,n,cs=;cin>>t;
     while(t--)
     {
         sam.init();
         scanf("%d%s",&n,sa);
         for(int i=;i<n;i++)
         {
             scanf("%s",sb);
             for(char *p=sb;*p;p++)  sam.add(*p-'a');
             sam.last=;
         }
         printf("Case #%d: ",cs++);
         sam.go(sa);
         printf("\n");
     }
     return ;
 }