(模板) hdu-5685-Problem A(扩展欧几里得)

Problem Description

度熊手上有一本字典存储了大量的单词，有一次，他把所有单词组成了一个很长很长的字符串。现在麻烦来了，他忘记了原来的字符串都是什么，神奇的是他竟然记得原来那些字符串的哈希值。一个字符串的哈希值，由以下公式计算得到：

H(s)=∏i≤len(s)i=1(Si−28) (mod 9973)

Si代表 S[i] 字符的 ASCII 码。

请帮助度熊计算大字符串中任意一段的哈希值是多少。

Input

多组测试数据，每组测试数据第一行是一个正整数N，代表询问的次数，第二行一个字符串，代表题目中的大字符串，接下来N行，每行包含两个正整数a和b，代表询问的起始位置以及终止位置。

1≤N≤1,000

1≤len(string)≤100,000

1≤a,b≤len(string)

Output

对于每一个询问，输出一个整数值，代表大字符串从 a 位到 b 位的子串的哈希值。

Sample Input

2 ACMlove2015 1 11 8 10 1 testMessage 1 1

Sample Output

6891 9240 88

题解 ：

        暴力会超时，所以需要把所有积依次累乘存在一个数组中（类似于前缀和形式），然后得到哈希值=（大前缀积/小前缀积）%9973；又因为除法不满足分配律，所以需要用到逆元。用扩展欧几里德求逆元！！

#include<stdio.h>
#include<string.h>
#include<math.h>
typedef long long ll;
ll M=9973;

void exgcd(ll a,ll b,ll &d,ll &x,ll &y){
if(b==0){ d=a; x=1; y=0; return;}
else{
exgcd(b,a%b,d,y,x);
y-=(a/b)*x;
}
}
ll fun(ll a,ll b){
ll d,x,y;
exgcd(a,b,d,x,y);
if(d==1){
return (x+b)%b;
}
else
 return -1;
}
int main(){
int n,a,b;
    char s[100005];
    ll sum[100005];                 // 保存前缀积的数组
while(scanf("%d",&n)!=EOF){
// getchar();
scanf("%s",s+1);
memset(sum,0,sizeof(sum));
int l=strlen(s+1);
sum[0]=1;
for(int i=1;i<=l;i++){           // 当初因为 i<=strlen(s+1); 一直超时；小细节！！！
sum[i]=(ll)(sum[i-1]*(s[i]-28)%M);
}
for(int i=0;i<n;i++){
scanf("%d %d",&a,&b);
   long long h=fun(sum[a-1],M);        // 得到sum[a-1] (小前缀积) 关于M的逆元
printf("%lld\n",(h*sum[b])%M);
   }
}
return 0;
}