筛法求欧拉函数（C++）

题目

给定一个正整数 n，求 1∼n 中每个数的欧拉函数之和。

输入格式

共一行，包含一个整数 n。

输出格式

共一行，包含一个整数，表示 1∼n 中每个数的欧拉函数之和。

数据范围

1≤n≤106

输入样例：

6

输出样例：

12

思路

 

i mod pj != 0 时如下图

 

 

代码

#include<iostream>

using namespace std;

const int N = 1e6 + 10;

typedef long long LL;

int primes[N], cnt;
int phi[N];//这个数的欧拉函数 
bool st[N];

LL get_eulers(int n)
{
	phi[1] = 1;
	
	for(int i = 2; i <= n; i ++)
	{
		if(!st[i])
		{
			primes[cnt ++] = i;
			phi[i] = i - 1;//质数的欧拉函数 
		}
		
		for(int j = 0; primes[j] <= n / i; j ++)
		{
			st[primes[j] * i] = true;
			
			if(i % primes[j] == 0)
			{
				phi[primes[j] * i] = phi[i] * primes[j];
				break;
			}
			phi[primes[j] * i] = phi[i] * (primes[j] - 1);
		}
	}
	
	LL res = 0;
	for(int i = 1; i <= n; i ++) res += phi[i];
	return res;
}

int main()
{
	int n;
	cin >> n;
	
	cout << get_eulers(n) << endl;
	
	return 0;
}