【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;
}