Codeforces#538 C. Trailing Loves (or L'oeufs?)（div 2）[number theory] b进制下n阶乘中尾零个数

传送门

给定两个整数 n 和 b （十进制表示），你的任务是计算b进制下n的阶乘尾部零的个数。

常规的解法在n较小的情况下，可以直接计算出n!的值，将其因式分解，计算出因子中b的个数，则是最终的解

但本题考虑到n数量达到了10的18次方(long long 级别)， 因此若直接计算n的阶乘一定会溢出。因此考虑用质因数分解法来缩小数的规模。

 

首先 对于 n!（） 可以表示为  n! = α*(b^t) ①

其中，t 即为 b 进制下 n! 中尾零个数。

接着考虑以下，

整数k (k >= 2) 都可将其分解为如下形式：

k = p1^q1+ p2^q2 + ... + pm^qm;

因此对于正整数 A = n!， A =  p1^q1* p2^q2 * ... * pm^qm;

对于正整数 B = b (进制),  B = p1^t1 * p2^t2 * ... * pn^tn; (此处的n与阶乘中的n不同，只是用来代指分解出质因数的个数)

即将此处 B代入①式可得 n! 的质因数分解后的表示

即         B = α( p1^t1 * p2^t2 * ... * pn^tn )^t

又因为 A = p1^q1* p2^q2 * ... * pn^qn

令 x = p1^t1 * p2^t2 * ... * pn^tn

此时的任务即计算A最多能整除多少个x， 答案即t的最大值，即尾零的个数

此时，取最大的x时，对应的pi^(ti * t)应不大于 pi^qi 即， ti * t <= qi

因此 t = qi / ti;

同时，为了能保证A 整除 B， 全局最大的 t 应为 min(qi / ti) (1 ≤ i ≤ n);

 

由于n!数目极大， 不容易分解质因数，因此

求 n! 下 素因子 p 的幂 用累除法，因为存在推论：

n!  下 p 的幂 = [ n/p ] + [ n/(p^2) ] + [ n/(p^3) ] + ... +  [ n/(p^k) ]( p^x 中x的个数)

 

总结：

本题考查数论知识，其中主要考查

(1) 将大整数通过质因数分解来减小规模

(2) 将大整数分解质因数的算法(累除法)

(3) 程序中注意处理数组越界问题

#include <iostream>
#include <algorithm>
#include <cmath>
#include <vector>
#include <map>
using namespace std;
typedef long long ll;
vector<ll> prime;
map<ll, int> numprimeB ;
void divideB(ll b) {   //分解质因数算法
	ll len = sqrt(double(b));
	for (ll i = 2; i <= len; i++) { 
		if (b % i == 0)
			prime.push_back(i);
		while (b % i == 0) {
			numprimeB [i]++;
			b /= i;
		}
	}
	if (b != 1) { prime.push_back(b); numprimeB [b]++ ;}
}
ll divideN(ll n, ll p) { //计算 n阶乘中 大小为p的质因子的个数x 即 p^x中x的个数
	ll cnt = 0;
	while (n) {
		cnt += n / p;
		n /= p;
	}
	return cnt;
}
int main() {
	ll n;
	ll b;
	cin >> n >> b;
	divideB(b);
	ll ans = 1LL << 62;
	for (int i = 0; i < prime.size(); i++) {
		ans = min(ans , divideN(n, prime[i])/ numprimeB[prime[i]]);
	}
	cout << ans << endl;
	//system("pause");
}