【bzoj1570】[JSOI2008]Blue Mary的旅行 动态加边网络流

题目描述

在一段时间之后,网络公司终于有了一定的知名度,也开始收到一些订单,其中最大的一宗来自B市。Blue Mary决定亲自去签下这份订单。为了节省旅行经费,他的某个金融顾问建议只购买U航空公司的机票。U航空公司的所有航班每天都只有一班,并且都是上午出发当天下午到达的,所以他们每人每天只能坐一班飞机。经过调查,他们得到了U航空公司经营的所有航班的详细信息,这包括每一航班的出发地,目的地以及最多能买到的某一天出发的票数。(注意: 对于一个确定的航班,无论是哪一天,他们最多能买到的那一天出发的票数都是相同的。) Blue Mary注意到他们一定可以只乘坐U航空公司的航班就从A市到达B市,但是,由于每一航班能买到的票的数量的限制,他们所有人可能不能在同一天到达B市。所以现在Blue Mary需要你的帮助,设计一个旅行方案使得最后到达B市的人的到达时间最早。

输入

第一行包含3个正整数N,M和T。题目中会出现的所有城市分别编号为1,2,…,N,其中城市A编号一定为1,城市B编号一定为N. U公司一共有M条(单向)航班。而连Blue Mary在内,公司一共有T个人要从A市前往B市。 以下M行,每行包含3个正整数X,Y,Z, 表示U公司的每一条航班的出发地,目的地以及Blue Mary最多能够买到的这一航班某一天出发的票数。(即:无论是哪一天,Blue Mary最多只能买到Z张U航空公司的从城市X出发到城市Y的机票。) 输入保证从一个城市到另一个城市的单向航班最多只有一个。

输出

仅有一行,包含一个正整数,表示最后到达B市的人的最早到达时间。假设他们第一次乘飞机的那一天是第一天。

样例输入

3 3 5
1 2 1
2 3 5
3 1 4

样例输出

6



题解

动态加边网络流

由于题目中限定一人每天只能坐一次飞机,所以可以将飞机的运行时间看作一天,出发时间看作第0天。

这样就可以拆点,将每个除1和n以外的点拆成totaltime+1个,看作0~totaltime+1层。注意点1和点n不需要拆,即1和n无论在多少层,编号始终为1和n。

首先加S->1,容量为T,代表T个人。

然后枚举时间t,设点x在第k层所拆的点为(x,k),则加边(i,t-1)->(i,t),容量为inf;对于原图中的边x->y(z),加边(x,t-1)->(y,t),容量为z。

这样动态加边后把n看做汇点,跑dinic,判断一下是否满流即可。

#include <cstdio>
#include <cstring>
#include <queue>
#define N 10010
#define M 1000010
#define inf 0x3f3f3f3f
#define pos(i , j) (i == 1 || i == n) ? i : (n - 2) * (j) + i + 1
using namespace std;
queue<int> q;
int head[N] , to[M] , val[M] , next[M] , cnt = 1 , s , t , dis[N] , dx[M] , dy[M] , dz[M];
void add(int x , int y , int z)
{
	to[++cnt] = y , val[cnt] = z , next[cnt] = head[x] , head[x] = cnt;
	to[++cnt] = x , val[cnt] = 0 , next[cnt] = head[y] , head[y] = cnt;
}
bool bfs()
{
	int x , i;
	memset(dis , 0 , sizeof(dis));
	while(!q.empty()) q.pop();
	dis[s] = 1 , q.push(s);
	while(!q.empty())
	{
		x = q.front() , q.pop();
		for(i = head[x] ; i ; i = next[i])
		{
			if(val[i] && !dis[to[i]])
			{
				dis[to[i]] = dis[x] + 1;
				if(to[i] == t) return 1;
				q.push(to[i]);
			}
		}
	}
	return 0;
}
int dinic(int x , int low)
{
	if(x == t) return low;
	int temp = low , i , k;
	for(i = head[x] ; i ; i = next[i])
	{
		if(val[i] && dis[to[i]] == dis[x] + 1)
		{
			k = dinic(to[i] , min(temp , val[i]));
			if(!k) dis[to[i]] = 0;
			val[i] -= k , val[i ^ 1] += k;
			if(!(temp -= k)) break;
		}
	}
	return low - temp;
}
int main()
{
	int n , m , k , i , j , ans = 0;
	scanf("%d%d%d" , &n , &m , &k);
	s = 0 , t = n , add(s , 1 , k);
	for(i = 1 ; i <= m ; i ++ ) scanf("%d%d%d" , &dx[i] , &dy[i] , &dz[i]);
	for(i = 1 ; i <= n + k ; i ++ )
	{
		for(j = 2 ; j < n ; j ++ ) add(pos(j , i - 1) , pos(j , i) , inf);
		for(j = 1 ; j <= m ; j ++ ) add(pos(dx[j] , i - 1) , pos(dy[j] , i) , dz[j]);
		while(bfs()) ans += dinic(s , inf);
		if(ans == k)
		{
			printf("%d\n" , i);
			return 0;
		}
	}
	return 0;
}
时间: 2024-10-17 18:35:25

【bzoj1570】[JSOI2008]Blue Mary的旅行 动态加边网络流的相关文章

BZOJ1570 [JSOI2008]Blue Mary的旅行

建分层图,每一层表示一天的情况 从S向第0层的1号点连边,每层的n向T连INF的边 枚举天数,每多一天就多建一层然后跑最大流,如果当前流量大于人数则输出答案 由于路径长度不会超过n,因此tot个人走这条路径总天数不会超过tot + n,故只需要建tot + n层即可 1 /************************************************************** 2 Problem: 1570 3 User: rausen 4 Language: C++ 5 R

[JSOI2008]Blue Mary的旅行

[**[JSOI2008]Blue Mary的旅行**](https://www.luogu.org/problem/P4400) #### 题解: 看见这些人流动的过程,我们很容易想到网络流,看见最早到达时间,~~我们很容易想到费用流~~. 我们可以看出题目中所说的z对应着网络流中的容量,但是一天过后又重新可以使用这一条边,我们显然不能直接跑网络流,对于一个人,它的贡献为$1$,所以这道题可以转化为一个最大流的模型. 考虑整个问题,朝贪心的方面去想,要尽量的选择一条$s$到$t$的最短路径,选

【bzoj1507】 JSOI2008—Blue Mary的旅行

http://www.lydsy.com/JudgeOnline/problem.php?id=1570 (题目链接) 题意 给出$m$个航班,每天只能做一次飞机,有$T$人从起点到终点,问最晚到达的人最早什么时候到. Solution 枚举答案分层建图最大流判断即可.之前的流量不要清空. 细节 ? 代码 // bzoj1570 #include<algorithm> #include<iostream> #include<cstdlib> #include<cs

BZOJ 1570 JSOI2008 Blue Mary的旅行 网络流

题目大意:给定一张有向图,每条边每天最多经过有限次,一个人每天只能经过一条边,T个人从1号点出发,问多少天之后能到达n点 将图分层,每一天分作一层,每一层的点向下一层连边 从源点向第0层的1号点连边 每层的n向T连INF的边 从1开始枚举天数,每多一天就多建一层然后跑最大流,如果当前T个人已经能到达点n则输出答案 由于1~n的路径长度不会超过n,因此T个人排队走这条路径总天数不会超过T+n 故只需要建n+T层即可出解 点数O(n^2+nT) 边数O(mn+mT) 都不是很大 可以跑出来 #inc

BZOJ 1567: [JSOI2008]Blue Mary的战役地图

二次联通门 : BZOJ 1567: [JSOI2008]Blue Mary的战役地图 /* BZOJ 1567: [JSOI2008]Blue Mary的战役地图 社会我栋哥 人怂P话多 暴力能A题 正解能WA0 */ #include <cstdio> #include <iostream> #define rg register inline void read (int &n) { rg char c = getchar (); for (n = 0; !isdig

BZOJ 1568: [JSOI2008]Blue Mary开公司(超哥线段树)

1568: [JSOI2008]Blue Mary开公司 Time Limit: 15 Sec  Memory Limit: 162 MBSubmit: 1080  Solved: 379[Submit][Status][Discuss] Description Input 第一行 :一个整数N ,表示方案和询问的总数. 接下来N行,每行开头一个单词“Query”或“Project”. 若单词为Query,则后接一个整数T,表示Blue Mary询问第T天的最大收益. 若单词为Project,则

数据结构(线段树):BZOJ 1568 [JSOI2008]Blue Mary开公司

1568: [JSOI2008]Blue Mary开公司 Time Limit: 15 Sec  Memory Limit: 162 MBSubmit: 602  Solved: 214[Submit][Status][Discuss] Description Input 第 一行 :一个整数N ,表示方案和询问的总数. 接下来N行,每行开头一个单词“Query”或“Project”. 若单词为Query,则后接一个整数T,表示Blue Mary询问第T天的最大收益. 若单词为Project,则

bzoj千题计划219:bzoj1568: [JSOI2008]Blue Mary开公司

http://www.lydsy.com/JudgeOnline/problem.php?id=1568 写多了就觉着水了... #include<cstdio> #include<iostream> #include<algorithm> using namespace std; #define N 100001 double a[N<<2],b[N<<2]; bool have[N<<2]; long long ans; void

1568: [JSOI2008]Blue Mary开公司

1568: [JSOI2008]Blue Mary开公司 题目描述 传送门 题目分析 简单分析可以发现就是不停给出了\(n\)条直线,要求每次给出一条直线后求出所有直线在横坐标为\(x\)时\(y\)的最大值. 李超树裸题. 不知道李超树的可以移步百度. 是代码呢 #include <bits/stdc++.h> using namespace std; #define ls rt<<1 #define rs rt<<1|1 #define mid ((l+r)/2)