求欧拉函数
2015-08-22 22:22
369 查看
在数论,对正整数n,欧拉函数是少于或等于n的数中与n互质的数的数目。此函数以其首名研究者欧拉命名,它又称为Euler's totient function、φ函数、欧拉商数等。 例如φ(8)=4,因为1,3,5,7均和8互质。 从欧拉函数引伸出来在环论方面的事实和拉格朗日定理构成了欧拉定理的证明。
φ函数的值 通式:φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数,x是不为0的整数。φ(1)=1(唯一和1互质的数(小于等于1)就是1本身)。 (注意:每种质因数只一个。比如12=2*2*3那么φ(12)=12*(1-1/2)*(1-1/3)=4若n是质数p的k次幂,φ(n)=p^k-p^(k-1)=(p-1)p^(k-1),因为除了p的倍数外,其他数都跟n互质。设n为正整数,以 φ(n)表示不超过n且与n互素的正整数的个数,称为n的欧拉函数值,这里函数φ:N→N,n→φ(n)称为欧拉函数。欧拉函数是积性函数——若m,n互质,φ(mn)=φ(m)φ(n)。特殊性质:当n为奇数时,φ(2n)=φ(n), 证明与上述类似。
求1~N之间的欧拉函数值O(nloglogn)
#include<bits/stdc++.h> using namespace std; typedef long long LL; LL N; LL Phi[5000000]; inline void get_Phi(){ for(int i=2;i<=N;i++) Phi[i]=0; Phi[1]=1; for(int i=2;i<=N;i++){ if(!Phi[i]){ for(int j=i;j<=N;j+=i){ if(!Phi[j]) Phi[j]=j; Phi[j]=Phi[j]/i*(i-1); } } } } int main(){ scanf("%d",&N); get_Phi(); for(int i=1;i<=N;i++){ cout<<Phi[i]<<" "; if(i%10==0){ cout<<endl; } } return 0; }
但是这个复杂度还不够优秀,网上有O(n)的筛法:
/*线性筛O(n)时间复杂度内筛出maxn内欧拉函数值*/ int m[maxn],phi[maxn],p[maxn],pt;//m[i]是i的最小素因数,p是素数,pt是素数个数 int make(){ phi[1]=1; int N=maxn; int k; phi[1]=1; for(int i=2;i<N;i++){ if(!m[i])//i是素数 p[pt++]=m[i]=i,phi[i]=i-1; for(int j=0;j<pt&&(k=p[j]*i)<N;j++) { m[k]=p[j]; if(m[i]==p[j])//为了保证以后的数不被再筛,要break { phi[k]=phi[i]*p[j]; /*这里的phi[k]与phi[i]后面的∏(p[i]-1)/p[i]都一样(m[i]==p[j]) 只差一个p[j],就可以保证∏(p[i]-1)/p[i]前面也一样了*/ break; } else phi[k]=phi[i]*(p[j]-1);//积性函数性质,f(i*k)=f(i)*f(k) } } }
还可以单独求欧拉函数的O(sqrt(n)):
int euler_Phi(int n){ int m=(int)sqrt(n+0.5); int ans=n; for(int i=2;i<=m;i++){ if(n%i==0){ ans=ans/i*(i-1); while(n%i==0) n/=i; } } if(n>1) ans=ans/n*(n-1); return ans; }
相关文章推荐
- DLX by java (example : hdu 2295)
- 关于数组实现队列的一些设计
- hdu1465 不容易系列之一
- 多校连萌15-8-12#A
- Unity3D开发——安卓版的虚拟按键实现
- tcp和udp
- 浪潮天梭k1主机并非全是基于安腾cpu的
- TSP+Floyd BestCoder Round #52 (div.2) 1002 Victor and Machine
- 引用和指针的区别和联系
- Objective-C 【继承、变量修饰符(私有变量/方法)、description方法】
- Delphi2010中DataSnap技术网摘
- Objective-C 【继承、变量修饰符(私有变量/方法)、description方法】
- LeetCode 29 Divide Two Integers 除法的原理
- R语言_非参数检验
- IOS弹出视图preferredContentSize
- ACM学习历程—HDU5418 Victor and World(动态规划 && 状压)
- 启动tomcat时出现The specified JRE installation does not exist 如何解决?
- 面向对象设计模拟简单ATM系统
- CCActionGird源码分析
- 第91讲:Akka第一个案例动手实战架构设计