筛法求欧拉函数

问题

https://www.acwing.com/problem/content/description/876/

分析

在前面求一个数所有的质数时，使用到了线性筛质数，时间复杂度 $O(n)=n$ ，通过线性筛质数将质数筛出来的同时，计算每个数的欧拉函数。

质数的欧拉函数 $\phi (n)=n\ast (1-\frac{1}{n})=n-1$

在线性筛质数的过程中，$i\%primes[j]==0$ 时，意味着 primes[j] 是 i 的最大质因子，在计算 $\phi (i\ast primes[j])$ 时只需要将 $primes[j]$ 乘到 $\phi (i)$ 即可，因为 $(1-\frac{1}{primes[j]})$ 在计算 $\phi (i)$ 时已经计算过了，只需要补上一个 $primes[j]$ 的乘积即可

代码

#include<iostream>

using namespace std ;

const int N = 1e6 + 10 ;

typedef long long LL ;

int n, cnt ;
int primes[N], st[N], phi[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] = 1 ;
			if(i % primes[j] == 0)
			{
				phi[primes[j] * i] = primes[j] * phi[i] ;
				break ;
			}
			// phi[primes[j] * i] = phi[primes[j]] * phi[i] 因为 primes[j] 是质数，所以可化为 (primes[j] - 1) * phi[i] 
			phi[primes[j] * i] = (primes[j] - 1) * phi[i] ;
		}
	}
	
	LL res = 0 ;
	for(int i : phi) res += i ;
	
	return res ;
}

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