bzoj2440 [中山市选2011]完全平方数莫比乌斯函数应用

最新推荐文章于 2025-11-26 13:52:00 发布

转载最新推荐文章于 2025-11-26 13:52:00 发布 · 96 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/mybing/p/8327602.html

文章标签：

#c/c++

本文提供了一道 ACM 题目 2440 的解题思路及 C++ 实现代码。该题目通过二分查找将问题进行转化，并运用了容斥原理与莫比乌斯函数来解决。文章包含了完整的代码实现及解析。

先用二分把问题转化一下，然后利用容斥定理和莫比乌斯函数搞一搞就好了。

/**************************************************************
    Problem: 2440
    User: mybing
    Language: C++
    Result: Accepted
    Time:4908 ms
    Memory:2556 kb
****************************************************************/
 
#include<bits/stdc++.h>
using namespace std;
int f[100010],tot=0,prime[100010],mu[100010],T,n;
 
bool not_prime[100010];
bool pd(long long p,long long w)
{
    long long sum=0;
    for(long long i=1;i*i<=p;i++)
    {
        sum+=p/(i*i)*mu[i];
    }
    if(sum>=w)return true;
    else return false;
}
void work(int p)
{
    long long l=1,r=3000000000;
    while(l+1<r)
    {
        long long mid=(l+r)>>1;
        if(pd(mid,(long long)p))r=mid;
        else l=mid;
    }
    if(pd(l,n))printf("%lld\n",l);
    else printf("%lld\n",r);
}
int main()
{
    memset(not_prime,true,sizeof(not_prime));
    mu[1]=1;
    for(int i=2;i<=100000;i++)
    {
        if(not_prime[i])
        {
            prime[++tot]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=tot&&i*prime[j]<=100000;j++)
        {
            not_prime[i*prime[j]]=false;
            if(i%prime[j]==0)
            {
                not_prime[prime[j]*i]=false;
                mu[prime[j]*i]=0;
                break;
            }
            mu[i*prime[j]]=-1*mu[i];
        }
    }
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        work(n);
    }
    return 0;
}

转载于:https://www.cnblogs.com/mybing/p/8327602.html