题目描述
令g(n)表示n能表示成几种不同的完全k次方数(k>1),求。
例如,,所以g(64)=3。
分析:
1.首先将题目进行转化。转化为枚举底数分别求贡献,则答案就是先枚举一个次数k(2<=k<=log n),再求和。
2.于是我们就可以预处理出前
105
项,剩下的每一项用代替,于是就是要计算。
#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;
}
^_^