博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
洛谷P2568 GCD(莫比乌斯反演)
阅读量:6852 次
发布时间:2019-06-26

本文共 884 字,大约阅读时间需要 2 分钟。

 

这题和一样……不过是n和m相同而已……

所以虽然正解是欧拉函数然而直接改改就行了所以懒得再码一遍了2333

不过这题卡空间,记得mu开short,vis开bool

1 //minamoto 2 #include
3 #define ll long long 4 const int N=1e7+5; 5 int p[1000005],n,m;short mu[N];bool vis[N];ll sum[N],ans; 6 void init(int n){ 7 mu[1]=1; 8 for(int i=2;i<=n;++i){ 9 if(!vis[i]) mu[i]=-1,p[++m]=i;10 for(int j=1;j<=m&&p[j]*i<=n;++j){11 vis[i*p[j]]=1;12 if(i%p[j]==0) break;13 mu[i*p[j]]=-mu[i];14 }15 }16 for(int j=1;j<=m;++j)17 for(int i=1;i*p[j]<=n;++i)18 sum[i*p[j]]+=mu[i];19 for(int i=1;i<=n;++i) sum[i]+=sum[i-1];20 }21 int main(){22 scanf("%d",&n),init(n);23 for(int l=1,r;l<=n;l=r+1){24 r=n/(n/l);25 ans+=(sum[r]-sum[l-1])*(n/l)*(n/l);26 }27 printf("%lld\n",ans);28 return 0;29 }

 

转载于:https://www.cnblogs.com/bztMinamoto/p/9712748.html

你可能感兴趣的文章
Extjs4.1.x 框架搭建 采用Application动态按需加载MVC各模块-完美解决(二)
查看>>
Jquery控制图片最大宽度
查看>>
[philosophy]空间
查看>>
myeclipse6.0添加maven
查看>>
【css】纯 css 制作带三角的边框
查看>>
【WinForm】创建自定义控件
查看>>
c++重新抛出异常
查看>>
dedecms 安装完成后登录后台出现空白
查看>>
HBase vs. BigTable Comparison
查看>>
JS:2.1,流程控制(if,switch)高级
查看>>
android native c 的so调试
查看>>
PHP 单一入口框架设计简析
查看>>
boost.asio系列——io_service
查看>>
Android 编程下如何获取有 Internet 访问权限的应用
查看>>
云计算里AWS和Azure的探究(6) - Amazon Simple Storage Service 和 Microsoft Azure Blob Storage
查看>>
流行趋势:27个引领时尚的网页设计作品【上篇】
查看>>
(visual)c++ 内存分配
查看>>
反射 简单工厂模式
查看>>
Linux学习之CentOS(一)--CentOS6.4环境搭建
查看>>
SQL Server 2008 R2 附加数据库错误
查看>>