数学问题合集

求n！首位数

---

使用斯特林公式得到近似结果：

$$\lim_{n\to+\infty} {n!\over{\sqrt {2 \pi n}} ({n\over e})^n} = 1$$

之后用10的幂计算取首位即可，注意数字较小时特判，代码如下：

#include <stdio.h>
#include <math.h>

const double PI =3.14159265358979;
const double E =2.718281828;

int main()
{
    int n,fn;
    double log_n_fac;
    while (scanf("%d", &n) != EOF){
        log_n_fac = 0.5 * log10(2 * PI *(double)n) + (double)n * log10((double)n / E);
        log_n_fac -=(int)log_n_fac;
        fn = pow(10, log_n_fac);//Stirling's approximation
        switch(n){
            case 0:printf("1\n");break;
            case 1:printf("1\n");break;
            case 2:printf("2\n");break;
            case 3:printf("6\n");break;
            case 7:printf("5\n");break;
            case 8:printf("4\n");break;
            default:printf("%d\n", fn);
        }
    }
    return 0;
}

求n^n首位

---

方法类似，代码：

#include <stdio.h>
#include <math.h>

int main()
{
    int n;
    scanf("%d",&n);
    while(n != 0){
        printf("%d\n",(int)pow(10,n*log10(n)-(int)(n*log10(n))));
        scanf("%d",&n);
    }
    return 0;
}

整数质因子分解

---

从小到大枚举因数，如果这个因数不是素数（如15），那么之前枚举过的素数已经消耗掉了这个因数（如3和5），因此n%15不会为0。所以这题并不需要筛选素数。

#include <stdio.h>

int n;

void solve(){
    int i;
    int m = n;
    for (i = 2; i <= n; i++){
        int cnt = 0;
        if (m % i) continue;
        while (m % i == 0){
            m /= i;
            cnt++;
        }
        printf("(%d,%d)", i, cnt);
        if (m == 1) break;
    }
    printf("\n");
}

int main(void)
{
    int t,i;
    scanf("%d", &t);
    for (i = 0; i < t; i++){
        scanf("%d", &n);
        solve();
    }
    return 0;
}

n！右端的0的个数

---

令 n!=$5^m 2 ^p x$， m,p,x为非负整数。同时由于n!中2的倍数远多于5的倍数，即p>m，所以m即为答案。要得到m，首先找从1到n中5的倍数，即n/5，而25=$5^2$，一个数就含有两个5在内，因此第二轮我们找从1到n/5中5的倍数。由于第一轮中已经被拿掉了一个5，第二轮25就只需再拿出一个5即可，n/25。第三轮则n/125，第i论${n \over 5^i}$，以此类推：

#include <stdio.h>

int main()
{
    int t,i,n,m,z;
    scanf("%d", &t);
    for (i = 0; i < t; i++){
        scanf("%d", &n);
        m = 5;z = 0;
        while (n >= m){
            z += n / m;
            m *= 5;
        }
        printf("case #%d:\n%d\n", i, z);
    }
    return 0;
}