SPOJ 220 Relevant Phrases of Annihilation（后缀数组）

最新推荐文章于 2020-01-25 17:48:19 发布

转载最新推荐文章于 2020-01-25 17:48:19 发布 · 83 阅读

本文提供了一种解决SPOJ PHRASES问题的有效算法，通过连接字符串并使用特殊字符隔开，利用后缀数组和高度数组进行二分查找，以找到在每个给定字符串中至少出现两次且不重叠的最长子串。

题目链接：http://www.spoj.pl/problems/PHRASES/

题意：给定n个字符串，求在每个字符串中至少出现两次且不重叠的最长字串。

思路：先将n个字符串连起来，中间用不相同的且没有出现在字符串中的字符隔开，求sa和height数组。然后二分答案，将后缀分组。判断的时候，要看是否有一组后缀在每个原来的字符串中至少出现两次，并且在每个原来的字符串中，后缀的起始位置的最大值与最小值之差是否不小于当前答案。




int r[N],sa[N],wa[N],wb[N],wd[N],rank[N],h[N];


int cmp(int *r,int a,int b,int L)
{
    return r[a]==r[b]&&r[a+L]==r[b+L];
}


void da(int *r,int *sa,int n,int m)
{
    int i,j,p,*x=wa,*y=wb,*t;
    FOR0(i,m) wd[i]=0;
    FOR0(i,n) wd[x[i]=r[i]]++;
    FOR1(i,m-1) wd[i]+=wd[i-1];
    FORL0(i,n-1) sa[--wd[x[i]]]=i;
    for(j=1,p=1;p<n;j<<=1,m=p)
    {
        p=0;
        FOR(i,n-j,n-1) y[p++]=i;
        FOR0(i,n) if(sa[i]>=j) y[p++]=sa[i]-j;
        FOR0(i,m) wd[i]=0;
        FOR0(i,n) wd[x[i]]++;
        FOR1(i,m-1) wd[i]+=wd[i-1];
        FORL0(i,n-1) sa[--wd[x[y[i]]]]=y[i];
        t=x;x=y;y=t;p=1;x[sa[0]]=0;
        FOR1(i,n-1) x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
    }
}



void calHeight(int *r,int *sa,int n)
{
    int i,j,k=0;
    FOR1(i,n) rank[sa[i]]=i;
    FOR0(i,n)
    {
        if(k) k--;
        j=sa[rank[i]-1];
        while(i+k<n&&j+k<n&&r[i+k]==r[j+k]) k++;
        h[rank[i]]=k;
    }
}

char s[N];
int n,b[N],M;


int OK(int mid)
{
    int i,j,k,t,p[15][2],L,R;
    FOR1(i,M)
    {
        L=i;
        while(L<=M&&h[L]<mid) L++;
        R=L;
        while(R<=M&&h[R]>=mid) R++;
        FOR0(j,n) p[j][0]=INF,p[j][1]=-INF;
        FOR(j,L-1,R-1)
        {
            k=b[sa[j]];
            if(k==n) continue;
            p[k][0]=min(p[k][0],sa[j]);
            p[k][1]=max(p[k][1],sa[j]);
        }
        FOR0(j,n)
        {
            if(p[j][0]==INF) break;
            if(p[j][1]-p[j][0]<mid) break;
        }
        if(j>=n) return 1;
        i=R;
    }
    return 0;
}

void deal()
{
    int low=1,high=10000,mid;
    while(low<=high)
    {
        mid=(low+high)>>1;
        if(OK(mid)) low=mid+1;
        else high=mid-1;
    }
    if(OK(low)) PR(low);
    else PR(high);
}

int main()
{
    int C;
    RD(C);
    while(C--)
    {
        RD(n);
        int L=0,i,j;
        FOR0(i,n)
        {
            RD(s);
            for(j=0;s[j];j++) r[L]=s[j]+100,b[L++]=i;
            if(i<n-1) r[L]=i+1,b[L++]=n;
        }
        r[L]=0;
        da(r,sa,L+1,250);
        calHeight(r,sa,L);
        M=L;
        deal();
    }
    return 0;
}