HDOJ 3605 Escape 解题报告 二分图匹配 dinic 网络流
2017-05-30 11:01
399 查看
【Problem Description】
2012 If this is the end of the world how to do? I do not know how. But now scientists have found that some stars, who can live, but some people do not fit to live some of the planet. Now scientists want your help, is to determine what all of people can live in these planets.
题目大意:
这道题就是给你n个人,m个星球,这n个人对m个星期有不同的意向。然而,每个星球都有自己的承受能力。我们想知道能否全部满足这n个人的意向。
【解题报告】
我们不难想出这个网络流怎么建:
源点S连接每个人,边权为1,每个人同星球相连,边权为1,每个星球通汇点T相连,边权即承受能力。最后再跑一边dinic。
然而人们发现,这样一来铁定会超时(点太多了)。所以,我们采取了一个缩点的策略,将去同一个星球的人缩成了一个点,这样一来,每个人同星球的边权自然就不是1,而是人数。
下面我们看看代码:
以上
2012 If this is the end of the world how to do? I do not know how. But now scientists have found that some stars, who can live, but some people do not fit to live some of the planet. Now scientists want your help, is to determine what all of people can live in these planets.
题目大意:
这道题就是给你n个人,m个星球,这n个人对m个星期有不同的意向。然而,每个星球都有自己的承受能力。我们想知道能否全部满足这n个人的意向。
【解题报告】
我们不难想出这个网络流怎么建:
源点S连接每个人,边权为1,每个人同星球相连,边权为1,每个星球通汇点T相连,边权即承受能力。最后再跑一边dinic。
然而人们发现,这样一来铁定会超时(点太多了)。所以,我们采取了一个缩点的策略,将去同一个星球的人缩成了一个点,这样一来,每个人同星球的边权自然就不是1,而是人数。
下面我们看看代码:
#include<cstdio> #include<cstring> #include<string> #include<algorithm> #include<vector> #include<queue> #include<cmath> using namespace std; #define Del(a,b) memset(a,b,sizeof(a)) const int N=1210; const int inf=0x3f3f3f3f; const double esp=1e-9; int n,m; struct Node { int from,to,cap,flow; }; vector<int> v ; int sum ; vector<Node> e; int vis ; //构建层次图 int cur ; int MIN(int x,int y) { return x>y?y:x; } void add_Node(int from,int to,int cap) { Node tmp1,tmp2; tmp1.from=from,tmp1.to=to,tmp1.cap=cap,tmp1.flow=0; e.push_back(tmp1); tmp2.from=to,tmp2.to=from,tmp2.cap=0,tmp2.flow=0; e.push_back(tmp2); int tmp=e.size(); v[from].push_back(tmp-2); v[to].push_back(tmp-1); } bool bfs(int s,int t) { Del(vis,-1); queue<int> q; q.push(s); vis[s] = 0; while(!q.empty()) { int x=q.front(); q.pop(); int i; for(i=0;i<v[x].size();i++) { Node tmp = e[v[x][i]]; if(vis[tmp.to]<0 && tmp.cap>tmp.flow) //第二个条件保证 { vis[tmp.to]=vis[x]+1; q.push(tmp.to); } } } if(vis[t]>0) return true; return false; } int dfs(int o,int f,int t) { if(o==t || f==0) //优化 return f; int a = 0,ans=0; int &i = cur[o]; for(;i<v[o].size();i++) //注意前面 ’&‘,很重要的优化 { Node &tmp = e[v[o][i]]; if(vis[tmp.to]==(vis[o]+1) && (a = dfs(tmp.to,MIN(f,tmp.cap-tmp.flow),t))>0) { tmp.flow+=a; e[v[o][i]^1].flow-=a; //存图方式 ans+=a; f-=a; if(f==0) //注意优化 break; } } return ans; //优化 } int dinci(int s,int t) { int ans=0; while(bfs(s,t)) { Del(cur,0); int tm=dfs(s,inf,t); ans+=tm; } return ans; } int main() { int n,m; while(~scanf("%d%d",&n,&m)) { int i,j; memset(sum,0,sizeof(sum)); int s=0,x,t=m+(1<<m)+3,nn=(1<<m)+2; for(i=1;i<=n;i++) { int tmp=0; for(j=0;j<m;j++) { scanf("%d",&x); if(x) tmp+=(1<<j); } sum[tmp]++; } for(i=0;i<(1<<m);i++) { if(sum[i]) { add_Node(s,i+1,sum[i]); for(j=0;j<m;j++) { if(i&(1<<j)) add_Node(i+1,nn+j,sum[i]); } } } for(i=0;i<m;i++) { scanf("%d",&x); add_Node(nn+i,t,x); } int ans=dinci(s,t); //printf("%d\n",ans); if(ans==n) printf("YES\n"); else printf("NO\n"); for(i=0;i<=t;i++) v[i].clear(); e.clear(); } return 0; }
以上
相关文章推荐
- [BZOJ1391]解题报告|网络流的又一类建图&Dinic的若干优化
- BZOJ 3894 文理分科 解题报告 最小割 网络流 DINIC
- HDOJ-1164-Eddy's research I 解题报告
- hdoj 2000 ASCII码排序 用java写的解题报告
- ACM--HDOJ2031解题报告
- ACM--HDOJ1201解题报告
- 网络流改进SAP算法模版、HDU 1532 Drainage Ditches(解题报告)
- HDU 3572 网络流最大流 解题报告
- ZOJ 3348 网络流最大流 解题报告
- 解题报告-HDOJ-1233(最小生成树——kruskal)
- HDOJ2072单词数 解题报告
- BZOJ 1822 计算几何+网络流+二分答案 解题报告
- [BZOJ3275]Number解题报告|网络流
- [置顶] HDOJ 1728 逃离迷宫 (BFS )解题报告
- hdoj2007解题报告
- HDOJ-1905-Pseudoprime numbers 解题报告
- hdoj2012解题报告
- ACM-HDoj暑假竞赛(7)-1009解题报告(
- LA-4287 & HDOJ-2767 Proving Equivalences 解题报告
- 线性规划与网络流24题解题报告