码迷,mamicode.com
首页 > 其他好文 > 详细

费马小定理 求乘法逆元

时间:2020-01-17 21:04:37      阅读:114      评论:0      收藏:0      [点我收藏+]

标签:快速幂   font   putc   pow   div   时间   复杂度   return   char   

//P3811 【模板】乘法逆元
#include<bits/stdc++.h>
using namespace std;
inline void write(long long X)
{
    if(X<0) {X=~(X-1); putchar(-);}
    if(X>9) write(X/10);
    putchar(X%10+0);
}
long long qpow(long long n,long long m,long long mod)
{
       long long ans=1;
    while(m)
    {
        if(m&1) ans=ans*n%mod;
        n=n*n%mod;
        m>>=1;
    }
    return ans%mod;
}
int main()
{
    long long n,p;
    scanf("%lld%lld",&n,&p);
    for(long long i=1;i<=n;i++)
        write(qpow(i,p-2,p)),printf("\n");
    return 0;
}

这一道题是洛谷P3811 【模板】乘法逆元

当然这一道题用费马小定理还是过不去的

不过可以证明这一做法的正确性

 

首先我们要保证题目给出的p是质数

所谓费马小定理就是

a^(p-1) ≡ 1 (mod p)

稍加化简就可以看出左边可以化为a*a^(p-2)

把那个单独的a挪到右边

所以a^(p-2)=a-1

也就是a^(p-2)就是a在模p意义下的乘法逆元

 

然而我们可以用快速幂来进行加速

这样一来算出一个乘法逆元的时间复杂度就是logn

 

但是一定要注意 必须p是负数!

 

上线性筛

 

//P3811 【模板】乘法逆元
#include<bits/stdc++.h>
using namespace std;
inline void write(long long X)
{
    if(X<0) {X=~(X-1); putchar(-);}
    if(X>9) write(X/10);
    putchar(X%10+0);
}
long long qpow(long long n,long long m,long long mod)
{
       long long ans=1;
    while(m)
    {
        if(m&1) ans=ans*n%mod;
        n=n*n%mod;
        m>>=1;
    }
    return ans%mod;
}
long long inv[3000005];
int main()
{
    long long n,p;
    scanf("%lld%lld",&n,&p);
//    for(long long i=1;i<=n;i++)
//        write(qpow(i,p-2,p)),printf("\n");
    inv[1]=1;write(1),printf("\n");
    for(int i=2;i<=n;i++)
    inv[i]=(p-p/i)*inv[p%i]%p,write(inv[i]),printf("\n");//线性筛 
    return 0;
}

 

还有一个需要注意的地方

就是那个取模运算是很耗时间的

一开始我把代码改成线性筛之后 

在输出 就是那个write的括号里又增加了一个%p

事实上是没有必要的

结果就因为这个多余的取模 y以至于我还是T一个点,,

所以千万不要多加太多%取模

费马小定理 求乘法逆元

标签:快速幂   font   putc   pow   div   时间   复杂度   return   char   

原文地址:https://www.cnblogs.com/akioi/p/12207274.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!