bzoj2320: 最多重复子串

2320: 最多重复子串

Time Limit: 40 Sec  Memory Limit: 128 MB
Submit: 246  Solved: 66
[Submit][Status][Discuss]

Description

一个字符串P的重复数定义为最大的整数R，使得P可以分为R段连续且相同的子串。比方说，“ababab”的重复数为3，“ababa”的重复数为1。

Your Task

对于给定的串S，找出S的一个子串K使得K的重复数最大。

 

Input

第一行T表示数据组数

对于每组数据，一行中一个仅包含小写字母的字符串S

 

Output

对于每组数据，在一行中输出K，如果有多个解，输出字典序最小的那一个

 

Sample Input

2
ccabababc
daabbccaa

Sample Output

ababab
aa

HINT

 

100%：T≤10，S的长度不超过100000

 

题解

……我就来贴个代码吧……orz下claris……

#include<cstdio>
#include<cstring>
#include<algorithm>
typedef unsigned long long ll;
const int N=100010,base=127;
int ans,k,st,ed,n,m,time,i,l;char s[N];ll f[N],pow[N];
using namespace std;
ll hash(int l,int r){return f[r]-f[l-1]*pow[r-l+1];}
int lcs(int x,int y)
{
    int l=0,r=min(x,y);
    while(l<r){
        int mid=(l+r)/2;
        if(hash(x-mid,x)==hash(y-mid,y))l=mid+1;else r=mid;
    }
    return l;
}
int lcp(int x,int y)
{
    int l=0,r=min(n-x+1,n-y+1);
    while(l<r){
        int mid=(l+r)/2;
        int a=hash(x,x+mid);
        int b=hash(y,y+mid);
        if(hash(x,x+mid)==hash(y,y+mid))l=mid+1;else r=mid;
    }
    return l;
}
void up(int x,int y)
{
    if(k>ans){ans=k;st=x;ed=y;return;}
    int l0=ed-st+1,l1=y-x+1,t=lcp(x,st);
    if(t==min(l0,l1)){
        if(t==l0)return;
        st=x;ed=y;return;
    }
    if(s[st+t]>s[x+t]){st=x;ed=y;}
}
int main()
{
    scanf("%d",&time);
    for(i=1,pow[0]=1;i<=N;i++)pow[i]=pow[i-1]*base;
    while(time--)
    {
        scanf("%s",s+1);n=strlen(s+1);
        for(ans=st=ed=i=1;i<=n;i++)
            if(s[i]<s[st])st=ed=i;
        for(f[0]=0,i=1;i<=n;i++)f[i]=f[i-1]*base+s[i];
        for(int i=1;i<=n/2;i++)
            for(int j=1;j<=n;j+=i){
                int x=j+i;if(x>n)break;
                if(s[j]!=s[x])continue;
                int a=lcp(j,x),b=lcs(j,x);
                k=(a+b-1)/i+1;
                l=k*i;
                if(k>=ans)for(int y=j-b+1;y<=i+j+a-l;y++)up(y,y+l-1);
            }
        for(int i=st;i<=ed;i++)putchar(s[i]);puts("");
    }
    return 0;
}