BZOJ2795: [Poi2012]A Horrible Poem
题解：
首先判断长度len是否是l-r的循环节可以做到O（1）
只需判断 l-l+len-1 和 r-len+1-r 是否相同即可，这一步可以用hash
然后循环节一定是总长度的约数。
这样可以做到n√n，但是会T。
一个鬼畜的优化是 枚举在每块内每个字母的个数，这个个数必须是总个数的约数，然后我们取一下gcd，再枚举约数。
实测可过。
另一种做法可以戳这里：http://blog.csdn.net/zeyu_king/article/details/41989279
代码：
1 #include<cstdio> 2 #include<cstdlib> 3 #include<cmath> 4 #include<cstring> 5 #include<algorithm> 6 #include<iostream> 7 #include<vector> 8 #include<map> 9 #include<set> 10 #include<queue> 11 #include<string> 12 #define inf 1000000000 13 #define maxn 600000 14 #define maxm 500+100 15 #define eps 1e-10 16 #define ull unsigned long long 17 #define pa pair<int,int> 18 #define for0(i,n) for(int i=0;i<=(n);i++) 19 #define for1(i,n) for(int i=1;i<=(n);i++) 20 #define for2(i,x,y) for(int i=(x);i<=(y);i++) 21 #define for3(i,x,y) for(int i=(x);i>=(y);i--) 22 #define mod 1000000007 23 using namespace std; 24 inline int read() 25 { 26 int x=0,f=1;char ch=getchar(); 27 while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} 28 while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();} 29 return x*f; 30 } 31 int n,m,a[maxn],c[maxn],s[maxn][30]; 32 ull b[maxn],power[maxn]; 33 const ull base=13131; 34 inline ull get(int x,int y){return b[x]-b[x+y]*power[y];} 35 inline int gcd(int x,int y){return y?gcd(y,x%y):x;} 36 inline int solve() 37 { 38 int x=read(),y=read(),z=y-x+1,len=y-x+1,t=0; 39 for1(i,26)z=gcd(z,s[y][i]-s[x-1][i]); 40 if(z==1)return len; 41 for(int i=1;i*i<=z;i++)if(z%i==0) 42 { 43 int tmp=len/z*i; 44 if(get(x,len-tmp)==get(y-(len-tmp)+1,len-tmp))return tmp; 45 c[++t]=len/i; 46 } 47 for3(i,t,1)if(get(x,len-c[i])==get(y-(len-c[i])+1,len-c[i]))return c[i]; 48 return len; 49 } 50 int main() 51 { 52 freopen("input.txt","r",stdin); 53 freopen("output.txt","w",stdout); 54 n=read();n=0; 55 char ch=getchar(); 56 while(ch<'a'||ch>'z')ch=getchar(); 57 for(;ch>='a'&&ch<='z';ch=getchar())a[++n]=ch-'a'+1; 58 for3(i,n,1)b[i]=b[i+1]*base+a[i]; 59 for1(i,n)for1(j,26)s[i][j]=s[i-1][j]+(a[i]==j); 60 power[0]=1; 61 for1(i,n)power[i]=power[i-1]*base; 62 m=read(); 63 while(m--)printf("%d\n",solve()); 64 return 0; 65 }
2795: [Poi2012]A Horrible Poem
Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 121 Solved: 72
[Submit][Status]
Description
给出一个由小写英文字母组成的字符串S，再给出q个询问，要求回答S某个子串的最短循环节。
如果字符串B是字符串A的循环节，那么A可以由B重复若干次得到。
Input
第一行一个正整数n (n<=500,000)，表示S的长度。
第二行n个小写英文字母，表示字符串S。
第三行一个正整数q (q<=2,000,000)，表示询问个数。
下面q行每行两个正整数a,b (1<=a<=b<=n)，表示询问字符串S[a..b]的最短循环节长度。
Output
依次输出q行正整数，第i行的正整数对应第i个询问的答案。
Sample Input
8
aaabcabc
3
1 3
3 8
4 8
aaabcabc
3
1 3
3 8
4 8
Sample Output
1
3
5
3
5
