题干:
描述
很久很久以前,森林里住着一群兔子。有一天,兔子们想要研究自己的 DNA 序列。我们首先选取一个好长好长的 DNA 序列(小兔子是外星生物,DNA 序列可能包含 26 个小写英文字母),然后我们每次选择两个区间,询问如果用两个区间里的 DNA 序列分别生产出来两只兔子,这两个兔子是否一模一样。注意两个兔子一模一样只可能是他们的 DNA 序列一模一样。
输入格式
第一行一个 DNA 字符串 S。
接下来一个数字 m,表示 m 次询问。
接下来 m 行,每行四个数字 l1, r1, l2, r2,分别表示此次询问的两个区间,注意字符串的位置从1开始编号。
其中 1 ≤ length(S), m ≤ 1000000
输出格式
对于每次询问,输出一行表示结果。如果两只兔子完全相同输出 Yes,否则输出 No(注意大小写)
样例输入
aabbaabb
3
1 3 5 7
1 3 6 8
1 2 1 2
样例输出
Yes
No
Yes
来源
罗翔宇,北京大学2014年数据结构与算法A(实验班)期末考试
解题报告:
这题让他自然溢出或者加个模数都可以,,毕竟数据水,,但是你的模数大于1e10就不行了、、(那就相当于两个模数了啊)
AC代码:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<queue>
#include<map>
#include<vector>
#include<set>
#include<string>
#include<cmath>
#include<cstring>
#define ll long long
#define ull unsigned long long
#define pb push_back
#define pm make_pair
#define fi first
#define se second
using namespace std;
const int MAX = 2e6 + 5;
char s[MAX];
ull Hash[MAX];
ull p = 31,p2=131;
ull mod = 1e9 + 7;
ull qpow(ull a,ull k) {
ull res = 1;
while(k) {
if(k&1) res =(a*res)%mod;
k>>=1;
a=(a*a)%mod;
}
return res%mod;
}
int main()
{
cin>>s+1;
int len = strlen(s+1);
for(int i = 1; i<=len; i++) {
Hash[i] = (Hash[i-1]*p + s[i]-'a'+1)%mod;
}
int m,l1,l2,r1,r2;
cin>>m;
while(m--) {
scanf("%d%d%d%d",&l1,&r1,&l2,&r2);
ull tmp1 = (Hash[r1] + mod - (Hash[l1-1] * qpow(p,r1-l1+1))%mod)%mod;
ull tmp2 = (Hash[r2] + mod - (Hash[l2-1] * qpow(p,r2-l2+1))%mod)%mod;
if(tmp1 == tmp2) puts("Yes");
else puts("No");
}
return 0 ;
}
总结:
刚开始直接用mod=1e13,1e16,1e18各试了一次,,都wa了,,还以为是有Hash Crash了,,还弄了个Hash2数组用不同的种子p2去哈希了一次,两次Hash判断如果均相同才输出Yes。。。但是发现不是这里的问题,就是模数大打了。
其二,其实这题可以不用快速幂的,因为你都是用的p,所以可以先预处理一下次幂。