AT_abc114_d 题解

Gavin0331

已于 2024-02-05 23:16:03 修改

阅读量376

点赞数 7

文章标签：算法数学建模

于 2024-02-05 22:23:42 首次发布

本文链接：https://blog.csdn.net/Gavin0331/article/details/136048891

版权

思路

题意是求 $N!$ 的因数中有多少个七五数。

把 $75$ 分解一下，即 $75=3 \times 5^2$ 。

如果要让因数 $a$ 有 $75$ 个因数，则有如下几种情况：

1. $a=p_1^{\alpha_1}p_2^{\alpha_2}p_3^{\alpha_3}$ ，且 $(\alpha_1+1)(\alpha_2+1)(\alpha_3+1)=75$ ，因为 $75=3 \times 5^2$ ，所以 $\alpha_1,\alpha_2,\alpha_3$ 为 $2,4,4$ 。

2. $a=p_1^{\alpha_1}p_2^{\alpha_2}$ ，且 $(\alpha_1+1)(\alpha_2+1)=75=3 \times25=5 \times15$ ，所以 $\alpha_1,\alpha_2$ 为 $2,24$ 或者 $4,14$ 。

3. $a=p_1^{\alpha_1}$ ，可知 $\alpha_1=74$ 。

把 $N!$ 分解之后每个因数枚举一下即可。

$code$ ：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 2e5 + 10;
int a[N];
int n;

bool is_prime(int k) {
    if (k <= 1) {
        return 0;
    }
    for (int i = 2; i * i <= k; i++) {
        if (k % i == 0) {
            return 0;
        }
    }
    return 1;
}

int main() {
    cin >> n;
    for (int i = 1; i <= n; i++) {
        if (is_prime(i)) {
            int tmp = 1, cnt = 0;
            while (tmp * i <= n) {
                tmp *= i;
                cnt += n / tmp;
            }
            a[cnt]++;
        }
    }
    for (int i = 100; i >= 1; i--) a[i] += a[i + 1];
    cout << a[74] + a[4] * (a[4] - 1) / 2 * (a[2] - 2) + a[24] * (a[2] - 1) + a[14] * (a[4] - 1);
    return 0;
}