首页 > 其他分享 >ABC180F

ABC180F

时间:2022-11-11 19:55:47浏览次数:76  
标签:int res 1ll ABC180F add fac mod

由于点的度数最大为 \(2\),于是这张图由链,孤立点,大小至少为 \(2\) 的环组成,为了方便把孤立点也看成链。

考虑容斥掉第三个条件,最大连通块大小恰好为 \(L\) 的方案数即为最大连通块大小至多为 \(L\) 的方案数减去大小至多为 \(L-1\) 的方案数。

考虑 DP,设 \(f(i,j)\) 表示 \(i\) 个点,\(j\) 条边的方案数,每次考虑加进去一条链或一个环,但是会计算重复,因为加入的顺序没有影响,借鉴 ABC226F 的思想,强制选剩下的点中编号最小的点即可。因为是图,于是不仅要考虑环的旋转同构,还要考虑翻转同构,但是注意大小为 \(1\) 的链和大小为 \(2\) 的环是不用考虑翻转同构的。

时间复杂度 \(\mathcal O(NML)\)。

Code:

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 305, mod = 1e9 + 7, inv2 = (mod + 1) / 2;
int n, m, L;
int fac[N], inv[N];
int f[N][N];

void add(int &a, int b) {
	a += b;
	if (a >= mod) a -= mod;
}

int qpow(int x, int y) {
	int res = 1;
	while (y) {
		if (y & 1) res = 1ll * res * x % mod;
		x = 1ll * x * x % mod;
		y >>= 1;
	}
	return res;
}

void init(int n) {
	fac[0] = 1;
	for (int i = 1; i <= n; ++i) fac[i] = 1ll * fac[i - 1] * i % mod;
	inv[n] = qpow(fac[n], mod - 2);
	for (int i = n - 1; ~i; --i) inv[i] = 1ll * inv[i + 1] * (i + 1) % mod;
}

int C(int n, int m) {
	if (n < 0 || m < 0 || n < m) return 0;
	return 1ll * fac[n] * inv[n - m] % mod * inv[m] % mod;
}

int solve(int lim) {
	memset(f, 0, sizeof f), f[0][0] = 1;
	for (int i = 0; i <= n; ++i)
		for (int j = 0; j <= m; ++j) {
			for (int k = 1; k <= lim; ++k) {
				if (i + k > n || j + k - 1 > m) continue;
				add(f[i + k][j + k - 1], 1ll * f[i][j] * C(n - i - 1, k - 1) % mod * fac[k] % mod * (k > 1 ? inv2 : 1) % mod);
			}
			for (int k = 2; k <= lim; ++k) {
				if (i + k > n || j + k > m) continue;
				add(f[i + k][j + k], 1ll * f[i][j] * C(n - i - 1, k - 1) % mod * fac[k - 1] % mod * (k > 2 ? inv2 : 1) % mod);
			}
		}
	return f[n][m];
}

int main() {
	scanf("%d%d%d", &n, &m, &L);
	init(n);
	printf("%d", (solve(L) - solve(L - 1) + mod) % mod);
	return 0;
}

标签:int,res,1ll,ABC180F,add,fac,mod
From: https://www.cnblogs.com/Kobe303/p/16881583.html

相关文章