阶乘浅析poj1150 3406 zoj1222 2358

原地址;http://blog.csdn.net/shiwei408/article/details/8739433

阶乘问题分为几类:

1.求阶乘末尾0的个数,,直接除以5,累加即可。

2.求阶乘的结果一共有多少位,stirling公式:n!≈sqrt(2*PI*n) * (n/e)^n,直接取以10为底的对数,整数部分即为位数。http://poj.org 第1423题

3.求阶乘的最后非零位,这类问题比较复杂,专题中我们着重讨论这个问题

 

首先看POJ1150

题目大意:求n的m排列的最后非零位。

题目分析:n的m排列即P(n,m)=n!/(n-m)!,所以这题是求两个阶乘商的最后非零位。我们处理阶乘时一般是逐个处理。首先看普遍性的对于一个数n的阶乘,我们如何处理它的末尾非零位。

10的因子是2和5,这两个数不属于模10的缩系,所以我们在处理数时要剔除掉这两个因子,剔除所有数中含5和2这两个因子。

首先我们用f(n,x)表示小于等于n的数中,末尾位x的数的个数。n我们可以分解为1,3,5,7,9...  2,4,6,8,10...

因为我们将偶数中的2因子剔除了,所以2,4,6,8,10...实际上又变成了1,3,5,7,9...

所以我们可以得到递推式f(n,x)=f(n/2,x)+g(n,x)g(n,x)表示的是小于等于n的奇数中尾数为x的数的个数。

那么g(n,x)怎么求呢,首先n/10必定是其中一部分其次要看n>=x是否成立,成立就加1,不成立就不加。最后一点要注意,我们在奇数中还有汗5因子的数,所以要剔除掉5这个因子,比如5 15 25 35 45... 剔除了5之后又变成了1 3 5 7 9...所以g(n,x)=n/10+(n%10>=x)+g(n/5,x)。

有了这两个递推式再加上求5因子和2因子个数的递归除法,我们很快就能得到2因子个数,5因子个数,末尾位3 7 9因子的个数。下面我们要列个循环mod2[4]={6,2,4,8},mod3[4]={1,3,9,7},mod7[4]={1,7,9,3},mod9[4]={1,9,1,9}表示的是各个因子自己相乘时尾数的循环节,这里要特别注意一点,由于2非模10的缩系,所以2^0=1,这一点要格外注意,所以只有它循环节不是从0开始,0要单独处理。因为n-m的阶乘必定是n的阶乘的子集,所以我们可以直接用因子数相减即可。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
int cnt[8];
int mod2[4]={6,2,4,8};
int mod3[4]={1,3,9,7};
int mod7[4]={1,7,9,3};
int mod9[4]={1,9,1,9};
int getpow(int n,int k)
{
	int sum=0;
	while(n)
	{
		sum+=n/k;
		n/=k;
	}
	return sum;
}
int g(int n,int k)
{
	if(n==0)return 0;
	return n/10+(n%10>=k)+g(n/5,k); 
}
int get(int n,int k)
{
	if(n==0)return 0;
	return get(n/2,k)+g(n,k);
}
int main()
{
	int n,m;
	while(~scanf("%d%d",&n,&m))
	{
		int res=1;
		memset(cnt,0,sizeof(cnt));
		cnt[2]=getpow(n,2)-getpow(n-m,2);
		cnt[5]=getpow(n,5)-getpow(n-m,5);
		cnt[3]=get(n,3)-get(n-m,3);
		cnt[7]=get(n,7)-get(n-m,7);
		cnt[9]=get(n,9)-get(n-m,9);
		if(cnt[2]>cnt[5]){res*=mod2[(cnt[2]-cnt[5])%4];}
		else if(cnt[2]<cnt[5]){res*=5;}
		res=res*mod3[cnt[3]%4]*mod7[cnt[7]%4]*mod9[cnt[9]%4];
		printf("%d\n",res%10);
	}
	return 0;
}

再来看POJ3406

求n!/(m!*(n-m)!)的最后非零位

与上题一样,但是由于m的阶乘和n-m的阶乘的乘积并不是n的阶乘的子集,但是我们可以通过观察发现,只有n的阶乘中3的个数会少于m的阶乘和n-m的阶乘3的数目之和,但是我们不要担心,这个时候我们可以通过将n的阶乘中的9的个数分解得到。还有一种方法是再多加一个循环模(我表示不理解)

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
int cnt[8];
int mod2[4]={6,2,4,8};
int mod3[4]={1,3,9,7};
int mod7[4]={1,7,9,3};
int mod9[4]={1,9,1,9};
int getpow(int n,int k)
{
	int sum=0;
	while(n)
	{
		sum+=n/k;
		n/=k;
	}
	return sum;
}
int g(int n,int k)
{
	if(n==0)return 0;
	return n/10+(n%10>=k)+g(n/5,k); 
}
int get(int n,int k)
{
	if(n==0)return 0;
	return get(n/2,k)+g(n,k);
}
int main()
{
	int n,m;
	while(~scanf("%d%d",&n,&m))
	{
		int res=1;
		memset(cnt,0,sizeof(cnt));
		cnt[2]=getpow(n,2)-getpow(n-m,2)-getpow(m,2);
		cnt[5]=getpow(n,5)-getpow(n-m,5)-getpow(m,5);
		cnt[3]=get(n,3)-get(n-m,3)-get(m,3);
		cnt[7]=get(n,7)-get(n-m,7)-get(m,7);
		cnt[9]=get(n,9)-get(n-m,9)-get(m,9);
		cnt[3]+=cnt[9]*2;cnt[9]=0;
		if(cnt[2]>cnt[5]){res*=mod2[(cnt[2]-cnt[5])%4];}
		else if(cnt[2]<cnt[5]){res*=5;}
		res=res*mod3[cnt[3]%4]*mod7[cnt[7]%4]*mod9[cnt[9]%4];
		printf("%d\n",res%10);
	}
	return 0;
}

最后看一个BT题ZOJ1222

这题在很多OJ上都有,只不过数量级很小,可以用上面的方法轻松ac,但是这题据说输入是一个字符串,可能有100+位,所以面对如此高精度大数,要用新的方法来解决。

思路不变,还是以2和5为核心,把每个数最后的尾数乘起来,只是由于数的极速扩大,我们不能用原来的递归求解了,要用到一个循环节:int mod[10]={1,1,2,6,4,4,4,8,4,6};这个表示从0乘到9(将5换为1)的非0尾数。当n<5时,直接输出mod[n].我们要去掉尾数0,就要去掉相同数量的2和5,我们发现2的个数是远远多于5的个数,并且因子2的数目是最多的,所以最后得到的结果一定是一个偶数,于是存在这样一个特殊的尾数除法:2/2=6,4/2=2,6/2=8,8/2=4。同时在除以2的个数时,还是由于2的个数比较多的原因,是存在一个除法的循环节的,循环节的长度为4.并且这10个尾数相乘为6,当n>10时,我们可以得到一个公式:

此时我们就可以应用这个公式进行求解了,值得注意的是n是一个非常大的数,在除以5的时候我们可以利用高精度除法。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
char str[1001];
int mod[10]={1,1,2,6,4,4,4,8,4,6};
int a[1001];
int getAns(int cnt)
{
	if(cnt==1&&a[0]<5)return mod[a[0]];
	int res=0;
	int w=a[0];
	for(int i=cnt-1;i>=0;--i)
	{
		res=res*10+a[i];
		a[i]=res/5;
		res=res%5;
	}
	if(a[cnt-1]==0)--cnt;
	int m=(a[1]*10+a[0])%4;
	int tmp=getAns(cnt)*mod[w]*6%10;
	while(m--)
	{
		if(tmp==2)tmp=6;
		else if(tmp==6)tmp=8;
		else tmp/=2;
	}
	return tmp;

}
int main()
{
	while(~scanf("%s",str))
	{
		int len=strlen(str);
		memset(a,0,sizeof(a));
		for(int i=0;i<len;++i)
		{
			a[i]=str[len-1-i]-'0';
		}
		int ans=getAns(len);
		printf("%d\n",ans);
	}
	return 0; 
}

4.数字n是否可以表示为若干个不同的阶乘的和

实际上这个问题跟数论无关,是一道搜索的问题http://acm.zju.edu.cn的第2358题,由于n比较小,只可能是0!~9!的数字和,所以可以直接用dfs搜索。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <map>
using namespace std;
int a[10]={1,1,2,6,24,120,720,5040,40320,362880};
int sum;
bool vis[11];
map<int,int>mymap;
void dfs(int n)
{
if(n==10)
{
if(mymap.find(sum)==mymap.end())mymap[sum]=1;
return;
}
sum+=a[n];
dfs(n+1);
sum-=a[n];
dfs(n+1);
}
int main()
{
memset(vis,false,sizeof(vis));
sum=0;
mymap.clear();
dfs(0);
int n; 
while(~scanf("%d",&n))

if(n<0)break;
if(n==0||mymap.find(n)==mymap.end())puts("NO");
else puts("YES");
}
return 0;
}

你可能感兴趣的:(算法,数学,ZOJ,poj)