hdu6217BBP Formula

题目

题解

BBP算法
题解

标程：

#include
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;i1000;i++) v+=powf(16,n-i)/(i*8+k);
    return v;
}
int main(){
    scanf("%d",&T);
    while (T--){
        scanf("%d",&n);
        n--;//求小数点后第n位就要把小数点向右移n-1位
        v=bbp(1)*4-bbp(4)*2-bbp(5)-bbp(6);
        v-=(int)v;
        if (v<0) v+=1;//bbp()中有取模，计算有减法，所以有负数
        v*=16;
        printf("Case #%d: %d ",++Case,n+1);
        out((int)v);
    }
}