热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

php如何实现n的阶乘

这篇文章主要介绍php如何实现n的阶乘,文中介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们一定要看完!php实现n的阶乘的方法:1、通过普

这篇文章主要介绍php如何实现n的阶乘,文中介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们一定要看完!

php实现n的阶乘的方法:1、通过普通递归实现,代码如“function fact(int $n): int{...}”;2、通过普通循环实现,代码如“while ($num <= $n) {$result = $result...}”。

本文操作环境:Windows10系统、PHP 7.2.15版,DELL G3电脑

php怎么实现n的阶乘?

据本人了解,阶乘的实现方法一般可以分为三种,通常意义下的递归和循环各算一种,还有一大类通过一些巧妙的数学方法减少运算次数(尤其是乘法运算次数),进而优化计算效率。

如果要考虑到高精度、大整数的阶乘,对于 PHP 语言而言,情况会更复杂一些,比如使用 BCMath 扩展提供的一些方法时,显式的数字与字符串转换操作比较频繁。

本文在只考虑 n 为整数的情况下,分别尝试实现上述的几种情况,每种情况给出可用的代码示例,并在文末附上几种方法的综合对比情况。

普通递归实现

首先是普通递归实现,根据递归的通用公式 fact(n) = n * fact(n-1) 很容易写出阶乘的计算代码。普通递归实现的优点在于代码比较简洁,和通用公式一样的过程使得代码容易理解。缺点则在于由于需要频繁地调用自身,需要大量的入栈出栈操作,整体的计算效率不高(见文末表格)。

function fact(int $n): int
{
    if ($n == 0) {
        return 1;
    }
    return $n * fact($n - 1);
}

普通循环实现

普通循环实现有些「动态规划」的味道,但由于中间态变量使用频率低,不需要额外存储空间,所以要比一般的动态规划算法简单。普通递归方法是自顶向下(由 n 到 1)的计算过程,而普通循环是自底向上进行计算。

因此相对而言,代码没有上述方法直观,但由于少了频繁的入栈出栈过程,计算效率会高一些(见文末表格)。

function fact(int $n): int
{
    $result = 1;
    $num = 1;
    while ($num <= $n) {
        $result = $result * $num;
        $num = $num + 1;
    }
    return $result;
}

自行实现的大整数阶乘

由于 PHP 中 int 类型的范围限制,上述两种方法最多只能精确计算到 20 的阶乘。如果只是考虑到 20 的阶乘的情况,那么用查表法实现会更快:事先计算好 0-20 的阶乘并存储到一个数组中,需要用时查询一次便可。

为了能够适应大数的阶乘,得到精确的计算结果,本文基于「普通循环方法」进行改进,使用数组存储计算结果中的每一位(由低到高位),通过相乘进位的方式依次计算每一位的结果。

不言而喻,本方法的优点在于可以适用于高精度的大数阶乘场合,缺点就是对于小数阶乘而言,计算过程复杂且速度慢。

function fact(int $n): array
{
    $result = [1];
    $num = 1;
    while ($num <= $n) {
        $carry = 0;
        for ($index = 0; $index < count($result); $index++) {
            $tmp = $result[$index] * $num + $carry;
            $result[$index] = $tmp % 10;
            $carry = floor($tmp / 10);
        }
        while ($carry > 0) {
            $result[] = $carry % 10;
            $carry = floor($carry / 10);
        }
        $num = $num + 1;
    }
    return $result;
}

BCMath 扩展方法

BCMath 是 PHP 的一个数学扩展,用于处理字符串表示的数字(任意大小和精度)的数值计算。由于是使用 C 语言实现的扩展,计算速度会比上述自行实现的快。

在本人的笔记本上,同样是计算 2000 的阶乘,自行实现的需要平均 0.5-0.6 秒,使用 BCMath 耗时 0.18-0.19 秒。该方法的缺点主要在于需要安装相应的扩展,属于非代码层面的改动,对于环境管理升级不便的应用而言,可实践性有待商榷。

function fact(int $n): string
{
    $result = &#39;1&#39;;
    $num = &#39;1&#39;;
    while ($num <= $n) {
        $result = bcmul($result, $num);
        $num = bcadd($num, &#39;1&#39;);
    }
    return $result;
}

优化算法

在本文开头有提到,优化算法尝试尽可能地减少运算次数(尤其是乘法的运算次数)来实现快速阶乘。考虑到对于小整数阶乘而言,最快的算法应该是查表法,时间复杂度为 O(1),所以本小节主要针对大整数的精确阶乘进行讨论和测试。

据了解,目前阶乘优化比较常见的是通过 n! = C(n, n/2) * (n/2)! * (n/2)! 式子进行复杂度优化,而该式子中的亮点主要在于 C(n, n/2) 的优化。考虑到大整数情况下,PHP 语言实现 C(n, n/2) 的效率不高,而且实现的代码可读性比较差(频繁的数字与字符串的显式转换),所以本文用的是另外一种比较巧妙的方法。

乘法的计算速度通常要低于加减法运算,通过减少乘法的运算次数可以提高整体运算速度。通过数学归纳可以发现,对于 n 的阶乘,可以依次求出比 (n/2)^2 小 1、1+3、1+3+5... 的数值,再依次相乘得到目标值。

该算法的优点在于计算速度较快,而缺点就是实现过程不直观、不易理解。经测试,以下代码计算 2000 的阶乘平均时间为 0.11 秒,大约是普通循环方法的一半耗时。

除了这种方法优化,也有看到其它的类似的思路,比如对 1...n 中的数反复检验是否被 2 整除,记录下被 2 整除的次数 x,并尝试归纳出共同的奇数相乘式,最后乘以 2^x 得到结果。

function fact(int $n): string
{
    $middleSquare = pow(floor($n / 2), 2);
    $result = $n & 1 == 1 ? 2 * $middleSquare * $n : 2 * $middleSquare;
    $result = (string)$result;
    for ($num = 1; $num < $n - 2; $num = $num + 2) {
        $middleSquare = $middleSquare - $num;
        $result = bcmul($result, (string)$middleSquare);
    }
    return $result;
}

综合对比

本文中提到的方法是按照由劣到优的顺序,因此,下列表格中每一行中提到优劣势,主要是和其上一两种方法对比。

表格中「测试耗时」一列的测试环境为个人笔记本,硬件配置为 Dell/i5-8250U/16GB RAM/256GB SSD Disk,软件配置为 Win 10/PHP 7.2.15。

php如何实现n的阶乘

以上是“php如何实现n的阶乘”这篇文章的所有内容,感谢各位的阅读!希望分享的内容对大家有帮助,更多相关知识,欢迎关注编程笔记行业资讯频道!


推荐阅读
  • 简单数据结构模板
    一,STL1&amp;gt;STL中数据结构常见操作a.assign(b.begin(),b.begin()+3);b为向量,将b的0~2个元素构成的向量赋给aa.as ... [详细]
  • C语言解题——指针解析(牛客网题目)
    指针,是C语言中的一个重要概念及其特点,也是掌握C语言比较困难的部分。指针也就是内存地址,指针变量是用来存放内存地址的变量,在同一CPU构架下,不同类型的指针变量所占用的存储单元长度是相同的,而 ... [详细]
  • 新手指南何为PHP语言_PHP教程:PHP语言经过长时间的发展,很多用户都很了解PHP语言了,这里我发表一下个人理解,和大家讨论讨论。比较关注网络的朋友,特别是关注网站建设技术的朋 ... [详细]
  • 人脸识别中的损失函数
    本文主要是针对人脸识别中的各种loss进行总结。 背景对于分类问题,我们常用的lossfunction是softmax,表示为:,当然有softmax肯定也有hardmax:,so ... [详细]
  • MFC中的文件读写操作
    由于要写个游戏工具方便编程,要用到文件类的操作,看到这篇文章不错,特转载计算机室如何管理自身所存放着的大量的信息的呢?win ... [详细]
  • CTF编程入门第一题 循环
    循环分值:10循环分值:10来源:北邮天枢战队难度:易参与人数:1478人GetFlag:467人答题人数:523人解题通过率:89%来源:北邮天枢战队难度:易参与人数:1478人 ... [详细]
  • 不懂性能优化,再强的计算机也白玩
    不懂性能优化,再强的计算机也白玩-Python的优秀有目共睹,不过说的性能,还真比不了Java、C、Go,有没有提升性能的技巧或方法呢?今天我们一起学习下提升Python性能的方式 ... [详细]
  • NPDP第五章 工具与度量
    什么是创意开发阶段作用:早起产生初步的产品概念;中期用于解决实施问题;后期用于规划上市。创意开发是创造、发展和传达新创意的创造性过程&# ... [详细]
  • 目录一.python循环1.wuile与else连用2.死循环3while的镶嵌二.for循环1.fuo和whlie循环2.循环字符串3.循环字符4.循环元祖5.循环元组6.循环集 ... [详细]
  • 设计模式——行为型模式(二) ... [详细]
  • 在实际开发中,现在安卓端和后台之间的数据交互,一般都是用JSON来传递数据信息。JSON大家一般都比较熟悉。我这边就以实际项目中的后台传过来的情况和大家分析下及如何处理。比如后台返 ... [详细]
  • 最近想用js做一个简单的计算器,不过网上的例子好像大部分都是直接从左到右挨个计算,就好像1+2*5,就会先计算1+2,再计算3*5,并没有实现运算符的优先级,这里找到了一种方法实现,来总结一下。不过这 ... [详细]
  • 九宫格计算. ... [详细]
  • 前端微服务二
    为了解决庞大的一整块后端服务带来的变更与扩展方面的限制,出现了微服务架构(Microservices):微服务是面向服务架构(SOA)的一种变体,把应用程序设计成一系列松耦合的细粒 ... [详细]
  • 导读:今天编程笔记来给各位分享关于php动态扩展怎么加载的相关内容,如果能碰巧解决你现在面临的问题,别忘了关注本站,现在开始吧!本文目录一览: ... [详细]
author-avatar
追梦的青春灬_176
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有