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

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

题意：

$给定正整数N，计算：$

$$G=\sum _{i=1}^{i<N}\sum _{j=i+1}^{j\le N}GCD(i,j)$$

$其中，GCD(i,j)表示i和j的最大公约数。$

输入：

$多组测试数据，读到0为止。$

$每组包括一个正整数N(1<N<4000001)。$

输出：

$一个正整数，表示答案。$

Sample Input

10
100
200000
0

Sample Output

67
13015
143295493160

---

分析：

$$\sum _{i=1}^{i<N}\sum _{j=i+1}^{j\le N}GCD(i,j)=\sum _{j=1}^{j\le N}\sum _{i=1}^{j-1}GCD(i,j)$$

$即对所有的j\in [1,N]，求[1,j-1]中的所有数与j的最大公约数的和。$

$记GCD(i,j)=g，则:$

$$\sum _{j=1}^{j\le N}\sum _{i=1}^{j-1}GCD(i,j)=\sum _{g=1}^{N}g\times \sum _{j=1}^{j\le N}\sum _{i=1}^{j-1}[GCD(i,j)=g]=\sum _{g=1}^{N}g\times \sum _{j=1}^{j\le \lfloor \frac{N}{g}\rfloor }\sum _{i=1}^{j-1}[GCD(i,j)=1]$$

$问题转化为，对每一个g\in [1,N]，统计[1,j-1]中与j互质的数的个数，j\in [1,N]。$

$由此发现，我们需要预处理[1,\lfloor \frac{N}{g}\rfloor ]的欧拉函数。$

$答案即：$

$$\sum _{g=1}^{N}g\times \sum _{j=1}^{j\le \lfloor \frac{N}{g}\rfloor }\varphi (j)$$

$最后我们枚举g，计算上式。$

优化：

$①、预处理\varphi (i)的前缀和。$

$②、由于我们枚举g，\lfloor \frac{N}{g}\rfloor 的值可能重复，这里利用整除分块优化。$

$故还需预处理前缀和sum[x]={\sum }_{i=1}^{x}i。$

代码：

#include<cstring>
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>

#define ll long long

using namespace std;

const int N=4e6+10;

int primes[N], cnt;
bool st[N];
ll f[N], sum[N], phi[N];

void get_prime(int n)
{
    phi[1]=0;
    for(int i=2;i<=n;i++)
    {
        if(!st[i]) 
        {
            primes[cnt++]=i;
            phi[i]=i-1;
        }
        for(int j=0;primes[j]*i<=n;j++)
        {
            int p=primes[j];
            st[p*i]=true;
            if(i%p==0)
            {
                phi[p*i]=phi[i]*p;
                break;
            }
            else phi[p*i]=(p-1)*phi[i];
        }
    }
    for(int i=1;i<=n;i++) phi[i]+=phi[i-1];
    for(int i=1;i<=n;i++) sum[i]=sum[i-1]+i;
}

int main()
{  
    get_prime(N-1);

    int n;
    while(~scanf("%d",&n),n)
    {
        ll res=0;
        for(int i=1,j; i<=n; i=j+1)
        {
            j=n/(n/i);
            res+=phi[n/i]*(sum[j]-sum[i-1]);
        }
        printf("%lld\n",res);
    }
    
    return 0;
}