概念:素数(质数),因数只有1和它自身的自然数(>=2,0和1不是素数)。
1.朴素做法(略)
->改进版:因为因子是成对存在的(8的因子2 4/4 2),i循环到sqrt(num)即可。
2.埃式筛:
原理:素数的倍数不是素数,对于每个素数i(从小到大滚动),筛掉i*j(j=2,3..)
代码:
int prime[MAXN];
void solve(){
memset(prime,1,sizeof(prime));
prime[0]=prime[1]=0;
for(int i=2;i<=MAXN;i++)
if(prime[i])
for(int j=2;i*j<=MAXN;j++)//j从2开始
prime[i*j]=0;
}
//小优化: j可以从i*i开始(感觉无所谓)
即 for(int j=i*i;j<=MAXN;j+=i)
prime[j]=0;
复杂度: O(nloglogn) (不会证明。。)
问题:对于一个合数,有可能会被筛多次。比如:10=2*5,5*2,由此引出优化版--欧拉筛。
3.欧拉筛(线性筛)
int n,cnt,prime[MAXN];//prime[]用来存放素数
bool isprime[MAXN];
void PrimeSolve(){
memset(isprime,1,sizeof(isprime));//初始全设为1
isprime[0]=isprime[1]=0;
for(int i=2;i<=n;i++){
if(isprime[i])//遇到素数,存进去
prime[++cnt]=i;
for(int j=1;j<=cnt&&i*prime[j]<=n;j++){//内层循环j遍历prime[]
isprime[i*prime[j]]=0;
if(i%prime[j]==0)
break;//***使得一个合数不会被重复筛(关键)***
}
}
}
原理:在筛素数时,保证任一合数被这样筛去:
“最小质因数 × 最大因数 = 这个合数” (排除因子是1和本身的情况)
(另一个尽可能小的质数去筛掉一个合数,并且这个合数不会被再次筛一遍)
if(i%prime[j]==0) break;
当prime[j]是i的一个因子时,假设 某合数=i*prime[j],那么可以写成i= k*prime[j],此时如果我们不跳出循环,继续消掉 i*prime[j+1] (也就是k * prime[j] *prime[j+1] ),那么很显然,当i=k*prime[j+1]时,会重复消去k*prime[j+1]*prime[j],并且很显然prime[j]是更小的素因子。
例如:i= 6,prime[j]=3,正常消去 18,prime[j+1]=5,i= 3 * 2(k=2),那么我们这里会消去6*5=30,但是当i=10时,3*10=30,并且很显然3是最小的素因子,因此我们需要跳出循环。
->说实话感觉理解不够透彻,还是需要手推一下辅助理解。
时间复杂度:O(n),不然为什么叫线性筛awa。