线性筛的代码和普通筛法真的很像, 唯一不同的就是线性筛是用一个 自然数 乘以 若干 质数 去筛, 而普通筛法是用一个质数乘以 若干个自然数去筛, 那它们的不同之处在哪里呢?
普通的筛法显然是有优化的余地的, 因为一个合数会被筛掉多遍(应该是它的质因数个数遍), 每找到一个素数时都要从它一直循环到N, 而且这是不可避免的, 下面看一下欧拉筛法的代码:
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#define N 100
using namespace std;
bool issu[N + 5];
int prim[N], pt;
int main()
{
for(int i = 2; i <= N; i ++){
if(!issu[i])prim[++ pt] = i;
for(int j = 1; j <= pt && i * prim[j] <= N; j ++){
issu[i * prim[j]] = 1;
if(i % prim[j] == 0)break;
}
}
for(int i = 1; i <= pt; i ++)printf("%d ", prim[i]);
system("pause");
return 0;
}
计算Phi:
bool com[MAXN];
int primes, prime[MAXN], phi[MAXN];
phi[1] = 1;
for (int i = 2; i <= n; ++i)
{
if (!com[i])
{
prime[primes++] = i;
phi[i] = i-1;
}
for (int j = 0; j < primes && i*prime[j] <= n; ++j)
{
com[i*prime[j]] = true;
if (i % prime[j])
phi[i*prime[j]] = phi[i]*(prime[j]-1);
else
{ phi[i*prime[j]] = phi[i]*prime[j]; break; }
}
}
贾志鹏的证明:
设p 为一个合数(A)的最小的质因数, p'为这个合数的另一个质因数,
k = A / p; k' = A / p'
易证 k' 是p的倍数。
p 可以筛掉这个 合数, 但是p' 筛不掉。
考虑循环到k' 时, 这时j循环到p的时候就会被break掉了, 所以就走不到p' 所以A 没有被p'筛掉。
由于这样做一个合数只有可能被其最小素数筛掉, 所以时间是严格O(n)的。
我上午学的东西下午比赛就用到了这是什么心态。。不过因为上午懒了, 只把那个PDF看了一半, 害的下午自己YY了好久怎么筛出 phi , 呜呜呜~~~~(>_<)~~~~
积性函数:
考虑一个定义域为ℕ+的函数f,对于任意两个互质的正整数a,b,均满足 f(ab)=f(a) f(b)则函数f被称为积性函数。
下面看欧拉函数的性质: