uva11426(gcd + 欧拉函数)

题意:

给出一个n,求

 

思路:

设f(n)=gcd(1,n)+gcd(2,n)+……+gcd(n-1,n)。


gcd(x,n)=i是n的约数(x<n),按照这个约数进行分类。设满足gcd(x,n)=i的约束有g(n,i)个,则有f(n)=sum(i*g(n,i))。


而gcd(x,n)=i等价于gcd(x/i,n/i)=1,因此g(n,i)等价于p(n/i).p(x)为欧拉函数。



#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define ll long long

const int N = 4000001;
ll sum[N];
ll f[N];
ll p[N];
int n;

void init(){
	memset(p, 0, sizeof(p));
	p[1] = 1;
	for(int i = 2; i < N; i++) {
		if(!p[i]) {
			for(int j = i; j < N; j += i) {
				if(!p[j])
					p[j] = j;
				p[j] = p[j] / i * (i - 1);
			}
		}	
	}
	memset(sum, 0, sizeof(sum));
	for(int i = 1; i < N; i++) {
		for(int j = i * 2; j < N; j += i) {
			sum[j] += i * p[j / i];
		}
	}
	f[2] = sum[2];
	for(int i = 3; i < N; i++)
		f[i] = f[i - 1] + sum[i];
}
int main() {
	init();
	while(scanf("%d", &n) && n)	{
		printf("%lld\n", f[n]);
	}
}





评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值