题目大意:
有一个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;
}