给定r,n求[1,r]内与n互素的个数有多少个?(容斥定理+位运算)

题目

给定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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值