博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
欧拉定理,欧拉降幂
阅读量:5290 次
发布时间:2019-06-14

本文共 3249 字,大约阅读时间需要 10 分钟。

在数论中,
欧拉定理,(也称费马-
欧拉定理)是一个关于同余的性质定理。了解欧拉定理之前先来看一下费马小定理:

a是不能被质数p整除的正整数,则有a^(p-1) ≡ 1 (mod p)

欧拉给出了推广形式

若n,a为正整数且互质,则,其中φ(n)表示小于等于n的数中与n互质的数的数目。可以看出费马小定理是欧拉定理的一种特殊情况。

证明(好文章):https://zhuanlan.zhihu.com/p/24902174

欧拉函数板子:

#include
using namespace std;typedef long long LL;const int maxn = 1e6+5;LL n;LL get_Euler(LL x){ LL res = x; ///初始值 for(LL i = 2LL; i * i <= x; ++i) { if(x % i == 0) { res = res / i * (i - 1); ///先除后乘,避免数据过大 while(x % i == 0) x /= i; } } if(x > 1LL) res = res / x * (x - 1); ///若x大于1,则剩下的x必为素因子 return res;}int main(){ while(cin >> n) { cout << get_Euler(n) << endl; ///求n的互质数的个数 // cout << n * get_Euler(n) / 2 << endl; ///求n的所有互质数之和 } return 0;} //预处理打表写法:/* #include
#include
#include
#include
#include
#include
using namespace std;typedef long long LL;const int maxn = 1e6+5;int n, phi[maxn];void phi_table() { phi[0] = 0, phi[1] = 1; ///1的欧拉函数值为1,唯一与1互质的数 for(int i = 2; i < maxn; ++i) phi[i] = i; //先初始化为其本身 for(int i = 2; i < maxn; ++i) { if(phi[i] == i) { ///如果欧拉函数值仍为其本身,说明i为素数 for(int j = i; j < maxn; j += i) ///把i的欧拉函数值改变,同时也把能被素因子i整除的数的欧拉函数值改变 phi[j] = phi[j] / i * (i - 1); } }}int main(){ phi_table(); int m; while(cin >> n>>m) { int sum=0; //cout << phi[n] << endl; for(int i=n;i<=m;i++) sum+=phi[i]; cout<
<

 

题:https://vjudge.net/problem/FZU-1759

#include
#include
#include
#include
using namespace std;const int M=1e6+6;char s[M];typedef long long ll;ll ph(ll x){ ll res = x; ///初始值 for(ll i = 2LL; i * i <= x; ++i) { if(x % i == 0) { res = res / i * (i - 1); ///先除后乘,避免数据过大 while(x % i == 0) x /= i; } } if(x > 1ll) res = res / x * (x - 1); ///若x大于1,则剩下的x必为素因子 return res;}ll ksm(ll a,ll b,ll mod){ ll t=1; while(b){ if(b&1) t=(t*a)%mod; b>>=1; a=(a*a)%mod; } return t;}int main(){ ll a,m; while(~scanf("%I64d %s %I64d",&a,s,&m)){ ll ans=0; int len=strlen(s); ll p=ph(m); for(int i=0;i
View Code

经典题:https://www.lydsy.com/JudgeOnline/problem.php?id=3884

指数是无穷的,但是模数是有限的,从不断减小p去考虑。

因为次数是无穷的,所以B肯定>=ph(),所以采用第二条式子去递归

#include
#include
#include
#include
using namespace std;typedef long long ll;ll ph(ll x){ ll ans=x; for(ll i=2ll;i*i<=x;i++){ if(x%i==0){ ans=ans/i*(i-1); while(x%i==0) x/=i; } } if(x>1ll) ans=ans/x*(x-1); return ans;}ll ksm(ll a,ll b,ll mod){ ll t=1; while(b){ if(b&1) t=(t*a)%mod; b>>=1; a=(a*a)%mod; } return t;}ll dfs(ll p){ if(p==1) return 0; ll k=ph(p); return ksm(2,dfs(k)+k,p);}int main(){ int t; scanf("%d",&t); while(t--){ ll p; scanf("%lld",&p); printf("%lld\n",dfs(p)); } return 0;}
View Code

 

转载于:https://www.cnblogs.com/starve/p/11443764.html

你可能感兴趣的文章
linux下编译openjdk8
查看>>
【python】--迭代器生成器装饰器
查看>>
Pow(x, n)
查看>>
安卓当中的线程和每秒刷一次
查看>>
MySQL Proxy
查看>>
关于Vue的组件的通用性问题
查看>>
随机颜色值
查看>>
每日一库:Modernizr.js,es5-shim.js,es5-safe.js
查看>>
目录相关的操作
查看>>
解决虚拟机vmware安装64位系统“此主机支持 Intel VT-x,但 Intel VT-x 处于禁用状态”的问题...
查看>>
C++----练习--引用头文件
查看>>
11.基本包装类型
查看>>
ajax连接服务器框架
查看>>
wpf样式绑定 行为绑定 事件关联 路由事件实例
查看>>
利用maven管理项目之POM文件配置
查看>>
用HttpCombiner来减少js和css的请问次数
查看>>
FUSE-用户空间文件系统
查看>>
将tiff文件转化为jpg文件并保存
查看>>
ubuntu 16.04 开机脚本
查看>>
 VS2012 C#调用C++ dll
查看>>