容斥原理部分例题及板子

最新推荐文章于 2024-08-10 18:27:05 发布

zsmj_sy

最新推荐文章于 2024-08-10 18:27:05 发布

阅读量530

点赞数

分类专栏：组合数学数学数论容斥

本文链接：https://blog.csdn.net/neuq_zsmj/article/details/82946746

版权

本文介绍了几个利用高精度计算和容斥原理解决的数论问题，包括求解一定范围内满足特定指数形式的数的数量、求解最大指数之和以及涉及质因子和最大公约数的问题。通过分析题意，转换思路，运用数学方法和算法，如高精度开根号、容斥原理中的奇加偶减策略，来解决这些竞赛题目。

摘要由CSDN通过智能技术生成

高精度开k次根号


ll pow_mul(ll x,ll k)
{
    ll ans=1;
    while(k)
    {
        if(k&1)
        {
            double judge=1.0*INF/ans;
            if(x>judge)return -1;
            ans*=x;
        }
        k>>=1;
        if(x>T&&k>0)return -1;
        x*=x;
    }
    return ans;
}

ll find(ll x,ll k)
{
    ll r=(ll)pow(x,1.0/k);
    ll p=pow_mul(r,k);
    if(p==x)return r;
    if(p>x||p==-1)r--;
    else
    {
        ll tmp=pow_mul(r+1,k);
        if(tmp!=-1&&tmp<=x)r++;
    }
    return r;
}

hdu 2204

https://vjudge.net/problem/HDU-2204

题意：求1到n中满足a=m^k的数量(m>=2) (n<=1e18)

假如k为合数，我们一定可以把它转为 m'^k'（k'为素数）的形式，所以我们只要考虑k为素数的情况

2^60>1e18,所以满足要求且在范围内的k最大不超过60。对于一个给定的n，我们求出最大的k，将所有小于等于k的素数全部放入vector中。当n和k给定时，所有满足a=m^k的数量就是m的最大值，也就是pow(n,1/k)，这个地方需要注意精度。

而对于一个数a，最多只有一组m^k是合法的，例如16=2^4=4^2，但k要取最大只能是2^4

所以这个地方需要容斥，奇加偶减

ll n;
vector<int>p;
bool check(int a)
{
    for(int i=2;i*i<=a;i++)
    {
        if(a%i==0)
        {
            return 0;
        }
    }
    return 1;
}

void get_prime(ll n)
{
    p.clear();
    for (int i = 2; i <= 65; i++)
    {
        if(!check(i))continue;
        if ((1ll << i )> n)
        {
            break;
        }
        p.push_back(i);
    }
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);cout.tie(0);

    while(cin>>n)
    {
        ll ans=0;
        get_prime(n);
        int len=p.size();
        for(int i=1;i<(1ll<<len);i++)
        {
            ll mul=1;
            int flag=0;
            for(int j=0;j<len;j++)
            {
                if(i&(1<<j))
                {
                    mul*=p[j];
                    flag++;
                }
            }
            ll now=(ll)(pow((double)n,(1.0/mul))+eps);
            if(flag%2)
            {
                ans+=now;
            }
            else
            {
                ans-=now;
            }
        }
        cout<<ans<<endl;
    }

    return 0;
}

hdu 3208

https://vjudge.net/problem/HDU-3208

题意：求[a-b]中每个数可以被分解为m^k（k取最大）的最大指数之和

可以求1-n的最大指数之和，然后减一下。

设一个num数组，num[i]表示指数为i的数的数量，最后遍历乘i就得到了总数。

num[i]的求法和上一题很像，1-n中满足a=m^i的数量就是确定了n和i之后pow(n,1/k)这个地方要高精度开k次根。<