dtoj3571. 乌(mi)

我有一段代码:

 

输入N,求G(N)的值.


sol

i=1nj=1n/ih[i×j]+=[gcd(i,j)==1]

di=1n/dj=1n/idh[ijd2]+=μ(d)
由于我们求h
calc(x)=i=1nj=1n/i1
di=1n/dj=1n/idμ(d)calc(n/d2)

前面d的枚举是sqrt(n)的 ,calc整数分块
那么效率是n+n/4+n/9...
感觉是nln(n)

复制代码
#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define maxn 100005
#define mod 998244353
using namespace std;
int n,N,pri[maxn],flag[maxn],tot,mu[maxn],ans;
int calc(int x){
    int sum=0;
    for(int i=1,nex;i<=x;i=nex+1){
        nex=x/(x/i);
        sum=(sum+1LL*(nex-i+1)*(x/i)%mod)%mod;
    }
    return sum;
}
int main(){
    cin>>n;N=sqrt(n)+1;mu[1]=1;
    for(int i=2;i<=N;i++){
        if(!flag[i]){pri[++tot]=i;mu[i]=-1;}
        for(int j=1;j<=tot&&i*pri[j]<=N;j++){
            flag[i*pri[j]]=1;
            if(i%pri[j]==0){mu[i*pri[j]]=0;break;}
            mu[i*pri[j]]=-mu[i];
        }
    }
    for(int i=1;i<=N;i++){
        if(mu[i]!=0){
            ans=(ans+mu[i]*calc(n/(i*i)))%mod;
        }
    }
    ans=(ans+mod)%mod;
    cout<<ans<<endl;
    return 0;
}
View Code
复制代码

 

posted @   liankewei123456  阅读(168)  评论(0编辑  收藏  举报
编辑推荐:
· DeepSeek 解答了困扰我五年的技术问题
· 为什么说在企业级应用开发中,后端往往是效率杀手?
· 用 C# 插值字符串处理器写一个 sscanf
· Java 中堆内存和栈内存上的数据分布和特点
· 开发中对象命名的一点思考
阅读排行:
· DeepSeek 解答了困扰我五年的技术问题。时代确实变了!
· PPT革命!DeepSeek+Kimi=N小时工作5分钟完成?
· What?废柴, 还在本地部署DeepSeek吗?Are you kidding?
· 赶AI大潮:在VSCode中使用DeepSeek及近百种模型的极简方法
· DeepSeek企业级部署实战指南:从服务器选型到Dify私有化落地
历史上的今天:
2019-02-18 [ZJOI2006]书架
点击右上角即可分享
微信分享提示