您的位置:首页 > 其它

BZOJ2795: [Poi2012]A Horrible Poem

2016-06-22 18:52 357 查看
题目大意:给一个字符串和多个询问,每次询问区间内最短循环节长度

假设我们找到一个长度,然后变成判定问题
如何判定呢?假如当前枚举长度为x,字符串总长为L,只需判断两点:
1.x是不是L的约数
2.这个字符串移除前x个字符和移除后x个字符之后的字符串是否完全匹配
因为满足条件2时,可以证明gcd(x,L)一定是一个合法答案,又因为x是L的约数,所以x就是一个合法答案

怎么判断字符串是否相等呢?可以hash,hash虚了可以用双hash或者后缀数组+RMQ
怎么枚举长度呢?直接找26个字母出现次数和总长度的公约数,然后根号L的时间复杂度就可以出解啦!

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 500010
using namespace std;
char s
;
int pre
[27];
int gcd(int a,int b)
{
if(!a) return b;
return gcd(b%a,a);
}
long long h
,bas=233,mod=1e9+7,cf
;
int geth(int x,int l)
{
return ((h[x+l-1]-h[x-1]*cf[l])%mod+mod)%mod;
}
bool judge(int x,int y,int l)
{
int L=y-x+1;
if(l==y-x+1) return true;
if(geth(x,L-l)==geth(x+l,L-l)) return true;
return false;
}
int main()
{
int n;
scanf("%d",&n);
scanf("%s",s+1);
int q,x,y,i,j;
scanf("%d",&q);
int maxl;
cf[0]=1;
for(i=1;i<=n;i++)
{
cf[i]=cf[i-1]*bas%mod;
h[i]=(h[i-1]*bas+s[i]-96)%mod;
for(j=1;j<=26;j++)
pre[i][j]=pre[i-1][j];
pre[i][s[i]-96]++;
}
while(q--)
{
scanf("%d%d",&x,&y);
maxl=y-x+1;
for(i=1;i<=26;i++)
maxl=gcd(maxl,pre[y][i]-pre[x-1][i]);
int ans=0,L=y-x+1;
for(i=1;i<=sqrt(maxl);i++)
if(maxl%i==0)
{
if(judge(x,y,L/(maxl/i))) {ans=L/(maxl/i);break;}
else if(judge(x,y,L/i)) ans=L/i;
}
printf("%d\n",ans);
}
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: