阶乘浅析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;
}