阶乘相关的算法及其C++实现

       有关阶乘的算法,不外乎两个方面:一是高精度计算;二是与数论相关。

  一、 高精度计算阶乘

  这实际上是最没有技术含量的问题,但是又会经常用到,所以还是得编写,优化它的计算。

  首先看小于等于12的阶乘计算(计算结果不会超出32位范围):

  int factorial(int n) {

  if (n == 1 || n == 0)

  return 1;

  return factorial(n-1)*n;

  }

  这个递归程序简单明了,非常直观,然而一旦n > 12,则超过32位int型的范围出现错误结果,所以上面这个递归程序仅适合n <= 12的阶乘计算,为了计算较大n的阶乘,需要将高精度乘法算法纳入到阶乘计算中来,高精度乘法过程可以如下简单的描述:(其中A * B = C,A[0], B[0], C[0]分别存储长度)

  for (i = 1; i <= A[0]; i++)

  for (j = 1; j <= B[0]; j++) {

  C[i+j-1] += A[i]*B[j];        // 当前i+j-1位对应项 + A[i] * B[j]

  C[i+j] += C[i+j-1]/10;       // 它的后一位 + 它的商(进位)

  C[i+j-1] %= 10;                // 它再取余即可

  }

  C[0] = A[0] + B[0];

  while (C[0] > 1 && C[C[0]] == 0) C[0]--;   // 去头0,获得实际C的长度

  有了这个高精度乘法之后,计算阶乘就可以简单的迭代进行:

  for (i = 2; i <= n; i++) {

  将i转换成字符数组;

  执行高精度乘法:将上一次结果乘上i

  }

  二、 与数论有关

  由于阶乘到后面越来越大,巧妙的利用数论求得一些有趣的数字(数值)等成为阶乘算法的设计点,下面给出几道相关的问题与分析:

  (1)   计算阶乘末尾第一个非0数字:

  这是一个比较经典的问题,比较复杂的算法是利用一个艰难的数学公式,可惜我不会,从网上的资料学习中,整理出下面这个简单易懂的算法:

  观察n!,可以发现在乘的过程中,对于任意 n > 1,n!的末尾第一个非0数字都是偶数。我们只需保留最后一位非零数。当要乘的数中含有因数5时,我们可以把所有的因数5都当作8来乘。这是因为:

  …x2*5=…10(舍)或…60,最后一位非零数为6。而恰好2*8=16,末位为6。

  …x4*5=…70(舍)或…20,最后一位非零数为2。而恰好4*8=32,末位为2。

  …x6*5=…30(舍)或…80,最后一位非零数为8。而恰好6*8=48,末位为8。

  …x8*5=…90(舍)或…40,最后一位非零数为4。而恰好8*8=64,末位为4。

  (对于n > 1时,最后一位不会出现 1, 7, 3, 9,而永远是2, 4, 6, 8的循环出现)

  因此,在迭代作乘法时,主要就是计算因子5的数量,同时可见因子5的个数以4为循环节(即只需要取它的数量对4取模)。那么对于不同情况下的因子5的数量,可以通过res[5][4] = {{0,0,0,0}, {2,6,8,4}, {4,2,6,8}, {6,8,4,2}, {8,4,2,6}}来得到,使用nonzero[i]表示i的阶乘的最后一位,那么:

  如果t是偶数,则直接乘:nonzero[i] = (nonzero[i-1]*t)%10。

  否则nonzero[i] = res[((nonzero[i-1]*t)%10)/2][five];

  其中t是除掉所有因子5的结果,five为因子5数量对4的模。

  (2)。 阶乘末尾有多少个0

  分析发现,实际上形成末尾0,就是因子5的数量,而计算1~n之间包含一个因子i的个数的简单算法就是:

  cnt = 0; while (n) { n /= i; cnt += n; }

  因此,直接将i换成5,就可以得到因子5的数量,也即n!末尾0的数量。

  (3)。 返回阶乘左边的第二个数字

  简单算法:用实数乘,超过100就除以10,最后取个位即可。因为整数部分的个位就是阶乘结果左边的第二个数字。相关题目:

  (4)。 判断数值 m 是否可以整除 n!

  算法:使用素因子判断法

  A. 首先直接输出两种特殊情况:

  m == 0 则 0肯定不会整除n!;

  n >= m 则 m肯定可以整除n!;

  B. 那么就只剩最后一种情况:m > n,我们从m的最小素因子取起,设素因子为i那么可以   求得m的素因子i的个数 nums1;再检查闭区间 i ~ n 之间的数,一共包含多少个素因子i,就可以简单的利用上面(2)中所介绍的数学公式进行计算得到nums2。如果nums2 < nums1,就表示1 ~ n中包含素因子的数量 < 除数m包含素因子i的数量,那么m必然不能整除n!,置ok = false。

  C. 最后:如果 !ok or m > n or m == 0 则不能整除;否则可以整除

  (5)。数字N能否表示成若干个不相同的阶乘的和:

  这里可以选择的阶乘为:0! ~ 9!,实际上这一题与数论无关,与搜索有关。

  分析,由于可供选择的阶乘数量较少,直接可以利用DFS搜索来做:

  A. 首先将0 ~ 9的阶乘作一个表A[10];再设置一个可以组成“和”的数组ans[N]。

  B. 深度优先搜索方法:

  search(n) {

  for(i = n; i <= 9; i++) {

  sum += A[i];      //求和

  如果sum在ans数组中不存在,则将sum插入到ans[]数组中

  search(n+1);

  sum -= A[i];       //回溯

  }

  }

  C. 最后对于输入n,就在ans数组中查找是否存在n,如果存在,则表示n可以表示成不同的阶乘和,否则不行。

时间: 2024-10-27 03:51:59

阶乘相关的算法及其C++实现的相关文章

精度计算-大数阶乘-ACM常用算法

语法  int result=factorial(int n); 参数: n: n 的阶乘 返回值: 阶乘结果的位数 注意:     本程序直接输出n!的结果,需要返回结果请保留long a[]  inlclude <math.h> 源程序:   int factorial(int n) { long a[10000]; int i,j,l,c,m=0,w; a[0]=1; for(i=1;i<=n;i++) { c=0; for(j=0;j<=m;j++) { a[j]=a[j]

求随机数相关的算法或思路。

问题描述 小弟的问题是可以描述为:现在有编号为1-20的20类不同的苹果,每次从这20类只取一个,(注意每种类型苹果个数不限)允许类型重复,现在取上10次,要求取到的苹果的所有类型为一确定数,如8种,或者在一个区间上,如为5-8种.注意,是在取了10次的基础上.然后还有一个随机数问题:在上面的基础上,我现在已经有了10个苹果,他们的编号允许存在重复.现在分3次拿走完,每次拿若干个,怎样确定拿走的的苹果的编号? 解决方案 解决方案二:importjava.util.*;publicclassApp

HDOJ(HDU) 2212 DFS(阶乘相关、)

Problem Description A DFS(digital factorial sum) number is found by summing the factorial of every digit of a positive integer. For example ,consider the positive integer 145 = 1!+4!+5!, so it's a DFS number. Now you should find out all the DFS numbe

图像分割相关算法

前言 图像分割(Image Segmentation)指的是将数字图像细分为多个图像子区域的过程,在OpenCv中实现 了三种跟图像分割相关的算法,它们分别是:分水岭分割算法.金字塔分割算法 以及均值漂移分割算法.它们的使用过程都很简单,下面的文章权且用于记录, 并使该系列保持完整吧. 分水岭分割算法 分水岭分割算法需要您或者先前算法提供标记,该标记用于指定哪些大致区 域是目标,哪些大致区域是背景等等:分水岭分割算法的分割效果严重依赖于提 供的标记.OpenCv中的函数cvWatershed实现

浅谈算法和数据结构 十二 无向图相关算法基础

从这篇文章开始介绍图相关的算法,这也是Algorithms在线课程第二部分的第一次课程笔记. 图的应用很广泛,也有很多非常有用的算法,当然也有很多待解决的问题,根据性质,图可以分为无向图和有向图.本文先介绍无向图,后文再介绍有向图. 之所以要研究图,是因为图在生活中应用比较广泛: 无向图 图是若干个顶点(Vertices)和边(Edges)相互连接组成的.边仅由两个顶点连接,并且没有方向的图称为无向图. 在研究图之前,有一些定义需要明确,下图中表示了图的一些基本属性的含义,这里就不多说明. 图的

排序算法的相关问题。

问题描述 排序算法的相关问题. while 循环条件里的 i > 0 有什么意义吗 解决方案 保证数组下标合法... 解决方案二: 快速排序算法相关排序算法的相关介绍 解决方案三: 快速排序算法相关排序算法的相关介绍 解决方案四: 循环的结束条件,i-- 当i<=0时结束循环 解决方案五: 循环的结束条件,i-- 当i<=0时结束循环 解决方案六: i>0: 比如这个数组2345,要插入的数是1,那么如果不加i>0会导致程序会无意义的再进入循环一次,我觉的这么做,除了坐标合法

千万不要被阶乘吓倒_C 语言

阶乘(Factorial)是个很有意思的函数,但是不少人都比较怕它,我们来看看两个与阶乘相关的问题: 1. 给定一个整数N,那么N的阶乘N!末尾有多少个0呢?例如:N=10,N!=3 628 800,N!的末尾有两个0.2.求N!的二进制表示中最低位1的位置. 有些人碰到这样的题目会想:是不是要完整计算出N!的值?如果溢出怎么办?事实上,如果我们从"哪些数相乘能得到10"这个角度来考虑,问题就变得简单了.首先考虑,如果N!= K×10^M,且K不能被10整除,那么N!末尾有M个0.再考

消息摘要算法-SHA算法实现

简述 SHA算法分为一代和二代,这个是美国安全局发布的一系列的密码散列算法. 基于如下两点安全理由: 1.从消息摘要反推输入信息,从计算理论上来说很困难 2.想找到两组不同的消息对应到相应的消息摘要,从计算理论上来说一额很困难.对任何消息的变动,都有很高几率产生消息摘要迥异   不过这个算法已经不安全了,被中国的教授破解了.这个技术通常用来做数字签名.常见的就是B2B,B2C网站的支付.数据post到支付中心的时候,是发送一些加密串,同时带上签名.而这个签名可以被伪造.所以就不安全了.   常见

SQL Server 中的时间算法

DECLARE @Date DATETIME SET @Date=GETDATE() --前一天,给定日期的前一天 SELECT DATEADD(DAY,-1,@Date) AS '前一天' --后一天,给定日期的后一天 SELECT DATEADD(DAY,1,@Date) AS '后一天' GO --月初,计算给定日期所在月的第一天 --这个计算的技巧是先计算当前日期到"1900-01-01"的时间间隔数,然后把它加到"1900-01-01"上来获得特殊的日期,