hdu 5685 Problem A (逆元)

题目

题意：H(s)=∏i≤len(s)i=1(Si−28) (mod 9973)，求一个字符串 子串(a 位到 b 位的)的哈希值。这个公式便是求字符串哈希值的公式，（字符的哈希值 = 字符的ASCII码 - 28），字符串的哈希值等于字符的哈希值的乘积（ ∏ 这个就是累乘符号 ）。

化简过后的题意就是求一段序列中的区间乘，由于询问次数比较多，直接求乘会超时。

就比如如下代码：

#include<stdio.h>char s[100010];int main(){    int T,a,b;    while(~scanf("%d",&T))    {        int ans = 1;        scanf("%s",s+1);        while(T--){                ans =1;            scanf("%d%d",&a,&b);            for(int i=a;i<=b;i++)                ans  = ans * (s[i] - 28 ) % 9973;                 printf("%d\n",ans);        }    }    return 0;}

所以用前缀乘的方法，例如：0到b 区间的哈希值，除以 0到a 区间的哈希值，得到的就是 a 到 b 的哈希值。

若用preMulit[i]表示前i个序列的前缀乘。要求[l,r]区间内的的数字全乘起来，那么用preMulit[r] / preMulit[l-1] 即可，考虑到取模的问题，用到逆元。

那么我就要问了，为什么考虑到取模的问题，就要用到逆元呢？而且什么是逆元呢？

（看链接）

接下来就是解题代码了

代码一：

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;typedef long long ll;const ll maxn=100005;const ll mod=9973;ll sum[maxn],inv[maxn],re[maxn];char s[maxn];int main(){    inv[1]=1;    for(int i=2;i<maxn;++i)    {        inv[i]=inv[mod%i]*(mod-mod/i)%mod;    }    int n;    //freopen("shuju.txt","r",stdin);    while(~scanf("%d",&n))    {        scanf("%s",s+1);        sum[0]=re[0]=1;        for(int i=1;s[i]!=0;++i)        {            sum[i]=(sum[i-1]*(s[i]-28))%mod;            re[i]=inv[sum[i]];        }        for(int i=0;i<n;++i)        {            ll a,b;            scanf("%lld%lld",&a,&b);            ll tp=re[a-1];            printf("%I64d\n",(sum[b]*tp)%mod);        }    }    return 0;}

代码二：（通过extgcd直接求）

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;typedef long long ll;const ll maxn=100005;const ll mod=9973;ll sum[maxn];char s[maxn];void extgcd(ll a,ll b,ll &x,ll &y){    if(!b)    {        x=1;y=0;        return;    }    extgcd(b,a%b,y,x);    y-=(a/b)*x;}ll inv(ll a,ll n){    ll d,x,y;    extgcd(a,n,x,y);    return (x+n)%n;}int main(){    int n;    //freopen("shuju.txt","r",stdin);    while(~scanf("%d",&n))    {        scanf("%s",s+1);        sum[0]=1;        for(int i=1;s[i]!=0;++i)        {            sum[i]=(sum[i-1]*(s[i]-28))%mod;        }        for(int i=0;i<n;++i)        {            ll a,b;            scanf("%lld%lld",&a,&b);            ll tp=inv(sum[a-1],mod);            printf("%I64d\n",(sum[b]*tp)%mod);        }    }    return 0;}