呵呵!刚刚在chinaunix注册了,这是第一篇日志,大部分借鉴的是hhccbb的算法。写的不好的,希望大牛们谅解。现在转入正题。
hdu1042要求求n!,而n的值最大到10000。属于规模型解题范围,10000!的位数为35670。
(一)首先,了解n!的阶乘范围:
(二)n!(n在两位数时)的基本求法:
{//a[0]存放n!的位数,a[1],...,a[a[0]]逆序存放每位数字。
//若于高精度连用,注意将a[1],...,a[a[0]]转置。
/*for(i=1,k=a[0];i<k;i++,k--) {temp=a[i];a[i]=a[k];a[k]=temp;}*/
int i,j,k=1;
a[1]=1;
for(i=2;i<maxsize;i++) a[i]=0;
for(i=2;i<=n;i++)
{
for(j=1;j<=k;j++) a[j]*=i;
for(j=1;j<=k;j++)
{
if(a[j]>=10)
{if(j==k) k++;//k计位数。
a[j+1]+=a[j]/10;
a[j]%=10;
}
}//for j
}//for i
a[0]=k;
}
(三)n!(n很大时,处于三四位时)的求解:
假设要求10000!(35670位数),若采用标题(二)的算法,利用每个数组元素来存放每位数字的话,显然造成了空间的浪费,而循环也使得算法(二)对10000!不适合。
所以可以考虑用每个数组元素来存放n!的每个10位数字,只需要把算法(二)的a[]类型设置为__int64型,将进制扩大为10^10。
hdu1042的代码设计如下:
#include<stdio.h>
//#include<math.h>
#define I 10000000000
#define maxsize 35670
//10000! 位数:35660
__int64 a[maxsize];
void fac2(int n,__int64 *a)
{int i,j;__int64 k;
if(n==0||n==1)
{a[0]=1;a[1]=1;return;
}
a[0]=1;a[1]=1;k=1;
for(i=2;i<maxsize;i++) a[i]=0;
for(i=2;i<=n;i++)
{
for(j=1;j<=k;j++) if(a[j]) a[j]*=i;
for(j=1;j<=k;j++)
{
if(a[j]>=I)
{
if(j==k) k++;
a[j+1]+=a[j]/I;
a[j]%=I;
}
}//for j
}//for i
a[0]=k;
}
int main()
{int n;__int64 i,j;
while(scanf("%d",&n)!=EOF)
//for(n=3;n<=50;n++)
{fac2(n,a);
//printf("%d!=/n",n);
printf("%I64d",a[a[0]]);
for(i=a[0]-1;i>=1;i--)
{
//若a[i]不足9位数,前面补零。
if(a[i]<10) printf("000000000");
else if(a[i]<100) printf("00000000");
else if(a[i]<1000) printf("0000000");
else if(a[i]<10000) printf("000000");
else if(a[i]<100000) printf("00000");
else if(a[i]<1000000) printf("0000");
else if(a[i]<10000000) printf("000");
else if(a[i]<100000000) printf("00");
else if(a[i]<1000000000) printf("0");
printf("%I64d",a[i]);
}
printf("/n");
}
return 1;
}注意n!的输出,由于除a[a[0]]外前面的每个数组元素位数不一定为10位。所以这些元素输出时,需要前面补零凑够10位。
if(a[i]<10) printf("000000000");
else if(a[i]<100) printf("00000000");
else if(a[i]<1000) printf("0000000");
else if(a[i]<10000) printf("000000");
else if(a[i]<100000) printf("00000");
else if(a[i]<1000000) printf("0000");
else if(a[i]<10000000) printf("000");
else if(a[i]<100000000) printf("00");
else if(a[i]<1000000000) printf("0");
此部分的设计,也可以先求出a[i]的位数,然后根据位数补足10位。
for(j=1;j<=10-digit(a[i]);j++) printf("0");
printf("%I64d",a[i]);//若a[i]不足10位数,前面补零。
digit()定义如下:
int digit(__int64 a)
{int num=0;
while(a) {a/=10;num++;}
return num;
}