[AcWing 197] 阶乘分解

2022-08-08 00:04:26 阅读：166 来源： 互联网

标签：frac int cdot 197 次方倍数阶乘质数 AcWing

点击查看代码

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;

const int N = 1e6 + 10;

int n;
vector<int> primes;
bool st[N];

void get_primes(int x)
{
    for (int i = 2; i <= x; i ++) {
        if (!st[i])
            primes.push_back(i);
        for (auto p : primes) {
            if (p * i > x)
                break;
            st[p * i] = true;
            if (i % p == 0)
                break;
        }
    }
}

void solve()
{
    cin >> n;
    get_primes(n);
    for (auto p : primes) {
        int s = 0;
        for (int i = n; i; i /= p)
            s += i / p;
        cout << p << ' ' << s << endl;
    }
}

signed main()
{
    ios::sync_with_stdio(false);
    cin.tie(nullptr);

    solve();

    return 0;
}

首先用线性筛法筛出所有小于等于 $n$ 的质数，对于一个质数 $p$，可以求出在 $n!$ 中出现多少次：
① $p$ 的倍数在 $1$ ~ $n$ 中出现的次数，即 $\frac{n}{p}$
② $p^2$ 的倍数在 $1$ ~ $n$ 中出现的次数，即 $\frac{n}{p^2}$
$\cdots$
③ $p^k$ 的倍数在 $1$ ~ $n$ 中出现的次数，即 $\frac{n}{p^k}$
由于求的是次方数，$p$ 的倍数是 $1$ 次方，$p^2$ 的倍数是 $2$ 次方，$\cdots$，$p^k$ 的倍数是 $k$ 次方，对于 $p^2$，在计算 $p$ 的倍数时会加 $1$ 次（$p^2 = p \cdot p$），在计算 $p^2$ 的倍数时会加 $1$ 次（$p^2 = 1 \cdot p^2$），总共加了 $2$ 次，刚好是次方数，对于 $p^3$，在计算 $p$ 的倍数时会加 $1$ 次（$p^3 = p^2 \cdot p$），在计算 $p^2$ 的倍数时会加 $1$ 次（$p^2 = p \cdot p^2$），在计算 $p^3$ 的倍数时会加 $1$ 次（$p^3 = 1 \cdot p^3$），总共加了 $3$ 次，刚好是次方数
即对于 $p^k$，会在 $p, p^2, \cdots ,p^k$ 各加 $1$ 次，总共 $k$ 次，刚好是次方数
故 $p$ 的次方数为 $s = \frac{n}{p} + \frac{n}{p^2} + \cdots + \frac{n}{p^k} $
任意质数必定至少出现一次，因为筛出来的是满足小于等于 $n$ 的质数

标签：frac,int,cdot,197,次方,倍数,阶乘,质数,AcWing
来源： https://www.cnblogs.com/wKingYu/p/16560266.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

[AcWing 197] 阶乘分解