【JZOJ 5272】 神奇的重复序列

Description

1<=|S|<=3000

Analysis

挖掘一下题目性质：
1.答案满足二分性（不过我并没有想到有什么卵用，如果有人想到较优复杂度的算法麻烦留言告知）
2.最关键的一点：考虑枚举两个相同子串的结尾位置 i 和 j（i < j），然后对于一个字符，考虑和它相同的字符有哪些。设 k=j-i，对于那些下标对 k 取模相同的位置，字符一定都是一样的，不同的互不影响。
然后如果枚举k，再枚举i或j，可以通过维护一些东西来做到O(n^2)复杂度

Code

#include<cstdio>#include<cstring>#include<algorithm>#define fo(i,a,b) for(int i=a;i<=b;i++)#define fd(i,b,a) for(int i=b;i>=a;i--)#define max(x,y) ((x)>(y)?(x):(y))#define min(x,y) ((x)<(y)?(x):(y))#define mset(a,x) memset(a,x,sizeof(a))using namespace std;typedef long long ll;const int N=3e3+5;char s[N];int n,p,l,r,k,use,ans,cnt[N][N],f[N][N],mx[N];int ins(int x){    x+=k;    int id=x%k,&y=cnt[id][s[x]];    --f[id][y],++f[id][++y];    if(y>mx[id]) mx[id]=y;else ++use;    return use;}int del(int x){    int id=x%k,&y=cnt[id][s[x]];    if(!--f[id][y] && y==mx[id]) mx[id]--;else --use;    ++f[id][--y];    return use;}int main(){    freopen("repeat.in","r",stdin);    //freopen("repeat.out","w",stdout);    scanf("%d\n%s",&p,s+1);n=strlen(s+1);    for(k=1;k<n;k++)    {        fo(i,1,n) f[i][0]=cnt[i][0]=0;        mset(mx,0);use=0;        fo(i,1,k) ins(i-k);        for(l=1,r=0;r+k<n;)        {            if(ins(++r)>p)                while(del(l++)>p);            ans=max(ans,r-l+1);        }        fo(i,l,r+k) del(i);    }    printf("%d",ans);    return 0;}