#include
#include
#include
struct list{ int value;list*next;list*pre;};
/*对于输入的n想办法昼精确地估计出n!所占的位数.就能确定数组元素的个数,
可以将n!表示成10的次幂,即n!=10^M(10的M次方)则不小于M的最小整数就是n!的位数，
对该式两边取对数，有=log10^n!即：
M = log10^1+log10^2+log10^3...+log10^n
循环求和,就能算得M值，该M是n!的精确位数*/
int GetBitNum(int n)
{
double sum = 1.0;
for(int i=1; i<=n; i++)
{
sum += log10((long double)i);
}
return int(sum);
}
/*创建链表，头结点不存储数据，每个节点表示一位
低位位在前，高位在后，如9302:
listhead--->2---->0---->3---->9;初始化为1*/
void CreateList(list**L,int n)
{
list*p;
list*prelist;
*L=(list*)(malloc(sizeof(list)));
prelist=*L;
(*L)->next=NULL;
for(int i=0;i {
p=(list*)malloc(sizeof(list));
p->value=0;
p->next=(*L)->next;
(*L)->next=p;
prelist->pre=p;
prelist=p;
}
(*L)->next->value=1;
(*L)->next->pre=*L;
}
/*定义链表和整数的乘法,模仿进位进行计算*/
void mul(list*L,int num)
{
int n=0;
L=L->next;
for(;L!=NULL;L=L->next)
{
n += (num*L->value);
L->value = n % 10;
n /= 10;
}
}
/*定义阶乘*/
void fac(list*L,int n)
{
for(int i=1;i<=n;i++)
{
mul(L,i);
}
}
/*输出链表*/
void echo(list*L,int n)
{
list*p=L->pre;
for(int i=0;ipre)
{
printf("%d",p->value);
if((i+1)%4==0)printf(" ");
}
}
int main()
{
int bitnum,n;
list*plist;
char flag='y';
int in;
while(flag=='y'||flag=='Y')
{
printf("输入一个1-100的数:");
fflush(stdin);
in=scanf("%d",&n);
if(!(n>=1&&n<=100)||in!=1)
{
printf("输入无效，请重新输入!\n");
continue;
}
bitnum=GetBitNum(n);
CreateList(&plist,bitnum);
fac(plist,n);
echo(plist,bitnum);
printf("\n是否再次继续运算(Y/N): ");
fflush(stdin);
scanf("%c",&flag);
}
}
