poj 1423 Big Number

    xiaoxiao2025-09-17  33

    第一次是在杭电碰到这个题目(1018),当时刚学java的大数,对于阶乘的问题当时都是一股脑的用java去做,而在这里就不行了,或者说没必要用java大数。(用java大数要么打表的时候超内存,要么不打表超时)

    题意:1<=n<=10^7,对应一个n输出数n!的位数。

    预备知识:

    对于一个数n,(log10(n)+1)再取整就是数n的位数。 例如: log10(1000)=3.000000,log10(1000)+1=4.000000再取整就是4,1000是个4位数。 log10(9999)=3.999957,log10(9999)+1=4.999957再取整也是4,9999是个4位数。

    法一(打表法):

    //Code: #include<iostream> #include<cstdio> #include<cmath> const int maxn =10000000; int a[maxn+5]; int main() { //freopen("in.txt","r",stdin); int t; double ans=0; for(int i=1; i<=maxn; i++) { ans+=log10(i*1.0); a[i]=floor(ans)+1; } scanf("%d",&t); while(t--) { int n; scanf("%d",&n); printf("%d\n",a[n]); } return 0; } 法二:

    斯特灵公式是一条用来取n阶乘近似值数学公式。一般来说,当n很大的时候,n阶乘的计算量十分大,所以斯特灵公式十分好用,而且,即使在n很小的时候,斯特灵公式的取值已经十分准确。                           -----引自维基百科

    由预备知识可知log10(n!)+1便是答案

    但是不能直接带斯特林公式,应该带入log10(n!)化简,不然(n/e)^n会爆炸的。

    依据公式把

    log10(n!)+1

    化成

    log10(sqrt(2*acos(-1.0)*n))+n*(log10(n/exp(1.0)))+1

    就好了

    #include<iostream> #include<cstdio> #include<cmath> int main() { //freopen("in.txt","r",stdin); int t; scanf("%d",&t); while(t--) { int n; scanf("%d",&n); if(n==1){puts("1");continue;} double ans=log10(sqrt(2*acos(-1.0)*n))+n*(log10(n/exp(1.0)))+1; printf("%d\n",(int)ans); } return 0; }

    注意这个公式对n==1并不成立,所以n==1要特判一下。

    转载请注明原文地址: https://ju.6miu.com/read-1302745.html
    最新回复(0)