高精度问题--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>>n)
{
memset(f,0,sizeof(f));
f[0]=1;
count=1;
for (int i=1;i<=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("%d",f[k]);
for(int j=k-1;j>=0;j--)
printf("%05d",f[j]);
cout<<endl;
}
return 0;
}
还可以用斯特林函数大约求出n！的位数
</count></stdio></string></iostream>