【BZOJ-1458】士兵占领 最大流

1458: 士兵占领

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 782  Solved: 456
[Submit][Status][Discuss]

Description

有一个M * N的棋盘，有的格子是障碍。现在你要选择一些格子来放置一些士兵，一个格子里最多可以放置一个士兵，障碍格里不能放置士兵。我们称这些士兵占领了整个棋盘当满足第i行至少放置了Li个士兵, 第j列至少放置了Cj个士兵。现在你的任务是要求使用最少个数的士兵来占领整个棋盘。

Input

第一行两个数M, N, K分别表示棋盘的行数，列数以及障碍的个数。 第二行有M个数表示Li。 第三行有N个数表示Ci。 接下来有K行，每行两个数X, Y表示(X, Y)这个格子是障碍。

Output

输出一个数表示最少需要使用的士兵个数。如果无论放置多少个士兵都没有办法占领整个棋盘，输出”JIONG!” (不含引号)

Sample Input

4 4 4
1 1 1 1
0 1 0 3
1 4
2 2
3 3
4 3

Sample Output

4
数据范围
M, N <= 100, 0 <= K <= M * N

HINT

Source

Solution

仰慕黄学长....

把题目转化一下，使用最少的，转换成删去最多的，那么可用最大流求解

对于某行或某列，如果可以放的个数小于必须放的个数，那么直接JIONG!

那么对于源向各行连边，容量为可以放的格子数 – 需求的格子数

对于各列向汇连边，容量同上

从可放置的点的行连至列，容量为1

Code

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
int read()
{
    int x=,f=; char ch=getchar();
    while (ch<'' || ch>'') {if (ch=='-') f=-; ch=getchar();}
    while (ch>='' && ch<='') {x=x*+ch-''; ch=getchar();}
    return x*f;
}
#define maxm 20000
#define maxn 500
int n,m,k;int L[],C[];
int zt[][];
int hang[],lie[],ans,tot,S,T;
struct Edgenode{int next,to,cap;}edge[maxm];
int head[maxn],cnt=;
void add(int u,int v,int w)
{cnt++;edge[cnt].to=v;edge[cnt].next=head[u];head[u]=cnt;edge[cnt].cap=w;}
void insert(int u,int v,int w)
{add(u,v,w);add(v,u,);}
//
int q[maxn],dis[maxn],cur[maxn];
bool bfs()
{
    for (int i=S; i<=T; i++) dis[i]=-;
    int he=,ta=; q[]=S; dis[S]=;
    while (he<ta)
        {
            int now=q[he++];
            for (int i=head[now]; i; i=edge[i].next)
                if (edge[i].cap && dis[edge[i].to]==-)
                    {
                        dis[edge[i].to]=dis[now]+;
                        q[ta++]=edge[i].to;
                    }
        }
    return dis[T]!=-;
}
int dfs(int loc,int low)
{
    if (loc==T) return low;
    int w,used=;
    for (int i=cur[loc]; i; i=edge[i].next)
        if (edge[i].cap && dis[edge[i].to]==dis[loc]+)
            {
                w=dfs(edge[i].to,min(low-used,edge[i].cap));
                edge[i].cap-=w; edge[i^].cap+=w;
                used+=w; if (edge[i].cap) cur[loc]=i;
                if (used==low) return low;
            }
    if (!used) dis[loc]=-;
    return used;
}
#define inf 0x7fffffff
int dinic()
{
    int tmp=;
    while (bfs())
        {
            for (int i=S; i<=T; i++) cur[i]=head[i];
            tmp+=dfs(S,inf);
        }
    return tmp;
}
void make()
{
    S=,T=n+m+;
    for (int i=; i<=m; i++)
        insert(S,i,n-L[i]-hang[i]);
    for (int i=; i<=n; i++)
        insert(i+m,T,m-C[i]-lie[T]);
    for (int i=; i<=m; i++)
        for (int j=; j<=n; j++)
            if (!zt[i][j]) insert(i,j+m,);
}
int main()
{
    m=read(),n=read(),k=read();
    for (int i=; i<=m; i++) L[i]=read();
    for (int i=; i<=n; i++) C[i]=read();
    tot=n*m-k;
    for (int x,y,i=; i<=k; i++)
        x=read(),y=read(),hang[x]++,lie[y]++,zt[x][y]=;
    for (int i=; i<=m; i++)
        if (hang[i]>n-L[i]) {puts("JIONG!");return ;}
    for (int i=; i<=n; i++)
        if (lie[i]>m-C[i]) {puts("JIONG!");return ;}
    make();
    ans=dinic();
    printf("%d\n",tot-ans);
    return ;
}

好厉害...一开始没想到...看到这种网格图，就想黑白染色QAQ