其他分享
首页 > 其他分享> > [AcWing 197] 阶乘分解

[AcWing 197] 阶乘分解

作者:互联网

image


点击查看代码
#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;
}

  1. 首先用线性筛法筛出所有小于等于 \(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} $
  2. 任意质数必定至少出现一次,因为筛出来的是满足小于等于 \(n\) 的质数

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