NYOJ-ACM-数的长度
来源:互联网 发布:sql with as等价 编辑:程序博客网 时间:2024/05/29 09:11
数的长度
时间限制: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
思路:
这道题的难点就是要用斯特林公式来求解n!的位数。下面是求解的公式和注意。
求n!的位数:
利用斯特林(Stirling)公式的进行求解。下面是推导得到的公式:
res=(long)( (log10(sqrt(4.0*acos(0.0)*n)) + n*(log10(n)-log10(exp(1.0)))) + 1 );
当n=1的时候,上面的公式不适用,所以要单独处理n=1的情况!
这种方法速度很快就可以得到结果。下面是我的代码,可以参考。
res=(long)( (log10(sqrt(4.0*acos(0.0)*n)) + n*(log10(n)-log10(exp(1.0)))) + 1 );
当n=1的时候,上面的公式不适用,所以要单独处理n=1的情况!
这种方法速度很快就可以得到结果。下面是我的代码,可以参考。
代码:
import java.util.Scanner;public class Main {public static void main(String[] args) {// TODO Auto-generated method stubScanner sc = new Scanner(System.in);int a = sc.nextInt();long[] b = new long[a];//定义一个long型的存放数据的数组long[] sum = new long[a];//定义一个long型的存放结果的数组for(int i=0;i<a;i++){b[i] = sc.nextInt();}for(int i=0;i<a;i++){//斯特林(Stirling)公式if(b[i] == 1){sum[i] = 1;}else{sum[i]=(long)( (Math.log10(Math.sqrt(4.0*Math.acos(0.0)*b[i])) + b[i]*(Math.log10(b[i])-Math.log10(Math.exp(1.0)))) + 1 );}}for(int i=0;i<a;i++){//输出结果System.out.println(sum[i]);}}}
0 0
- NYOJ-ACM-数的长度
- NYOJ - 数的长度
- NYOJ-数的长度
- NYOJ - 数的长度
- NYOJ 数的长度
- acm-数的长度
- NYOJ 69 数的长度
- nyoj-69-数的长度
- NYOJ 69 数的长度
- NYOJ 69 数的长度
- NYOJ 69【数的长度】
- NYOj 69 数的长度
- NYOJ-69数的长度
- NYOJ:69-数的长度
- NYOJ 69 数的长度
- NYOJ 69 数的长度
- NYOJ 题目69 数的长度
- nyoj 题目69 数的长度
- Linux Ubuntu 16.04 初次安装使用总结zzz
- Searching for Patterns の Set 1 (Naive Pattern Searching)
- 项目中的web.xml基础理解
- 使用DBUtils+DBCP2+pool2+logging第三方工具类,创建连接池
- poj1661 最短路
- NYOJ-ACM-数的长度
- 路由
- Golang一键启动程序(windows平台 部署 微服务)
- [Unity框架]PureMVC在unity中的简单使用
- 导入项目到SVN server 中
- 十进制和二进制
- 建议1 正确操作字符串
- matplotlib的基本用法(十)——绘制Image
- 技术转型产品学习(交互设计)——第十二天