来一发莫比乌斯反演的题解。

\[\begin{split} ans & =\sum^{n-1}_{i=1} \sum^{n-1}_{j=1} [\gcd(i,j)=1]\\ & = \sum^{n-1}_{i=1} \sum^{n-1}_{j=1} \sum_{d|\gcd(i,j)} \mu(d)\\ & = \sum^n_{d=1}\mu(d)\sum^{\frac{n-1}{d}}_{i=1}\sum^{\frac{n-1}{d}}_{j=1}\\ & =\sum^n_{d=1}\mu(d)\times\left(\dfrac{n-1}{d}\right)^2 \end{split} \]

\(\dfrac{n-1}{g}\) 可以用数论分块求出。

//P2158 [SDOI2008] 仪仗队
#include <cstdio>
#include <cstdlib>
using namespace std;
const int MAXN=40005;
int prime[MAXN],cnt,mu[MAXN],s[MAXN],n,ans;
bool flag[MAXN];
void sieve()
{
    mu[1]=1;
    for(int i=2;i<=4e4;i++)
    {
        if(!flag[i])
        {
            prime[++cnt]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=cnt&&i*prime[j]<=4e4;j++)
        {
            flag[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                break;
            }
            mu[i*prime[j]]=-mu[i];
        }
    }
    return;
}
int main()
{
    sieve();
    scanf("%d",&n);
    if(!--n)
    {
        printf("0\n");
        exit(0);
    }
    for(int i=1;i<=n;i++)
    {
        s[i]=s[i-1]+mu[i];
    }
    for(int l=1,r;l<=n;l=r+1)
    {
        r=n/(n/l);
        ans+=(s[r]-s[l-1])*(n/l)*(n/r);
    }
    printf("%d\n",ans+2);
    return 0;
}
/*
 * 洛谷
 * https://www.luogu.com.cn/problem/P2158
 * C++17 -O0
 * 2022.10.12
 */

原文地址：http://www.cnblogs.com/2020gyk080/p/16784088.html

1. 本站所有资源来源于用户上传和网络，如有侵权请邮件联系站长！ 2. 分享目的仅供大家学习和交流，请务用于商业用途! 3. 如果你也有好源码或者教程，可以到用户中心发布，分享有积分奖励和额外收入！ 4. 本站提供的源码、模板、插件等等其他资源，都不包含技术服务请大家谅解！ 5. 如有链接无法下载、失效或广告，请联系管理员处理！ 6. 本站资源售价只是赞助，收取费用仅维持本站的日常运营所需！ 7. 如遇到加密压缩包，默认解压密码为"gltf",如遇到无法解压的请联系管理员！ 8. 因为资源和程序源码均为可复制品，所以不支持任何理由的退款兑现，请斟酌后支付下载 声明：如果标题没有注明"已测试"或者"测试可用"等字样的资源源码均未经过站长测试.特别注意没有标注的源码不保证任何可用性