HDU 6217 BBP Formula (数学)
题目链接: HDU 7217
题意:
题目给你可以计算 \(π\) 的公式:
\(\pi = \sum_{k=0}^{\infty}[\frac{1}{16^k}(\frac{4}{8k+1})-(\frac{2}{8k+4})-(\frac{1}{8k+5})-(\frac{1}{8k+6})]\)
告诉你可以求十六进制下的小数点后 \(π\) 的第 \(n\) 位,而不用计算前 \(n-1\) 项。
十六进制表示下,问你 \(π\) 的小数点后的第 \(n\) 位是多少 $ (1 ≤ n ≤ 100000)$ 。
Paper链接:
BBP Paper http://www.experimentalmath.info/bbp-codes/bbp-alg.pdf
简要题解:
其实看上面的 \(Paper\) 就知道怎么做了。
我简单解析一下。
把公式的第一项拿出来分析:
\(\sum_{k=0}^{\infty}\frac{1}{16^k}(\frac{4}{8k+1}) = \sum_{k=0}^{\infty}(\frac{4}{16^k(8k+1)}) = \sum_{k=0}^{\infty}(\frac{1}{16^k(8k+1)})\).
把公式拆分:
\(\sum_{k=0}^{\infty}(\frac{1}{16^k(8k+1)}) = \sum_{k=0}^{\infty}(\frac{1}{16^k(8k+1)}) = \sum_{k=0}^{n}(\frac{1}{16^k(8k+1)}) + \sum_{k=n + 1}^{\infty}(\frac{1}{16^k(8k+1)})\).
拆分之后我们就可以得到第 \(n\) 位。
将式子乘上 \(16^n\) ,使得小数点往后移动 \(n\) 位。
\([\sum_{k=0}^{n}(\frac{1}{16^k(8k+1)}) + \sum_{k=n + 1}^{\infty}(\frac{1}{16^k(8k+1)})]*16^{n}==> \sum_{k=0}^{n}(\frac{16^{n-k}}{(8k+1)}) + \sum_{k=n + 1}^{\infty}(\frac{16^{n-k}}{(8k+1)})\).
前一项 \(\sum_{k=0}^{n}(\frac{16^{n-k}}{(8k+1)})\) 为了避免高精度,可以化成 \(\sum_{k=0}^{n}(\frac{16^{n-k} mod (8k+1)}{(8k+1)})\).
后一项 \(\sum_{k=n + 1}^{\infty}(\frac{16^{n-k}}{(8k+1)})\) 就不用简化了,将 \(\infty\) 取够一定范围就可以了。
令 \(S_1 = \sum_{k=0}^{n}(\frac{16^{n-k}}{(8k+1)}) + \sum_{k=n + 1}^{\infty}(\frac{16^{n-k}}{(8k+1)})\).
那么,答案就是 \(4S_1 - 2S_2 - S_3 - S_4\)的小数部分。因为得到的只是小数部分,所以再乘以 \(16\) 后,得到的整数部分转化成十六进制就可以啦。
时间复杂度:\(O(nlogn)\)
所以,我是不是可以出一道关于计算二进制表示下的 \(log2\) 的题 ???
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
char print(int x)
{
if(x>=0 && x<=9)return x + '0';
return x+55;
}
ll qpower(ll a, ll b, ll mod)
{
ll res = 1;
while(b)
{
if(b & 1) res = a * res % mod;
b >>= 1;
a = a * a % mod;
}
return res;
}
double bbp(int n,ll k,ll b)
{
double res = 0;
for(int i=0;i<=n;i++)
{
res += (qpower(16,n-i,8*i+b) * 1.0/(8*i+b));
}
for(int i = n + 1;i <= n + 1000 + 1;i++)
{
res += (powf(16,n-i)* 1.0/(8*i+b));
}
return k * res;
}
int main()
{
int t,n;
cin>>t;
int cas = 1;
while(t--)
{
double ans = 0;
cin>>n;
n--;
ans = bbp(n,4,1) - bbp(n,2,4) - bbp(n,1,5) - bbp(n,1,6);
// cout<<"ans="<<ans<<endl;
ans = ans - (int)ans;
if(ans<0)ans+=1;
ans*=16;
char c ;
c = print(ans);
printf("Case #%d: %d %c\n",cas++,n+1,c);
}
return 0;
}
HDU 6217 BBP Formula (数学)的更多相关文章
- hdu 6217 A BBP Formula 公式题
题意 已知公式:$\pi=\sum_{k=0}^{\infty}\left[\frac{1}{16^{k}}\left(\frac{4}{8 k+1}-\frac{2}{8 k+4}-\frac{1} ...
- HDU 4342History repeat itself 数学
C - History repeat itself Time Limit:1000MS Memory Limit:32768KB Description Tom took the D ...
- HDU 4816 Bathysphere(数学)(2013 Asia Regional Changchun)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4816 Problem Description The Bathysphere is a spheric ...
- HDU 5584 LCM Walk 数学
LCM Walk Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=5584 ...
- HDU 4336 Card Collector 数学期望(容斥原理)
题目地址: http://acm.hdu.edu.cn/showproblem.php?pid=4336 题意简单,直接用容斥原理即可 AC代码: #include <iostream> ...
- HDU 5570 balls 期望 数学
balls Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=5570 De ...
- hdu 4710 Balls Rearrangement (数学思维)
意甲冠军:那是, 从数0-n小球进入相应的i%a箱号.然后买一个新的盒子. 今天的总合伙人b一个盒子,Bob试图把球i%b箱号. 求复位的最小成本. 每次移动的花费为y - x ,即移动前后盒子编号 ...
- HDU 4790 Just Random 数学
链接:pid=4790">http://acm.hdu.edu.cn/showproblem.php?pid=4790 意:从[a.b]中随机找出一个数字x,从[c.d]中随机找出一个 ...
- HDU 1018-Big Number(数学)
Big Number Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) Total ...
随机推荐
- Python中import和from的一些事。。。
摘自python学习手册, 用于记录. 客户端可以执行import或from语句.如果模块还没有加载,这两个语句会去搜索.编译以及执行模块文件程序.主要差别在于,import会读取整个模块,所以必须进 ...
- oracle查询字段大于指定长度的数据
select * from MES_MACHINE_RECORD t where length(t.bar_code2)<10 ;
- css————获取样式的各种方法
元素.style.样式:只能获取行间样式,css中的样式不能获取,且获得的内容是字符串. 元素.offsetWidth:可以获取无论css样式还是行间样式,但只能获得width,height,pare ...
- CODEVS——T1052 地鼠游戏
http://codevs.cn/problem/1052/ 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题解 查看运行结果 题目描述 D ...
- c++中六种构造函数的实现以及9中情况下,构造函数的调用过程
六种构造函数的实现代码例如以下: #include<iostream> using namespace std; //c++中六种默认的构造函数 class Test { public: ...
- 论Nim中的 proc 和 method
在Nim中.proc 是定义过程的keyword.method 是定义方法的keyword.它们之间根本的差别是proc定义的过程是静态绑定.method定义的方法是动态绑定.谈到静态绑定.动态绑定又 ...
- crm2011处理save事件和获取当前窗口信息
//防止保存记录,eContext:当前上下文对象 function My_PreventSaveFunction(eContext) { eContext.getEventArgs().pr ...
- android framework 02
Android底层开发1.安装Ubuntu系统2.Ubuntu配置开发环境: sudo apt-get install git-core gnupg flex bison gperf zip sudo ...
- WHU 1470 Join in tasks 水题
http://acm.whu.edu.cn/land/problem/detail?problem_id=1470 大概是给你一个队列,每次移动队头的数到队尾并减1,如果本身这个数为1就删去. 然后a ...
- watch---周期性的方式执行给定的指令
watch命令以周期性的方式执行给定的指令,指令输出以全屏方式显示. 选项 -n:指定指令执行的间隔时间(秒): -d:高亮显示指令输出信息不同之处: -t:不显示标题.