码迷,mamicode.com
首页 > 编程语言 > 详细

poj1423---求一个大数的位数方法,我猜网站上统计输入字符少于多少位的那个算法

时间:2015-05-02 00:57:46      阅读:210      评论:0      收藏:0      [点我收藏+]

标签:

法一:对一个数求它的对数,+1取整为其位数

问题转化为int (log10(N!)+1),对数性质log10(N!)=log10(N)+log10(N-1)+...+log10(1)

/*用log10求位数*/

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

int main()
{
    int tim,N;
    scanf("%d",&tim);
    while(tim--)
    {
        int i;
        double NumOfDigit=1;
        scanf("%d",&N);
        for(i=N;i>=1;i--)
        {
            NumOfDigit+=log10(i);
        }
        printf("%d\n",(int)NumOfDigit);
    }
}

  当n偏大的时候,时间长,TLE

法二:Stirling公式

log(n!) = log10(sqrt(2*pi*n)) + n*log10(n/e);

/*用斯特灵求位数*/

#include<stdio.h>
#include<math.h>
#define e 2.718281828459045
#define pi 3.141592653589793239
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n;
        scanf("%d",&N);
        double num_digit;
        num_digit=log10(sqrt(2*pi*N)) + N*log10(N/e);
        printf("%d\n",(int)num_digit+1);
    }
    return 0;
}

  WA两次原因:

 num_digit=log10(sqrt(2*pi*N)) + N*log10(N/e)+1;
        printf("%d\n",(int)num_digit);

  当N=1,num_digit=0.几,因为当n=1时,

log10(sqrt(2*pi*N)) + N*log10(N/e)=-0.03

  最后值就为0

总结:用斯特林公式求位数时,考虑到N=1,加1放在取了整之后

(int)(log10(sqrt(2*pi*N)) + N*log10(N/e))=0
加1放在取了整之后

int(3.1+1)=4
int(3.1)+1=4
int(3)+1=4
int(3+1)=4
int(-0.03)+1=1
int(-0.03+1)=0

poj1423---求一个大数的位数方法,我猜网站上统计输入字符少于多少位的那个算法

标签:

原文地址:http://www.cnblogs.com/gabygoole/p/4471397.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!