1. 程式人生 > >[bzoj 1458] 士兵占領(網絡流dinic)

[bzoj 1458] 士兵占領(網絡流dinic)

def 一行 algo body aps 一次 需要 sta include

1458: 士兵占領

Time Limit: 10 Sec Memory Limit: 64 MB
Submit: 784 Solved: 458
[Submit][Status][Discuss]

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
可以考慮拆點建網絡流圖 由源點向所有可以放置士兵的位置連一條容量為2的邊, 再由每個位置向它所屬的行和列的點連邊容量為1, 然後由所有的行和列向匯點連邊,容量為該行(該列)所需的士兵個數。 然後套dinic版子。。。。話說那玩意真的難寫5555

考慮 為什麽源點向放置士兵的位置連邊為2?

我是這麽想的,因為一個位置不僅需要向他所在的列連邊還需要向他所在的行連邊,但是行和列是互不影響的,所以只能分別補充,因為補充的容量為1,所以每個點只能算一次。

因為這樣連邊每個點被算了兩次 所以算出來之後ans要除以2 剛入cnblogs,第一篇題解...... 代碼有點醜,將就看下23333 技術分享圖片
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
using
namespace std; #define inf 2147483647 void read(int &x) { char ch; bool ok; for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch==-) ok=1; for(x=0; isdigit(ch); x=x*10+ch-0,ch=getchar()); if(ok) x=-x; } struct edge{ int w,v,nxt; }e[1000000]; int head[30000],tot,level[30000],n,m,k; void ad(int u,int v,int w) { e[tot].v=v; e[tot].w=w; e[tot].nxt=head[u]; head[u]=tot++; } void add(int u,int v,int w){ad(u,v,w),ad(v,u,0);} int bfs(int s,int t) { memset(level,-1,sizeof(level)); queue<int> q; q.push(s);level[s]=0; while(!q.empty()) { int u=q.front();q.pop(); for(int i=head[u];~i;i=e[i].nxt) if(e[i].w>0&&level[e[i].v]<0) level[e[i].v]=level[u]+1,q.push(e[i].v); } return level[t]>0; } int dfs(int u,int t,int f) { if(u==t) return f; for(int i=head[u];~i;i=e[i].nxt) if(level[e[i].v]>level[u]&&e[i].w>0) { int d=dfs(e[i].v,t,min(f,e[i].w)); if(d>0){e[i].w-=d,e[i^1].w+=d;return d;} } level[u]=-1; return 0; } int solve(int s,int t) { int flow=0,f; while(bfs(s,t)){while(f=dfs(s,t,inf)) flow+=f;} return flow; } bool mp[104][104]; int main() { memset(mp,true,sizeof(mp)); read(m),read(n),read(k); memset(head,-1,sizeof(head)); int x,y,z; int s=0,t=n*m+m+n+1;int ans=0; for(int i=1;i<=m;i++)read(x),add(i+n+n*m,t,x); for(int i=1;i<=n;i++)read(x),add(i+n*m,t,x); for(int i=1;i<=k;i++) read(x),read(y),mp[x][y]=false; for(int i=1;i<=m;i++) for(int j=1;j<=n;j++) if(mp[i][j]) add(s,(i-1)*n+j,2),add((i-1)*n+j,i+n+n*m,1),add((i-1)*n+j,j+n*m,1); printf("%d\n",solve(s,t)/2); }
View Code

[bzoj 1458] 士兵占領(網絡流dinic)