#include<stdlib.h>
#include<iostream.h>
#include<iomanip.h>

const int N=1000; 

int compute(unsigned int *s,int n)//s用來存儲一次的計算結果,n為本次計算的乘數,函數返回結果中有效數據的節數
{
    unsigned long p;  //暫時存放一節的結果
	unsigned long more=0;//一次乘法的進位
	int i;  //循環變量

	static int m=1;//存放結果中的節數,三位一節
	static int q=0;//存放結果中最后為零的節數

	//計算本次的結果
    for(i=q;i<m;i++)
	{
		p=(long)s[i]*(long)n+more;
		more=p/N;
		s[i]=p-more*N;
	}

	//計算結果中最后為零的節數
	while(s[q]==0) q++;

	//處理最高位
	for(more=p/N;more;)
	{
		p=more;
		more=p/N;
		s[i++]=p-more*N;
		m++;    //有效節數增一
	}

	return m;  //返回有效節數
}

void main()
{
	unsigned int *s; //存放結果
	int i;           //循環變量
	int m;			 //存放節數
	int n;           //求n的階乘

	cout<<"請輸入一個正整數<0----10000>:";
	cin>>n;

	if(n<0) 
	{
		cout<<"輸入數據錯誤!"<<endl;
		return;
	}
	else if(n==0)
	{
		cout<<"0!==1"<<endl;
		return;
	}

	//初始化s
	s=(unsigned int *)malloc(n*sizeof(s)*10);
    for(i=1;i<n;i++) s[i]=0;
	s[0]=1;

	//計算結果
	for(i=2;i<=n;i++) m=compute(s,i);

	while(s[m]==0) m--;//過濾掉前面的零

	//輸出結果
	cout<<n<<"!=="<<s[m--];//最高一節
    
	int num=1;
	for(i=m;i>=0;i--) 
	{   
		num++;
		cout<<','<<setw(3)<<setfill('0')<<s[i];
	    if((num%10)==0) cout<<endl;
	}
	cout<<endl;

	//釋放空間
	free(s);
}