uva 11426 欧拉函数的运用

2019-08-01 10:45:02 浏览数 (1)

题意: 求 n 以内gcd( i , j ) i < j 的和的个数。

n 的欧拉函数值 表示 小于 n 与 n 互质的数的个数。

对于互质的两个数 a,b (a<b<=n), gcd(a,b) = 1 是显然的。

对 b 而言互质的数有 fai(b) = m个 ,那么b 的贡献就是 m。如果是 2*b呢? 那贡献就是2*m

因为将与 b 互质的m个数都乘2, 就得到 m 个 gcd(2*a,2*b) = 2的数

但2*b不仅仅有 2*m个,因为它也可以通过别的数相乘得到 。那么只需要求出每个数的欧拉函数,

再对 i * j 下标累加 j * fai(i) 个贡献就行,j从1倍一直到小于maxn。

求出的是 i 以内的gcd()和,前 i 项gcd() 和累加就行。

代码语言:javascript复制
//uva 11426 欧拉函数的运用 
#include <iostream>
#include <algorithm>
#include <cstring>
#define ll long long
using namespace std;
const int maxn = 4000051;
bool is_prime[maxn];
int prime[maxn],phi[maxn];
ll a[maxn];
void getphi(){
	memset(is_prime,0,sizeof(is_prime));
	phi[1] = 1;
	for(int i=2;i<maxn;i  ){
		if(!is_prime[i]){
			prime[  prime[0]] = i;
			phi[i] = i-1;
		}
		for(int j=1;j<=prime[0]&&i*prime[j]<maxn;j  ){
			is_prime[prime[j]*i] = 1;
			if(i%prime[j]==0){
				phi[i*prime[j]] = phi[i]*prime[j];
			}else{
				phi[i*prime[j]] = phi[i]*(prime[j]-1);
			}
		}
	}
}
void solve(){
	getphi();
	for(int i=2;i<maxn;i  ){
		for(int j=1;j*i<maxn;j  ){
			a[j*i] =j*phi[i];
		}
	}
	for(int i=1;i<maxn;i  )
		a[i] =a[i-1];
}
int main() {
	solve();
	int n;
	while(scanf("%d",&n)&&n){
		printf("%lldn",a[n]);
	}
	return 0;
}

0 人点赞