其他分享
首页 > 其他分享> > 【HDU2204】Eddy's爱好

【HDU2204】Eddy's爱好

作者:互联网

题目大意:求从 1 到 N 中共有多少个数可以表示成 \(M^K,K \gt 1\)。\(N \le 1e18\)

题解:
发现 N 很大,若直接枚举 M 的话有 1e9 级别的数据量,肯定超时,因此考虑枚举幂次。发现对于幂次为 k 的符合条件的数有 N 开 K 次方下取整个,同时注意到 k 的取值范围最大为 60,因为 2 的 60 次方为 1e18 级别。因此考虑从小到大进行枚举幂次即可,但是发现有些数字会产生重复,如:\((2^3)^2=(2^2)^3=2^6\),即:同一个数字被计入了三次贡献,因此涉及到了容斥原理,即:若幂次的质因子分解为奇数个时累加贡献,反之减去贡献,因此只需对 1-60 中的素数进行容斥操作即可,即:问题转化成了多重集合的组合问题。

代码如下

#include <bits/stdc++.h>
using namespace std;
const int maxn=61;
typedef long long LL;

LL n;
int prime[]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61};

void solve(){
    LL ans=0;
    int ub=0;
    while((1LL<<prime[ub+1])<=n)++ub;
    for(int i=1;i<1<<ub;i++){
        LL ret=1;
        int cnt=0;
        for(int j=0;j<ub;j++)
            if(i>>j&1){
                ++cnt,ret*=prime[j];
            }
        if(cnt&1)ans+=(LL)(pow(n,1.0/ret)+1e-8);
        else ans-=(LL)(pow(n,1.0/ret)+1e-8);
    }
    printf("%lld\n",ans+1);
}
int main(){
    while(scanf("%lld",&n)!=EOF){
        solve();
    }   
    return 0;
} 

标签:幂次,int,HDU2204,爱好,ret,60,Eddy,ans,LL
来源: https://www.cnblogs.com/wzj-xhjbk/p/10987589.html