一道清华上机题目(求n! 高精度)

//writen by sunboy 
//2007.08.30
//说明:输入n(0<=n<=1000)  输出n!
#include <iostream>
#include <string>
#include <cstring>
#include<sstream> 
using namespace std;
string hp_mul(string &a,string &b);   //返回a*b ,字符串乘法
void change(int a[],string &str)      //将string转换为数组表示方式
...{
   int i,len=str.length ();
   for(i=len-1;i>=0;i--)
       a[len-i-1]=str[i]-'0';
}
string itos(int i)   //将int 转换成string 
...{ 
  stringstream s; 
  s << i; 
  return s.str(); 
} 
int main()
...{
   string f1,f2,f3;
   int n,i;
   while(cin>>n)
   ...{
       if(n==0)
       ...{
           cout<<0<<endl;
           continue;
       }
      f1="1";
      for(i=2;i<=n;i++)
        f1=hp_mul(f1,itos(i));
      cout<<f1<<endl;
  }
   return 0;
}

string hp_mul(string &a,string &b) //返回a*b ,字符串乘法
...{
   int i,j,x[5120],y[5120],z[5120];
   memset(x,0,sizeof(x));
   memset(y,0,sizeof(y));
   memset(z,0,sizeof(z));
   int len1=a.length ();
   int len2=b.length ();
   change(x,a);
   change(y,b);
   for(i=0;i<len1;i++)
       for(j=0;j<len2;j++)
           z[i+j]+=x[i]*y[j];
   for(i=0;i<=len1+len2;i++)
           if(z[i]>=10)
           ...{ z[i+1]+=z[i]/10;
            z[i]%=10;
           }
   for(i=len1+len2+2;i>=0;i--)
       if(z[i])break;
        string res="";
               
   for(j=i;j>=0;j--)
       res+=z[j]+'0';
   return res;
}