uva 11107 Life Forms

题意：给你N个串，求一个串在大于等于N/2的模板串中连续出现。如果有多解按字典序最小输出。

白书模板题。二分答案+合并模板串成一个新串，扫秒新串的height数组。

考查后缀数组+LCP

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 using namespace std;

 const int maxn =  *  + ;

 struct SuffixArray
 {
     int s[maxn];      // 原始字符数组（最后一个字符应必须是0，而前面的字符必须非0）
     int sa[maxn];     // 后缀数组
     int rank[maxn];   // 名次数组. rank[0]一定是n-1，即最后一个字符
     int height[maxn]; // height数组
     int t[maxn], t2[maxn], c[maxn]; // 辅助数组
     int n; // 字符个数

     void clear()
     {
         n = ;
         memset(sa, , sizeof(sa));
     }

     // m为最大字符值加1。调用之前需设置好s和n
     void build_sa(int m)
     {
         int i, *x = t, *y = t2;
         for(i = ; i < m; i++) c[i] = ;
         for(i = ; i < n; i++) c[x[i] = s[i]]++;
         for(i = ; i < m; i++) c[i] += c[i-];
         for(i = n-; i >= ; i--) sa[--c[x[i]]] = i;
         for(int k = ; k <= n; k <<= )
         {
             int p = ;
             for(i = n-k; i < n; i++) y[p++] = i;
             for(i = ; i < n; i++) if(sa[i] >= k) y[p++] = sa[i]-k;
             for(i = ; i < m; i++) c[i] = ;
             for(i = ; i < n; i++) c[x[y[i]]]++;
             for(i = ; i < m; i++) c[i] += c[i-];
             for(i = n-; i >= ; i--) sa[--c[x[y[i]]]] = y[i];
             swap(x, y);
             p = ;
             x[sa[]] = ;
             for(i = ; i < n; i++)
                 x[sa[i]] = y[sa[i-]]==y[sa[i]] && y[sa[i-]+k]==y[sa[i]+k] ? p- : p++;
             if(p >= n) break;
             m = p;
         }
     }

     void build_height()
     {
         int i, j, k = ;
         for(i = ; i < n; i++) rank[sa[i]] = i;
         for(i = ; i < n; i++)
         {
             if(k) k--;
             int j = sa[rank[i]-];
             while(s[i+k] == s[j+k]) k++;
             height[rank[i]] = k;
         }
     }
 };

 const int maxc =  + ; // 串的个数
 const int maxl =  + ; // 每个串的长度

 SuffixArray sa;
 int n;
 char word[maxl];
 int idx[maxn];
 int flag[maxc];

 // 子串[L,R) 是否符合要求
 bool good(int L, int R)
 {
     memset(flag, , sizeof(flag));
     if(R - L <= n/) return false;
     int cnt = ;
     for(int i = L; i < R; i++)
     {
         int x = idx[sa.sa[i]];
         if(x != n && !flag[x])
         {
             flag[x] = ;
             cnt++;
         }
     }
     return cnt > n/;
 }

 void print_sub(int L, int R)
 {
     for(int i = L; i < R; i++)
         printf("%c", sa.s[i] -  + 'a');
     printf("\n");
 }

 bool print_solutions(int len, bool print)
 {
     int L = ;
     for(int R = ; R <= sa.n; R++)
     {
         if(R == sa.n || sa.height[R] < len)   // 新开一段
         {
             if(good(L, R))
             {
                 if(print) print_sub(sa.sa[L], sa.sa[L] + len);
                 else
                     return true;
             }
             L = R;
         }
     }
     return false;
 }

 void solve(int maxlen)
 {
     if(!print_solutions(, false))
         printf("?\n");
     else
     {
         int L = , R = maxlen, M;
         while(L < R)
         {
             M = L + (R-L+)/;
             if(print_solutions(M, false)) L = M;
             else R = M-;
         }
         print_solutions(L, true);
     }
 }

 // 给字符串加上一个字符，属于字符串i
 void add(int ch, int i)
 {
     idx[sa.n] = i;
     sa.s[sa.n++] = ch;
 }

 int main()
 {
     int kase = ;
     while(scanf("%d", &n) ==  && n)
     {
         if(kase++ > ) printf("\n");
         int maxlen = ;
         sa.clear();
         for(int i = ; i < n; i++)
         {
             scanf("%s", word);
             int sz = strlen(word);
             maxlen = max(maxlen, sz);
             for(int j = ; j < sz; j++)
                 add(word[j] - 'a' + , i);
             add( + i, n); // 结束字符
         }
         add(, n);

         if(n == ) printf("%s\n", word);
         else
         {
             sa.build_sa( + n);
             sa.build_height();
             solve(maxlen);
         }
     }
     return ;
 }