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;
}