2016"百度之星" - 资格赛(Astar Round1)
http://bestcoder.hdu.edu.cn/contests/contest_showproblem.php?cid=690&pid=1001
Problem Description
度熊手上有一本字典存储了大量的单词,有一次,他把所有单词组成了一个很长很长的字符串。现在麻烦来了,他忘记了原来的字符串都是什么,神奇的是他竟然记得原来那些字符串的哈希值。一个字符串的哈希值,由以下公式计算得到:(公式略,为连乘求模) 请帮助度熊计算大字符串中任意一段的哈希值是多少。
Input
多组测试数据,每组测试数据第一行是一个正整数NN,代表询问的次数,第二行一个字符串,代表题目中的大字符串,接下来NN行,每行包含两个正整数aa和bb,代表询问的起始位置以及终止位置。
1≤N≤1,000
1≤len(string)≤100,000
1≤a,b≤len(string)
Output
对于每一个询问,输出一个整数值,代表大字符串从 aa 位到 bb 位的子串的哈希值。
Sample Input
Copy
2
ACMlove2015
1 11
8 10
1
testMessage
1 1
Sample Output
6891
9240
88
分析
我AC这题的代码采用的是线段树,但看到别人用的方法更好,这边介绍别人方法。
首先是前缀积:
H[0] = 1;
for(int i = 1 ;i<=len;i++)
{
H[i]=H[i-1]*(Hstr[i-1]-28)%mods;
}
然后是
逆元:
为何可以用逆元 http://blog.csdn.net/cqlf__/article/details/7953039 摘取博客部分:
比如: (8/2)%5 我们求a*b*c*d*e*f*g..../z 前面乘积部分LL存不下所以要一边mod一边乘。最后处理到除z时,不一定能除尽。比如前面那个例子,8%5=3,3除不尽2就乘以2%5的逆元在%5
2%5的逆元=2^(5-2)=8 这是计算逆元的一种方法,后面讲。还有一直哦你方法是扩展欧几里德算法也是后面详细讲。
(3*8)%5=4=4%5
===============
在计算(a/b)%Mod时,往往需要先计算b%Mod的逆元p(b有逆元的条件是gcd(b,Mod)==1,显然素数肯定有逆元),然后由(a*p)%Mod得结果c。这里b的逆元p满足(b*p)%Mod=1。先来简单证明一下:
(a/b)%Mod=c; (b*p)%Mod=1; ==》 (a/b)*(b*p) %Mod=c; ==》 (a*p)%Mod=c;
从上面可以看出结论的正确性,当然这里b需要是a的因子。接下来就需要知道根据b和Mod,我们怎么计算逆元p了。扩展欧几里德算法,大家应该都知道,就是已知a、b,求一组解(x,y)使得a*x+b*y=1。这里求得的x即为a%b的逆元,y为b%a的逆元(想想为什么?把方程两边都模上b或a看看)。调用ExtGcd(b,Mod,x,y),x即为b%Mod的逆元p。
求b%Mod的逆元p还有另外一种方法,即p=b^(Mod-2)%Mod,因为b^(Mod-1)%Mod=1(这里需要Mod为素数)。
逆元详解 http://blog.csdn.net/acdreamers/article/details/8220787
重点(1/a)(mod m)=a^(m-2)。注意要像本题中这样使用逆元,要满足素数条件。
对于a*b*c/(x*y*z),等于(a*x的逆元)*(b*y的逆元)*(c*z的逆元),求大数组合的时候可以利用这点。
最后是快速幂取模,不要求素数条件:
原理为:
a为底,b为幂
可以把b按二进制展开为:b = p(n)*2^n + p(n-1)*2^(n-1) +…+ p(1)*2 + p(0)
其中p(i) (0<=i<=n)为 0 或 1
这样 a^b = a^ (p(n)*2^n + p(n-1)*2^(n-1) +...+ p(1)*2 + p(0))
= a^(p(n)*2^n) * a^(p(n-1)*2^(n-1)) *...* a^(p(1)*2)*a^p(0)
LL mod_pow(LL x,LL n,LL mod) { // x是底数,n是幂数,mod是取余数
LL res = 1;
x%=mod;
while(n>0){
if(n & 1){
res = res * x % mod;
}
x = x * x % mod;
n>>=1;
}
return res;
}
《剑指offer》p93 面试题11数值的整数次方也介绍了快速幂方法,代码为递归形式。
另一种求逆元方法:http://www.w2bc.com/article/137005
inv[1] = 1;
for (int i = 2; i<MAXN; i++)
inv[i] = inv[MOD%i]*(MOD-MOD/i)%MOD;
代码
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;
const int MAXN = 1e5 + 5;
int H[MAXN];
char Hstr[MAXN];
int N,l,r;
const int mods = 9973;
typedef long long LL;
//快速幂取余
LL mod_pow(LL x,LL n,LL mod) {
LL res = 1;
//x%=mod;
while(n>0){
if(n & 1){
res = res * x % mod;
}
x = x * x % mod;
n>>=1;
}
return res;
}
int main ()
{
freopen("in.txt","r",stdin);
while(scanf("%d",&N)!=EOF){
scanf("%s",Hstr);
int len =strlen(Hstr);
H[0] = 1;
for(int i = 1 ;i<=len;i++){
H[i]=H[i-1]*(Hstr[i-1]-28)%mods;
}
while(N--){
scanf("%d%d",&l,&r);
printf("%I64d\n",(LL)H[r]*mod_pow(H[l-1],mods-2,mods)%mods);
}
}
}