Acwing_888（求组合数Ⅳ）

笨笨sg

于 2024-03-15 22:24:23 发布

阅读量358

点赞数 7

分类专栏： Acwing算法基础课文章标签：组合数高精乘

本文链接：https://blog.csdn.net/a131529/article/details/136751481

版权

Acwing算法基础课专栏收录该内容

77 篇文章 0 订阅

订阅专栏

本文介绍了如何使用C++解决AcWing平台上的888题，涉及高精度计算、素数分解以及组合数的计算，通过计算特定素数在给定范围内的出现次数来解决问题。

摘要由CSDN通过智能技术生成

原题链接：

888. 求组合数 IV - AcWing题库

题解：

思路：

具体推导过程参见：AcWing 888. 求组合数 IV(高精度-素数组合) - AcWing

代码：

#include<bits/stdc++.h>
using namespace std;
const int N = 5010;

int primes[N], cnt;
int sum[N];
bool st[N];

void get_primes(int n)
{
	for (int i = 2;i <= n;i++)
	{
		if (!st[i])primes[cnt++] = i;
		for (int j = 0;primes[j] * i <= n;j++)
		{
			st[primes[j] * i] = true;
			if (i % primes[j] == 0)break;//==0每次漏
		}
	}
}
// 对p的各个<=a的次数算整除下取整倍数
int get(int n, int p)
{
	int res = 0;
	while (n)
	{
		res += n / p;
		n /= p;
	}
	return res;
}
//高精度乘
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;
	}
	// while(C.size()>1 && C.back()==0) C.pop_back();//考虑b==0时才有pop多余的0 ,题目中b>=1, b!=0不需要这行
	return c;
}

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(a - b, p) - get(b, p);//是a-b不是b-a
	}

	vector<int> res;
	res.push_back(1);

	for (int i = 0; i < cnt; i++)
		for (int j = 0; j < sum[i]; j++)//primes[i]的次数
			res = mul(res, primes[i]);

	for (int i = res.size() - 1; i >= 0; i--) printf("%d", res[i]);
	puts("");
}