UVa 11426题解 GCD - Extreme (II) 题解

…我怎么又在做水题

---

题目传送门

---

翻译：

这道题是说给一个n然后求所有满足$1<=i<j<=n$的数对$(i,j)$,求他们的$gcd$之和

---

题解：

我们来推结论

我们要求的是$gcd$之和，如果我们把$gcd(1,n),gcd(2,n),gcd(3,n)...gcd(n-1,n)$的和命为$f(n)$

如果令$S(n)$为$n$对于题目中要求的答案，那么$S(n)=f(1)+f(2)+...+f(n)$

显然我们可以递推，$S(n)=S(n-1)+f(n)$

所以现在我们只需要求出$f(n)$

对于一个$n$，我们令$g(n,i)$表示和$n$的最大公约数为$i$的数有多少个

那么$f(n)=\sum_{i=1}^ni*g(n,i)*[i|n]$

所以现在任务就是求$g(n,i)$

注意到$gcd(x,n)==i$的充要条件是$gcd(\frac{x}{i},\frac{n}{i})==1$

所以我们只需要去求和$\frac{n}{i}$互质的数有多少个就可以了，对于一个$i$，和$\frac{n}{i}$互质的有多少个数

就有多少个$\frac{x}{i}$也就有多少个$x$

所以我们在对于每个$i$枚举的时候，都找$i$的倍数，使得$i$是$n$的因数，然后更新累加枚举的$n$的答案

然后先预处理一个欧拉函数值出来就可以了

---

欧拉筛筛素数求欧拉函数代码：

void get_phi(){
    phi[1]=1;
    for(register int i=2;i<=MAXN;i++){
        if(!is_not_prime[i]){
            prime[++prime[0]]=i;
            phi[i]=i-1;
        }
        for(register int j=1;j<=prime[0]&&i*prime[j]<=MAXN;j++){
            is_not_prime[i*prime[j]]=1;
            if(i%prime[j]==0){
                phi[i*prime[j]]=prime[j]*phi[i];//important
                break;
            }else{
                phi[i*prime[j]]=phi[i]*(prime[j]-1);//gcd(i,prime[j])==1
            }//gcd(a,b)==1 -> phi[a*b]==phi[a]*phi[b]
        }
    }
}

完整代码：

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int MAXN=4000000;
long long S[MAXN+10],prime[MAXN],f[MAXN+10],phi[MAXN+10],is_not_prime[MAXN];
void get_phi(){
    phi[1]=1;
    for(register int i=2;i<=MAXN;i++){
        if(!is_not_prime[i]){
            prime[++prime[0]]=i;
            phi[i]=i-1;
        }
        for(register int j=1;j<=prime[0]&&i*prime[j]<=MAXN;j++){
            is_not_prime[i*prime[j]]=1;
            if(i%prime[j]==0){
                phi[i*prime[j]]=prime[j]*phi[i];
                break;
            }else{
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
            }
        }
    }
}
int main(){
    get_phi();
    for(register int i=1;i<=MAXN;i++)
        for(register int n=i*2;n<=MAXN;n+=i) f[n]+=i*phi[n/i];
    S[2]=f[2];
    for(register int n=3;n<=MAXN;n++)S[n]=S[n-1]+f[n];
    int n;
    while(scanf("%d",&n)==1&&n) printf("%lld\n",S[n]);
    return 0;
}