POJ 3281 (最大流+匹配+拆点)
题目链接:http://poj.org/problem?id=3281
题目大意:有一些牛,一堆食物,一堆饮料。一头牛要吃一份食物喝一份饮料才算满足,而且牛对某些食物和饮料才有好感,问最多有多少头牛是满足的。
解题思路:
没有费用的匹配最大流题。
我一开始是这么考虑的,S->牛->食物->饮料->T,cap都是1,啊哈,多简单。
这样是不对的。因为牛比较挑剔,假设牛1喜欢食物1,而不喜欢饮料1,那么食物1和饮料1之间连不连边呢?
不连吧,其它牛怎么办?
连吧,牛1怎么办?
果断不能这么建图。于是改成S->食物->牛->饮料->T,这样就解决了模棱两可的食物与饮料之间的联系。
但是新的问题又来了,假设有食物1->牛1->饮料1,那么又会有食物2->牛1->饮料2,这样一头牛就吞了多份食物和饮料。
解决方案是拆点,把牛插成牛->牛',cap=1,这样一头牛只会被用一次了。
最终的建图方案:S->食物->牛->牛’->饮料->T。
#include "cstdio"
#include "vector"
#include "cstring"
#include "queue"
using namespace std;
#define maxn 405
#define inf 100000000
struct Edge
{
int from,to,cap,flow;
Edge(int FROM,int TO,int CAP,int FLOW):from(FROM),to(TO),cap(CAP),flow(FLOW) {}
};
int d[maxn],p[maxn],gap[maxn],cur[maxn];
bool vis[maxn];
vector<int> G[maxn],food[],drink[];
vector<Edge> edges;
void addedge(int from,int to,int cap)
{
edges.push_back(Edge(from,to,cap,));
edges.push_back(Edge(to,from,,));
int m=edges.size();
G[from].push_back(m-);
G[to].push_back(m-);
}
void bfs(int s,int t)
{
memset(vis,false,sizeof(vis));
memset(d,,sizeof(d));
memset(p,,sizeof(p));
d[t]=;vis[t]=true;
queue<int> Q;Q.push(t);
while(!Q.empty())
{
int u=Q.front();Q.pop();
for(int v=;v<G[u].size();v++)
{
Edge e=edges[G[u][v]^];
if(!vis[e.from]&&e.cap>e.flow)
{
vis[e.from]=true;
d[e.from]=d[u]+;
Q.push(e.from);
}
}
}
}
int augment(int s,int t)
{
int x=t,a=inf;
while(x!=s)
{
Edge e=edges[p[x]];
a=min(a,e.cap-e.flow);
x=e.from;
}
x=t;
while(x!=s)
{
edges[p[x]].flow+=a;
edges[p[x]^].flow-=a;
x=edges[p[x]].from;
}
return a;
}
int maxflow(int s,int t)
{
int flow=,u=s;
bfs(s,t);
memset(gap,,sizeof(gap));
memset(cur,,sizeof(cur));
for(int i=;i<=t;i++) gap[d[i]]++;
while(d[s]<t+)
{
if(u==t)
{
flow+=augment(s,t);
u=s;
}
bool flag=false;
for(int v=cur[u];v<G[u].size();v++) //Advance
{
Edge e=edges[G[u][v]];
if(e.cap>e.flow&&d[u]==d[e.to]+)
{
flag=true;
p[e.to]=G[u][v];
cur[u]=v;
u=e.to;
break;
}
}
if(!flag) //Retreat
{
int m=t+;
for(int v=;v<G[u].size();v++)
{
Edge e=edges[G[u][v]];
if(e.cap>e.flow) m=min(m,d[e.to]);
}
if(--gap[d[u]]==) break;
gap[d[u]=m+]++;
cur[u]=;
if(u!=s) u=edges[p[u]].from;
}
}
return flow;
}
int main()
{
int N,M,K,f,d,t;
while(scanf("%d%d%d",&N,&M,&K)!=EOF)
{
for(int i=;i<=M;i++) addedge(,i,); //S-food
for(int i=;i<=K;i++) addedge(i+M+N*,K+M+N*+,); //drink-T
for(int i=;i<=N;i++)
{
addedge(i+M,i+M+N,); //cow-cow'
scanf("%d%d",&f,&d);
for(int j=; j<=f; j++)
{
scanf("%d",&t);
addedge(t,i+M,); //food-cow
}
for(int j=; j<=d; j++)
{
scanf("%d",&t);
addedge(i+M+N,t+M+*N,); //cow'-drink
}
}
printf("%d\n",maxflow(,*N+M+K+));
}
}
13456393 | neopenx | 3281 | Accepted | 320K | 16MS | C++ | 3222B | 2014-09-19 00:28:49 |
POJ 3281 (最大流+匹配+拆点)的更多相关文章
- poj 3281 最大流+建图
很巧妙的思想 转自:http://www.cnblogs.com/kuangbin/archive/2012/08/21/2649850.html 本题能够想到用最大流做,那真的是太绝了.建模的方法很 ...
- POJ 3281 Dining(网络流拆点)
[题目链接] http://poj.org/problem?id=3281 [题目大意] 给出一些食物,一些饮料,每头牛只喜欢一些种类的食物和饮料, 但是每头牛最多只能得到一种饮料和食物,问可以最多满 ...
- poj 3281 最大流建图
题目链接:http://poj.org/problem?id=3281 #include <cstdio> #include <cmath> #include <algo ...
- POJ 3281 最大流
Dining Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 17251 Accepted: 7643 Descripti ...
- [poj 3281]最大流+建图很巧妙
题目链接:http://poj.org/problem?id=3281 看了kuangbin大佬的思路,还用着kuangbin板子orz http://www.cnblogs.com/kuangb ...
- poj 3281 Dining(网络流+拆点)
Dining Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 20052 Accepted: 8915 Descripti ...
- POJ 3281 网络流 拆点保证本身只匹配一对食物和饮料
如何建图? 最开始的问题就是,怎么表示一只牛有了食物和饮料呢? 后来发现可以先将食物与牛匹配,牛再去和饮料匹配,实际上这就构成了三个层次. 起点到食物层边的容量是1,食物层到奶牛层容量是1,奶牛层到饮 ...
- POJ 3281 Dining(最大流)
POJ 3281 Dining id=3281" target="_blank" style="">题目链接 题意:n个牛.每一个牛有一些喜欢的 ...
- poj 3281 Dining 网络流-最大流-建图的题
题意很简单:JOHN是一个农场主养了一些奶牛,神奇的是这些个奶牛有不同的品味,只喜欢吃某些食物,喝某些饮料,傻傻的John做了很多食物和饮料,但她不知道可以最多喂饱多少牛,(喂饱当然是有吃有喝才会饱) ...
随机推荐
- HDU 2897
Problem Description 当日遇到月,于是有了明.当我遇到了你,便成了侣. 那天,日月相会,我见到了你.而且,大地失去了光辉,你我是否成侣?这注定是个凄美的故事.(以上是废话) 小t和所 ...
- 使用Discuz关键词服务器实现PHP中文分词
不同于使用自己的服务器进行分词,Discuz!在线中文分词服务是基于API返回分词结果的.在项目中,我们只需要一个函数即可方便地进行分词.关键词提取.以下是根据Discuz!在线分词服务API写的函数 ...
- shell脚本检测局域网内存活主机
<1> d211 admin # for i in {3..254} ; do ping -c 1 192.168.1.$i &>/dev/null && e ...
- EtherCAT报文寻址
EtherCAT通信通过主站发送EtherCAT数据帧读写从站设备的内部存储区实现. 一个EtherCAT网段相当于一个以太网设备,主站首先通过以太网数据帧头的MAC地址寻址到网段,然后使用Ether ...
- zpf 命名规则
2014年8月19日 18:48:39 所有控制器都要继承main类,main类是一个入口类,他里边根据请求初始化了一些变量,也初始化了一些系统变量等等,这些变量和函数可以被控制器类直接使用 控制器类 ...
- MongoDB副本集学习(一):概述和环境搭建
MongoDB副本集概述 以下图片摘自MongoDB官方文档:http://docs.mongodb.org/manual/core/replication-introduction/ Primary ...
- Linux下的ip命令,除了ifconfig还有很多
linux的ip命令和ifconfig类似,但前者功能更强大,并旨在取代后者.使用ip命令,只需一个命令,你就能很轻松地执行一些网络管理任务.ifconfig是net-tools中已被废弃使用的一个命 ...
- app后端架构设计(转)
(1)Restful设计原则 Restful风格:RESTfu设计原则,它被Roy Felding提出(在他的”基于网络的软件架构“论文中第五章).而REST的核心原则是将你的API拆分为逻辑上的资源 ...
- jquery easy ui 1.3.4 快速入门(1)
什么是easyui jQuery EasyUI是一组基于jQuery的UI插件集合,而jQuery EasyUI的目标就是帮助web开发者更轻松的打造出功能丰富并且美观的UI界面.开发者不需要编写复杂 ...
- javascript二叉树基本功能实现
都是常用的功能. 删除是最复杂的.. <!DOCTYPE html> <html lang="en"> <head> <meta char ...