大数阶乘的位数

利用log10的性质

我们知道整数n的位数的计算方法为：log10(n)+1
故n!的位数为log10(n!)+1
即：log10(n)+log10(n-1)+…+log10(1)+1

#include<iostream>
using namespace std;

double sum[1000];
int main()
{
    sum[1]=1;
    for(int i=2; i<1000; i++)
        sum[i]=sum[i-1]+log10(i);
     for(int i=1;i<16;i++){
        cout<<i<<" "<<sum[i]<<endl;
    }
｝

斯特林（Stirling）公式

用公式求n！就是log10((2*PI*n)^1/2*(n/e)^n)+1
即 1/2*log10(2*PI*n)+n*log10(n/e)+1

#include<iostream>
using namespace std;
#define PI 3.141592654
#define E 2.71828182846
int l(int n)
{
    int s=1;
    if(n>3)
        s=log10(2*PI*n)/2+n*log10(n/E)+1;
    return s;
}

实验

通过这里看出阶乘n>6后，位数是独立的。就可以做很多事情了。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
using namespace std;

typedef long long LL;
const int mod=1e9+7;

double sum[1000];
int main()
{
    sum[1]=1;
    for(int i=2; i<1000; i++)
        sum[i]=sum[i-1]+log10(i);
    LL q=1;
    for(int i=1;i<16;i++){
        q*=i;
        cout<<i<<" "<<sum[i];
        cout<<" "<<q<<endl;
    }
    return 0;
}