LightOj 1370欧拉函数快速筛法
2016-08-16 16:42
435 查看
欧拉函数
首先介绍下什么是欧拉函数吧,欧拉函数phi(x)代表小于等于x的数中和x互质的数的个数(小于显然只对1成立), 比如说小于等于9的数中与9互质的有1,2,4,5,7,8,则phi(9)=6.求phi(x)得公式由欧拉给出(神一般的男人,几何学,数论,统计学,物理学,统计学,据说连金融学都有以欧拉命名的公式,膜拜…….)欧拉指出,以phi(x)表示小于等于x的数中与x互质的数的个数可以这样得到:
phi(x) = x * (1 - 1/p1) * (1-1/p2) * (1 - 1 / p3)…..*(1 - 1/pn)
其中p1,p2,p3….pn是x的所有质因数,每个质因数只使用一次,使用上述公式,则phi(9)=9*(1-1/3)=6.
如果只是单纯求解某一个数的欧拉值,这是很简单的,对这个数分解质因数求解即可,如果x不是很大,甚至你直接用gcd来求也行,但是,如果我们要求解1-100万的所有数的欧拉函数值,上述的做法便太耗时间了,做法并不难,效仿素数筛法,我们有对应欧拉筛法.
欧拉筛法的原理如下,我们观察phi(x)的求解式子,首先另phi(x)=x,然后仿照素数筛法,如果x能被2,3,5,7….这些素数筛到,则执行phi(x)=phi(x)*(1-1/prime)
代码如下,一看就明白
//欧拉筛法,快速筛出1-100万的欧拉函数值 const int maxn = 1000000; int phi[maxn+1]; bool isPrime[maxn+1]; void Eular() { for(int i=1;i<=maxn;i++) phi[i]=i; memset(isPrime,true,sizeof(isPrime)); isPrime[0]=isPrime[1]=false; phi[1]=0; for(int i=2;i<=maxn;i++) { if(isPrime[i]) { for(int j=i;j<=maxn;j+=i) { isPrime[j]=false; phi[j] -= phi[j]/i; } } } }
好吧,我也是今天才学会这个筛法,素数筛法还有一种更快的接近O(n)的筛法,叫做线性素数筛法,对应的,欧拉函数筛法也有一个线性的,有兴趣自己google吧.
今天写的一道用上此知识点的题目:
LightOj 1370
题目链接
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; const int maxn = 1100000; int phi[maxn+1]; bool isPrime[maxn+1]; int mp[maxn+1]; //欧拉筛法,快速筛出1-100万的欧拉函数值 void Eular() { for(int i=1;i<=maxn;i++) phi[i]=i; memset(isPrime,true,sizeof(isPrime)); isPrime[0]=isPrime[1]=false; phi[1]=0; for(int i=2;i<=maxn;i++) { if(isPrime[i]) { for(int j=i;j<=maxn;j+=i) { isPrime[j]=false; phi[j] -= phi[j]/i; } } } } void getMap() { memset(mp,0,sizeof(mp)); for(int i=1;i<=maxn;i++) // the length of bamboo { for(int j=phi[i];j>=0 && mp[j]==0;j--) mp[j]=i; } } int main() { int T,cas=1; scanf("%d" ,&T); Eular(); getMap(); while(T--) { int n,x; scanf("%d",&n); long long sum=0; for(int i=0;i<n;i++) { scanf("%d" ,&x); sum+=mp[x]; } printf("Case %d: %lld Xukha\n",cas++,sum); } return 0; }
相关文章推荐
- LightOJ - 1370 欧拉函数的应用
- 【欧拉函数打表】LightOJ - 1370 Bi-shoe and Phi-shoe
- LightOJ 1370(欧拉函数)
- 【LightOJ】1370 - Bi-shoe and Phi-shoe(欧拉函数,素数打表)
- LightOJ - 1370 Bi-shoe and Phi-shoe【欧拉函数的性质】
- LightOJ1370 - Bi-shoe and Phi-shoe(欧拉函数+打表)
- 【欧拉函数】Bi-shoe and Phi-shoe【LightOJ1370】------长篇阅读专场D题
- 数论快速入门(同余、扩展欧几里德、中国剩余定理、大素数测定和整数分解、素数三种筛法、欧拉函数以及各种模板)
- lightoj1370——Bi-shoe and Phi-shoe(欧拉函数应用)
- LightOJ - 1370 - Bi-shoe and Phi-shoe【欧拉函数预处理】
- lightoj 1370 - Bi-shoe and Phi-shoe(欧拉函数)
- lightoj1370欧拉函数/素数筛
- LightOJ 1370 Bi-shoe and Phi-shoe(欧拉函数)
- LightOJ - 1370【欧拉筛法*经典】
- lightoj 1370 Bi-shoe and Phi-shoe 【欧拉函数应用】
- lightoj-1370 Bi-shoe and Phi-shoe( 欧拉筛 求 欧拉函数)
- 质数快速筛法;欧拉函数线性打表;因子个数,因子和打表
- LightOJ 1370 Bi-shoe and Phi-shoe【欧拉函数 && 质数】
- LightOJ 1370 Bi-shoe and Phi-shoe 欧拉函数
- LightOJ 1370-Bi-shoe and Phi-shoe(欧拉函数)