资格赛 A题

画图说话:

前缀积推导

分析:9937 为质数,费马小定理+逆元,除法变乘法,快速幂取余

费马小定理

公式推导:(b/a)mod p=> (b*a的逆元)mod p 根据费马小定理,a^(p-1) =1 mod p <=> a^(p-2) a= 1 mod p <=> a的逆元为 a^(p-2),所以(b/a)mod p <=>(ba^(p-2)) mod p

用到的板子:快速幂取余、前缀积

快速幂取余

// x是底数,n是幂数,mod是取余数
LL mod_pow(LL x,LL n,LL mod)
{
    LL res = 1;
    while(n>0)
    {
        if(n & 1)
        {
            res = res * x % mod;
        }
        x = x * x % mod;
        n >>= 1;
    }
    return res;
}

前缀积

H[0] = 1;
for(int i = 1 ;i<=len;i++)
    {
          H[i]=H[i-1]*(Hstr[i-1]-28)%mods;
    }

AC代码

#include <cstdio>
#include <cstring>
#include <algorithm>
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;
    while(n>0)
    {
        if(n & 1)
        {
            res = res * x % mod;
        }
        x = x * x % mod;
        n >> 1;
    }
    return res;
}

int main ()
{
    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);
            if(l>r)
            {
                swap(l,r);
            }
            printf("%I64D\n",(LL)H[r]*mod_pow(H[l-1],mods-2,mods)%mods);
        }
    }
}
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 归去来兮。 1.1 说明 本篇为《挑战程序设计竞赛(第2版)》[http://www.ituring.com.cn...
    尤汐Yogy阅读 14,507评论 0 160
  • MD5的全称是Message-Digest Algorithm 5,在90年代初由MIT的计算机科学实验室和RSA...
    没能唱给你的歌曲阅读 999评论 2 6
  • 欧几里得算法:递归版本: 迭代版本: 扩展欧几里德算法:基本算法:对于不完全为 0 的非负整数 a,b,gcd(a...
    Gitfan阅读 690评论 0 0
  • 有人说: 母校就是你一天骂它八遍,却不许别人骂一遍的地方。 因为那里: 有我们在课堂上的争论 有我们在球场上的奔跑...
    田觅阅读 712评论 11 11
  • 长大成熟了就不会再纠结了。原因有很多: 1.看的多了,知道然并卵。 2.无所谓了,皮糙肉厚的。 3.定型了,没那么...
    饶实阅读 331评论 0 1