求组合数4

题目

思路

这次ab数字都很大，需要用到高精度乘法，为了方便计算，采用一下计算步骤

1.分解质因数（Cab），每个质数p的次数等于a中p的次数减去b中p的次数减去（a - b）中p的次数(因为是除法),利用get函数求出次数

2.高精度乘法算出

​​​​​​​

代码

#include<iostream>
#include<vector>

using namespace std;

const int N = 5010;

bool st[N];
int primes[N], sum[N], cnt;//分解质因数，质数和其指数 

void get_primes(int a)
{
	for(int i = 2; i <= a; i ++)
	{
		if(!st[i]) primes[cnt ++] = i;	
		for(int j = 0; primes[j] <= a / i; j ++)
		{
			st[primes[j] * i] = true;
			if(i % primes[j] == 0) break;
		}	
	} 
}

vector<int> mul(vector<int> a, int b)
{
	vector<int> c;
	int t = 0;
	for(int i = 0; i < a.size(); i ++)
	{
		t += a[i] * b;
		c.push_back(t % 10);
		t /= 10; 
	}
	
	while(t)
	{
		c.push_back(t % 10);
		t /= 10;
	}
	
	return c;
}

int get(int a, int p)
{
	int res = 0;
	while(a)
	{
		res += a / p;
		a /= p;
	}
	return res;
}

int main()
{
	int a, b;
	cin >> a >> b;
	
	get_primes(a);
	
	//质因数的指数 
	for(int i = 0; i < cnt; i ++)
	{
		int p = primes[i];
		sum[i] = get(a, p) - get(b, p) - get(a - b, p);
	}
	
	vector<int> res;
	res.push_back(1);
	
	for(int i = 0; i < cnt; i ++)
		for(int j = 0; j < sum[i]; j ++)
			res = mul(res, primes[i]);
			
	for(int i = res.size() - 1; i >= 0; i --) cout << res[i]; 
	
	return 0;
}