[UVALive-8198] Keep it covered

题目链接：
https://vjudge.net/contest/248768#problem/K

题目大意：

给你一个n*m的格子图， 初始时一些格子上有点’o’， 其余格子是空的。 要求用4种拼图(可以旋转), 来铺满原图， 点’o’只能用拼图1， 要保证线段的连通性。 问是否存在解。
$(n, m \leq 20)$

题目思路：
乍一看可能想到插头dp或者什么别的， 其实本质是一个二分图。
考虑一条线段就是给相邻格子提供度数， 我们维护度数平衡这个性质。
考虑点’o’格子， 只能用1号拼图， 即他能向相邻格子提供1的度数。
考虑空格子， 只能用2，3，4号拼图， 这三种拼图算上旋转各种组合， 即该格子能向所有相邻格子一共提供2的度数。
根据棋盘图的性质， 黑白染色后建立二分图。
点’o’格子初始流量为1， 空格子初始流量为2， 相邻格子之间的流量设为2， 跑最大流， check一下源点和汇点的边都满流即有解。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <bitset>
#include <map>

#define pi pair<int, int>
#define fi first
#define se second
#define mp make_pair
#define ll long long
#define ls (x << 1)
#define rs ((x << 1) | 1)
#define mid ((l + r) >> 1)
#define pb push_back

const int N = (int)2020 + 10;

using namespace std;

int gi(){
    char c = getchar(); int ret = 0;
    while (!isdigit(c)) c = getchar();
    while (isdigit(c)){
        ret = ret * 10 + c - '0';
        c = getchar();
    }
    return ret;
}

int n, m, k, s, t; char G[N][N];
int cnt, lst[N * N], nxt[N * N], to[N * N], f[N * N];
void add(int u, int v, int c){
    nxt[++ cnt] = lst[u]; lst[u] = cnt; to[cnt] = v; f[cnt] = c;
    nxt[++ cnt] = lst[v]; lst[v] = cnt; to[cnt] = u; f[cnt] = 0;
}

int mv[] = {1, 0, 1};
bool over(int x, int y){
    return x < 1 || y < 1 || x > n || y > m;
}
int id(int i, int j) {return (i - 1) * m + j;} 

const int inf = 1 <<30;
int head, tail, que[N * N], d[N * N];
bool bfs(){
    head = 0;
    que[tail = 1] = s;
    for (int i = 1; i <= t; i ++) d[i] = 0;
    d[s] = 1;

    while (head < tail){
        int u = que[++ head];
        for (int j = lst[u]; j; j = nxt[j]){
            int v = to[j];
            if (!f[j] || d[v]) continue;
            d[v] = d[u] + 1;
            que[++ tail] = v;
        }
    }

    return d[t];
}
int dfs(int u, int flow){
    if (u == t) return flow;
    int a, ret = 0;
    for (int j = lst[u]; j; j = nxt[j]){
        int v = to[j];
        if (f[j] && flow && d[v] == d[u] + 1 && (a = dfs(v, min(flow, f[j])))){
            flow -= a, ret += a;
            f[j] -= a, f[j ^ 1] += a;
        }
    }
    return ret;
}

int main()
{
    while (scanf("%d %d", &n, &m) != EOF){
        cnt = 1;
        s = n * m + 1, t = n * m + 2;
        for (int i = 1; i <= t; i ++) lst[i] = 0;

        for (int i = 1; i <= n; i ++) scanf("%s", G[i] + 1);
        for (int i = 1; i <= n; i ++)
            for (int j = 1; j <= m; j ++)
                for (int k = 0; k < 2; k ++){
                    int x = i + mv[k], y = j + mv[k + 1];
                    if (over(x, y)) continue;
                    if ((i + j) & 1) add(id(i, j), id(x, y), 2);
                    else add(id(x, y), id(i, j), 2);
                }

        int A = 0, B = 0;
        for (int i = 1; i <= n; i ++)
            for (int j = 1; j <= m; j ++)
                if ((i + j) & 1) add(s, id(i, j), 2 - (G[i][j] == 'o')), A += 2 - (G[i][j] == 'o');
                else add(id(i, j), t, 2 - (G[i][j] == 'o')), B += 2 - (G[i][j] == 'o');

        if (A != B) {puts("N"); continue;}

        int ans = 0, res;
        while (bfs()){
            while (res = dfs(s, inf)) ans += res; 
        }

        puts(ans == A ? "Y" : "N");
    }

    return 0;
}