http://www.bnuoj.com/v3/problem_show.php?pid=19844
输入正整数n,求gcd(1,2)+gcd(1,3)+...+gcd(1,n) + gcd(2,3)+...+g(2,n) + ...... + gcd(n-1,n)
即所有满足1<=i<j<=n的数对(i,j)对应的gcd(i,j)之和,
n=10答案为67 , n=100答案为13015。
分析:设f(n)=gcd(1,n)+gcd(2,n)+...+gcd(n-1,n),则所求答案S(n)=f(2)+f(3)+...+f(n),
可以看出S(n)=S(n-1)+f(n)
解析f(n):因为gcd(x,n)=i一定是n的约数,而且很多gcd(x1,n)=gcd(x2,n)=...=i,有多少这样的gcd呢?
gcd(x,n)=i <==> gcd(x/i,n/i)=1 ,所以满足条件的x/i有phi[n/i]个,即满足gcd(x,n)=i的x有phi[n/i]个。
所以f(n)=∑(i*phi[n/i]) | i是n的约数
由此我们可以想到类似素数筛法的方式预处理出f数组。
p之后S(n)=S(n-1)+f(n)预处理出S数组即可
#include<cstring>
#include<cstdio>
#define MT(x,i) memset(x,i,sizeof(x))
using namespace std;
typedef long long LL;
const int MAX = 4000000;
int phi[MAX+5];
void phi_table(int n)///欧拉函数表
{
for(int i=2;i<=n;i++) phi[i]=0;
phi[1]=1;
for(int i=2;i<=n;i++) if(!phi[i])
for(int j=i;j<=n;j+=i){
if(!phi[j]) phi[j]=j;
phi[j]=phi[j]/i*(i-1);
}
}
LL S[MAX+5],f[MAX+5];
void init()
{
phi_table(MAX);
///初始化f
MT(f,0);
for(int i=1;i<=MAX;i++)
for(int n=i*2;n<=MAX;n+=i) f[n]+=(LL)i*phi[n/i];
///初始化S
S[2]=f[2];
for(int n=3;n<=MAX;n++) S[n]=S[n-1]+f[n];
}
int main()
{
init();
int n;
while(scanf("%d",&n)&&n){
printf("%lld\n",S[n]);
}
return 0;
}