【数论】线性筛素数，线性筛欧拉函数，求前N个数的约数个数

先来最基本的线性筛素数，以后的算法其实都是基于这个最基本的算法：

#include<stdio.h> #include<string.h> #define M 10000000 int prime[M/3]; bool flag[M]; void get_prime() { int i,j,k; memset(flag,false,sizeof(flag)); k=0; for(i=2;i<M;i++){ if(!flag[i]) prime[k++]=i; for(j=0;j<k&&i*prime[j]<M;j++){ flag[i*prime[j]]=true; if(i%prime[j]==0) break; } } } int main() {}

利用了每个合数必有一个最小素因子，每个合数仅被它的最小素因子筛去正好一次，所以是线性时间。
代码中体现在： if(i%prime[j]==0) break;
----------------------------------------------------------------------- 我是低调的分割线 ------------------------------------------------------------------------------------------
然后可以利用这种线性筛法求欧拉函数，需要用到以下几个性质：
//(1) 若(N%a==0 && (N/a)%a==0) 则有:E(N)=E(N/a)*a;
//(2) 若(N%a==0 && (N/a