拓展欧拉定理及其应用

Posted aininot260

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了拓展欧拉定理及其应用相关的知识,希望对你有一定的参考价值。

回顾一下上节的欧拉定理:

其化简的形式为:

a^φ(m)≡1(mod m)

就有:

a^x≡a^(x%φ(m))≡a^(x%φ(m)+φ(m))(mod m)

看一道题:

Given A,B,C, You should quickly calculate the result of A^B mod C. (1<=A,C<=1000000000,1<=B<=10^1000000)

这个B是非常大的一个指数,那么幂就更大了,需要降幂

若gcd(a,c)=1,那么使用欧拉定理即可:a^b≡a^(b%φ(m))(mod m)

若gcd(a,m)>1,且b>φ(m),则有“求幂大法”

a^b≡a^(b%φ(m)+φ(m))(mod m)

这就是刚才的那个式子,(当b<=φ(m)时直接用快速幂即可)

 1 //A^B %C=A^( B%phi(C)+phi(C) ) %C
 2 #include<cstdio>
 3 using namespace std;
 4 int c;
 5 long long a,nb;
 6 char tb[1000005];
 7 int phi(int x)
 8 {
 9     int ans=x;
10     for(int i=2;i*i<=x;i++)
11     {
12         if(x%i==0)
13         {
14             ans=(ans/i)*(i-1);
15             while(x%i==0) x=x/i;
16         }
17     }
18     if(x!=1) ans=(ans/x)*(x-1);
19     return ans;
20 }
21 long long qpow(long long m,long long n,long long k)
22 {
23     long long ans=1;
24     while(n!=0)
25     {
26         if(n&1) ans=(ans*m)%k;
27         n=(n>>1);
28         m=(m*m)%k;
29     }
30     return ans;
31 }
32 int main()
33 {
34     while(scanf("%I64d%s%d",&a,tb,&c)!=EOF)
35     {
36         int PHI=phi(c);
37         long long res=0;
38         for(int i=0;tb[i];i++)
39         {
40             res=(res*10+tb[i]-0);
41             if(res>c) break;
42         }
43         if(res<=PHI) printf("%I64d
",qpow(a,res,c));
44         else
45         {
46             res=0;
47             for(int i=0;tb[i];i++) res=(res*10+tb[i]-0)%PHI;
48             printf("%I64d
",qpow(a,res+PHI,c));
49         }
50     }
51     return 0;
52 }

拓展欧拉定理,肯定试用范围最广喽

以上是关于拓展欧拉定理及其应用的主要内容,如果未能解决你的问题,请参考以下文章

BSGS(Baby-Step-Giant-Step)算法及其应用

欧拉定理及其应用

初等数论及其应用——欧拉函数

Luogu4139 上帝与集合的正确用法 拓展欧拉定理

暑假集训Day1 B(拓展欧拉定理)

[六省联考2017]相逢是问候(线段树+拓展欧拉定理)