2023牛客暑期多校训练营3 B Auspiciousness

#include <bits/stdc++.h>
#define LL long long
#define int LL
#define INF 0x3f3f3f3f
#define PII pair<int,int>
#define kanm7 ios::sync_with_stdio(false), cin.tie(0), cout.tie(0)
#define endl '\n'
using namespace std;
const int N = 310;
int dp[N][N][2];
int inv[N * 2], c[N][N];
int n, mod;

void solve() {
	cin >> n >> mod;
	for(int i = 0; i <= n; i ++) {
		for(int j = 0; j <= n; j ++) 
			{
				dp[i][j][0] = dp[i][j][1] = 0;
			}
	}
	memset(inv, 0, (n * 2 + 5) * 4);
	memset(c, 0, ((n + 3) * (n + 3)) * 4);
	inv[0] = 1;
	for(int i = 1; i <= n * 2; i ++) inv[i] = inv[i - 1] * i % mod;
	for(int i = 0; i <= n; i ++) {
		for(int j = 0; j <= i; j ++) {
			if(!j) c[i][j] = 1;
			else c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % mod;
		}
	}
	int ans = inv[n * 2];
	dp[0][0][0] = 1, dp[0][0][1] = 1;
	for(int i = 0; i <= n; i ++) {
		for(int j = 0; j <= n; j ++) {
			for(int k = 1; k <= i; k ++) {
				dp[i][j][0] += dp[i - k][j][1] * c[n - (i - k)][k];
				dp[i][j][0] %= mod;
			}
			for(int k = 1; k <= j; k ++) {
				dp[i][j][1] += dp[i][j - k][0] * c[n - (j - k)][k];
				dp[i][j][1] %= mod;
			}
			if(i == j && (i == 0 || i == n)) continue;
            // 这里dp[i][j][0] + dp[i][j][1]记录的是可以拿第i + j + 1张
牌的所有情况,但是
            当i = j = 0的情况,在初始的时候
已经算上了,当 i = j = n的时候,已经拿完所有的牌了。
			ans += (dp[i][j][0] + dp[i][j][1]) * inv[n * 2 - i - j];
			ans %= mod;
		}
	}
	cout << ans << endl;
}

signed main() {
	kanm7;
	int T; cin >> T;
	while(T --) {
		solve();
	}
	return 0;
}

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

kanm7

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值