[Luogu_MayR1 D] 核心密码B

题目描述

令g(n)表示n能表示成几种不同的完全k次方数（k>1），求。

例如，，所以g(64)=3。

分析：
1.首先将题目进行转化。转化为枚举底数分别求贡献，则答案就是先枚举一个次数k（2<=k<=log n），再求和。
2.于是我们就可以预处理出前$10^5$ 项，剩下的每一项用代替，于是就是要计算。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn=100010;
typedef long double ld;
typedef long long LL;
inline ld Pow(ld x,int y){
    ld ret=1;
    while(y){
        if(y&1) ret=ret*x;
        x=x*x;
        y>>=1;
    }
    return ret;
}
ld sum[maxn][65];
int MAX=100000;
void init(){
    for(int k=2;k<=60;k++){
        int Max=pow(1000000000000000000,1.0/k);
        for(int i=2;i<=MAX;i++){
            sum[i][k]=sum[i-1][k];
            if(i<=Max) sum[i][k]+=1.0/pow(i*1.0,k);
        }
    }
}
ld calc(ld st,ld ed,LL k){
    return ((1.0/Pow(ed,k-1))-(1.0/Pow(st,k-1)))/(1.0-1.0*k);
}
LL n;
ld ans;
int main(){
    init();
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%lld",&n);
        ans=0.0;
        if(n<4){
            printf("0\n");
            continue;
        }
        int k=1;
        for(LL i=4;i<=n;i<<=1){
            k++;
            LL tmp=pow(n,1.0/k);
            if(tmp<=1LL*MAX){
                ans+=sum[tmp][k];continue;
            }
            ans+=sum[MAX][k];
            int st=MAX+1;if(st>tmp) continue;
            ans+=calc(1.0*st-0.5,1.0*tmp+0.5,k);
        }
        printf("%.15lf\n",(double)ans);
    }
    return 0;
}

^_^