问题引入和简单分析
对于质数,大家应该都不陌生,现在有一个问题,如果给定我们正整数 n n n和 p p p,如何求出 n ! n! n!里面有多少个质因子 p p p
- 可以仔细考虑一下,根据唯一分解定理,任何一个大于1的正整数都能够写成若干个质数乘积的形式,所以如果 n ! ∣ p n!|p n!∣p,我们可以一直除下去直到不能整除,这时候能除多少次就有多少个质因子 p p p
- 上述思路是容易的,但是不能忽略一个重要问题,C++能够表示的数据并不是无限大的, n ! n! n!又是一个非常大的数,那么在不能直接得到 n ! n! n!的时候,能不能直接求出这个答案呢?
求解
因为 n ! = 1 × 2 × 3 × . . . × n n!=1\times2\times3\times...\times n n!=1×2×3×...×n,可以分别求 1 , 2 , 3 , . . . n 1,2,3,...n 1,2,3,...n的质因子 p p p的数量再相加,这样时间复杂度大约是 O ( n l o g ( n ) ) O(nlog(n)) O(nlog(n)),原理显然,是一种比较暴力的方法
- 考虑优化,设 n ! n! n!以内质因子 p p p的数量为 a n s ans ans,可以这样想, n ! n! n!以内可能是 p p p倍数的情况有哪些呢?显然有 p , 2 p , 3 p , 4 p , . . . , k p p,2p,3p,4p,...,kp p,2p,3p,4p,...,kp
- 如果 k < p k\lt p k<p,那么显然 a n s = k = ⌊ n p ⌋ ans=k=\lfloor\frac{n}{p}\rfloor ans=k=⌊pn⌋
- 如果 k = p k=p k=p,由于最后一项 p 2 p^2 p2有两个 p p p,所以 a n s = ⌊ n p ⌋ + 1 ans=\lfloor\frac{n}{p}\rfloor+1 ans=⌊pn⌋+1
- 如果 p < k < p 2 p\lt k\lt p^2 p<k<p2, n ! n! n!以内可能是 p p p的倍数的情况变为 p , 2 p , 3 p , . . . , p 2 , p 2 + p , p 2 + 2 p , . . . , 2 p 2 , . . . , 3 p 2 , . . . , k p p,2p,3p,...,p^2,p^2+p,p^2+2p,...,2p^2,...,3p^2,...,kp p,2p,3p,...,p2,p2+p,p2+2p,...,2p2,...,3p2,...,kp这个时候 a n s = k + ⌊ n p 2 ⌋ ans=k+\lfloor\frac{n}{p^2}\rfloor ans=k+⌊p2n⌋。解释一下这个式子, k k k的意思是这些一共有 k k k项, k = ⌊ n p ⌋ k=\lfloor\frac{n}{p}\rfloor k=⌊pn⌋,但是这里面还有 p 2 , 2 p 2 , 3 p 2 . . . p^2,2p^2,3p^2... p2,2p2,3p2...这样的数,这样的数一共有 ⌊ n p 2 ⌋ \lfloor\frac{n}{p^2}\rfloor ⌊p2n⌋个,这两部分和是 a n s ans ans
- 同理也可以推出
k
k
k更大的情况,解法同理,略去,总结成一个式子就是
a
n
s
=
⌊
n
p
⌋
+
⌊
n
p
2
⌋
+
⌊
n
p
3
⌋
+
.
.
.
+
0
ans=\lfloor\frac{n}{p}\rfloor+\lfloor\frac{n}{p^2}\rfloor+\lfloor\frac{n}{p^3}\rfloor+...+0
ans=⌊pn⌋+⌊p2n⌋+⌊p3n⌋+...+0
也就是
int solve(int n, int p){//计算n!所含质因子数量
int ans = 0;
while(n){
ans += n / p;
n /= p;
}
return ans;
}