题目
给定r,n求[1,r]内与n互素的个数有多少个?
解题思路
正面想太过复杂,所以从反面来求,那么问题就变成[1,r]中有多少个和n不互质的的数。而要解决这个问题就可以把n的质数都枚举出来,在通过这些质数倍数个数和,减去重叠部分(容斥定理)就可以了。而我写这篇博客的目的是下面代码对位运算的应用让我吃惊。。。祭奠这次发现新大陆。
代码部分
LL Solve(LL n,LL r)
{
vector<LL> p;
for(LL i=2; i*i<=n; i++)
{
if(n%i==0)
{
p.push_back(i);
while(n%i==0) n/=i;
}
}
if(n>1)
p.push_back(n);
LL ans=0;
for(LL msk=1; msk<(1<<p.size()); msk++) ///用2进制存储每个质数是否被使用
{
LL multi=1,bits=0;
for(LL i=0; i<p.size(); i++)
{
if(msk&(1<<i)) //判断第几个因子目前被用到
{
++bits;
multi*=p[i];
}
}
LL cur=r/multi; ///容斥定理
if(bits&1) ans+=cur;
else ans-=cur;
}
return r-ans;
}