N!(N的阶乘)最末位非0的求解方法

问题是求关于N!的最后一位非0位, 如3!=6,最后一位非0位为6, 5!=120, 最后一位非0位为2.怎么样快速的求出最后一位非0位呢?

最朴素的想法就是先求出N!的结果,再求出结果的最后一位非0位.当N比较小时,是可以承受的,但是当N达到一定规模的时候,时间,空间都不会太理想.这里需要一些技巧.既然是求最后一位非0位,我们就可以先除去所有对结果没有影响的数,如10的倍数.于是先把N!因子分解得到形如2^a*5^b*c.这个时候我们去掉一个b个5因子和b个2因子,最后一位非0位是不变的.(N!中2的因子一定不会比5的因子少).

于是我们的要求的结果就变为(2^(a-b)*c)%10.由(a*b)%10=((a%10)*(b%10))%10我们可以得((2^(a-b)%10)*(c%10))%10,由于c不会产生未位为0,故只保留c的最未位即可.于是可将c转化为1,3,7,9因子的相乘得到的结果的最未位(因为1,3,7,9因子相乘不会产生最未非0位,故去掉高位不会对结果产生影响,同时1*n=n可以去掉1的因子).

2,3,7,9因子规律如下:

2^1=2, 2^2=4, 2^3=8, 2^4=16->(6), 2^5=32->(2)

3^0=1, 3^1=3, 3^2=9, 3^3=27->(7), 3^4=81->(1)

7^0=1, 7^1=7, 7^2=49->(9), 7^3=343->(3), 7^4=2401->(1)

9^0=1, 9^1=9, 9^2=81->(1), 9^3=729->(9), 9^4=6561->(1)

它们都是以4为循环周期的.于是我们只要求出2, 5, 3, 7, 9因子的个数即可.

首先我们求2,5因子在N!中的个数.2的因子的每个偶数到少有1个,同时将数列中每个数/2,其中的偶数还有一个2因子.直至n=1或n=0结束.5因子求法相同.代码如下:

int getFactor_2_5(int n, int f){
    int ret=0;
    while(n>0){
        ret+=n/f;
        n/=f;
    }
    return ret;
}

3,7,9因子的个数有多少呢?对于1,2,3,4......n-1,n来说,未尾以3,7,9结束的数的个数为

n/10+(n%103f?1:0),(f=3,7, 9).

同时我们对于对于奇数数列/5可以得到一个新的数列也有3,7,9因子,对于偶数数列/2也可以得到新的数列也有3,7,9的因子,将所有的3,7,9因子相加即可得到总的3,7,9因子的个数.得到3,7,9因子的个数后,我们可以将其全部转化为因子3的个数.因为9=3*3(3^2), 7=(3*3*3(3^3))%10,设f3, f7, f9为3, 7, 9因子的个数,全部转化为因子3的个数为f3+2*f9+3*f7.

于是我们可以用递归同时求2,3,5,7,9因子的个数,代码如下:

void getFactor(int n){
    if(n==0)
        return;
    for(int m=n; m>0; m/=5){
        int t=m/10, r=m%10;
        f3+=t+(r>=3);
        f5+=t+(r>=5);
        f7+=t+(r>=7);
        f9+=t+(r>=9);
    }
    f2+=n/2;
    getFactor(n/2);
}

可用两个数组表示循环:

  1. int p2[4]={6, 2, 4, 8};
  2. int p3[4]={1, 3, 9, 7};

故结果为(1):当2, 5因子个数相同时,只与3因子相关,结果为p3[f3%4]%10;

(2):当2因子大于5因子时,结果同时与3因子和2因子相关,为(p2[f2%4]*p3[f3%4])%10.

其实通过N!的最未非0位的方法我们可以求排列组合数NPM,C(N,M)的最未非0位,用上面的各因子个数减去下面的各因子个数就是结果的各因子个数.只是此时需要注意的是5的因子可能会比2的因子多.当5的因子比2的因子多时,未位一定为5.其余情况与上面相同.

时间: 2024-11-09 04:51:13

N!(N的阶乘)最末位非0的求解方法的相关文章

算法-计算阶乘n!末尾0的个数

算法逻辑转载自计算阶乘n!末尾0的个数: 问题描述    给定参数n(n为正整数),请计算n的阶乘n!末尾所含有"0"的个数.    例如,5!=120,其末尾所含有的"0"的个数为1:10!= 3628800,其末尾所含有的"0"的个数为2:20!= 2432902008176640000,其末尾所含有的"0"的个数为4. 计算公式    这里先给出其计算公式,后面给出推导过程.    令f(x)表示正整数x末尾所含有的&q

NYOJ1026 阶乘末尾非0 【模板】

阶乘末尾非0 时间限制:2000 ms  |  内存限制:65535 KB 难度:3 描述 我们的问题很是简单,n!末尾非0数是几? 比如n=5的时候,n!=120,那么n!末尾非0数是2. 输入 多组数据, 每组数据占一行,每行一个整数0<=n<=10^1000 输出 n!末尾非0数. 样例输入 5 样例输出 2 直接用的网上的模板 /*==================================================*| 阶乘最后非零位,复杂度 O(nlogn) \*==

微信临时二维码时为32位非0整型

生成带参数的临时二维码时,需要指定场景值ID(scene_id), 由于项目中需要通过场景值来详细记录二维码的来源,故想对这个场景值id做出一个统一的生成机制. 微信文档中的表述如下: scene_id:场景值ID,临时二维码时为32位非0整型,永久二维码时最大值为100000(目前参数只支持1–100000). 然后团队中的人就以为这里的32位非0整型,是指长度为32位的数字就可以了. 结果当我们指定scene_id为 12345678901234567890123456789012 时, 我

计算阶乘n!末尾0的个数

一.问题描述 给定一个正整数n,请计算n的阶乘n!末尾所含有“0”的个数.例如: 5!=120,其末尾所含有的“0”的个数为1: 10!= 3628800,其末尾所含有的“0”的个数为2: 20!= 2432902008176640000,其末尾所含有的“0”的个数为4. 二.算法分析 此类问题很显然属于数学问题,一定要找到其中的本质规律才能得到正确的数学模型. 两个大数字相乘,都可以拆分成多个质数相乘,而质数相乘结果尾数为0的,只可能是2*5.如果想到了这一点,那么就可以进一步想到:两个数相乘

N阶乘尾部的0个数

N阶乘尾部的0个数 描述 设计一个算法,计算出n阶乘中尾部零的个数 思路: 1.1 * 2 * 3 * ... * n --> 1 * 2 * 3 * (2 * 2) * 5 * (2 * 3) * 7 (2 2 * 2) * (3 * 3) (2 5) * ...化成质数相乘,只有2 * 5才可能得到结果有尾数中有0 2.因为2的个数是比5多的,求0的个数问题就转化成了求5的个数的问题 3.5 * 5 * 5 * 5 * 5 * ... * 5有n个5 ; 得到有n个5:有几个 4.- ; -

Myeclipse 2015 stable 2.0 完美破解方法

2015-08-21 以前写了一篇<Myeclipse 2015 stable 1.0 完美破解方法>,现 在跟新一下Myeclipse 2015 stable 2.0 破解方法,此方法来源于贴吧的一篇文章:http://tieba.baidu.com/p/3809637196?pn=1 :在此先感谢一下作者的付出.下面是我在MAC下破解的图文过程(Linux和windows其实大同小异).破解方法和我写得<Myeclipse 2015 stable 1.0 完美破解方法>基本一样

nodejs+express4.0+mongodb安装方法 for Linux, Mac

废话不多说 1:下载nodejs包 下载地址如下:http://www.nodejs.org/download/ 下载source code版本需要解压后到其目录执行./configure,然后make && make install; 而binaries的tar.gz解压后即可. 2:下载mongodb 下载地址如下:http://www.mongodb.org/downloads 而binaries的tar.gz解压后即可. 3:配置环境变量 cd ~:此目录下有个.bash_prof

PHP Warning: Module &#39;modulename&#39; already loaded in Unknown on line 0 的解决方法

今天无间断服务加载php-fpm时,爆出了一个错误:PHP Warning:  Module 'xhprof' already loaded in Unknown on line 0 <br /> 从英文可以看出PHP的扩展xhprof已经加载了,所以就查看php.ini,搜索xhprof.so发现有两个,删除其中一个,然后无间断服务加载php-fpm就解决了该问题. PHP Warning: Module 'modulename' already loaded in Unknown on l

Axure7.0汉化方法

下载汉化包 AxureRP7CN_汉化包.rar 首先退出正在运行中的 Axure (如果您正在使用). 将 汉化包.rar 文件解压, 得到 lang 文件夹,  然后将其复制到 Axure 安装目录.然后重新打开Axure,就OK Axure7.0汉化方法,布布扣,bubuko.com