Processing math: 100%

洛谷 - P1390 - 公约数的和 - 莫比乌斯反演 - 欧拉函数

https://www.luogu.org/problemnew/show/P1390

求 $\sum\limits_{i=1}{n}\sum\limits_{j=1} gcd(i,j) $

不会,看题解:

类似求gcd为p的求法:

$ f(n) = \sum\limits_{i=1}{n}\sum\limits_{j=1} gcd(i,j) =\sum\limits_{i=1}^{N} d \sum\limits_{i=1}{n}\sum\limits_{j=1} [gcd(i,j)==d] $

提出 d :
f(n)=Ni=1dndi=1mdj=1[gcd(i,j)==1]

d|nμ(d)=[n==1] 替换,反演:
f(n)=Ni=1dNk=1μ(k)nkdmkd

T=kd :
f(n)=NT=1d|Tdμ(Td)nTmT

提出 T
f(n)=NT=1nTmTd|Tdμ(Td)

因为:
d|nμ(d)d=φ(n)n

f(n)=NT=1nTmTφ(T)

#include<bits/stdc++.h>
using namespace std;
#define ll long long

#define N 2000005
int phi[N],pri[N],cntpri=0;
bool notpri[N];

void sieve_phi(int n)
{
    notpri[1]=phi[1]=1;
    for (int i=2;i<=n;i++)
    {
        if (!notpri[i]) pri[++cntpri]=i,phi[i]=i-1;
        for (int j=1;j<=cntpri&&i*pri[j]<=n;j++)
        {
            notpri[i*pri[j]]=1;
            if (i%pri[j]) phi[i*pri[j]]=phi[i]*phi[pri[j]];
            else {phi[i*pri[j]]=phi[i]*pri[j];break;}
        }
    }
}

int main(){
    int n;
    cin>>n;
    sieve_phi(n);
    ll ans=0;
    for(int i=1;i<=n;i++){
        ans+=1ll*phi[i]*(n/i)*(n/i);
    }
    cout<<(ans-(1ll*(1+n)*n)/2)/2<<endl;
}


另一种奇怪的做法:
$ f(n) = \sum\limits_{d=1}^{n} d \sum\limits_{i=1}{n}\sum\limits_{j=1} [gcd(i,j)==d] $

提d:
$ \sum\limits_{d=1}^{n} d \sum\limits_{i=1}{\frac{n}{d}}\sum\limits_{j=1} [gcd(i,j)==1] $

后面是欧拉函数的定义:
nd=1dndi=2φ(i)

这里有个bug是因为1和1互质但是1和1相同,所以要去掉 φ(1)

#include<bits/stdc++.h>
using namespace std;
#define ll long long

#define N 2000000+5

int phi[N],pri[N],cntpri=0;
bool notpri[N];

ll prefix[N];

void sieve_phi(int n) {
    notpri[1]=phi[1]=1;
    prefix[0]=0;
    prefix[1]=1;
    for(int i=2; i<=n; i++) {
        if(!notpri[i])
            pri[++cntpri]=i,phi[i]=i-1;
        for(int j=1; j<=cntpri&&i*pri[j]<=n; j++) {
            notpri[i*pri[j]]=1;
            if(i%pri[j])
                phi[i*pri[j]]=phi[i]*phi[pri[j]];
            else {
                phi[i*pri[j]]=phi[i]*pri[j];
                break;
            }
        }
        prefix[i]=prefix[i-1]+phi[i];
    }
}

ll solve(ll n){
    ll ans=0;
    for(int d=1;d<=n;d++){
        ans+=d*((prefix[n/d])-1);
    }
    return ans;
}

int main() {
    sieve_phi(2000000+1);
    int n;
    while(cin>>n) {
        ll ans=solve(n);
        cout<<ans<<endl;
    }
}

posted @   韵意  阅读(107)  评论(0编辑  收藏  举报
(评论功能已被禁用)
编辑推荐:
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
· 现代计算机视觉入门之:什么是图片特征编码
· .NET 9 new features-C#13新的锁类型和语义
· Linux系统下SQL Server数据库镜像配置全流程详解
阅读排行:
· Sdcb Chats 技术博客:数据库 ID 选型的曲折之路 - 从 Guid 到自增 ID,再到
· 语音处理 开源项目 EchoSharp
· 《HelloGitHub》第 106 期
· mysql8.0无备份通过idb文件恢复数据过程、idb文件修复和tablespace id不一致处
· 使用 Dify + LLM 构建精确任务处理应用
点击右上角即可分享
微信分享提示