杭电１０４２

这是一道典型的大数题，涉及到大数加法，大数乘法，大数分成字数的乘法方法。这个比较麻烦。而且容易超时，和超过内存，所以今天就不用这种方法，但是后面还是会贴上一个代码。我们首先想想两个十进制的数字相乘（两个数不大于１０）那么结果进位的情况最多也是进到百位置上，那么两个数（不大于１００）最多是１００００，那么不大于１００００的话，最差的是１００００×１００００；在这里这个题中是不大于１００００，那么我们用一个字符窜承载的话，应该需要４位（１００００的时候会进位不会保存）。那么这就引出了今天的解法。数组存储法。我们的想法是将结果没４位分成一组，然后将这些数字直接用整数形式的数组存储，那么再进位上我们必须以１００００为基本（结果进位必须是大于１００００的数字），原因就是上面那几句，如果不理解，你可以试试如果以１０００为基本看看，结果会发生什么样的事情。想打这一点我们自然而然就形成一个想法，这个想法就是将１００００看成是１０然后直接按照平常的乘法运算，这样的话，有点类似二进制压缩的思想。

注意：当我们数组中有一个是１６００，这时候需要乘以１３结果是１０８００会向高位进１０８００／１００００＝１；而这一位就会变成８００.本来的假设他的上一位是１０，那么乘以１３再加上１就是１３１，原本的数字是１３１０８００，但是你的这时候两个数组变量合起来就是１３１８００中间少了一位０，这时候我们需要补上，那么什么情况补上呢？答案是什么情况都需要，我们必须保证每一个数组的数字都是四位数，这就是在ｃ里面的printf("%04d",a[i]);他的意识是结果以四位数字显示，不够的在前面加上０，这个时候就可以保证分割后的数字和以前的是一样的了。

下面是一个AC的代码：

#include<iostream>
using namespace std;
#include<cstdlib>
const int Max=10000;
int a[Max];
int main()
{
   int N;
   int i,j,k;
   int t;
   
   while(cin>>N)
   {
  memset(a,0,sizeof(a));
  a[0]=1;
        t=0;
      for(i=1;i<=N;i++)
 {
 k=0;
 for(j=0;j<=t;j++)
 {
 a[j]=a[j]*i+k;
 k=a[j]/10000;
 a[j]=a[j]%10000;
 }
 if(k!=0)
 {
 t++;
 a[t]=k;
 }
 }
 for(i=t;i>=0;i--)
 {
 if(i==t)
 cout<<a[i];
 else 
 printf("%04d",a[i]);
 }
 cout<<endl;
   }
   return 0;
}

再给大家一个用大数做的代码：（注：这个是自己写的代码，试了几十个例子都没有问题但是ＡＣ不了，贴这个代码旨在说一下方法）：

#include<iostream>
using namespace std;
#include<cstdlib>
#include<string>
string plus(string s1,string s2,int p)
{
int a1,b1;
int i,j;
a1=s1.size();
b1=s2.size();
int num=0;
if(a1>=(b1+1))
{
for(i=a1-1-p,j=b1-1;j>=0;i--,j--)
{
// cout<<"hao4"<<endl;
s1[i]+=s2[j]-'0';
}
//cout<<"hao5"<<endl;
for(i=a1-1;i>0;i--)
if(s1[i]>='9')
s1[i-1]+=1;
if(s1[0]>'9')
{ 
s1[0]-=10;
s1='1'+10;
}
return s1;
}
else 
{
for(i=a1-2,j=b1-1;i>=0;i--,j--)
{
//cout<<"hao3"<<endl;
s2[j]+=s1[i]-'0';
}
if(s1!="0")
{
if(p>a1)
{
for(i=p-1;i>=a1;i--)
s2=s2+'0';
for(i=a1-1;i>=0;i--)
s2=s2+s1[i];
}
else for(i=a1-1;i>=0;i--)
s2=s2+s1[i];
}
for(i=b1;i>0;i--)
if(s2[i]>'9')
s2[i-1]+=1;
if(s2[0]>'9')
{
s2[0]-=10;
s2='1'+s2;
}
return s2;
}
}
string ss(string s,int n)
{
string s1,s2;
s2="0";
int i;
//cout<<s2.size()<<endl;
int a[6];
memset(a,0,sizeof(a));
int j=0;
int k=0;
while(n)
{
a[k++]=n%10;
n/=10;
}
int num=0;
int sum=0;
int m=s.size();
for(j=0;j<k;j++)
{
s1=s;
num=0;
for(i=m-1;i>=0;i--)
{ 
/// cout<<"hao2"<<endl;
sum=(s1[i]-'0')*a[j];
s1[i]=(s1[i]-'0')*a[j]%10+num+'0';
if(s1[i]>'9')
{
s1[i]-=10;
num=sum/10+1;
}
else 
num=sum/10;
}
if(num!=0)
{

s1='1'+s1;
s1[0]+=(num-1);
if(s1[0]>'9')
{
s1[0]-=10;
s1='1'+s1;
}
}
s2=plus(s2,s1,j);
}
return s2;
}
int main()
{
int N;
int i;

while(cin>>N)
{
string s="1";
for(i=1;i<=N;i++)
{
//cout<<s<<"hao6"<<endl;
s=ss(s,i);
//  cout<<"hao1"<<endl;
}
if(N==0)
cout<<"0"<<endl;
else 
cout<<s<<endl;
}
return 0;
}