数论 组合素数

组合素数 - nefu119

题意：求解C(2n,n)恰好被p整除多少次

由组合数公式可以将其转化为：
(2 * n)!/(n! * n!)

则将题目转化为求解分子关于q的幂，减去分子关于q的幂
用公式

• n!的素因子分解中素数p的幂为[n/p]+[n/p^2]+[n/p^3]...0 (p^t <= n)
• 我们首先计算的是2*n可以被P整除的个数，但是包含了(n!*n!)的个数，及包含了n!的两个个数
• 代码如下

#include<iostream>
#include<vector>
using namespace std;

int pr(int n,int p)
{
	int r = 2 * n / p;
	int sum = 0;
	int q = p;
	while (r)
	{
		sum += (2 * n / p - 2*(n / p));
		p *= q;
		r = 2 * r / p;
	}

	return sum;
}
void main()
{
	cout << pr(2, 3) << endl;
	system("pause");

}