poj 3294 后缀数组求至少在k个串中出现的公共子串

做这道题伤心了，思路是对的，后缀数组的经典用法嘛。可是就是一直OLT，后来认真的思考了一下程序的逻辑，发现二分的时候，L应该从1开始，否则会有bug。

每次二分都写不好，真的伤心了 T_T。本以为这就ok了，结果还是WA。真的崩溃，后来没办法了，看了discuss，发现说G++肯能过，我就抱着试试的态度，提交上去结果RE。

我就发现是我有一个数组开小了。瞬间想哭，debug这么久，C++为什么就不提示RE呢。表示无语。

说了这么多废话，是因为真的想吐槽一下。

这个题的思路应该都知道，就是后缀数组把各个串连接起来，求height数组，然后二分分组，求公共前缀大于k的后缀有多少个，并且看它们是否分别在原来不同的串中。

关键就是实现细节了。其实认真的话，这道题是可以1A的。由于有一百多个串，中间加特殊字符就比较尴尬了，所以可以考虑把字符转化成整形来做，这样就可以很方便的加分隔符了。注意二分的时候，L=1，否则可能有bug ，反正我每次写二分都会出问题，表示很伤心也很郁闷。然后就是标记每个字符是在原来的那个串中，belong[j]=i;

大于等于k，进行分组的时候，要看各个后缀所在的串是否已经被标记，这样的话，要用到好多次memset，可能会很耗时。建议用不同的数来标记，这样应该会省点时间吧。具体的看代码，judge()里的appear，来进行不同标记。如果不懂，看一眼代码就知道什么意思了~

代码如下：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define MAXD 110010
#define MAXL 2111
using namespace std;
int r[MAXD], Rank[MAXD], height[MAXD];
int sa[MAXD], wa[MAXD], wb[MAXD], WS[MAXD], wv[MAXD];
int cmp(int *p, int x, int y, int l)
{
return p[x] == p[y] && p[x + l] == p[y + l];
}
void da(int n, int m)      //n为字符串长度，m为字符的最大值
{
int i, j, p, *x = wa, *y = wb, *t;
for(i = 0; i < m; i ++)
WS[i] = 0;
for(i = 0; i < n; i ++)
++ WS[x[i] = r[i]];
for(i = 1; i < m; i ++)
WS[i] += WS[i - 1];
for(i = n - 1; i >= 0; i --)
sa[-- WS[x[i]]] = i;
for(p = 1, j = 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 ++)
WS[i] = 0;
for(i = 0; i < n; i ++)
++ WS[wv[i]];
for(i = 1; i < m; i ++)
WS[i] += WS[i - 1];
for(i = n - 1; i >= 0; i --)
sa[-- WS[wv[i]]] = y[i];
for(t = x, x = y, y = t, x[sa[0]] = 0, p = 1, i = 1; i < n; i ++)
x[sa[i]] = cmp(y, sa[i - 1], sa[i], j) ? p - 1: p ++;
}
}
void calheight(int n)  //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 ++);
}
char str[1111];
int n,half,totlen,L,R;
int belong[MAXD],used[111];
int start[1110];
int length;
void input()
{
int i,j=0,k=0,len;
for(i=1;i<=n;i++)
{
gets(str);
len=strlen(str);
for(j=0;j<len;j++)
{
r[k]=str[j];
belong[k]=i;
k++;
}
r[k]=300+i;
belong[k++]=i;
}
r[k]=0;
totlen=k;
}
int judge(int k)
{
int i,j,num=0,t=0;
int appear=1;
memset(used,0,sizeof(used));
for(i=1;i<=totlen;i++)
{
if(height[i]>=k)
{
int th=belong[sa[i]];
if(used!=appear)
{
used=appear;
num++;
if(num==half-1)
{
start[t++]=sa[i];
length=k;
}
}
}
else
{
num=0;
appear++;
used[belong[sa[i]]]=appear;
}
}
return t;
}
int main()
{
int i,j;
while(scanf("%d",&n)!=EOF)
{
getchar();
if(n==0) break;
input();
memset(height,0,sizeof(height));
da(totlen+1,2000);
calheight(totlen);
half=n/2+1;
L=1;R=1005;
int T=0,t;
length=0;
while(L<=R)
{
int mid=(L+R)/2;
t=judge(mid);
if(t) {L=mid+1;T=t;}
else R=mid-1;
}
if(T==0) { printf("\?\n\n"); continue; }
for(i=0;i<T;i++)
{
for(j=start[i];j<start[i]+length;j++)
printf("%c",char(r[j]));
printf("\n");
}
printf("\n");
}
return 0;
}