[TC_SRM_466]DrawingBlackCrosses
[TC_SRM_466]DrawingBlackCrosses
试题描述
\(n \times m\)(\(n, m \le 20\))的棋盘
其中至多有 \(8\) 个格子为黑色,其他格子为白色
每次可以选一个白格子把它所在的行、列包括它本身变成黑色
求把棋盘全变成黑色的操作方案数
输入
传给你一个 string[] 类型的参数
输出
返回一个整数表示答案
输入示例
{"B..B", "B.B.", "...B", "BB.B", "...."}
输出示例
Returns: 324
数据规模及约定
见“试题描述”
题解
由于黑格子不超过 \(8\) 个,并且每次操作都是整行整列变黑,所以可以任意交换行列,不会影响最终结果。
于是我们可以把它交换成左上角最多 \(8 \times 8\) 区域内有黑格子,剩下一个反 \(L\) 型的全白区域。
对于全白区域,我们只需要关心它有几行几列就行了;而对于有黑格子的地方就需要用状压。
于是令 \(f(s_x, s_y, i, j)\) 表示对于有黑格子的区域行覆盖的集合是 \(s_x\),列覆盖集合为 \(s_y\),全白区域覆盖了 \(i\) 行 \(j\) 列,然后转移显然。
最后累计答案时需要开个数组模拟一下确定哪些状态是最终状态。(好像也可以看哪些状态没有后继状态吧,这样快一点,但是我不管了。这题搞了我一个上午,弄得我生活不能自理)
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <algorithm>
#include <vector>
using namespace std;
#define rep(i, s, t) for(int i = (s); i <= (t); i++)
#define dwn(i, s, t) for(int i = (s); i >= (t); i--)
const int maxspace = 12010101, maxn = 21, MOD = 1000000007;
#define LL long long
int n, m, bn, bm, sn, sm, _x[maxn*maxn], _y[maxn*maxn], cb, xnum[maxn*maxn], ynum[maxn*maxn];
char Map[maxn][maxn];
bool g[maxn][maxn], tg[maxn][maxn];
int F[maxspace];
#define f(si, sj, i, j) F[(si)*(sm+1)*(n-bn+1)*(m-bm+1)+(sj)*(n-bn+1)*(m-bm+1)+(i)*(m-bm+1)+(j)]
class DrawingBlackCrosses {
public:
int count(vector <string> field) {
n = field.size(); m = field[0].length();
rep(i, 0, n - 1) strcpy(Map[i], field[i].c_str());
rep(i, 0, n - 1) rep(j, 0, m - 1) if(Map[i][j] == 'B') {
_x[cb] = i; _y[cb++] = j;
xnum[bn++] = i; ynum[bm++] = j;
}
sort(xnum, xnum + bn); sort(ynum, ynum + bm);
bn = unique(xnum, xnum + bn) - xnum; bm = unique(ynum, ynum + bm) - ynum;
sn = (1 << bn) - 1; sm = (1 << bm) - 1;
rep(i, 0, cb - 1) {
_x[i] = lower_bound(xnum, xnum + bn, _x[i]) - xnum;
_y[i] = lower_bound(ynum, ynum + bm, _y[i]) - ynum;
g[_x[i]][_y[i]] = 1;
}
int sx = 0, sy = 0, fullx = 0, fully = 0; // all_black
rep(i, 0, n - 1) {
bool ab = 1;
rep(j, 0, m - 1) if(!g[i][j]){ ab = 0; break; }
if(ab && i < bn) sx |= 1 << i;
ab = 1;
rep(j, 0, bm - 1) if(!g[i][j]){ ab = 0; break; }
if(i < bn) fullx |= (int)ab << i;
}
rep(j, 0, m - 1) {
bool ab = 1;
rep(i, 0, n - 1) if(!g[i][j]){ ab = 0; break; }
if(ab && j < bm) sy |= 1 << j;
ab = 1;
rep(i, 0, bn - 1) if(!g[i][j]){ ab = 0; break; }
if(j < bm) fully |= (int)ab << j;
}
f(sx, sy, n - bn, m - bm) = 1;
rep(si, 0, sn) rep(sj, 0, sm) dwn(i, n - bn, 0) dwn(j, m - bm, 0) if(f(si, sj, i, j)) {
int now = f(si, sj, i, j), tsi, tsj, ti, tj;
rep(ini, 0, bn) rep(inj, 0, bm) {
if(ini < bn) tsi = si | (1 << ini), ti = i;
else tsi = si, ti = i - 1;
if(inj < bm) tsj = sj | (1 << inj), tj = j;
else tsj = sj, tj = j - 1;
if(ini < bn && inj < bm && !g[ini][inj] && !(si >> ini & 1) && !(sj >> inj & 1)) (f(tsi, tsj, ti, tj) += now) %= MOD;
if(ini < bn && inj == bm && !(si >> ini & 1) && j > 0) (f(tsi, tsj, ti, tj) += (LL)now * j % MOD) %= MOD;
if(ini == bn && inj < bm && i > 0 && !(sj >> inj & 1)) (f(tsi, tsj, ti, tj) += (LL)now * i % MOD) %= MOD;
if(ini == bn && inj == bm && i > 0 && j > 0) (f(tsi, tsj, ti, tj) += (LL)now * i % MOD * j % MOD) %= MOD;
}
}
int ans = 0;
rep(si, 0, sn) rep(sj, 0, sm) rep(i, 0, n - bn) rep(j, 0, m - bm) if(f(si, sj, i, j)) {
memcpy(tg, g, sizeof(g));
rep(x, 0, bn - 1) if(si >> x & 1)
rep(y, 0, m - 1) tg[x][y] = 1;
rep(y, 0, bm - 1) if(sj >> y & 1)
rep(x, 0, n - 1) tg[x][y] = 1;
dwn(x, n - 1, n - (n - bn - i)) rep(y, 0, m - 1) tg[x][y] = 1;
dwn(y, m - 1, m - (m - bm - j)) rep(x, 0, n - 1) tg[x][y] = 1;
bool all1 = 1;
rep(x, 0, n - 1) rep(y, 0, m - 1) if(!tg[x][y]){ all1 = 0; break; }
if(all1) (ans += f(si, sj, i, j)) %= MOD;
}
return ans;
}
};
[TC_SRM_466]DrawingBlackCrosses的更多相关文章
随机推荐
- 进入Windows之前发出警告
实现效果: 知识运用: 通过注册表中HKLM:\SOFTWARE\Microsoft\Windows NT\CurrentVersion\Winlogon\子键下的LegalNoticeCaption ...
- 制作带复选框的ListView控件
实现效果: 知识运用 ListView控件的GridLines //设置是否在ListView控件中显示网格线 public bool GridLines{get;set} 和CheckBoxes ...
- 【转】mongoDB 学习笔记纯干货(mongoose、增删改查、聚合、索引、连接、备份与恢复、监控等等)
mongoDB 学习笔记纯干货(mongoose.增删改查.聚合.索引.连接.备份与恢复.监控等等) http://www.cnblogs.com/bxm0927/p/7159556.html
- CentOS6.5-DHCP配置
原文链接主机配置: 1.查看主机名:hostname 2.修改主机名:vi /etc/sysconfig/network NETWORKING=yes HOSTNAME=master.303kfy.n ...
- configParser模块详谈
前言 使用配置文件来灵活的配置一些参数是一件很常见的事情,配置文件的解析并不复杂,在python里更是如此,在官方发布的库中就包含有做这件事情的库,那就是configParser configPars ...
- UML类图关系模式(C++代码说明)
在UML类图中的关系模式主要有以下几种: 泛化(Generalization), 实现(Realization), 关联(Association), 聚合(Aggregation), 依赖(Depe ...
- Network of Schools POJ - 1236 (强联通)
一些学校连接到了一个计算机网络.网络中的学校间有如下约定:每个学校维护一个列表,当该学校收到软件或信息后将会转发给列表中的所有学校(也就是接收方列表).需要注意的是如果B学校在A学校的接收方列表中,A ...
- nova boot添加volume_type参数支持
早前由于添加了全SSD的高性能Ceph集群,区别于现有的HDD集群,在OpenStack端需要能够选择使用两种集群.Cinder配置多Ceph后端的文档早已整理,整理文件夹时发现这篇为nova boo ...
- PowerShell批量配置VM端点
我们可以通过PowerShell脚本批量添加VM端点.请您参考以下方案. 准备工作 – PowerShell连接China Azure 1. 从官网下载页面,下载并安装Windows Azure Po ...
- P1194 买礼物(建模)
P1194 买礼物 题目描述 又到了一年一度的明明生日了,明明想要买B样东西,巧的是,这B样东西价格都是A元. 但是,商店老板说最近有促销活动,也就是: 如果你买了第I样东西,再买第J样,那么就可以只 ...