2017.9.19 Gcd 思考记录

一开始想错了

有时候还是要用欧拉函数的。。

用反演的话需要提取质数，比较麻烦，

先利用初中奥数分成gcd和两个质数，然后剩下的就是互质了

设i=gcd*a , j=gcd*b , b>a所以

ans= phi（2/2）+phi（4/2）+phi（6/2）+...phi（n/2）

+ phi（3/3）+phi（6/3）+phi（9/3）+...phi（n/3）

...

前缀和即可

码：

#include<iostream>
#include<cstdio>
using namespace std;
#define N 10000007
long long su[N],phi[N],n,i,ans,tot,lin,j;
bool he[N];
void eular(int n)
{phi[1]=1;for(i=2;i<=n;i++){if(!he[i])  {su[++tot]=i;phi[i]=i-1;            }for(j=1;j<=tot&&su[j]*i<=n;j++){he[su[j]*i]=1;if(i%su[j]==0){phi[i*su[j]]=phi[i]*su[j];break;}else phi[i*su[j]]=phi[su[j]]*phi[i]; }   }   for(i=2;i<=n;i++)phi[i]+=phi[i-1];
}
int main()
{scanf("%lld",&n);eular(n);
for(i=1;i<=tot;i++)
{//cout<<su[i];ans+=(phi[n/su[i]]*2-1);
}printf("%lld",ans);
}

2017.9.19 Gcd 思考记录相关推荐

2017.7.19 Race 思考记录
此题没有AC!! 但对拍了几个小时了,没有发现问题..但就是过不了如果有路过的dalao可以帮忙指导一下..感激不尽码(秒WA): #include<iostream> #in ...
2017.4.19 数列思考记录
这个题还是挺好的,虽然搞了半天这个题特殊性质挺多的,抓住一个就好办了首先对于后面的次方数,前面的数的和肯定小于后面的一个次方数所以没一个次方数的前面和后面一个都是独立的可以考虑递推 ...
2017.10.5 高速公路思考记录
..一开始死活不知道样例8/3是怎么算出来的.. 后来才发现右端点要-1.. 这个题说是期望,但因为每个点的概率相同,所以直接求和然后除就可以了维护每个区间的答案,听说直接维护答案很难写,,感觉还可 ...
2017.9.5 组合数学思考记录
这题风格很正统.所以我感觉如果自己能不靠题解把这个题做出来,那应该就学得不错了.. 于是在历尽千辛万苦之后,终于1A了.(事后证明结论题的结论强行无视乱搞也是有可能AC的) xym&ISA ...
2017.9.4 Nim 思考记录
只要知道nim的结论这题就是链剖板子题.1A 注意bzoj 的换行符是 \r!! \r!! \r!! 就说一下链剖的易错点吧: 1.记录父节点.sz.深度 2.当 ...
2017.6.6 数表思考记录
这个题卡在最后一步了. 最原始的式子是: 这个d|i 且 d|j 看起来非常和善,是不是可以反演? 但d要求同时能整除i和j,怎么办呢? 由于i和j的最大公因数一定包含所有因数,所以我们可以对gc ...
2017.10.28 排序思考记录
这个题有一种套路,就是大小关系转化成01串,这样就变成了二分检验问题,, 就是把排序变成区间修改,然后单点查询.. 把所有比他小的赋成0,比他大的赋成1 然后判断要求位是0还是1来判断答案与当前值的大 ...
2017.10.28 压缩思考记录
这个题首先要知道M和R不是栈的搭配,是贪心的搭配,,就是向左找到第一个M.. 如 abababab cdcdcdcd abababab cdcdcdcd = abRRcdcdcdcdR != Ma ...
2017.10.6 单词思考记录
少有的会做的ac自动机题 .这个题是匹配模式串. 所以就不需要匹配的过程,只需要建出fail树,然后直接拓扑由于每次加进去一个单词,这个单词的每一个字母都有可能作为最后一个字母影响他的失配,所以词路 ...

2017.9.19 Gcd 思考记录

2017.9.19 Gcd 思考记录相关推荐

最新文章

热门文章