GCD - Extreme (II) (UVA - 11426)欧拉函数

题目链接

一道感觉挺有意思的题,自己想了好久没想出来,最后看了其他人的博客,才发现是这样哦。

如果设f(x)为x和比他小的数的gcd之和 s(n)=f(n)+f(n-1)+......+f(1);

s(n)为所要求计算的C;

f(12)=gcd(1,12)+gcd(2,12)+gcd(3,12)+gcd(4,12)+gcd(5,12)+gcd(6,12)+gcd(7,12)+gcd(8,12)+gcd(9,12)+gcd(10,12)+gcd(11,12)

       =1+2+3+4+1+6+1+4+3+2+1

其中1的数量为phi(12);

其中2的数量为2个,12=2*6.如果gcd(12,x)=2,则12/2和x/2是互素的。而x/2<6 且x/2还和6互素 那x的数量不就是phi(6)嘛。

这样的话f(i*j)+=i*phi(j);计算出所有的i和j就可以算出f(n),然后也就可以算出 s(n);

#include<iostream>
#include<cstring>
using namespace std;
const int maxn=4000010;
long long  euler[maxn]; 
long long  f[maxn];
void getEuler() { 
	memset(euler,0,sizeof(euler)); 
	euler[1] = 1; 
	for(int i = 2;i < maxn;i++) {
		if(!euler[i]) 
			for(int j = i;j < maxn; j += i) { 
				if(!euler[j]) euler[j] = j; 
				euler[j] = euler[j]/i*(i-1); 
			} 
		for(int j=1;j*i<maxn;j++){
			f[j*i]+=j*euler[i];
		}
	}
	for(int i=1;i<maxn;i++){
		f[i]+=f[i-1];
	} 
}
int main()
{
	getEuler();
	int n;
	while(cin>>n,n) cout<<f[n]<<endl;
	return 0;
}


  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值