线性筛法

贾志鹏写的线性筛 

线性筛的代码和普通筛法真的很像， 唯一不同的就是线性筛是用一个 自然数 乘以 若干 质数 去筛， 而普通筛法是用一个质数乘以 若干个自然数去筛， 那它们的不同之处在哪里呢？

普通的筛法显然是有优化的余地的， 因为一个合数会被筛掉多遍（应该是它的质因数个数遍）， 每找到一个素数时都要从它一直循环到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被称为积性函数。

下面看欧拉函数的性质：