阶乘算法

寄信人: Fll (沧海一声笑)
标  题: 阶乘算法
发信站: 郁金香BBS站 (2004年08月17日13:48:27 星期二)
来  源: 10.13.3.118

#include<stdio.h>
#include<iostream.h>
int fact(int n)
{
        int a[900];/*记录结果*/
        int flag=1;/*结果位数标志*/
        int data=n;/*计算过程中的被乘数*/
        int i=0;/*循环标志*/
        if((n<0)||(n>1042)){      /*只对0———1041之间的数进行计算*/
                printf("数字溢出!/n");
                return 0;
        }
        if(n>=1000) flag = 2;  /*如果输入的值大于2位，初始化头两个数组值*/
        for(i=0;i<900;i++)
            a[i]=0;
        a[0]=1;
        while((data>0)&&(flag<900)){   /*主运算*/
            for(i=0; i<flag; i++)/*将上次的结果和下一个被乘数进行一轮运算*/
                a[i]=a[i]*data;
            for(i=0; i<flag-1; i++){ /*调整结果*/
                a[i+1]+=a[i]/1000;
                a[i]=a[i]%1000;
            }
            if(a[flag-1]>1000){/*调整结果*/
                flag++;
                a[i+2]=a[i+1]/1000;
                a[i+1]=a[i+1]%1000;
            }
            data--;            /*被乘数自减*/
        }
        printf("%d != ",n);  /*结果的打印*/
        printf("%d",a[flag-1]);
        for(i=flag-2;i>=0;i--){          /*针对三位数字的头两位是否为0进行处
理*/
            if(a[i]<100) printf("0");
            if(a[i]<10)  printf("0");
            printf("%d",a[i]);
        }
        printf("/n");