c语言 10000以内的阶层,10000以内的阶乘（ZT）

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

#include

# include 

// 计算

# define N 1000

int cal(unsigned int *s,int n)

{

unsigned long p; // p是对每一位乘法中的值加上进位，如34*5,4*5是20,3*5的加上进位2是17

unsigned long k=0; // k是一次乘法中的进位，如10进制乘法中，34*5,4*5的进位是2,3*5的进位是1

int i;

static int m=1; // m是位数，表示s有总共有多少位数字，注意：是1000进制

static int b=0;

/* b用来记录后面的0，比如213,000,000,000，则b=3,后面的3个000不必再参与计算了 */

//

for(i=b;i

{

p=(long)s[i]*(long)n+k;

k=p/N;

s[i]=p-k*N;

}

// b是低位乘出来的000的数目,增加后加1

while(!s[b])

b++;

// 最高位的进位处理

for(k=p/N;k;)

{

p=k;

k=p/N;

s[i++]=p-k*N;

m++; // 进一次m加一次

}

return m;

}

int main(int argc,char**argv)

{

/* s是用来存计算结果的，以N为进位，这里N=1000，如s[0]=1,s[1]=21,s[2]=213,s[3]以上都为0，

则结果是：213,021,001

*/

unsigned int *s;

int i;

int m; // m是位数，表示s有总共有多少位数字，注意：是1000进制

int n; /* 求n! (0

scanf("%d",&n); // 输入n的值

s=(unsigned int *)malloc(n*sizeof(s)*10); // 数组开足够大，其实要不了这么大

// 附初值，S=1，即s[0]=1,以上都为0

for(i=1;i

s[i]=0;

s[0]=1;

// 循环，s乘以2、3、4,……,n

for(i=2;i<=n;i++)

m=cal(s,i);

// 输出：n!=

printf("\n%d!=",n);

m=n-1;

// 滤掉前面的0,前面说了，其实

while(!s[m])

m--;

// 输出第一位

printf("%d",s[m--]);

// 输出后面的位数

for(i=m;i>=0;i--)

printf(",%03d",s[i]);

printf("\n");

// 释放内存空间

free(s);

}