数论素数筛法以及用欧拉筛求欧拉函数
2018-02-01 17:45
211 查看
http://blog.csdn.net/nk_test/article/details/46242401来自一位大神博客
素数筛
最经典的是最经典的埃拉特斯特尼筛法。时间复杂度为O(n
loglog n)
最好的筛法——欧拉筛线性O(n) 因为不会重复筛
然后利用每个合数必有一个最小质因子,每个合数仅被它的最小质因子筛去正好一次,所以是线性的
在代码上体现为if(i%prime[j]==0) break;
prime数组中的素数是单调递增的,当i 能整除prime[j],那么i*prime[j+1]这个合数肯定被prime[j] 乘以某个数筛掉
因为i中含有prime[j], prime[j] 比 prime[j+1] 小。接下去的素数同理。所以不用筛下去了。
在满足i%prme[j]==0这个条件之前以及第一次满足改条件时,prime[j]必定是prime[j]*i的最小因子。
如果还不是很理解,可以手动模拟一下。
欧拉函数:在数论中 ,对于正整数n ,欧拉函数是少于或等于n的数中与n互质的数的数目
现在还理解不了这段代码!!!!FUCK
using namespace std;
const int MAXN=3000001;
int prime[MAXN];//保存素数
bool vis[MAXN];//初始化
int phi[MAXN];//欧拉函数
void Prime(int n)
{
int cnt=0;
memset(vis,0,sizeof(vis));
for(int i=2;i<n;i++)
{
if(!vis[i])
{
prime[cnt++]=i;
phi[i]=i-1;// if p is prime,then phi[i]=i-1
}
for(int j=0;j<cnt&&i*prime[j]<n;j++)
{
__int64 k=i*prime[j];
vis[k]=1;
if(i%prime[j]==0)//关键
{
phi[k]=phi[i]*prime[j];
break;
}
else
phi[k]=phi[i]*(prime[j]-1);
}
}
}
素数筛
最经典的是最经典的埃拉特斯特尼筛法。时间复杂度为O(n
loglog n)
int ans[MAXN]; void Prime(int n) { int cnt=0; memset(prime,1,sizeof(prime)); prime[0]=prime[1]=0; for(int i=2;i<n;i++) { if(vis[i]) { ans[cnt++]=i;//保存素数 for(int j==i*i;j<n;j+=i) // i*i开始进行了稍微的优化 prime[j]=0; //不是素数 } } return ; }备注:从i*i而不从i*2 开始 是因为i*3,i*2早被 2 3筛过了
最好的筛法——欧拉筛线性O(n) 因为不会重复筛
const int MAXN=3000001; int prime[MAXN];//保存素数 bool vis[MAXN];//初始化 int prime(int n) { int cnt=0; memset(vis,0,sizeof(vis)); for(int i=2;i<n;i++) { if(!vis[i]) prime[cnt++]=i; for(int j=0;j<cnt&&i*prime[j]<n;j++) { vis[i*prime[j]]=1; //min_prime[i*prime[j]] = prime[j]; 求每个数的最小质因子 if(i%prime[j]==0)//关键 break; } } return cnt;//返回小于n的素数的个数 }理解:任何合数都可以表示成 一系列素数的积
然后利用每个合数必有一个最小质因子,每个合数仅被它的最小质因子筛去正好一次,所以是线性的
在代码上体现为if(i%prime[j]==0) break;
prime数组中的素数是单调递增的,当i 能整除prime[j],那么i*prime[j+1]这个合数肯定被prime[j] 乘以某个数筛掉
因为i中含有prime[j], prime[j] 比 prime[j+1] 小。接下去的素数同理。所以不用筛下去了。
在满足i%prme[j]==0这个条件之前以及第一次满足改条件时,prime[j]必定是prime[j]*i的最小因子。
如果还不是很理解,可以手动模拟一下。
欧拉函数:在数论中 ,对于正整数n ,欧拉函数是少于或等于n的数中与n互质的数的数目
现在还理解不了这段代码!!!!FUCK
using namespace std;
const int MAXN=3000001;
int prime[MAXN];//保存素数
bool vis[MAXN];//初始化
int phi[MAXN];//欧拉函数
void Prime(int n)
{
int cnt=0;
memset(vis,0,sizeof(vis));
for(int i=2;i<n;i++)
{
if(!vis[i])
{
prime[cnt++]=i;
phi[i]=i-1;// if p is prime,then phi[i]=i-1
}
for(int j=0;j<cnt&&i*prime[j]<n;j++)
{
__int64 k=i*prime[j];
vis[k]=1;
if(i%prime[j]==0)//关键
{
phi[k]=phi[i]*prime[j];
break;
}
else
phi[k]=phi[i]*(prime[j]-1);
}
}
}
相关文章推荐
- 数论快速入门(同余、扩展欧几里德、中国剩余定理、大素数测定和整数分解、素数三种筛法、欧拉函数以及各种模板)
- 【数论】gcd|扩展gcd|素数筛法|快速幂|欧拉函数(各种模板)
- 【hihocoder 1298】 数论五·欧拉函数
- HDU 2588 GCD [欧拉函数]【数论】
- poj——2689(数论之素数筛法)
- hdu - GCD -欧拉函数 & GCD原理 - 数论
- Toxophily-数论以及二分三分
- ACM第一课---数论之欧拉函数
- ACM常用数论之欧拉函数
- HDU 4002 Find the maximum(数论-欧拉函数)
- luogu 1865 数论 线性素数筛法
- 洛谷 2158 数论 打表 欧拉函数
- UVa 11440 Help Tomisu (数论欧拉函数)
- Bzoj 2705: [SDOI2012]Longge的问题 欧拉函数,数论
- 数论(素数以及计算方程)
- BZOJ 4805 欧拉函数求和 数论
- 51Nod-1136 欧拉函数【数论】
- hihocoder 1298 : 数论五·欧拉函数
- 欧拉函数概念与应用拓展 - 数论
- BZOJ_P2818 GCD(数论GCD+欧拉函数线性筛)