欧拉函数

对正整数n，欧拉函数是小于n的正整数中与n互质的数的数目（其中φ(1)=1）例如φ(8)=4，因为1,3,5,7均和8互质
φ（n） = n * (1-1/p1) * (1-1/p2)…(1-1/pn)。其中p1,p2…pn是n的质因数
对公式的理解:
先理解一下这个定理：任意一个大于1的正整数，都可以写成一系列质数的积。
对于任意的一个pi，因为pi与n互质且是n的因子，那么1~n中有1/pi的部分是pi的倍数。对于另一个与n互质的pj，有1/pj的部分是pj的倍数。那么就有（1-pi）* (1-pj)的部分既不是pi的倍数，也不是pj的倍数。
所以对于1~n中所有与n互质的因子，就有(1-1/p1) * (1-1/p2)…(1-1/pn)的部分都不是他们的倍数。
那么因为对于任意一个大于1的正整数，都可以写成一系列质数的积，而这些数都不是他们的倍数，那么这些数与n一定没有除1以外的任何一个公因数，即就是与n互质。
代码实现

//任意一个大于1的正整数，都可以写成一系列质数的积。 
int eular( int n )
{
	int res = n , a = n;
	for( int i = 2 ; i * i <= n ; i++ )
	{
		if( a%i == 0 )
		{
			res -= res / i;
			while( a%i == 0 ) a /= i;   //把质因数全部除掉 
		}
	}
	if( a > 1 ) res -= res / a;    //n为质数的情况
	return res;
}

复杂度O(sqrt(n))

求1~n的欧拉函数的值
如果要求1~n的欧拉函数值，我们可以求每个的欧拉函数，复杂度为O(n^(3/2))
但其实我们可以利用欧拉筛来计算，复杂度为O(n)

/*
特性 :
1.若a为质数,phi[a]=a-1;
2.若a为质数,b mod a=0,phi[a*b]=phi[b]*a
3.若a,b互质,phi[a*b]=phi[a]*phi[b]
*/

#include <iostream>
#include <vector>
using namespace std;

bool check[10000005];
vector<int> prime;
int phi[1000005];

int n;

void euler()
{
	check[1] = true;
	phi[1] = 1;
	for (int i = 2; i <= n; i++)
	{
		if( !check[i] ) 
		{
			prime.push_back(i);
			phi[i] = i - 1;   //定理1 
		}
		for (int j = 0; j < prime.size(); j++)
		{
			if( i * prime[j] > n ) break;
			check[i*prime[j]] = true;
			if( i % prime[j] == 0 ) 
			{
				phi[i*prime[j]] = phi[i] * prime[j];     //定理2 
				break;
			}else phi[i*prime[j]] = phi[i] * phi[prime[j]];  //定理3 
		} 
	}
}

int main()
{
	cin >> n;
	euler(); 
	for (int i = 1; i <= n; i++)
	{
		cout << phi[i] << " ";
	}
	cout << endl;
	return 0;
}

欧拉函数的一个小应用

/*
求n以内与n互质的元素和，这个和为phi[n]*n/2
证明：与n互质的集合{x1,x2,...xi}，
由gcd(n,xi) = gcd(n-xi,n) = 1可得，
与n互质的一个新集合{n-x1,n-x2,...n-xi}。
将这两个集合相加除二即可得与n互质的元素和，即为phi[n]*n/2;
*/ 
#include <iostream>
using namespace std;

typedef long long ll;

ll mod = 1e9 + 7;

ll eular( int n )
{
	ll res = n , a = n;
	for( int i = 2 ; i * i <= n ; i++ )
	{
		if( a%i == 0 )
		{
			res -= res / i;
			while( a%i == 0 ) a /= i;   //把质因数全部除掉 
		}
	}
	if( a > 1 ) res -= res / a;    //n为质数的情况
	return res;
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	ll n;
	while( cin >> n && n != 0 )
	{
		ll ans = eular(n) * n / 2;
		ans %= mod;
		ll t = n * (n-1) / 2;
		t -= ans;
		ans = t % mod;
		cout << ans << '\n';
	}
	return 0;
}