编程之美： 第二章 数字之魅 2.2不要被阶乘吓到

/*
不要被阶乘吓到:
1给定一个整数N，那么N的阶乘等于N!，末尾有多少个0呢，例如N=10，N!=3628800，N!的末尾有两个0
2求N!的二进制表示中最低位为1的位置。

分析:
哪些数相乘能够得到10
如果N!= K*10^M,且K不能被10整除，那么N!末尾有M个0，再考虑对N!进行质数分解，N!=(2^X)*(3^Y)*(5^Z)...,由于10 = 1*5，所以M只跟X和Z有关，每一对2和5相乘
可以得到一个10，于是M= min(X,Z)。又因为X>=Z，因为能被2整除的概率比能被5整除的概率大得多。
所以M=Z

输入:
3
10
100
输出：
0
2

2
9

24
3

解法2:
公式Z = [N/5] + [N/5^2] + [N/5^3] + ...总存在一个k使得5^k大于N,使得[N/5^k] = 0.公式中[N/5]表示不大于N的数中5的倍数贡献一个5，[N/5^2]表示不大于N的
数中5^2的倍数再贡献一个5

问题2求N!的二进制表示中最低位为1的位置，例如N=3,N！= 6 = 0110的最低位1在第二位

把一个二进制数除以2:
判断最后一个二进制位是否为0，如果为0，则将次二进制数右移一位，即为商值；反之若为1，则说明这个二进制数是奇数，无法被2整除。
这个问题实际上等同求N!含有质因数2的个数。答案= N!含有质因数2的个数+1
*/

/*
关键:
1 		iRes = i;//注意这里是拿每一个数试除
		while(iRes % 5 == 0)//求质因数中5的指数即可
		{
			iRet++;
			iRes /= 5;
公式Z = [N/5] + [N/5^2] + [N/5^3] + ...总存在一个k使得5^k大于N,使得[N/5^k] = 0.公式中[N/5]表示不大于N的数中5的倍数贡献一个5，[N/5^2]表示不大于N的
数中5^2的倍数再贡献一个5
2 	while(n)
	{
		iRet += (n >> 1);
		n >>= 1;
	}
	return iRet + 1;//求最低位为1的位置，因为n-1个位置上为0，表示这个数可以被2^(n-1)整除，需要再加上1，例如8表示1000，最终为3+1 = 4
3判断一个数是否为2的方幂，n > 0 && ((n & (n-1)) == 0)
*/

#include <stdio.h>


int times_0(int n)
{
	int iRes;
	int iRet = 0;
	for(int i = 1 ; i <= n ; i++)
	{
		iRes = i;//注意这里是拿每一个数试除
		while(iRes % 5 == 0)//求质因数中5的指数即可
		{
			iRet++;
			iRes /= 5;
		}
	}
	return iRet;
}

int times_0(int n,int i)
{
	int iRet = 0;
	while(n)
	{
		iRet += n / 5;//参见公式Z= [N/5] + [N/5^2] + [N/5^3] + ....(直到存在5^K > N,使得[N/5^K]= 0)
		n /= 5;
	}
	return iRet;
}

int lowestOne(int n)
{
	int iRet = 0;
	while(n)
	{
		iRet += (n >> 1);
		n >>= 1;
	}
	return iRet + 1;//求最低位为1的位置，因为n-1个位置上为0，表示这个数可以被2^(n-1)整除，需要再加上1，例如8表示1000，最终为3+1 = 4
}

void process()
{
	int n;
	while(EOF != scanf("%d",&n))
	{
		//printf("%d\n",times_0(n,1));
		printf("%d\n",times_0(n));
		printf("%d\n",lowestOne(n));
	}
}

int main(int argc,char* argv[])
{
	process();
	getchar();
	return 0;
}