BZOJ1458: 士兵占领

题解：

裸上下界网络流。。。就当复习了。。。这个资料不错：http://hi.baidu.com/newfarking/item/b9780317201c9651f0090e64

代码：

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<string>
#define inf 1000000000
#define maxn 500+5
#define maxm 100000
#define eps 1e-10
#define ll long long
#define pa pair<int,int>
#define for0(i,n) for(int i=0;i<=(n);i++)
#define for1(i,n) for(int i=1;i<=(n);i++)
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
#define for4(i,x) for(int i=head[x],y;i;i=e[i].next)
#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=10*x+ch-'0';ch=getchar();}
return x*f;
}
int  n,m,s,t,ss,tt,tot=1,cnt,head[maxn],cur[maxn],h[maxn],a[maxn];
queue<int>q;
bool b[maxn][maxn];
struct edge{int go,next,v;}e[maxm];
inline void ins(int x,int y,int w)
{
e[++tot]=(edge){y,head[x],w};head[x]=tot;
e[++tot]=(edge){x,head[y],0};head[y]=tot;
}
inline void insert(int x,int y,int l,int r)
{
cnt+=l;
ins(ss,y,l);
ins(x,tt,l);
ins(x,y,r-l);
}
bool bfs()
{
for(int i=0;i<=n+m+3;i++)h[i]=-1;
q.push(s);h[s]=0;
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=head[x];i;i=e[i].next)
if(e[i].v&&h[e[i].go]==-1)
{
h[e[i].go]=h[x]+1;q.push(e[i].go);
}
}
return h[t]!=-1;
}
int dfs(int x,int f)
{
if(x==t) return f;
int tmp,used=0;
for(int i=cur[x];i;i=e[i].next)
if(e[i].v&&h[e[i].go]==h[x]+1)
{
tmp=dfs(e[i].go,min(e[i].v,f-used));
e[i].v-=tmp;if(e[i].v)cur[x]=i;
e[i^1].v+=tmp;used+=tmp;
if(used==f)return f;
}
if(!used) h[x]=-1;
return used;
}
int dinic()
{
int maxflow=0;
while(bfs())
{
for (int i=0;i<=n+m+3;i++)cur[i]=head[i];maxflow+=dfs(s,inf);
}
return maxflow;
}

int minflow()
{
s=ss;t=tt;
if(dinic()!=cnt){printf("JIONG\n");return 0;}
int ans=e[tot].v;
e[tot].v=e[tot^1].v=0;
s=n+m+1;t=0;
return ans-dinic();
}
int main()
{
freopen("input.txt","r",stdin);
freopen("output.txt","w",stdout);
n=read();m=read();int k=read();
s=0;t=n+m+1;ss=t+1;tt=t+2;
for1(i,n+m)a[i]=read();
for1(i,k){int x=read(),y=read();b[x][y]=1;}
for1(i,n)for1(j,m)if(!b[i][j])insert(i,j+n,0,1);
for1(i,n)insert(s,i,a[i],inf);
for2(i,n+1,n+m)insert(i,t,a[i],inf);
insert(t,s,0,inf);
printf("%d\n",minflow());
return 0;
}

View Code

1458: 士兵占领

Time Limit: 10 Sec Memory Limit: 64 MB
Submit: 569 Solved: 325
[Submit][Status]

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