hdu5685（乘法逆元）

最新推荐文章于 2021-02-24 00:30:01 发布

reallsp

最新推荐文章于 2021-02-24 00:30:01 发布

阅读量457

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/reallsp/article/details/77716367

版权

数论专栏收录该内容

6 篇文章 0 订阅

订阅专栏

hdu5685

题意：给定字符串，定义字符串哈希值的求法为：

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

求从给定区间的子串的哈希值。

思路：暴力求法超时，可以预处理存下（0，i）间的哈希值Hash[i][。则Hash[a,b]=Hash[b]/Hash[a-1],

注意运用乘法逆元解决精度问题。

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<map>
#include<string>
#include<vector>
#include<climits>
#include<queue>
#include<set>
#include<list>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
const int maxn = 100000;
int Hash[maxn+1];
ll exgcd(ll a,ll b,ll &x,ll &y)// 扩展欧几里得
{
    if(b==0)
    {
        x=1;y=0;
        return a;
    }
    ll temp=exgcd(b,a%b,y,x);
    y-=x*(a/b);
    return temp;
}
ll finv(ll a,ll m)// 求逆元
{
    ll x,y;
    ll g=exgcd(a,m,x,y);
    x=(x%m+m)%m;//
    return x;
}
int main()
{
    int n;
    int a,b;
    string s;
    while(cin>>n)
    {
        cin>>s;
        Hash[0]=1;
        for(int i=0;i<s.length();i++)
        Hash[i+1]=Hash[i]*(s[i]-28)%9973;
        while(n--)
        {
            cin>>a>>b;
            a=Hash[a-1];
            b=Hash[b];
            int temp=finv(a,9973);
            int ans=(b*temp)%9973;
            printf("%d\n",ans);
        }
    }
    return 0;
}