网络流24题 P2754 [CTSC1999]家园
思路
如图,建立分层图跑dinic
每次在残余网络里加边继续跑
跑到ans>=k时候的i就是答案
诶呀啊,忘记弄箭头了,最后一列是向上的箭头,不过聪明的你们应该没啥影响
代码
#include <bits/stdc++.h>
#define FOR(i,a,b) for(int i=a;i<=b;++i)
using namespace std;
const int N = 1e6 + 7;
const int inf = 0x3f3f3f3f;
int read() {
int x=0,f=1;char s=getchar();
for(;s>'9'||s<'0';s=getchar()) if(s=='-') f = -1;
for(;s>='0'&&s<='9';s=getchar()) x = x*10 + s - '0';
return x*f;
}
int n, m, k, h[100], S, T;
std::vector<int> G[100];
struct node {
int v, nxt, cap;
}e[N];
int head[N], tot = 1;
void add_edge(int u, int v, int cap) {
e[++tot].v = v;
e[tot].cap = cap;
e[tot].nxt = head[u];
head[u] = tot;
}
int dis[N];
queue <int> Q;
bool bfs() {
memset(dis, -1, sizeof(dis));
while(!Q.empty()) Q.pop();
dis[S] = 0;
Q.push(S);
while(!Q.empty()) {
int u = Q.front();
Q.pop();
for(int i = head[u]; i;i = e[i].nxt) {
int v = e[i].v;
if(dis[v] == -1 && e[i].cap) {
dis[v] = dis[u] + 1;
Q.push(v);
}
}
}
return dis[T] != -1;
}
int dfs(int now, int f) {
if(now == T) return f;
int rest = f;
for(int i = head[now]; i;i = e[i].nxt) {
int v = e[i].v;
if(dis[v] == dis[now] + 1 && e[i].cap && rest) {
int t = dfs(v, min(e[i].cap, rest));
if(!t) dis[v] = 0;
e[i].cap -= t;
e[i^1].cap += t;
rest -= t;
}
}
return f-rest;
}
namespace bcj {
int siz[N], f[N];
int find(int x) {
return x == f[x] ? x : f[x] = find(f[x]);
}
void uu(int x, int y) {
int fx = find(x), fy = find(y);
if(siz[fx] <= siz[fy]) {
f[fx] = fy;
if(siz[fx] == siz[fy]) siz[fy] ++;
} else f[fy] = fx;
}
}
int main() {
n = read() ,m = read(), k = read();
S = 5e5+1, T = 5e5+2;
FOR(i,1,1000) bcj::f[i] = i;
FOR(i, 1, m) {
h[i] = read();
int x = read();
FOR(j, 1, x) {
int tmp = read();
if(tmp==0) tmp = n+1;
if(tmp==-1) tmp = n+2;
G[i].push_back(tmp);
}
}
FOR(i,1,m) {
for(std::vector<int>::iterator it = G[i].begin() + 1; it != G[i].end(); ++ it) {
bcj::uu(*it, *(it-1));
}
}
if(bcj::find(n+1) != bcj::find(n+2)) {puts("0"); return 0;}
int ans = 0,i = 1;
add_edge(S, n + 1, inf);
add_edge(n + 1, S, 0);
add_edge(n + 2, T, inf);
add_edge(T, n + 2, 0);
for(;; ++i) {
// my -> my
FOR(j, 1, n + 1) {
add_edge(j + (i-1) * (n + 2), j + i * (n + 2), inf);
add_edge(j + i * (n + 2), j + (i-1) * (n + 2), 0);
}
add_edge((i + 1) * (n + 2), i * (n + 2), inf);
add_edge(i * (n + 2), (i + 1) * (n + 2), 0);
// go go go
FOR(j, 1, m) {
int a = i % G[j].size() - 1;
int b = i % G[j].size();
a= a<0 ? G[j].size()-1 : a;
a =G[j][a] + (i - 1) * (n + 2);
b =G[j][b] + i * (n + 2);
add_edge(a, b, h[j]);
add_edge(b, a, 0);
}
while(bfs()) ans += dfs(S,inf);
if(ans >= k) break;
}
printf("%d\n", i);
return 0;
}
网络流24题 P2754 [CTSC1999]家园的更多相关文章
- [CTSC1999][网络流24题] 星际转移
36. [CTSC1999][网络流24题] 星际转移 ★★★☆ 输入文件:home.in 输出文件:home.out 简单对比时间限制:1 s 内存限制:128 MB «问题描述: ...
- COGS743. [网络流24题] 最长k可重区间集
743. [网络流24题] 最长k可重区间集 ★★★ 输入文件:interv.in 输出文件:interv.out 简单对比时间限制:1 s 内存限制:128 MB «问题描述: «编 ...
- Cogs 14. [网络流24题] 搭配飞行员
这道题其实蛮好想的,因为分为正,副飞行员.所以就把正飞行员当作Boy,副飞行员当作Girl.然后做Hungry即可. #include<bits/stdc++.h> using names ...
- 【线性规划与网络流 24题】已完成(3道题因为某些奇怪的原因被抛弃了QAQ)
写在前面:SDOI2016 Round1滚粗后蒟蒻开始做网络流来自我拯救(2016-04-11再过几天就要考先修课,现在做网络流24题貌似没什么用←退役节奏) 做的题目将附上日期,见证我龟速刷题. 1 ...
- cogs_14_搭配飞行员_(二分图匹配+最大流,网络流24题#01)
描述 http://cojs.tk/cogs/problem/problem.php?pid=14 有一些正飞行员和副飞行员,给出每个正飞行员可以和哪些副飞行员一起飞.一架飞机上必须一正一副,求最多多 ...
- BZOJ_1221_ [HNOI2001]_软件开发(最小费用流,网络流24题#10)
描述 http://www.lydsy.com/JudgeOnline/problem.php?id=1221 n天,每天需要r个毛巾,用完以后可以洗,要么花fa洗a天,要么花fb洗b天,毛巾不够了可 ...
- 【网络流24题】 No.12 软件补丁问题(最小转移代价 最短路)
[题意] T 公司发现其研制的一个软件中有 n 个错误, 随即为该软件发放了一批共 m 个补丁程序. 每一个补丁程序都有其特定的适用环境, 某个补丁只有在软件中包含某些错误而同时又不包含另一些错误时才 ...
- 【网络流24题】No.1 搭配飞行员(飞行员配对方案问题)
[问题描述] 飞行大队有若干个来自各地的驾驶员,专门驾驶一种型号的飞机,这种飞机每架有两个驾驶员,需一个正驾驶员和一个副驾驶员.由于种种原因,例如相互配合的问题,有些驾驶员不能在同一架飞机上飞 ...
- CJOJ 1494 【网络流24题】 搭配飞行员(二分图最大匹配)
CJOJ 1494 [网络流24题] 搭配飞行员(二分图最大匹配) Description 飞行大队有若干个来自各地的驾驶员,专门驾驶一种型号的飞机,这种飞机每架有两个驾驶员,需一个正驾驶员和一个副驾 ...
随机推荐
- C#7.0新特性(VS2017可用)
分享一下其实2016年12月就已经公布了的C#7.0的新特性吧,虽然很早就出来了,但咱这IDE不支持啊.. 不过在昨天的VS2017中已经完美可以支持使用了. E文好的,移步官方介绍地址:https: ...
- 在lnmp1.3布置的web服务器上运行thinkphp3.2.3项目pathinfo路径模式
通过我的经历希望能给大家带来一些帮助: 我是在Linux系统上通过https://lnmp.org/install.html设置Nginx服务器,使用的是lnmp1.3版本,之后将一个thinkphp ...
- yii2 modules模块配置指南
在Yii2 中模块是可以无限级嵌套的,也就是说,模块可以包含另一个包含模块的模块,我们称前者为父模块,后者为子模块, 子模块必须在父模块的yiibaseModule::modules属性中申明,例如: ...
- QWidget 的 close 与 Qt::WA_DeleteOnClose
[1]close 与 Qt::WA_DeleteOnClose简介 1.1 Qt源码 /*! Closes this widget. Returns \c true if the widget was ...
- hive 常见时间日期函数的使用
1.时间戳函数 日期转时间戳:从1970-01-01 00:00:00 UTC到指定时间的秒数 获得当前时区的UNIX时间戳: select unix_timestamp(); 1533716607 ...
- JAVA中异常状况总结
之前在<会当凌绝顶>这本书中学到过对于异常处理的知识,当时也是根据书上的代码,自己进行编写大概知道是怎么回事儿,王老师给我们上了一节课之后,发现异常处理可以发挥很大的作用. 通过在网络上 ...
- Linux基础命令---ping
ping ping指令可以发送ICMP请求到目标地址,如果网络功能正常,目标主机会给出回应信息.ping使用ICMP协议强制发送ECHO_REQUEST报文到目标主机,从主机或网关获取ICMP ECH ...
- git 提交命令
git stash -u 占存本地版本 git commit git fetch 提交 git rebase git stash pop 将本地没有提交的代码暂存,然后切换到其他分支,然后再回到当前分 ...
- this和构造器的内存分析(***)
this的含义: 1.区分成员变量和局部变量: 2.谁调用这个this就指向谁: public class Person{ private String name; private int age; ...
- 2019/3/20 wen 流程控制