CF1237F Balanced Domino Placements

Solution

显然可以先考虑横着的骨牌，再考虑竖着的骨牌。但是思路卡在了选取横着的骨牌会对竖着的骨牌的相邻对数产生影响。

然而事实上我们只需要换一个统计顺序，先考虑横着的骨牌的列和竖着的骨牌的行，然后再考虑横着的骨牌的行和竖着的骨牌的列（因为单行/列没有双行/列的必须相邻的限制，因此可以简单地通过剩余行/列数统计）。

前面一项可以简单地dpdpdp得到，后面一项组合计数即可。

时间复杂度O(n2)O(n^2)O(n2)。

Code

#include <bits/stdc++.h>using namespace std;template<typename T> inline bool upmin(T &x, T y) { return y < x ? x = y, 1 : 0; }
template<typename T> inline bool upmax(T &x, T y) { return x < y ? x = y, 1 : 0; }#define MP(A,B) make_pair(A,B)
#define PB(A) push_back(A)
#define SIZE(A) ((int)A.size())
#define LEN(A) ((int)A.length())
#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define fi first
#define se secondtypedef long long ll;
typedef unsigned long long ull;
typedef long double lod;
typedef pair<int, int> PR;
typedef vector<int> VI; const lod eps = 1e-11;
const lod pi = acos(-1);
const int mods = 998244353;
const int oo = 1 << 30;
const ll loo = 1ll << 62;
const int MAXN = 4005;
const int INF = 0x3f3f3f3f;//1061109567
/*--------------------------------------------------------------------*/
inline int read() {int f = 1, x = 0; char c = getchar();while (c < '0' || c > '9') { if (c == '-') f = -1; c = getchar(); }while (c >= '0' && c <= '9') { x = (x << 3) + (x << 1) + (c ^ 48); c = getchar(); }return x * f;
}
int tagx[MAXN], tagy[MAXN], fn[MAXN][MAXN], fm[MAXN][MAXN], fac[MAXN], inv[MAXN];
int upd(int x , int y) { return x + y >= mods ? x + y - mods : x + y; }
int C(int x, int y) { return x < y ? 0 :1ll * fac[x] * inv[y] % mods * inv[x - y] % mods; }
int quick_pow(int x, int y) {int ret = 1;for(; y ; y >>= 1) {if (y & 1) ret = 1ll * ret * x % mods;x = 1ll * x * x % mods;}return ret;
}
void Init(int n) {fac[0] = 1;for (int i = 1; i <= n ; ++ i) fac[i] = 1ll * fac[i - 1] * i % mods;inv[n] = quick_pow(fac[n], mods - 2);for (int i = n - 1; i >= 0 ; -- i) inv[i] = 1ll * inv[i + 1] * (i + 1) % mods;
}
signed main() {int n = read(), m = read(), k = read();Init(max(n, m));for (int i = 1; i <= k ; ++ i) tagx[read()] = 1, tagy[read()] = 1, tagx[read()] = 1, tagy[read()] = 1;fn[0][0] = 1;for (int i = 1; i <= n ; ++ i)for (int j = 0; j * 2 <= i; ++ j) fn[i][j] = (i > 1 && !tagx[i] && !tagx[i - 1]) ? upd(fn[i - 1][j], fn[i - 2][j - 1]) : fn[i - 1][j];fm[0][0] = 1;for (int i = 1; i <= m ; ++ i)for (int j = 0; j * 2 <= i; ++ j) fm[i][j] = (i > 1 && !tagy[i] && !tagy[i - 1]) ? upd(fm[i - 1][j], fm[i - 2][j - 1]) : fm[i - 1][j];int rn = 0, rm = 0;for (int i = 1; i <= n ; ++ i) rn += !tagx[i];for (int i = 1; i <= m ; ++ i) rm += !tagy[i];int ans = 0;for (int i = 0; i * 2 <= rn ; ++ i)for (int j = 0; j * 2 <= rm; ++ j) ans = upd(ans, 1ll * fac[i] * fac[j] % mods * C(rm - j * 2, i) % mods * C(rn - i * 2, j) % mods * fn[n][i] % mods * fm[m][j] % mods);printf("%d\n", ans);return 0;
}

CF1237F Balanced Domino Placements(组合计数，dp)相关推荐

CF1237F Balanced Domino Placements（dp+组合计数）
CF1237F Balanced Domino Placements problem solution code problem 题目链接 solution 骨牌横着放会占用一行两列,骨牌竖着放会占用 ...
Educational Codeforces Round 81 (Rated for Div. 2) F.Good Contest \ 洛谷划艇组合计数dp
cf传送门 P3643 [APIO2016]划艇文章目录题意: 思路: 题意: aia_iai在[li,ri][l_i,r_i][li,ri]等概率随机选一个数,求aaa数组不增的概率. 思 ...
bzoj2111,P2606-[ZJOI2010]排列计数【Lucas,组合计数,dp】
正题题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=2111 https://www.luogu.org/problem/P2606 题 ...
Sumsets POJ - 2229（计数dp）
题意: 给一个数,是集合的总数和,集合元素只能为2的次幂数,问这样的集合有多少? 题目: Farmer John commanded his cows to search for different ...
数学知识——组合计数
组合计数文章目录组合计数概述动态规划牡牛和牝牛思路代码隔板法方程的解思路代码序列统计思路代码加法 & 乘法原理加法原理乘法原理车的摆放思路代码容斥原理 ...
解题报告（五）组合计数（ACM / OI）超高质量题解
繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量题解和代码,题目难度不一定按照题号排序,我会在每道题后面加上题目难度指数(1∼51 \sim 51∼5),以模板题难度 11 ...
解题报告（二）E、（BZOJ3513） [MUTC2013] idiots（生成函数 + FFT + 组合计数）
繁凡出品的全新系列:解题报告系列 -- 超高质量算法题单,配套我写的超高质量题解和代码,题目难度不一定按照题号排序,我会在每道题后面加上题目难度指数(1∼51 \sim 51∼5),以模板题难度 11 ...
UVA 10564 计数DP
也是经典的计数DP题,想练练手,故意不写记忆化搜索,改成递推,还是成功了嘞...不过很遗憾一开始WA了,原来是因为判断结束条件写个 n或s为0,应该要一起为0的,搞的我以为自己递推写挫了,又改了一下, ...
【Luogu4921】情侣？给我烧了！（组合计数）
[Luogu4921]情侣?给我烧了!(组合计数) 题面洛谷题解很有意思的一道题目. 直接容斥?怎么样都要一个平方复杂度了. 既然是恰好\(k\)对,那么我们直接来做: 首先枚举\(k\)对人出 ...

CF1237F Balanced Domino Placements(组合计数，dp)

CF1237F Balanced Domino Placements

Solution

Code

CF1237F Balanced Domino Placements(组合计数，dp)相关推荐

最新文章

热门文章