标签：get int res 网格 vector ans 卡特兰 size

原题链接

卡特兰数公式

变形一下

用分解质因数求组合数解

#include<iostream>
#include<cmath>
#include<vector>

using namespace std;

typedef long long LL;
const int N = 25050;
int primes[N], cnt;
int sum[N];
bool st[N];
int x, y;

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

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

void get_sum(int a, int b)
{
	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 > mul(vector<int> &a, int b)
{
	vector<int> c;
	int t = 0;
	for (int i = 0; i < a.size() || t; i++)
	{
		if (i < a.size())t += a[i] * b;
		c.push_back(t % 10);
		t /= 10;
	}
	return c;
}

void get_ans(vector<int> &res, int a, int b)
{
	get_sum(a, b);
	for (int i = 0; i < cnt; i++)
	{
		for (int j = 0; j < sum[i]; j++)
		{
			res = mul(res, primes[i]);
		}
	}
}
vector<int> sub(vector<int>& a, vector<int>& b) {
	
	vector<int> c;
	int t = 0;//进位
	for (int i = 0; i < a.size(); i++)
	{
		t = a[i] - t;
		if (i < b.size())t = t - b[i];
		c.push_back((t + 10) % 10);
		if (t < 0)t = 1;
		else t = 0;
	}
	while (c.size() > 1 && c.back() == 0)c.pop_back();//去除前导 0
	return c;
}

int main()
{
	cin >> x >> y;
	
	int a = x + y, b = x, c = y - 1;
	
	get_primes(a);
	vector<int> res, ans;
	res.push_back(1), ans.push_back(1);
		
	get_ans(res, a, b);
	get_ans(ans, a, c);
	
	ans = sub(res, ans);
	
	for (int i = ans.size() - 1; i >= 0; i--)cout << ans[i];
	cout << endl;
		
	return 0;
}

标签：get,int,res,网格,vector,ans,卡特兰,size
来源： https://blog.csdn.net/qq_56427460/article/details/122389947

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。