杭电1018

大数的阶乘

开始做的时候一直超限

#include<stdio.h>
#include<math.h>

int main(void)
{
    int a,j;
    double b=1.0;
    scanf("%d",&j);
    while(j--)
    {
        scanf("%d",&a);
        for(int i=1;i<=a;i++)
            b*=i;
        printf("%d\n",(int)log10(b)+1);
        b=1.0;
    }

    return 0;
}

不知道该怎么做～～

然后百度啊～～

看到了文库有一个文档，豁然开朗～～

1.位数=log10（n！）+1

   =log1+log2......+log(n)+1

#include<stdio.h>
#include<math.h>

int main(void)
{
    int a,j;
    double b=0;
    scanf("%d",&j);
    while(j--)
    {
        scanf("%d",&a);
        for(int i=1;i<=a;i++)
            b+=log10(i*1.0);
        printf("%d\n",(int)b+1);
        b=0;
    }

    return 0;
}

2.斯特灵公式（stirling公式）

lnN!=NlnN-N +0.5*ln( 2*N*pi)

（这是百科http://baike.baidu.com/link?url=3NdlL8WA1cTJlQkDs9paBb9PsHRjaLWyOXiaA9zOXYCVH2ETFo1F5n5B52IV3befUO8Ipx5BABrRcaZlJSnw6q#2)

ceil(a)//返回大于等于a的最小整数

为了避免e的麻烦不用log10，而采用log10（x)=l0gx/log10

注意：a=1时，log1=0，因而1要特殊判断

#include<stdio.h>
#include<math.h>

int main(void)
{
    int a,j;
    double b;
    scanf("%d",&j);
    while(j--)
    {
        scanf("%d",&a);
        b=ceil((a*log(a)-a+0.5*log(2*a*3.1415926))/log(10));
        if(a==1)
            b=1;
        printf("%d\n",(int)b);
    }

    return 0;
}