[Luogu P2891/POJ 3281/USACO07OPEN ]吃饭Dining
传送门:https://www.luogu.org/problemnew/show/P2891
题面
\
Solution
网络流
先引用一句真理:网络流最重要的就是建模
今天这道题让我深有体会
首先,观察数据范围,n=100,一般这种100-1000的图论题,很有可能是网络流.
那就直接从网络流的角度入手
考虑这样建模
建模要点如下:
1.建权值为1的边,保证每个食物和水仅用一次
2.没了
对以上的图求一个最大流,那不就是我们想要的最大的匹配数吗?
看起来是不是很OjbK?
其实不然,这样子一头牛有可能脚踏N条食物和水,但是题目要求一头牛只能吃喝一次
反例如下:
所以说,我们要对一头牛吃的东西做一个限制,保证其只流过1
怎么限制呢?
直接把一头牛拆成两头牛,中间连一条边就OK了嘛
如下图:
接下来就可以考虑如何给点编号了
我的编号方法很直接,很暴力
源点:1
食物: 2 ~ 2+f-1
牛 2+f ~ 2+f+2n -1
水: 2+f+2n ~ 2+f+2n+d
汇点:1000
如下图所示
按照以上方法,每种物品对应的点为:
食物i : 1+i
牛i : 1+f+i
牛i的右边的分点: 1+f+n+i
水i: 1+f+2*n+i
然后dinic直接求一波最大流就可以带走啦
DINIC教程传送门: http://www.cnblogs.com/SYCstudio/p/7260613.html (个人感觉写得很清楚)
Code
//Luogu P2891 [USACO07OPEN]吃饭Dining
//Apr,2ed,2018
//Dinic求最大流
#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
using namespace std;
long long read()
{
long long x=0,f=1; char c=getchar();
while(!isdigit(c)){if(c=='-') f=-1;c=getchar();}
while(isdigit(c)){x=x*10+c-'0';c=getchar();}
return x*f;
}
const int N=100+10;
const int M=N*10;
const int inf=0x3f3f3f3f;
struct line
{
int t,rev,w;
};
vector <line> e[M];
inline void AddLine(int s,int t)
{
line temp;
temp.t=t,temp.rev=e[t].size(),temp.w=1;
e[s].push_back(temp);
temp.t=s,temp.rev=e[s].size()-1,temp.w=0;
e[t].push_back(temp);
}
int n,f,d;
int dl[M],head,tail,depth[M];
bool visited[M];
bool bfs()
{
memset(visited,0,sizeof visited);
dl[1]=1,depth[1]=1;
visited[1]=true;
head=1,tail=2;
while(head<tail)
{
int now=dl[head],size=e[now].size();
for(int i=0;i<size;i++)
if(visited[e[now][i].t]==false and e[now][i].w>0)
{
visited[e[now][i].t]=true;
dl[tail++]=e[now][i].t;
depth[e[now][i].t]=depth[now]+1;
}
head++;
}
return visited[1000];
}
int dfs(int now,int f)
{
if(now==1000)
return f;
int size=e[now].size(),ans=0;
for(int i=0;i<size;i++)
if(depth[e[now][i].t]==depth[now]+1 and e[now][i].w>0)
{
int t_ans=dfs(e[now][i].t,min(f,e[now][i].w));
e[now][i].w-=t_ans;
e[e[now][i].t][e[now][i].rev].w+=t_ans;
f-=t_ans,ans+=t_ans;
if(f==0) break;
}
return ans;
}
int Dinic()
{
int ans=0;
while(bfs())
ans+=dfs(1,inf);
return ans;
}
int main()
{
n=read(),f=read(),d=read();
for(int i=1;i<=1000;i++)
e[i].reserve(16);
for(int i=1;i<=f;i++)
AddLine(1,1+i);
for(int i=1;i<=d;i++)
AddLine(1+f+2*n+i,1000);
for(int i=1;i<=n;i++)
{
int F=read(),D=read(),temp;
for(int j=1;j<=F;j++)
{
temp=read();
AddLine(1+temp,1+f+i);
}
for(int j=1;j<=D;j++)
{
temp=read();
AddLine(1+f+n+i,1+f+2*n+temp);
}
}
for(int i=1;i<=n;i++)
AddLine(1+f+i,1+f+n+i); printf("%d",Dinic());
return 0;
}
C++
后记
写的时候发现自己真的有点生疏了,这种比较靠记忆的算法还是久不久搞一下来恢复记忆为妙
网络流的时间复杂度真玄学
[Luogu P2891/POJ 3281/USACO07OPEN ]吃饭Dining的更多相关文章
- 「洛谷P2891」[USACO07OPEN]吃饭Dining 解题报告
P2891 [USACO07OPEN]吃饭Dining 题目描述 Cows are such finicky eaters. Each cow has a preference for certain ...
- P2891 [USACO07OPEN]吃饭Dining
漂亮小姐姐点击就送:https://www.luogu.org/problemnew/show/P2891 题目描述 Cows are such finicky eaters. Each cow ha ...
- P2891 [USACO07OPEN]吃饭Dining(最大流+拆点)
题目描述 Cows are such finicky eaters. Each cow has a preference for certain foods and drinks, and she w ...
- P2891 [USACO07OPEN]吃饭Dining 最大流
\(\color{#0066ff}{ 题目描述 }\) 有F种食物和D种饮料,每种食物或饮料只能供一头牛享用,且每头牛只享用一种食物和一种饮料.现在有n头牛,每头牛都有自己喜欢的食物种类列表和饮料种类 ...
- 洛谷P2891 [USACO07OPEN]吃饭Dining
题目描述 Cows are such finicky eaters. Each cow has a preference for certain foods and drinks, and she w ...
- 洛谷 P2891 [USACO07OPEN]吃饭Dining
裸的最大流. #include <cstdio> #include <cstring> #include <queue> const int MAXN = 4e3 ...
- [USACO07OPEN]吃饭Dining
嘟嘟嘟 这应该是网络流入门题之一了,跟教辅的组成这道题很像. 把每一只牛看成书,然后对牛拆点,因为每一只牛只要一份,食物和饮料分别看成练习册和答案. #include<cstdio> #i ...
- bzoj1711[USACO07OPEN]吃饭Dining
题意 有F种食物和D种饮料,每种食物或饮料只能供一头牛享用,且每头牛只享用一种食物和一种饮料.现在有n头牛,每头牛都有自己喜欢的食物种类列表和饮料种类列表,问最多能使几头牛同时享用到自己喜欢的食物和饮 ...
- poj 3281 Dining 网络流-最大流-建图的题
题意很简单:JOHN是一个农场主养了一些奶牛,神奇的是这些个奶牛有不同的品味,只喜欢吃某些食物,喝某些饮料,傻傻的John做了很多食物和饮料,但她不知道可以最多喂饱多少牛,(喂饱当然是有吃有喝才会饱) ...
随机推荐
- MySQL 8.0索引合并
简介 参考https://dev.mysql.com/doc/refman/8.0/en/index-merge-optimization.html#index-merge-intersection. ...
- Centos-系统内存信息-free
free 显示系统内存信息,包括物理内存.虚拟内存.共享内存和系统缓存 相关选项 -b 以字节byte为单位显示内存使用情况 -k 以k为单位显示内存使用情况 -m 以MB为单位显示内存使用情况 - ...
- python 系统设置
1. 设置python运行环境为utf-8 import sys #引用sys模块 reload(sys) #重新加载sys sys.setdefaultencoding("utf-8&qu ...
- 当安装、卸载件包时,出现依赖问题 error: Failed dependencies解决办法
error: Failed dependencies:-- 依赖关系非常复杂,当你试图先安装任何一个包时都会出现这样的依赖关系错误,这时候你就应该强制安装了,我认为只要你把服务或软件需要的包都装上,强 ...
- vs调试程序缺少 msvcp140d.dll 解决方法
简介一下吧: 如果只是为了解决问题请直接看第 7 点 ,谢谢. vs2013运行刚安装的opencv问题总结,尤其是电脑还很渣的情况下------花了我起码2天样子----很无奈 ...
- Serial.begin
串口波特率的设置:通常我们使用Serial.begin(speed)来完成串口的初始化,这种方式,只能配置串口的波特率. 使用Serial.begin(speed, config)可以配置数据位.校验 ...
- Keil ARm新建项目
一.新建一个工程 选好芯片后确认,完成创建 二.新建一个文件 保存为后缀名为*.c的文件 三.把文件添加进项目里面 四.测试 发现有警告 五.给项目添加特定的文件,去除警告或错误 现在保存项目的文件夹 ...
- 谈谈InnoDB中的B+树索引
索引类似于书的目录,他是帮助我们从大量数据中快速定位某一条或者某个范围数据的一种数据结构.有序数组,搜索树都可以被用作索引.MySQL中有三大索引,分别是B+树索引.Hash索引.全文索引.B+树索引 ...
- 关于IPA文件重签名后如何跟踪管理APP的技术探讨和实践演示
前言:开发iOS的朋友都知道,在功能开发完成后,我们就会用申请的苹果账号在后台做证书配置,然后提交到AppStore,但是也有部分APP我们不需要提交到AppStore,比如内部测试用的APP.定制给 ...
- Fiddler抓包工具 请求图标为一个锁的图标的设置
第一步,Fiddler抓包的数据 前面的都是一个锁的图标,的设置方法, 然后 点击打开 按此设置图一 在图二, 图三. 其他默认就好