欧拉 uva 11426 - GCD - Extreme (II)

关注

原创最新推荐文章于 2023-01-26 22:23:01 发布 456 阅读 AI 写同款· GEO 优化 ›

UVA11426 GCD - Extreme (II) (欧拉函数/莫比乌斯反演) UVA11426 GCD - Extreme (II) 题目描述 PDF 输入输出格式输入格式：输出格式：输入输出样例输入样例#1： 10 100 200000 0 输出样例#1： 67 13015 143295493160 Solution 这道题我用莫比乌斯反演和欧拉函数都写了一遍,发现欧拉函数比莫比乌斯反演优秀? 求所有\(gcd=k\)的数对的个数,记作\(f[k... 阅读详情

题意：（详见蓝皮书P125）

给定正整数n,求下列表达式的值：

G =
i<N j<=N
∑ ∑ GCD(i, j)
i=1 j=i+1

Input

The input file contains at most 100 lines of inputs. Each line contains an integer N (1 < N < 4000001).The meaning of N is given in the problem statement. Input is terminated by a line containing a singlezero.

Output

For each line of input produce one line of output. This line contains the value of G for the correspondingN. The value of G will fit in a 64-bit signed integer.

Sample Input

100

2000000

Sample Output

13015

143295493160

分析：

记f[n]= ∑gcd(i,n) (1<=i<=n) ,那么所求就是f[1]+f[2]+......+f[n]=sum[n];

那么怎么求f[n]呢?

设d=gcd（x,n）,显然d为n的因数，

那么 gcd(x/d,n/d)==1; （*）

根据phi的定义，（*）式有 phi(n/d)个解，那么这一种情况的解的和就是 d*phi(n/d);

这样用逐个分解因数的方法求一个f[n]的时间复杂度为n^(1/2)，显然会TLE。

可以模拟筛法求phi的过程求出f[n]。

注意使用 long long

代码如下：

// by spark 2016/5/21 Friday 23:50
#include<cstdio>
#include<iostream>
#define LL long long
using namespace std;
const int maxn=4000000+5;
int phi[maxn];
LL sum[maxn],f[maxn];
void phi_table(int n){
	int i,j;
	for(i=1;i<=n;i++)phi[i]=i;
	for(i=2;i<=n;i++){
		if(phi[i]==i)
			for(j=i;j<=n;j+=i)
				phi[j]=phi[j]/i*(i-1);
	}
}
void f_table(int n){
	int i,j;
	for(i=1;i<=n;i++)
		for(j=i+i;j<=n;j+=i)
			f[j]+=i*phi[j/i];
}
int main(){
	int n,i;
	phi_table(maxn);   //打出phi表 
	f_table(maxn);    //计算f 
	for(i=1;i<=maxn;i++)
		sum[i]=sum[i-1]+f[i];    //前缀和 
	while(true){
		cin>>n;
		if(!n) return 0;
		cout<<sum[n]<<endl;
	}
	return 0;
}