码迷,mamicode.com
首页 > 其他好文 > 详细

poj 1423 打表/斯特林公式

时间:2018-02-07 14:50:23      阅读:154      评论:0      收藏:0      [点我收藏+]

标签:主程   blog   防止   iostream   math   while   alt   exp   strong   

技术分享图片

对于n位数的计算,我们可以采用(int)log10(n) + 1的方法得到n的位数

第一种方法:

对于n!位数的计算,log10(n!) = log10(1) + log10(2) + ... + log10(n)

为防止直接暴力超时这部分运算可以打表等待主程序调用

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

const int MAXN = 1e7;
int ans[MAXN +1];

void action(int m)//打表计算n!位数,存在ans数组中 
{
    double d = 0;
    for(int i = 1;i<=m;i++)
    {
        d += log10(double(i));//累加log10(i) 
        ans[i] = (int)d + 1;//向下取整并+1 
    }
}

int main()
{
    int n,m;
    cin>>n;
    action(MAXN);
    while(n--)
    {
        cin>>m;
        cout<<ans[m]<<endl;
    }
    return 0;
}

第二种方法:

对于n!的计算,也可以用斯特林公式:

技术分享图片

然后直接计算(int)log10(n!) + 1

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

double pi = acos((double)-1);

int main()
{
    int n,m;
    cin>>n;
    while(n--)
    {
        cin>>m;
        cout<<(int)(log10(sqrt(2 * m * pi) )+ m * log10(m / exp((double)1)))+ 1<<endl;
    }
}

 

poj 1423 打表/斯特林公式

标签:主程   blog   防止   iostream   math   while   alt   exp   strong   

原文地址:https://www.cnblogs.com/ZZUGPY/p/8425636.html

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