POJ2406【Power Strings】（字符串哈希模板题）

题目地址：POJ2406【Power Strings】

题目描述：
Given two strings a and b we define ab to be their concatenation. For example, if a = “abc” and b = “def” then ab = “abcdef”. If we think of concatenation as multiplication, exponentiation by a non-negative integer is defined in the normal way: a^0 = “” (the empty string) and a^(n+1) = a*(a^n).

输入描述：
Each test case is a line of input representing s, a string of printable characters. The length of s will be at least 1 and will not exceed 1 million characters. A line containing a period follows the last test case.

输出描述：
For each s you should print the largest n such that s = a^n for some string a.

输入描述：

abcd
aaaa
ababab
.

输出样例：

1
4
3

注释：
This problem has huge input, use scanf instead of cin to avoid time limit exceed.

题意：
给出两个字符串a和b，定义ab是它们的串联。
例如：如果a=“abc”，b=“def”，则ab=“abcdef”。
如果把串联视为相乘，非负整数指数则定义为：a^0=" "（空串），而a^(n+1)=a* (a^n)。
输入
每个测试用例是在一行中给出一个可打印字符的字符串s。s的长
度将至少为1，并且不会超过1000000（一百万）个字符。在最后
一个测试用例后面，给出包含句点的一行。
输出
对于每个s，请输出最大的n，使得对某个字符串a，s=a^n。
提示
本题海量输入，为避免超时，请使用scanf替代cin。

思路：
设字符串s的长度为L=strlen (s+1) ，字符串s的下标从1开始。
首先计算出字符串s中每个前缀的hash函数值，
即h[i] = h[i-1] * base + str[i] - ‘a’ + 1 (1<=i<=L);
利用unsigned long long溢出取模
然后按照长度递增的顺序枚举s中可能存在的相邻子串。
若L%x=0,则说明s中可能存在长度为x且满足相乘关系的相邻子串，
即对于等长子串 s[1…x],s[x+1…2x],s[(n-1)*x+1…l]
如果h[x] = h[x+1…2x] = … = h[(n-1)*x+1…l]
其中，s[l…r]的hash值为 h[r] - h[l - 1] * k^(r - l + 1)
k的n次方可以用一个数组p来存，
将上式改写为: h[r] - h[l - 1] * p[r - l + 1]
则相乘关系成立，即s为连续n个子串a，s=a^n。
由于此时子串长度x是最小的，因此次幂n=L/x为最大，
n即为问题的解。

字符串哈希讲解

代码：

#include<iostream>
#include<string.h>
using namespace std;
typedef unsigned long long ULL;
const int N = 1000010,base = 131;
char str[N];
ULL h[N],p[N];
int len; 
ULL get(int l,int r)//获得 l 到 r 的 hash 值 
{
	return h[r] - h[l - 1] * p[r - l + 1];
}
bool check(int x)//若所有长度为 x 的相邻子串对应的散列函数值相等，则返回 true；否则返回 false
{
	for(int i = x*2 ;i <= len;i += x)
		if(h[x] != get(i-x+1,i))
			return false;
	return true;
}
int main()
{
	while(~scanf("%s",str + 1))
	{
		len = strlen(str + 1);
		if(len == 1&&str[1] == '.')
			return 0;		
		h[0] = 0,p[0] = 1;
		for(int i = 1;i <= len;i ++ )
		{
			h[i] = h[i-1] * base + str[i] - 'a' + 1; 
			p[i] = p[i-1] * base;
		}		
		for(int i = 1;i <= len;i ++ )
			if(len%i == 0 && check(i))
			{
				printf("%d\n",len/i);
				break;
			}
	}
	return 0;
}