阶乘

转载请注明出处:http://blog.csdn.net/ns_code/article/details/28335353

前言

    主要看两道有关阶乘的题目,从中能够看出一些规律来。

题目一

    N!末尾0的个数

    找末尾0出现的个数,那我们就要找产生0的乘数,即哪些数相乘会得到10。我们须要对N!进行质因数分解,因为10 = 2*5,因此0的个数至于N!中2和5出现的的对数有关,而能被2整除的数出现的频率比能被5整除的数要多得多,因此我们找出N!中质因数5出现的个数,即为N!末尾0的个数。

    方法一

    最直接的方法,就是计算1-N中每一个数的因式分解中5的个数,然后相加,代码例如以下:

int FactorialNum0_1(int n)
{
	int count = 0;
	int i;
	for(i=1;i<=n;i++)
	{
		int j = i;
		while(j%5 == 0)
		{
			count++;
			j /= 5;
		}
	}
	return count;
}
    这样的方法的时间复杂度为O(n)。

    方法二

    方法一中对不含质因数5的数也进行了推断,时间复杂度较高。

    另外一种方法要用到例如以下结论:

    N!中含有的质因数k的个数为:[N/k]+[N/k^2]+[N/k^3]+...(总会存在一个t,使得k^t>N,便有[N/k^t]=0)

    这个公式的证明并不难,自己尝试去理解下,当中[N/k]等于1,2,3...N中能被k整除的数的个数。

    用该方法写成的代码例如以下;

int FactorialNum0_2(int n)
{
	int count = 0;
	while(n)
	{
		count += n/5;
		n /= 5;
	}
	return count;
}
    该方法的时间复杂度为O(log5n)(这里是以5为底n的对数)。


题目二

    N!二进制表示中最低位1的位置

    比方3!为6,二进制为1010,那么最低位的1的位置为2,这里最左边的位置从1開始计算。

    方法一

    我们如果最低位的1的位置为第k位,则N!=2^(k-1)+...+2^t+...+2^p+...,这里t,p等意味着后面的第t+1位,第p+1位也为1,且p>t>k。这样非常easy看出来,N!中质因数2的个数为k-1个,也就是说,最低位1的位置即为N!中质因数2的个数加1,因此问题又转化为了求N!中质因数2的个数了,相同利用结论:

N!中含有的质因数k的个数为:[N/k]+[N/k^2]+[N/k^3]+...(总会存在一个t,使得k^t>N,便有[N/k^t]=0)

    这样写出的代码例如以下:

int Lowest1(int n)
{
	int count = 0;
	while(n)
	{
		count += n/2;
		n /= 2;
	}
	return count + 1;
}
    该方法时间复杂度为O(log2n)(这里是以2为底n的对数)。

    方法二

    用例如以下结论:N!中含有质因数2的个数,等于N减去N的二进制表示中1的个数,关于怎样求N的二进制表示中1的个数,參考我的这篇博文;http://blog.csdn.net/ns_code/article/details/25425577,最快的方法的操作步骤仅仅与二进制中1的个数相等,因此这样的方法的时间复杂度更好一些,尽管方法一已经非常快了,这样的方法的时间复杂度为O(k),当中k为N!中1的个数。这样的方法的代码不再给出。


posted @ 2014-06-20 09:19  mfrbuaa  阅读(714)  评论(0编辑  收藏  举报