[SDOI2012]Longge的问题

莫比乌斯反演???
好像不用也行。
因为gcd(i,N)一定是N的因数,枚举N的每个因数,找一下有几个数和N的gcd是这个数就行了。那么这样的话(i/gcd,N/gcd)=1。求的东西就转化成:
d|Nϕ(d)dd|Nϕ(d)d
因为N太大了,所以枚举因子,边算ϕϕ即可

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;
long long n,ans;
long long ph(long long x) {
    long long tp=sqrt(x),a=x;
    for(long long i=2ll;i<=tp;i++) {
        if(x%i==0) {
            a=a/i*(i-1ll);
            while(x%i==0) x/=i;
        }
    }
    if(x>1ll) a=a/x*(x-1ll);
    return a;
}
int main() {
    cin>>n;
    long long tp=sqrt(n);
    for(long long i=1ll;i<=tp;i++) {
        if(n%i==0) {
            ans+=(ph(n/i)*i)+(i*i!=n)*(ph(i)*(n/i));
        }
    }
    cout<<ans;
}
posted @   SWHsz  阅读(98)  评论(0编辑  收藏  举报
编辑推荐:
· 智能桌面机器人:用.NET IoT库控制舵机并多方法播放表情
· Linux glibc自带哈希表的用例及性能测试
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
阅读排行:
· 手把手教你在本地部署DeepSeek R1,搭建web-ui ,建议收藏!
· 新年开篇:在本地部署DeepSeek大模型实现联网增强的AI应用
· Janus Pro:DeepSeek 开源革新,多模态 AI 的未来
· 互联网不景气了那就玩玩嵌入式吧,用纯.NET开发并制作一个智能桌面机器人(三):用.NET IoT库
· 【非技术】说说2024年我都干了些啥
点击右上角即可分享
微信分享提示