CodeForces - 903F Clear The Matrix

博主的 BiBi 时间

能想象博主勤勤恳恳修完锅，交上去直接 $\mathtt{M}\mathtt{L}\mathtt{E}$ 的快乐？（为什么不直接用滚动数组是真的憨）

然后勤勤恳恳改成从 $0$ 开始。

Solution

我们把矩阵翻过来，改成 $n$ 行 $4$ 列。对于左上角在 $i$ 行的操作，最多只能影响到 $i+3$ 行。

所以就令 $f[i][a][b][c][d]$ 为前 $i-1$ 行为 ‘.’，$i$ 到 $i+3$ 行状态分别为 $a,b,c,d$ 的最小花费。（且达到这个目标状态的操作的左上角的行数不超过 $i$）

然后枚举每种操作，与现在状态与起来就行了。（我是令 ‘.’ 为 $1$，‘*’ 为 $0$）需要注意的是将 $n$ 到 $n+3$ 行的原始状态赋值为 $2^4-1$，因为虽然有些行的操作已经越出矩阵，这个操作仍然合法，仍有可能改变有效行的状态，而多出的行需要强制合法。

详见代码。

Code

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

const int lim = (1 << 4);

int n, w[4], g[1005], f[1002][16][16][16][16], s[4][4][4], num[4][4];
char str[4][1000];

int read() {
	int x = 0, f = 1; char s;
	while((s = getchar()) > '9' || s < '0') if(s == '-') f = -1;
	while(s >= '0' && s <= '9') x = (x << 1) + (x << 3) + (s ^ 48), s = getchar();
	return x * f;
}

int main() {
	int t;
	n = read();
	for(int i = 0; i < 4; ++ i) w[i] = read();
	for(int i = 0; i < 4; ++ i) scanf("%s", str[i]);
	for(int i = 0; i < n; ++ i) {
		t = 0;
		for(int j = 0; j < 4; ++ j) t <<= 1, t += (str[j][i] == '.');
		g[i] = t;
	}
	for(int i = n; i < n + 5; ++ i) g[i] = lim - 1;
	memset(f, 0x3f, sizeof f);
	f[0][g[0]][g[1]][g[2]][g[3]] = 0;
	for(int i = 0; i < 4; ++ i) 
		for(int j = 0; j < 4 - i; ++ j) {
			for(int k = 0; k < 4; ++ k)
				for(int u = 0; u < 4; ++ u)
					num[k][u] = 0;
			for(int k = 0; k <= i; ++ k)
				for(int u = j; u <= j + i; ++ u)
					num[k][u] = 1;
			for(int k = 0; k < 4; ++ k) {
				t = 0;
				for(int u = 0; u < 4; ++ u) t <<= 1, t += num[k][u];
				s[i][j][k] = t;
			}
		}
	for(int i = 0; i < n; ++ i) {
		for(int a = 0; a < lim; ++ a)//DP 时下标趋势从小变大，所以枚举从小到大
		for(int b = 0; b < lim; ++ b)
		for(int c = 0; c < lim; ++ c)	 	
		for(int d = 0; d < lim; ++ d) {
			if(f[i][a][b][c][d] >= 0x3f3f3f3f) continue;//玄学优化
			if(a == lim - 1) f[i + 1][b][c][d][g[i + 4]] = min(f[i + 1][b][c][d][g[i + 4]], f[i][a][b][c][d]);//如果第 i 行已经填完
			for(int j = 0; j < 4; ++ j)
				for(int k = 0; k < 4 - j; ++ k) {
					f[i][a | s[j][k][0]][b | s[j][k][1]][c | s[j][k][2]][d | s[j][k][3]] = min(f[i][a | s[j][k][0]][b | s[j][k][1]][c | s[j][k][2]][d | s[j][k][3]], f[i][a][b][c][d] + w[j]);
				}
		}
	}
	printf("%d\n", f[n - 1][lim - 1][lim - 1][lim - 1][lim - 1]);
	return 0;
}