求大数阶乘代码

[xuchaoroger]

求大数阶乘代码

[idontlikenickname]

// 3000!得结果超出了电脑能显示得范围，所以最好采用数组来记录每位
#include <ctype.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#define MAXN 0X7000

int a[MAXN];

int main(int argc,char *argv[])
{
int n,m,i,j,c,t;
printf("Enter n(n>=2) :");
while(1)
{
scanf("%d",&n);
if(n>=2&&n<=3276) break;
printf("must be 2<=n<=3276");
}
a[0]=1; m=1;
for(i=2;i<=n;i++)
{
for(c=0,j=0;j<m;j++)
{
t=a[j]*i+c;
a[j]=t%10;
c=t/10;
}
while(c)
{
a[m++]=c%10;
c/=10;
}
}
printf("%d!=\n",n);
for(i=m-1;i>=0;i--) putchar(a[i]+0x30);
printf("\npress any key to continue .\n");
getch();
return 0;
}

[idontlikenickname]

// 另外一个.
#include <stdio.h>

void main()
{
int data[40];
int digit;
int i,j,r,k;
int n;

for(i=1;i<=40;i++) /* 将数组初始值设为0 */
data[i]=0;

data[0]=1;
data[1]=1;
digit=1;

printf("Enter a number what you want to calculus : ");
scanf("%d",&n);

for(i=1;i<=n;i++)
{
for(j=1;j<=digit;j++) /* 每位上等乘上阶数 digit 决定有几位*/
data[j]*=i;
for(j=1;j<=digit;j++)
{
if(data[j]>10)
{
for(r=1;r<=digit;r++)
{
if(data[digit]>10)
digit++;
data[r+1]+=data[r]/10;
data[r]=data[r]%10;
}
}
}
printf("%d!=",i);
for(k=digit;k>0;k--)
printf("%d",data[k]);
printf("\n");
}
getch();
}