zoj3296 Connecting the Segments,后缀数组、dp

zoj3296 Connecting the Segments,后缀数组、dp

又是后缀数组。。嗯

求一个串最少能由几个回文串拼成。这里的拼接可以相互覆盖。

其实思路不难。先求出串中有哪些是回文串，然后就是经典的DP问题了，大家都应该懂的。
不过回文串最差可能有O(len^2)个。

思考一下，有一个重要的结论，如果一个回文串被另一个包含，那么这一个回文串是没有意义的，因为用另一个总是更好。
这样就简单了，对于每个点，我们求以这个点为中心的回文串最长的左右边界。
i从小到大做，将这些回文串排序，然后每次更新最大的右边界，每次f[i]只需要更新能拼到的最大的f[j]。

就搞定了。

要注意串的长度为奇数和偶数哦。。。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define maxn 101000
#define F(x) ((x)/3+((x)%3==1?0:tb))
#define G(x) ((x)<tb?(x)*3+1:((x)-tb)*3+2)

int r[maxn*3],sa[maxn*3];//DC3中r与sa数组的大小为3×maxn!!
int rank[maxn],height[maxn];
int wa[maxn],wb[maxn],wv[maxn],wd[maxn];

int c0(int *r,int a,int b) {return r[a]==r[b]&&r[a+1]==r[b+1]&&r[a+2]==r[b+2];}
int c12(int k,int *r,int a,int b){
     if(k==2) return r[a]<r[b]||r[a]==r[b]&&c12(1,r,a+1,b+1);
     else return r[a]<r[b]||r[a]==r[b]&&wv[a+1]<wv[b+1];
}
void sort(int *r,int *a,int *b,int n,int m){
    int i;
    for(i=0;i<n;i++) wv[i]=r[a[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--) b[--wd[wv[i]]]=a[i];
}
void dc3(int *r,int *sa,int n,int m){//n=r最后一个元素下标+1,strlen(s)+1
    int i,j,*rn=r+n,*san=sa+n,ta=0,tb=(n+1)/3,tbc=0,p;
    r[n]=r[n+1]=0;
    for(i=0;i<n;i++) if(i%3!=0) wa[tbc++]=i;
    sort(r+2,wa,wb,tbc,m);
    sort(r+1,wb,wa,tbc,m);
    sort(r,wa,wb,tbc,m);
    for(p=1,rn[F(wb[0])]=0,i=1;i<tbc;i++)
        rn[F(wb[i])]=c0(r,wb[i-1],wb[i])?p-1:p++;
    if(p<tbc) dc3(rn,san,tbc,p);
    else for(i=0;i<tbc;i++) san[rn[i]]=i;
    for(i=0;i<tbc;i++) if(san[i]<tb) wb[ta++]=san[i]*3;
    if(n%3==1) wb[ta++]=n-1;
    sort(r,wb,wa,ta,m);
    for(i=0;i<tbc;i++) wv[wb[i]=G(san[i])]=i;
    for(i=0,j=0,p=0;i<ta && j<tbc;p++)
        sa[p]=c12(wb[j]%3,r,wa[i],wb[j])?wa[i++]:wb[j++];
    for(;i<ta;p++) sa[p]=wa[i++];
    for(;j<tbc;p++) sa[p]=wb[j++];
}

void calheight(int *r,int *sa,int *height,int *rank,int n){//n==r的最后一个元素(添加的最小元素)，strlen(s);
    int i,j,k;
    for (i=1;i<=n;i++)rank[sa[i]]=i;
    for (i=0,height[0]=k=0;i<n;height[rank[i++]]=k)
         for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
}


//LCP----
int mi[maxn][17];
void lcpinit(int len,int *height){
     int i,j,m=-1,tt=len;
     while(tt) {m++;tt>>=1;}
     for(i=1;i<=len;++i)   mi[i][0]=height[i];
     for(j=1;j<=m;++j)
        for(i=1;i<=len;++i){
            if (i+(1<<(j-1))>len) mi[i][j]=mi[i][j-1];
            else mi[i][j]=min(mi[i][j-1],mi[(1<<(j-1))+i][j-1]);
        }
}

int lcpquery(int l,int r){//注意l=r需要特判
    l=rank[l];r=rank[r];
    if (l>r) swap(l,r);
    l++;
    int m=-1,tt=r+1-l;
    while(tt) {m+=1;tt>>=1;}
    return min(mi[l][m],mi[r+1-(1<<m)][m]);
}
//------------------

char s[maxn];
int len,l1;
int h[maxn][2],f[maxn];
pair<int,int> sv[maxn*2];

int apple(int l1,int tsv){
    int i,j;
    for(i=0;i<l1;++i){
        f[i]=i+1;
    }
    sort(sv+1,sv+tsv+1);
    int p=1,ma=0;
    while(p<=tsv&&sv[p].first==0) {
        ma=max(ma,sv[p].second);
        p++;
    }
    f[ma]=1;
    for(i=ma;i<l1;++i){
        while(p<=tsv&&sv[p].first<=i+1){
            ma=max(ma,sv[p].second);
            p++;
        }
        f[ma]=min(f[i]+1,f[ma]);
    }
    return f[l1-1];
}

int main(){
    //freopen("3296in.txt","r",stdin);
    int i,l1,tsv,ans;
    while(scanf("%s",&s)!=EOF){
        len=0;
        l1=strlen(s);
        for(i=0;i<l1;++i){
            r[len++]=s[i]-'a'+2;
        }
        r[len++]=1;//注意这两个分隔符不能相同
        for(i=l1-1;i>=0;--i){
            r[len++]=s[i]-'a'+2;
        }
        r[len]=0;

        dc3(r,sa,len+1,35);
        calheight(r,sa,height,rank,len);
        lcpinit(len,height);

        tsv=0;

        for(i=0;i<l1;++i){
            h[i][0]=lcpquery(i,l1+l1-i); //回文长度为奇数
            ++tsv;
            sv[tsv]=make_pair(i+1-h[i][0],i-1+h[i][0]);

            if (i!=0) {               //长度为偶数
                h[i][1]=lcpquery(i,l1+l1-i+1);
                ++tsv;
                sv[tsv]=make_pair(i-h[i][1],i+h[i][1]-1);
            }
        }
        ans=apple(l1,tsv);
        printf("%d\n",ans-1);
    }
    return 0;
}