【集训Day2】cinema ticket

SSL_GYX

于 2022-07-17 21:37:36 发布

阅读量139

点赞数

CC 4.0 BY-SA版权

分类专栏：数论/数学狗妈OJ 文章标签： c++

本文链接：https://blog.youkuaiyun.com/SSL_guyixin/article/details/125837393

狗妈OJ 同时被 2 个专栏收录

16 篇文章

订阅专栏

数论/数学

14 篇文章

订阅专栏

本文介绍了一种计算卡特兰数特定项的高精度算法，通过递推公式推导出答案，并使用C++实现，适用于需要精确计算的场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

cinema ticket

在这里插入图片描述

解题思路

易得递推式为 $f_{i,j}=f_{i-1,j}+f_{i,j-1}$ 。
已知当 $i == j$ ， $f_{i,j}$ 为卡特兰数的第 $i$ 项。
可以推得 $ans=(n+1−m)(n+m)!m!(n+1)!ans=\frac{(n+1-m)(n+m)!}{m!(n+1)!}$ 。
需要高精度。

code

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

int n,m;

struct num{
	int s[10010];
	num operator *(int x)
	{
		num t={0};
		int lst=0;
		for(int i=10000;i;i--)
		{
			t.s[i]=s[i]*x+lst;
			lst=t.s[i]/10;
			t.s[i]%=10;
		}
		return t;
	}
	num operator /(int x)
	{
		num t={0};
		int lst=0;
		for(int i=1;i<=10000;i++)
			lst=lst*10+s[i],t.s[i]=lst/x,lst%=x;
		return t;
	}
	num operator -(num x)
	{
		num t={0};
		for(int i=10000;i;i--)
		{
			if(s[i]>=x.s[i])
				t.s[i]=s[i]-x.s[i];
			else
				t.s[i]=s[i]+10-x.s[i],s[i-1]--;
		}
		return t;
	}
}a,b;

int main()
{
	cin>>n>>m;
	int i,j;
	b.s[10000]=1;
	for(i=1,j=n+m;i<m;i++,j--)
		b=b*j,b=b/i;
	a=b*j,a=a/i;
	num ans=a-b;
	i=0;
	while(!ans.s[i]) i++;
	while(i<=10000) cout<<ans.s[i++];
}