dtoj3571. 乌(mi)
我有一段代码:
输入N,求G(N)的值.
sol
由于我们求
前面d的枚举是的 ,calc整数分块
那么效率是
感觉是
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
#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; }
分类:
数学
, 数学---莫比乌斯反演
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· DeepSeek 解答了困扰我五年的技术问题
· 为什么说在企业级应用开发中,后端往往是效率杀手?
· 用 C# 插值字符串处理器写一个 sscanf
· Java 中堆内存和栈内存上的数据分布和特点
· 开发中对象命名的一点思考
· DeepSeek 解答了困扰我五年的技术问题。时代确实变了!
· PPT革命!DeepSeek+Kimi=N小时工作5分钟完成?
· What?废柴, 还在本地部署DeepSeek吗?Are you kidding?
· 赶AI大潮:在VSCode中使用DeepSeek及近百种模型的极简方法
· DeepSeek企业级部署实战指南:从服务器选型到Dify私有化落地
2019-02-18 [ZJOI2006]书架