使用BBP公式计算Pi的代码

php中文网

发布时间：2016-07-25 09:01:23

2927人浏览过

来源于php中文网

原创

BBP公式号称可以直接获取Pi的第N位结果，本来我以为这种获取的速度在任意一个位置都是相同的呢，但是从bellard的主页（http://bellard.org/pi/）下载一个改进的C程序测试后发现N越大，计算时间越长，我转换成php代码后，计算第1K位时会超过1分钟的时间，原C程序计算的时候也是很慢

这个程序一次只会输出9个数字（囧rz，不知道怎样改能多输出几位）

/**
* 圆周率计算(BBP)
* @author Moyo
* @url http://moyo.uuland.org/code/php-pi-calc/
* @version 1.0
* @date 2013.01.12
*/
class pi
{
public static function calc($__N__)
{
$n = (int)$__N__;
$av = $a = $vmax = $N = $num = $den = $k = $kq = $kq2 = $t = $v = $s = $i = 0;
$sum = 0.0;
$N = (int)(($n + 20) * log(10) / log(2));
$sum = 0;
for ($a = 3; $a {
$vmax = (int)(log(2 * $N) / log($a));
$av = 1;
for ($i = 0; $i {
$av = ($av * $a);
}
$s = 0;
$num = 1;
$den = 1;
$v = 0;
$kq = 1;
$kq2 = 1;
for ($k = 1; $k {
$t = $k;
if ($kq >= $a)
{
do
{
$t = (int)($t / $a);
$v --;
}
while (($t % $a) == 0);
$kq = 0;
}
$kq ++;
$num = self::mul_mod($num, $t, $av);
$t = (2 * $k -1);
if ($kq2 >= $a)
{
if ($kq2 == $a)
{
do
{
$t = (int)($t / $a);
$v ++;
}
while (($t % $a) == 0);
}
$kq2 -= $a;
}
$den = self::mul_mod($den, $t, $av);
$kq2 += 2;
if ($v > 0)
{
$t = self::inv_mod($den, $av);
$t = self::mul_mod($t, $num, $av);
$t = self::mul_mod($t, $k, $av);
for ($i = $v; $i {
$t = self::mul_mod($t, $a, $av);
}
$s += $t;
if ($s >= $av)
{
$s -= $av;
}
}
}
$t = self::pow_mod(10, ($n - 1), $av);
$s = self::mul_mod($s, $t, $av);
$sum = (double)fmod((double)$sum + (double)$s / (double)$av, 1.0);
}
return array(
'n' => $n,
'v' => sprintf('%09d', (int)($sum * 1e9))
);
}
private static function next_prime($n)
{
do
{
$n ++;
}
while (!self::is_prime($n));
return $n;
}
private static function is_prime($n)
{
$r = $i = 0;
if (($n % 2) == 0)
{
return 0;
}
$r = (int)(sqrt($n));
for ($i = 3; $i {
if (($n % $i) == 0)
{
return 0;
}
}
return 1;
}
private static function mul_mod($a, $b, $m)
{
return fmod((double)$a * (double)$b, $m);
}
private static function inv_mod($x, $y)
{
$q = $u = $v = $a = $c = $t = 0;
$u = $x;
$v = $y;
$c = 1;
$a = 0;
do
{
$q = (int)($v / $u);
$t = $c;
$c = $a - $q * $c;
$a = $t;
$t = $u;
$u = $v - $q * $u;
$v = $t;
}
while ($u != 0);
$a = $a % $y;
if ($a {
$a = $y + $a;
}
return $a;
}
private static function pow_mod($a, $b, $m)
{
$r = $aa = 0;
$r = 1;
$aa = $a;
while (1)
{
if ($b & 1)
{
$r = self::mul_mod($r, $aa, $m);
}
$b = $b >> 1;
if ($b == 0)
{
break;
}
$aa = self::mul_mod($aa, $aa, $m);
}
return $r;
}
}
?>

复制代码

PHP中日期范围交集计算与优化实践

怎么用php代码解密_用PHP代码实现多场景解密教程【技巧】

WordPress开发：在页脚动态插入短代码并实现样式控制

PHP Datepicker实现年龄验证：确保用户年龄不低于18岁

PHP 日期处理：高效获取月份内所有周的范围 (Carbon 实践)

本站声明：本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系admin@php.cn

上一篇：php imagick 将图像色彩CMKY转RGB 下一篇：php 获取中文字符串首字母

热门AI工具

DeepSeek

幻方量化公司旗下的开源大模型平台

AI大模型

开放平台

豆包大模型

字节跳动自主研发的一系列大型语言模型

AI大模型

通义千问

阿里巴巴推出的全能AI助手

AI大模型

腾讯元宝

腾讯混元平台推出的AI助手

文档处理

Excel 表格

文心一言

文心一言是百度开发的AI聊天机器人，通过对话可以生成各种形式的内容。

AI大模型

中文写作

讯飞写作

基于讯飞星火大模型的AI写作工具，可以快速生成新闻稿件、品宣文案、工作总结、心得体会等各种文文稿

中文写作

写作工具

即梦AI

一站式AI创作平台，免费AI图片和视频生成。

图片拼接

图画生成

ChatGPT

最最强大的AI聊天机器人程序，ChatGPT不单是聊天机器人，还能进行撰写邮件、视频脚本、文案、翻译、代码等任务。

AI大模型

中文写作

智谱清言 - 免费全能的AI助手

AI大模型

PDF 文档

相关专题

excel制作动态图表教程

本专题整合了excel制作动态图表相关教程，阅读专题下面的文章了解更多详细教程。

2025.12.29

freeok看剧入口合集

本专题整合了freeok看剧入口网址，阅读下面的文章了解更多网址。

2025.12.29

俄罗斯搜索引擎Yandex最新官方入口网址

Yandex官方入口网址是https://yandex.com；用户可通过网页端直连或移动端浏览器直接访问，无需登录即可使用搜索、图片、新闻、地图等全部基础功能，并支持多语种检索与静态资源精准筛选。本专题为大家提供相关的文章、下载、课程内容，供大家免费下载体验。

197

2025.12.29

python中def的用法大全

def关键字用于在Python中定义函数。其基本语法包括函数名、参数列表、文档字符串和返回值。使用def可以定义无参数、单参数、多参数、默认参数和可变参数的函数。本专题为大家提供相关的文章、下载、课程内容，供大家免费下载体验。

2025.12.29

python改成中文版教程大全

Python界面可通过以下方法改为中文版：修改系统语言环境：更改系统语言为“中文（简体）”。使用 IDE 修改：在 PyCharm 等 IDE 中更改语言设置为“中文”。使用 IDLE 修改：在 IDLE 中修改语言为“Chinese”。本专题为大家提供相关的文章、下载、课程内容，供大家免费下载体验。

2025.12.29

C++的Top K问题怎么解决

TopK问题可通过优先队列、partial_sort和nth_element解决：优先队列维护大小为K的堆，适合流式数据；partial_sort对前K个元素排序，适用于需有序结果且K较小的场景；nth_element基于快速选择，平均时间复杂度O(n)，效率最高但不保证前K内部有序。本专题为大家提供相关的文章、下载、课程内容，供大家免费下载体验。

2025.12.29