hdu3948（后缀数组）
题意：给一串字符，需要你求不相同的回文子串个数.......
同ural1297，链接：http://www.cnblogs.com/ziyi--caolu/archive/2013/06/09/3129877.html
#include<iostream>
#include<stdio.h>
#include<string.h>
using namespace std;
#define maxn 210000
#define max(x,y) x>y? x:y
#define min(x,y) x>y? y:x
int dp[maxn][32],log2[maxn];
int wa[maxn],wb[maxn],wsf[maxn],wv[maxn],sa[maxn];
int rank[maxn],height[maxn],s[maxn];
char str[maxn];
int cmp(int *r,int a,int b,int k)
{
return r[a]==r[b]&&r[a+k]==r[b+k];
}
void getsa(int *r,int *sa,int n,int m)
{
int i,j,p,*x=wa,*y=wb,*t;
for(i=0;i<m;i++) wsf[i]=0;
for(i=0;i<n;i++) wsf[x[i]=r[i]]++;
for(i=1;i<m;i++) wsf[i]+=wsf[i-1];
for(i=n-1;i>=0;i--) sa[--wsf[x[i]]]=i;
p=1;
j=1;
for(;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++) wsf[i]=0;
for(i=0;i<n;i++) wsf[wv[i]]++;
for(i=1;i<m;i++) wsf[i]+=wsf[i-1];
for(i=n-1;i>=0;i--) sa[--wsf[wv[i]]]=y[i];
t=x;
x=y;
y=t;
x[sa[0]]=0;
for(p=1,i=1;i<n;i++)
x[sa[i]]=cmp(y,sa[i-1],sa[i],j)? p-1:p++;
}
}
void getheight(int *r,int n)
{
int i,j,k=0;
for(i=1;i<=n;i++) rank[sa[i]]=i;
for(i=0;i<n;i++)
{
if(k)
k--;
else
k=0;
j=sa[rank[i]-1];
while(r[i+k]==r[j+k])
k++;
height[rank[i]]=k;
}
}
void rmq_init(int n)
{
for(int i=1;i<=n;i++) dp[i][0]=height[i];
for(int j=1;(1<<j)<=n;j++)
for(int i=1;i+(1<<j)-1<=n;i++)
dp[i][j]=dp[i][j-1]<dp[i+(1<<(j-1))][j-1]? dp[i][j-1]:dp[i+(1<<(j-1))][j-1];
}
int rmq(int left,int right)
{
int k=log2[right-left+1];
//printf(":%d %d %d\n",k,left,right);
return min(dp[left][k],dp[right-(1<<k)+1][k]);
}
int main()
{
int text,xxx=0;
scanf("%d",&text);
int num = 0, tot = 1;
while(tot < maxn)
{
for(int i = tot; i < tot * 2 && i < maxn; ++i)
log2[i] = num;
num++; tot *= 2;
}
while(text--)
{
scanf("%s",str);
int len=strlen(str),j=0;
for(int i=0;i<len;i++)
{
s[j++]=str[i]-'a'+1;
}
s[j++]=29;
for(int i=len-1;i>=0;i--)
{
s[j++]=str[i]-'a'+1;
}
s[j]=0;
int len1=j;
getsa(s,sa,len1+1,32);
getheight(s,len1);
rmq_init(len1);
int sum=0;
int cnt=0;
for(int i=1;i<len1;i++)
{
cnt=cnt>height[i] ? height[i]:cnt;
if(sa[i]<len)
{
int tmp=2*len-sa[i];
int pp,qq;
if(i<rank[tmp])
pp=i;
else
pp=rank[tmp];
if(i>rank[tmp])
qq=i;
else
qq=rank[tmp];
//int pp=min(i,rank[tmp]);
int k=rmq(pp+1,qq);
//printf("%d %d\n",i,rank[tmp]);
//printf("%d %d %d\n",pp+1,qq,tmp);
//printf("1:%d %d %d\n",i,rank[tmp],k);
if(k>cnt)
{
sum+=k-cnt;
cnt=k;
}
}
}
cnt=0;
//printf("\n\n");
for(int i=1;i<len1;i++)
{
cnt= (cnt>height[i]) ? height[i] : cnt;
if(!sa[i])
continue;
if(sa[i]<len)
{
int tmp=2*len-sa[i]+1;
int pp,qq;
if(i<rank[tmp])
pp=i;
else
pp=rank[tmp];
if(i>rank[tmp])
qq=i;
else
qq=rank[tmp];
int k=rmq(pp+1,qq);
printf("%d %d\n",pp+1,qq);
//printf("%d %d\n",min(i,rank[tmp])+1,max(i,rank[tmp]));
//printf("2:%d %d %d\n",i,rank[tmp],k);
if(k>cnt)
{
sum+=k-cnt;
cnt=k;
}
}
}
//printf("Case #%d: %d\n",++xxx,sum);
}
return 0;
}