2017.9.14 仪仗队 思考记录
2017-09-14 14:36
211 查看
似乎只要会了莫比乌斯反演就不用会phi了、
感觉反演更顺手、
码:
#include<iostream>
#include<cstdio>
using namespace std;
#define N 40005
int mu
,i,j,n,tot,ans,su
;
bool he
;
void eular()
{
mu[1]=1;
for(i=2;i<=n;i++)
{
if(!he[i])
{
su[++tot]=i;
mu[i]=-1;
}
for(j=1;j<=tot&&su[j]*i<=n;j++)
{
he[su[j]*i]=1;
if(i%su[j]==0)
{
mu[i*su[j]]=0;
break;
}else mu[i*su[j]]=-mu[i];
}
}
}
int main()
{
scanf("%d",&n);
if(n==1)
{
printf("0");
return 0;
}
n--;
eular();
for(i=1;i<=n;i++)
ans+=mu[i]*(n/i)*(n/i);
//cout<<ans;
ans+=2;
printf("%d",ans);
}
感觉反演更顺手、
码:
#include<iostream>
#include<cstdio>
using namespace std;
#define N 40005
int mu
,i,j,n,tot,ans,su
;
bool he
;
void eular()
{
mu[1]=1;
for(i=2;i<=n;i++)
{
if(!he[i])
{
su[++tot]=i;
mu[i]=-1;
}
for(j=1;j<=tot&&su[j]*i<=n;j++)
{
he[su[j]*i]=1;
if(i%su[j]==0)
{
mu[i*su[j]]=0;
break;
}else mu[i*su[j]]=-mu[i];
}
}
}
int main()
{
scanf("%d",&n);
if(n==1)
{
printf("0");
return 0;
}
n--;
eular();
for(i=1;i<=n;i++)
ans+=mu[i]*(n/i)*(n/i);
//cout<<ans;
ans+=2;
printf("%d",ans);
}
相关文章推荐
- 2017.9.14 棘手的操作 思考记录
- 2017.9.14 dispatching 思考记录
- 某系统日志记录的思考----第三篇:难点
- 2017.9.28 降雨量 思考记录
- 2017.4.19 细胞分裂 思考记录
- 2017.9.29 数三角形 思考记录
- 阅读 思考 感悟 记录 升华 幸福
- 2017.10.10 狼和羊的故事 思考记录
- 2017.7.18 树上距离=k 思考记录
- 2017.5.9 寻找道路 思考记录
- 写博文不仅仅是记录,更是思考、总结和提高
- 2017.5.11 道路修建 思考记录
- 2017.10.24 上升序列 思考记录
- 2017.9.9 股票交易 思考记录
- 代码中“抛异常”与“记录log"的思考
- 2010webRebuild 大会 广州站记录 & 后记 (在职前端人员的发展与思考)
- 【开卷有益】记录一次高并发下的死锁解决思考过程
- 2017.9.25 社交网络 思考记录
- 2017.4.22 进制转换 思考记录
- 2017.10.7 弹飞绵羊 思考记录