Description
在一段时间之后,网络公司终于有了一定的知名度,也开始收到一些订单,其中最大的一宗来自B市。Blue Mary决定亲自去签下这份订单。为了节省旅行经费,他的某个金融顾问建议只购买U航空公司的机票。U航空公司的所有航班每天都只有一班,并且都是上午出发当天下午到达的,所以他们每人每天只能坐一班飞机。经过调查,他们得到了U航空公司经营的所有航班的详细信息,这包括每一航班的出发地,目的地以及最多能买到的某一天出发的票数。(注意: 对于一个确定的航班,无论是哪一天,他们最多能买到的那一天出发的票数都是相同的。) Blue Mary注意到他们一定可以只乘坐U航空公司的航班就从A市到达B市,但是,由于每一航班能买到的票的数量的限制,他们所有人可能不能在同一天到达B市。所以现在Blue Mary需要你的帮助,设计一个旅行方案使得最后到达B市的人的到达时间最早。
2 <= N <= 50
1 <= M <= 2450
1 <= T <= 50
1 <= X,Y <= N
X != Y
1 <= Z <= 50
Solution
一开始还以为是动态加边呢。。
首先考虑拆点跑网络流,我们把第i个点拆成s个分别表示不同时刻
注意到答案是单调的,可以二分答案建图最大流
Code
#include <stdio.h>
#include <string.h>
#include <queue>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))
const int INF=0x3f3f3f3f;
const int N=500505;
const int E=5000005;
struct edge {int x,y,w,next;} e[E];
struct data {int fi,se,th;} p[N];
int ls[N],edCnt,n,m,T;
int dis[N],cur[N];
int read() {
int x=0,v=1; char ch=getchar();
for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
return x*v;
}
void add_edge(int x,int y,int w) {
e[++edCnt]=(edge) {x,y,w,ls[x]}; ls[x]=edCnt;
e[++edCnt]=(edge) {y,x,0,ls[y]}; ls[y]=edCnt;
}
bool bfs(int st,int ed) {
std:: queue <int> que;
que.push(st);
rep(i,st,ed) dis[i]=-1; dis[st]=1;
for (;!que.empty();) {
int now=que.front(); que.pop();
for (int i=ls[now];i;i=e[i].next) {
if (e[i].w>0&&dis[e[i].y]==-1) {
dis[e[i].y]=dis[now]+1;
if (e[i].y==ed) return true;
que.push(e[i].y);
}
}
}
return false;
}
int find(int now,int ed,int mn) {
if (now==ed||!mn) return mn;
int ret=0;
for (int &i=cur[now];i;i=e[i].next) {
if (e[i].w>0&&dis[now]+1==dis[e[i].y]) {
int d=find(e[i].y,ed,std:: min(mn-ret,e[i].w));
e[i].w-=d; e[i^1].w+=d; ret+=d;
if (ret==mn) break;
}
}
return ret;
}
int dinic(int st,int ed) {
int ret=0;
for (;bfs(st,ed);) {
rep(i,st,ed) cur[i]=ls[i];
ret+=find(st,ed,INF);
}
return ret;
}
bool check(int mid) {
edCnt=1; fill(ls,0);
rep(i,1,m) {
rep(j,1,mid-1) {
add_edge((p[i].fi-1)*mid+j,(p[i].se-1)*mid+j+1,p[i].th);
}
}
rep(i,1,n) rep(j,1,mid-1) {
add_edge((i-1)*mid+j,(i-1)*mid+j+1,INF);
}
int ret=dinic(1,n*mid);
return (ret>=T);
}
int main(void) {
n=read(),m=read(),T=read();
rep(i,1,m) {
p[i].fi=read(); p[i].se=read();
p[i].th=read();
}
int l=1,r=8000;
for (;l<=r;) {
int mid=(l+r)>>1;
if (check(mid)) r=mid-1;
else l=mid+1;
}
printf("%d\n", r);
return 0;
}