题目
题解
BBP算法
题解
标程:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
double v;
int n,T,Case;
void out(int x){
if (x<10) putchar(x+48);
else putchar(x-10+'A');
putchar('\n');
}
ll pow(ll x,ll y,ll M){
ll z=1;
for (;y;y>>=1,x=x*x%M)
if (y&1) z=z*x%M;
return z;
}
double bbp(ll k){
double v=0;
for (int i=0;i<=n;i++) v+=pow(16,n-i,i*8+k)*1.0/(i*8+k);
for (int i=n+1;i<n+1000;i++) v+=powf(16,n-i)/(i*8+k);
return v;
}
int main(){
scanf("%d",&T);
while (T--){
scanf("%d",&n);
n--;
v=bbp(1)*4-bbp(4)*2-bbp(5)-bbp(6);
v-=(int)v;
if (v<0) v+=1;
v*=16;
printf("Case #%d: %d ",++Case,n+1);
out((int)v);
}
}