POJ 1226

http://poj.org/problem?id=1226

题意：
给你一堆串，问你这些串中出现次数最多串的长度，反向出现的也算。

思路：因为问出现次数最多的串那么就要求这个串出现在每一个所给的串中，那么就是表明了需要出现N次，那么问题就简单明了了。就是N个串以及他们的反串中出现次数为K的串的最长长0度。那么我们只要 将height数组用一个数组分开就是要标记一下即可，之后二分答案，之后判断一下就好了。

代码很容易理解的

#include <stdio.h>
#include <iostream>
#include <string.h>
#include <algorithm>
#define maxs 2000010
using namespace std;
#define MME(i,j) memset(i,j,sizeof(i))
int s[maxs],sa[maxs],rank[maxs],height[maxs];
int wa[maxs],loc[maxs],wd[maxs],wb[maxs],wv[maxs];
char input[maxs];
int cmp(int *r,int a,int b,int k)
{
    return r[a]==r[b]&&r[a+k]==r[b+k];
}

void get_sa(int *r,int n,int m)
{
    int i,j,p,*x=wa,*y=wb;
    for(i=0;i<m;i++) wd[i]=0;
    for(i=0;i<n;i++) wd[x[i]=r[i]]++;
    for(i=1;i<m;i++) wd[i]+=wd[i-1];
    for(i=n-1;i>=0;i--) sa[--wd[x[i]]]=i;

    for(j=1,p=1;p<n;j*=2,m=p)
    {
        for(p=0,i=n-j;i<n;i++) y[p++]=i;
        for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;

        for(i=0;i<n;i++) wv[i]=x[y[i]];
        for(i=0;i<m;i++) wd[i]=0;

        for(i=0;i<n;i++) wd[wv[i]]++;
        for(i=1;i<m;i++) wd[i]+=wd[i-1];

        for(i=n-1;i>=0;i--) sa[--wd[wv[i]]]=y[i];

        for(swap(x,y),p=1,x[sa[0]]=0,i=1;i<n;i++){

            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
        }
    }
    /*  for(int i=0;i<n;i++)
        printf("SA[%d] is %d\n",i,sa[i]);*/
}

void build_height(int *r,int n)
{
    int i,j,k=0;
    for(i=1;i<=n;i++) rank[sa[i]]=i;
    for(i=0;i<n;height[rank[i++]]=k)
    {
        for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
    }
    /*for(int i=0;i<=n;i++)
        printf("h[%d] is %d\n",i,height[i]);*/
}

bool check(int mid,int len,int k)
{
    int ans=0;
    bool vis[1010];
    MME(vis,0);
    for(int i=2;i<=len;i++)
    {
        if(height[i]>=mid)
        {
            if(vis[loc[sa[i]]]==0)
            {
                ans+=1;
                vis[loc[sa[i]]]=1;
            }
            if(vis[loc[sa[i-1]]]==0)
                {
                    vis[loc[sa[i-1]]]=1;
                    ans++;
                }
            if(ans==k)
            return 1;
        }
        else
        {
            MME(vis,0);
            ans=0;
        }

    }
    return 0;
}

char rever[2*maxs];
int main()
{
    int n,t,len,templen=0;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        MME(s,0);
        char tempre[1010];
        //getchar();
        int spos=0,temp=140;
        for(int i=1;i<=n;i++)
        {
            scanf("%s",input);
            len=strlen(input);
           for(int k=0;k<len;k++)
           {

               loc[spos]=i;
               s[spos++]=input[k];
           }
           loc[spos]=temp;
           s[spos++]=temp++;

           for(int k=len-1;k>=0;k--)
           {
               loc[spos]=i;
               s[spos++]=input[k];
           }
          loc[spos]=temp;
          s[spos++]=temp++;
        }
     //   rever[spos]='\0';
        s[spos]=0;

        get_sa(s,spos+1,temp+1);
        build_height(s,spos);

        int left=1,right=strlen(input),mid,ans=0;
        while(right>=left)
        {
            mid=(left+right)/2;
            if(check(mid,spos,n))
            {
                ans=mid;
                left=mid+1;
            }
            else
                right=mid-1;
        }
        printf("%d\n",ans);
    }
    return 0;
}