前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >逆元模板

逆元模板

作者头像
attack
发布2018-04-13 15:37:11
6060
发布2018-04-13 15:37:11
举报

对于(a/b)%m==?

1.当m是素数的时候,根据费马小定理,直接输出b^(n-2)即可

2.否则,扩展欧几里得exgcd(b,m,x,y)

代码语言:javascript
复制
 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<cmath>
 5 using namespace std;
 6 int a,b,m;
 7 int x,y;
 8 int exgcd(int a,int b,int &x,int &y)
 9 {
10     if(b==0)
11     {
12         x=1;
13         y=0;
14         return a;
15     }
16     int r=exgcd(b,a%b,x,y);
17     int tmp=x;
18     x=y;
19     y=tmp-a/b*y;
20     return r;
21 }
22 int fastpow(int a,int p)
23 {
24     int base=a;int ans=1;
25     while(p!=0)
26     {
27         if(p%2==1)ans=ans*base;
28         base=base*base;
29         p=p/2;
30     }
31     return ans;
32 }
33 int main()
34 {
35     scanf("%d%d%d",&a,&b,&m);
36     for(int i=1;i<=sqrt(m);i++)
37     {
38         if(m%i==0)
39         {
40             int ans=exgcd(b,m,x,y);
41             printf("%d",(a*ans)%m);
42             return 0;    
43         }
44     }
45     printf("%d",fastpow(b,m-2));
46 }
本文参与 腾讯云自媒体分享计划,分享自作者个人站点/博客。
原始发表:2017-05-16 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档