codevs1022 覆盖[Hungary 二分图最大匹配]

codevs1022 覆盖

有一个N×M的单位方格中，其中有些方格是水塘，其他方格是陆地。如果要用1×2的矩阵区覆盖（覆盖过程不容许有任何部分重叠）这个陆地，那么最多可以覆盖多少陆地面积。

输入描述 Input Description

输入文件的第一行是两个整数N，M  (1<=N，M<=100)，第二行为一个整数K( K<=50)，接下来的K行，每行两个整数X，Y表示K个水塘的行列位置。(1<=X<=N，1<=Y<=M)。

输出描述 Output Description

输出所覆盖的最大面积块（1×2面积算一块）。

样例输入 Sample Input

4 4

6

1 1

1 4

2 2

4 1

4 2

4 4

样例输出 Sample Output

4

---

相邻的可以组成1*2建边，根据i+j奇偶分类就成了二分图，hungary求最大匹配（dinic当然也可以，本题规模没必要）

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
const int N=;
inline int read(){
    char c=getchar();int x=,f=;
    while(c<''||c>''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
    return x*f;
}
int n,m,k,x,y,g[N][N];
struct edge{
    int v,ne;
}e[N*N*];
int h[N*N],cnt=;
inline int id(int i,int j){return (i-)*n+j;}
inline void ins(int u,int v){
    cnt++;
    e[cnt].v=v;e[cnt].ne=h[u];h[u]=cnt;
}
void buildGraph(){
    for(int i=;i<=n;i++)
        for(int j=;j<=m;j++) if(!g[i][j]&&(i+j)&){
            if(i<n&&!g[i+][j]) ins(id(i,j),id(i+,j));
            if(j<m&&!g[i][j+]) ins(id(i,j),id(i,j+));
            if(i>&&!g[i-][j]) ins(id(i,j),id(i-,j));
            if(j>&&!g[i][j-]) ins(id(i,j),id(i,j-));
        }
}
int le[N*N],vis[N*N];
bool find(int u){
    for(int i=h[u];i;i=e[i].ne){
        int v=e[i].v;
        if(vis[v]) continue;
        vis[v]=;
        if(!le[v]||find(le[v])){
            le[v]=u;
            return true;
        }
    }
    return false;
}
int hungary(){
    int ans=,nn=id(n,m);
    for(int i=;i<=nn;i++){
        memset(vis,,sizeof(vis));
        if(find(i)) ans++;
    }
    return ans;
}
int main(int argc, const char * argv[]) {
    n=read();m=read();k=read();
    while(k--){x=read();y=read();g[x][y]=;}
    buildGraph();
    printf("%d",hungary());
    return ;
}