高精度问题--n的阶乘



高精度求阶乘

#include <iostream>
#define maxn 100000
using namespace std;

int main()
{
    int a[maxn],x,c;
    cin>>x;
    memset(a,0,sizeof(a));
    a[0]=1;
    for(int i=2;i<=x;i++)
    {
      c=0;
      for(int j=0;j<maxn;j++)
      {
        int s=a[j]*i+c;
        a[j]=s%10;
        c=s/10; 
      } 
    }
    for(int i=maxn-1;i>=0;i--)
    if(a[i])
    break;
    printf("%d",a[i]);
    for(int j=i-1;j>=0;j--)
    printf("%d",a[j]);
    printf("\n");
    return 0;
}

优化代码

#include <iostream>
#define maxn 100000
using namespace std;

int main()
{
    int a[maxn],x,c;
    cin>>x;
    memset(a,0,sizeof(a));
    a[0]=1;
    for(int i=2;i<=x;i++)
    {
      c=0;
      for(int j=0;j<maxn;j++)
      {
        int s=a[j]*i+c;
        a[j]=s%100000;
        c=s/100000; 
      } 
    }
    for(int i=maxn-1;i>=0;i--)
    if(a[i])
    break;
    printf("%d",a[i]);
    for(int j=i-1;j>=0;j--)
    printf("%05d",a[j]);
    printf("\n");
    return 0;
}

终极优化代码

#include <iostream>
#include <string h="">
#include <stdio h="">
using namespace std;
const int maxn=8000;
int f[maxn];
int main()
{
    long long n,k,c,count;
    while(cin&gt;&gt;n)
    {
     memset(f,0,sizeof(f));
     f[0]=1;
     count=1;
     for (int i=1;i&lt;=n;i++)
     {
         c=0;
         for(int j=0;j<count j="" int="" s="f[j]*i+c;" f="" j="" s="" 100000="" c="s/100000;" while="" c="" f="" count="" c="" 100000="" c="c/100000;" for="" k="count;k">=0;k--)
  if (f[k]) break;
  printf(&quot;%d&quot;,f[k]);
  for(int j=k-1;j&gt;=0;j--)
  printf(&quot;%05d&quot;,f[j]);
  cout&lt;&lt;endl;
    }
    return 0;
}
还可以用斯特林函数大约求出n！的位数
</count></stdio></string></iostream>