【bzoj1458】士兵占领 有源汇最小流

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

题解
在yzh大神的指导下，终于学会了最小流。
Orz YZH

S向每一行i连[li,inf]的边
每一行j向T连[cj,inf]的边
如果< i , j >不是障碍 i向j连[0,1]的边
跑最小流即可

给YZH神犇跪烂烂！

代码

#include<cstdio>
#include<iostream>
#include<queue>
#define ll long long
#define inf 1000000000
#define mod 1000000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
queue<int>q;
int S,T,SS,TT;
int n,cnt=1,ans,m,k;
int dis[205],in[205],cur[205];
int Head[205],Next[100005],ret[100005],len[100005];
int a[105],b[105],l[105],c[105];
bool f[105][105];
void ins(int u,int v,int w)
{
    ret[++cnt]=v;Next[cnt]=Head[u];Head[u]=cnt;len[cnt]=w;
}
void insert(int u,int v,int w)
{ins(u,v,w);ins(v,u,0);}
void build()
{
    for (int i=1;i<=n;i++)
    {
        insert(S,i,inf);
        in[S]-=l[i];in[i]+=l[i];
    }
    for (int i=1;i<=m;i++)
    {
        insert(i+n,T,inf);
        in[i+n]-=c[i];in[T]+=c[i];
    }
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            if (!f[i][j])insert(i,j+n,1);
    for (int i=1;i<=TT;i++)
        if (in[i]>0) insert(SS,i,in[i]);
        else insert(i,TT,-in[i]);
}
bool bfs(int x,int y)
{
    for (int i=1;i<=TT;i++)dis[i]=-1;
    dis[x]=0;q.push(x);
    while (!q.empty())
    {
        int now=q.front();q.pop();
        for (int i=Head[now];i;i=Next[i])if (dis[ret[i]]==-1&&len[i])
        {
            dis[ret[i]]=dis[now]+1;
            q.push(ret[i]);
        }
    }
    return dis[y]!=-1;
}
int dfs(int x,int y,int aug)
{
    if (x==y||aug==0) return aug;
    int flow=0;
    for (int i=cur[x];i;i=Next[i])if (dis[ret[i]]==dis[x]+1)
    {
        int t=dfs(ret[i],y,min(aug-flow,len[i]));flow+=t;
        len[i]-=t;len[i^1]+=t;if (len[i]) cur[x]=i;
        if (flow==aug) return flow; 
    } 
    if (flow==0) dis[x]=-1;
    return flow;
}
void dinic(int x,int y)
{
    while (bfs(x,y))
    {
        for (int i=1;i<=TT;i++) cur[i]=Head[i];
        dfs(x,y,inf<<1);
    }
}
int main()
{
    n=read();m=read();k=read();S=n+m+1;T=S+1;SS=T+1;TT=SS+1;
    for (int i=1;i<=n;i++) l[i]=read(),a[i]=m;
    for (int i=1;i<=m;i++) c[i]=read(),b[i]=n;
    for (int i=1;i<=k;i++)
    {
        int x=read(),y=read();
        f[x][y]=1;--a[x];--b[y];
    }
    for (int i=1;i<=n;i++) if (a[i]<l[i]) return puts("JIONG!"),0;
    for (int i=1;i<=m;i++) if (b[i]<c[i]) return puts("JIONG!"),0;
    build();
    dinic(SS,TT);
    insert(T,S,inf);
    dinic(SS,TT);
    printf("%d",len[cnt]);
    return 0;
}