zoukankan      html  css  js  c++  java
  • 求n^n和n!的最左边的数字

    n^n 的位数 k = [lg(n^n)]+1=[n*lg(n)]+1;

    最左边的数 x = n^n/10^(k-1);

    取对数:lg(x) = n * lg(n) - (k-1) = (n*lg(n) - [n*lg(n)]);

    最左边的数:[x] = [10^lg(x)] = [10^(n*lg(n) - [n*lg(n)])];

    #include<iostream>
    #include<cmath>
    #include<cstdio>
    using namespace std;
    int main() {
        int n, t;
        scanf("%d", &t);
        while (t--) {
            scanf("%d", &n);
            double a = n * log10(n);
            printf("%d
    ", (int)pow(10.0, a-(long long)a)); 
        }
    }
    Stirling公式:
    根据Stirling公式n! 与√(2πn) * (n/e)^n的值十分接近 
    n! 的位数 k = [lg(n!)]+1;
    最左边的数 x = n!/10^(k-1);
    取对数:lg(x) =  lg(n!) - (k-1) ;
    lg(n!)  = lg(√(2πn) * (n/e)^n) = 1/2*(lg(2π) + lg(n)) + n * (lg(n) - lg(e));

    代入[x] = [10^lg(x)] = [10^(lg(n!) - [lg(n!)])]即可。

    #include<iostream>
    #include<cmath>
    #include<cstdio>
    using namespace std;
    const double PI = 3.1415926;
    const double e = 2.718;
    int main() { int n, t; scanf("%d", &t); while (t--) { scanf("%d", &n); double b = (log10(2 * PI) + log10(n)) / 2 + n * (log10(n) - log10(e)); printf("%d ", (int)pow(10.0, b-(long long)b)); } }
  • 相关阅读:
    P4047 部落划分
    P1440 求m区间的最小值
    P2880 平衡的阵容
    P2700 逐个击破
    P2814 家谱 map模版
    P4403 秦腾与教学评估
    无油无糖低脂酸奶芒果蛋糕
    紫薯铜锣烧
    Spring In Action ③
    Spring In Action ②
  • 原文地址:https://www.cnblogs.com/a863886199/p/7798313.html
Copyright © 2011-2022 走看看