博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
poj-1423 NYOJ_69 数字长度 斯特林公式 对数应用
阅读量:5256 次
发布时间:2019-06-14

本文共 1408 字,大约阅读时间需要 4 分钟。

数的长度

时间限制:
3000 ms  |  内存限制:
65535 KB
难度:
1
描述

    N!阶乘是一个非常大的数,大家都知道计算公式是N!=N*(N-1)······*2*1.现在你的任务是计算出N!的位数有多少(十进制)?

输入
首行输入n,表示有多少组测试数据(n<10)
随后n行每行输入一组测试数据 N( 0 < N < 1000000 )
输出
对于每个数N,输出N!的(十进制)位数。
样例输入
31332000
样例输出
11130271 刚拿到这个题,因为是在数学分类里的,我想着是不是又公式呀!!一般阶乘的题都有特殊方法,首先想到的是用数组存 高位数 的方法,但是当大于5000就开始慢了,要超时。没办法,看了留言 板,说是 可以用对数 或者斯特林公式。 对数的方法: 可设想n!的结果是不大于10的M次幂的数,即n!<=10^M(10的M次方),则不小于M的最小整数就是 n!的位数,对 该式两边取对数,有 M =log10^n! 即:M = log10^1+log10^2+log10^3...+log10^n 循环求和,就能算得M值,该M是n!的精确位数。当n比较大的时候,这种方法方法需要花费很多的时间。
View Code
#include 
#include
int main(){int n,N,i;double sum;scanf("%d",&N);while (N--){scanf("%d",&n);sum=1;for(i=1;i<=n;i++)sum+=log10((double)i);printf("%d\n",(int)sum);}return 0;}

斯特林公式:

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

其中pi是圆周率,e是自然对数。
1时 要特殊处理:
View Code
#include 
#include
#include
#include
#include
#define e 2.718281828459045#define pi 3.141592653589793239using namespace std;int main (){ int cas,n; scanf("%d",&cas); while (cas --) { scanf("%d",&n); double t = log10(sqrt(2*pi*n)) + n * log10(n/e); printf ("%d\n",(int)t + 1); }return 0;}

 

转载于:https://www.cnblogs.com/zibuyu/archive/2013/03/11/2953486.html

你可能感兴趣的文章
Spring-hibernate整合
查看>>
exit和return的区别
查看>>
Django 相关
查看>>
Python(软件目录结构规范)
查看>>
codeforces水题100道 第二十二题 Codeforces Beta Round #89 (Div. 2) A. String Task (strings)
查看>>
c++||template
查看>>
条件断点 符号断点
查看>>
.net学习之继承、里氏替换原则LSP、虚方法、多态、抽象类、Equals方法、接口、装箱拆箱、字符串------(转)...
查看>>
python的多行注释
查看>>
连接Oracle需要jar包和javadoc文档的下载
查看>>
Dreamweaver cc新版本css单行显示
查看>>
【android】安卓的权限提示及版本相关
查看>>
Java基础教程——网络基础知识
查看>>
Kruskal基础最小生成树
查看>>
【hdu 1429】胜利大逃亡(续)
查看>>
关于收费软件
查看>>
javascript之Style物
查看>>
图的深度优先遍历
查看>>
C# 之 提高WebService性能大数据量网络传输处理
查看>>
Factory Design Pattern
查看>>