BZOJ1458 士兵占领

题目大意：

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

扯点别的：

做网络流的题最重要的就是建模，模型建好了代码都不是什么问题。往往做一道题卡很久，突然灵机一动就想出来正解，觉得自己之前贼傻缺这么容易的题也要想半天。

解析：

看到这道题，他让求最少士兵，但是我并不知道怎么求最小流，于是开始向某个定值减去最大流的方向考虑。那么对于这道题，假设每个士兵只对某行或者某列产生了贡献，那么设ans=∑Li+∑Cj，但是事实上有的士兵对于所在行和列都产生了贡献，我们应该最大化这个值，所以最大值再用ans减去它就可以了。这样建图的方式就显而易见了。

如图，所求出的最大流即为中间的边的最大匹配，对应了最大化产生两个贡献的士兵的值。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
#define inf 1e9
struct point
{
    int next;
    int to;
    int dis;
}e[1000001];
int n,m,k,num,ans,s,t;
int h[10001],head[10001],c[10001],l[10001],ll[10001],cc[10001];
int g[2010][2010];
void add(int from,int to,int dis)
{
    e[num].next=head[from];
    e[num].dis=dis;
    e[num].to=to;
    head[from]=num++;
}
int dfs(int x,int dis)
{
    if(x==t)
        return dis;
    int sum=0;
    for(int i=head[x];i!=0;i=e[i].next)
    {
        int to=e[i].to;
        if(h[to]==h[x]+1&&e[i].dis)
        {
            int diss=dfs(to,min(dis,e[i].dis));
            if(diss>0)
            {
                sum+=diss;
                dis-=diss;
                e[i].dis-=diss;
                e[i^1].dis+=diss;
                if(!sum)
                    break;
            }
        }
    }
    if(!sum)
        h[x]=-1;
    return sum;
}
bool bfs()
{
    memset(h,0,sizeof(h));
    h[s]=1;
    queue<int>q;
    q.push(s);
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        for(int i=head[x];i!=0;i=e[i].next)
        {
            int to=e[i].to;
            if(!h[to]&&e[i].dis)
            {
                q.push(to);
                h[to]=h[x]+1;
            }
        }
    }
    return h[t];
}
int dinic()
{
    int tot=0;
    while(bfs())
        tot+=dfs(s,inf);
    return tot;
}
int main()
{
    cin>>m>>n>>k;
    t=m+n+1;
    for(int i=1;i<=m;i++)
    {
        cin>>l[i];
        ans+=l[i];
    }
    for(int i=1+m;i<=n+m;i++)
    {
        cin>>c[i];
        ans+=c[i];
    }
    for(int i=1;i<=k;i++)
    {
        int x,y;
        cin>>x>>y;
        g[x][y+m]=1;
        ll[x]++;
        cc[y]++;
        if(ll[x]+l[x]>m||cc[y]+c[y]>n)
        {
            cout<<"JIONG";
            return 0;
        }
    }
    for(int i=1;i<=m;i++)
    {
        add(i,0,0);
        add(0,i,l[i]);
        for(int j=1;j<=n;j++)
            if(!g[i][j+m])
            {
                add(j+m,i,0);
                add(i,j+m,1);
            }
    }
    for(int i=m+1;i<=n+m;i++)
    {
        add(t,i,0);
        add(i,t,c[i]);
    }
    cout<<ans-dinic();
    return 0;
}