【BZOJ 2818】Gcd

最新推荐文章于 2021-11-15 16:50:24 发布

原创最新推荐文章于 2021-11-15 16:50:24 发布 · 904 阅读

0 ·

CC 4.0 BY-SA版权

数论专栏收录该内容

14 篇文章

订阅专栏

本文深入探讨了一个复杂的数学问题，通过巧妙地转换和应用欧拉函数，揭示了解决该问题的高效算法。利用线性筛法计算欧拉函数，并结合前缀和技巧，实现了O(n)的时间复杂度求解。此方法不仅简洁明了，而且在实际应用中展现出卓越的性能。

这题一开始我竟然想要用与能量采集差不多的思路去做= =（no zuo no die,why you try?）
有个显然的转化
$\sum_{x=1}^{n}\sum_{y=1}^{n} [gcd(x,y)==P]= \sum_{x=1}^{\lfloor n/P \rfloor}\sum_{y=1}^{\lfloor n/P \rfloor}[gcd(x,y)==1]=2*\sum_{i=2}^{\lfloor n/P \rfloor} \phi (i)+1$
然后我们线性筛出欧拉函数，在前缀和就可以 $O(n)$ 求出答案了
code：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
long long f[10000001];
int phi[10000001],p[700001];
long long ans=0;
long long s[10000001];
int a[10000001],n;
int main()
{
    long long i,t=0,j;
    scanf("%lld",&n);
    a[1]=1; phi[1]=1;
    for (i=2;i<=n;++i)
      {
        if (!a[i])
          {
            p[++t]=i;
            phi[i]=i-1;
          }
        for (j=1;j<=t;++j)
          {
            if (i*p[j]>n) break;
            a[i*p[j]]=1;
            if (i%p[j]==0)
              {
                phi[i*p[j]]=phi[i]*p[j];
                break;
              }
            else
              phi[i*p[j]]=phi[i]*(p[j]-1);
          }
      }
    s[1]=0;
    for (i=2;i<=n;++i) s[i]=s[i-1]+phi[i];
    for (i=1;i<=n;++i) f[i]=2*s[i]+1;
    for (i=1;i<=t;++i) ans=ans+f[n/p[i]];
    printf("%lld\n",ans);
}