ACM 阶乘数位数

描述

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

输入
首行输入n,表示有多少组测试数据(n<10)
随后n行每行输入一组测试数据 N( 0 < N < 1000000 )
输出
对于每个数N,输出N!的(十进制)位数。
样例输入
3
1
3
32000
样例输出
1
1
130271

/*    NYOJ69 阶乘数位长度
 *    方法一:
 *    可设想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比较大的时候,这种方法方法需要花费很多的时间。
 *
 *    方法二:
 *    利用斯特林(Stirling)公式的进行求解。下面是推导得到的公式:
 *    res=(long)( (log10(sqrt(4.0*acos(0.0)*n)) + n*(log10(n)-log10(exp(1.0)))) + 1 );
 *    当n=1的时候,上面的公式不适用,所以要单独处理n=1的情况!
 *    有关斯特林(Stirling)公式及其相关推导,这里就不进行详细描述,有兴趣的话可看这里。
 *    这种方法速度很快就可以得到结果。详细证明如下:
 *    http://episte.math.ntu.edu.tw/articles/mm/mm_17_2_05/index.html
*/
#include<iostream>
#include <cmath>
using namespace std;
int normal(double n)
{
    double x=0;
    while(n)
    {
        x +=log10(n);
        n--;
    }
    return (int)x+1;
}
long stirling(double n)
{
    long x=0;
    if( n ==1 )
        x = 1;
    else
    {
        x = (long)( (log10(sqrt(4.0*acos(0.0)*n)) + n*(log10(n)-log10(exp(1.0)))) + 1 );
    }
    return x;
}
int main()
{
    int n;
    cin>>n;
    while(n--)
    {
        int x;
        cin>>x;
        cout<<stirling(x)<<endl;
    }
    return 0;
}                

时间: 2024-10-14 14:32:30

ACM 阶乘数位数的相关文章

求一个数阶乘的位数

flyfish 2015-8-15 例如 7!=5040 ,7的阶乘结果是4位数(10进制) 求一个数的位数 1 循环方法 int get_digit_loop(int N) { int digit = 0; do { digit ++; } while ((N /= 10) > 0); return digit; } 2 递归方式 int get_digit_recursion(int N) { int digit = 0; digit = N < 10 ? 1 : 1 + get_digi

计算一个大数n的阶乘的位数宽度(十进制)转载

计算一个大数n的阶乘的位数宽度(十进制)(log i累加法 )转载 输入: 每行输入1个正整数n, (0<n<1000 000) 输出: 对于每个n,输出n!的(十进制)位数. 分析: 这道题采用蛮力法.根据定义,直接求解! 所谓n!的十进制位数,就是 log(n)+1, 根据数学公式有:n!=1*2*3*.....*n; lg(n!)=lg(2)+......lg(n); 代码: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22

大数阶乘的位数和精确值计算【转】

来源:http://www.cnblogs.com/stonehat/p/3603267.html 在此,顶礼膜拜一下原文作者呵呵 我们知道整数n的位数的计算方法为:log10(n)+1 故n!的位数为log10(n!)+1 如果要求出n!的具体值,对很大的n(例如n=1000000)来说,计算会很慢,如果仅仅是求阶乘的位数,可以用斯特林(Stirling)公式求解 斯特林(Stirling)公式: 于是求n!的位数就是求log10((2*PI*n)^1/2*(n/e)^n)+1 即 1/2*l

求一个数阶乘后位数问题

问题:In many applications very large integers numbers are required. Some of these applications are using keys for secure transmission of data, encryption, etc. In this problem you are given a number, you have to determine the number of digits in the fa

ACM——完数

完数 时间限制(普通/Java):1000MS/3000MS          运行内存限制:65536KByte 总提交:1930            测试通过:413 描述 自然数中,完数寥若晨星,请在从1到某个整数范围中打印出所有的完数来.所谓“完数”是指一个数恰好等于它的所有不同因子之和.例如,6是完数,因为6=1+2+3.而24不是完数,因为24≠1+2+3+4+6+8+12=36. 输入 输入数据中含有一些整数n(1<n<10000). 输出 对于每个整数n,输出所有不大于n的完

计算一个整数(N)的阶乘的位数

对于一个正整数N, 计算N! 的位数.例如N=4, 4!=24,那么位数就是2. 直接计算N!的数值,然后再去数位数,这个很难,因为N!很有可能超过int(32bit) 或long(64bit)的表达范围. 换一种思路,假设要求的位数为x, 那么一定满足 10^(x-1) <=N!<10^x.两边取10为底的对数,得到x-1<=log10(N!)<x.最终x 取 int(log10(N!))+1. log10(N!) = log10(N)+log10(N-1)+log10(N-2)

组合数(阶乘数质因子分解)

C(n,m)   表示组合数,n>=m>=0 以下适用范围: n<=1e6(or 1e7...) 爆龙龙的答案需取模,允许取合数模. 时间复杂度 线性筛略大一点点 大概还是nlon(n) C(n,m)=n!/(m!*(n-m)!) 举例说明一下为什么可以通过下面代码去计算阶乘数的所有质因子各自的数目. void factor_jc(int n)//分解n!的素数因子并记录个数 { //fac[i]=(prim[i]这个素数因子有几个) int i,up=n;for(i=1;i<=p

计算一个大数n的阶乘的位数宽度(十进制)(log i累加法 )

输入: 每行输入1个正整数n, (0<n<1000 000) 输出: 对于每个n,输出n!的(十进制)位数. 分析: 这道题采用蛮力法.根据定义,直接求解! 所谓n!的十进制位数,就是 log(n)+1, 根据数学公式有:n!=1*2*3*.....*n; lg(n!)=lg(2)+......lg(n); 代码: //输入一个数字n,请你计算该数的阶乘的十进制数的位数宽度 //比如:3!=6, 则宽度为1 //样例数据: //n=3 输出1 //n=32000 输出130271 //n=10

ACM 阶乘的0

阶乘的0 时间限制:3000 ms  |  内存限制:65535 KB 难度:3 描述 计算n!的十进制表示最后有多少个0 输入 第一行输入一个整数N表示测试数据的组数(1<=N<=100)每组测试数据占一行,都只有一个整数M(0<=M<=10000000) 输出 输出M的阶乘的十进制表示中最后0的个数比如5!=120则最后的0的个数为1 样例输入 6 3 60 100 1024 23456 8735373 样例输出 0 14 24 253 5861 2183837 令f(x)表示