筛法欧拉函数

Number数论

欧拉函数PHI

筛法欧拉函数

欧拉函数
在数论，对正整数n，欧拉函数是小于或等于n的正整数中与n互质的数的数目（因此φ(1)=1）。此函数以其首名研究者欧拉命名(Euler’s totient function)，它又称为Euler’s totient function、φ函数、欧拉商数等。 例如φ(8)=4，因为1,3,5,7均和8互质。 从欧拉函数引伸出来在环论方面的事实和拉格朗日定理构成了欧拉定理的证明。
百度百科链接: link.
欧拉定理是用来阐述素数模下，指数同余的性质。
对于欧拉定理还有几个相关的结论：
如果n为素数，则φ( n ) = n-1
因为对于n来说，小于等于n的正整数有n个，但是因为n是素数，n与其本身并不互质，所以有n-1个数跟它互质。
若n为质数p的k次方 那么φ( n ）=(p-1)p^k
n = p^k(p是质数），所以一个数要想跟n互质的话应该不是p的倍数，p的倍数有p * p, 2 * p, 3 * p……p^n-1 * p,共有p^n-1个，小于等于n的正整数有p^k个，所以φ( n ）= (p-1)p^k.

埃氏筛

时间复杂度为O(n loglog n)，如果一个数的因数除了1和他本身没有其他因数的数是素数，所以如果他是除了一和它本身的一个数的倍数的话那么他就是合数，所以我们就从2开始，如果这个数没有被访问过说明之前他没有素数因子，所以我们就把在2到n中他的倍数给标记，如果一个数被访问过，说明这个数是合数，可以直接跳过了，一直找下去就好。
第二重循环我们在找的时候是从i*i开始的，是因为i * (比i小的数)，我们之前找过了，就像2 * i, 3 * i ……所以不会被重复筛查。

#include<iostream>
#include<cstring>
using namespace std;
int visited[1000001];
int main(){
	int n;
	cin >> n;
	memset(visited,0,sizeof(visited));
	for(int i=2;i<=n;i++){
		if(visited[i]) continue;
		else{
			for(int j = i*i;j<=n;j+=i){
				visited[j] = 1;
			}
		}
	}
	for(int i=2;i<=n;i++){
		if(visited[i]==0) cout << i << " ";
	}
	return 0;
}

但是其中还有一个需要我们注意的地方就是就像30，我们在筛查2的倍数的术后筛查过了，但是筛查5的时候同时又被筛查了一遍，为了解决这种合数被重复筛查的情况我们就有了时间复杂度是O(n)的欧拉线性筛法。

欧拉线性筛法

为此我们需要知道的一些结论：
①每一个合数都可以分解为一个质数与另外一个数相乘（合数是指在大于1的整数中除了能被1和本身整除外，还能被其他数（0除外）整除的数，所以只要是合数就能够被分解，依次递归分解下去）；
②if(i%prime[j]==0) break;
如果这个数能够被prime整除，说明这个数可以分解为prime[j]和另外一个大于等于prime[j]的数相乘，另外的这个数又可以分解成为若干个大于等于prime[j]的数相乘，所以如果继续下去会导致计算判断重复，故在此处停止，继续向下执行，由此保证了O(n)的复杂度。
代码如下：

#include<iostream>
#include<cstring>
using namespace std;
int vis[10001];//判断是否访问过
int prime[100001];//用以存储找到的质数
int cnt = 0;
int main(){
	int n;
	cin >> n;
	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;
			if(i%prime[j]==0) break;
		}
	}
	for(int i=0;i<cnt;i++){
		cout << prime[i] << " ";
	}
	return 0;
}