【BZOJ1570】【JSOI2008】Blue Mary的旅行 最大流check
广告:
#include <stdio.h>
int main()
{
puts("转载请注明出处[vmurder]谢谢");
puts("网址:blog.csdn.net/vmurder/article/details/44099255");
}
题解:
枚举需要几天,然后每天新加一层点跑最大流看maxflow是否够数。
代码:
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define P 55
#define N 5010
#define M 1001000
#define inf 0x3f3f3f3f
/*
2 <= N <= 50
1 <= M <= 2450
1 <= T <= 50
1 <= X,Y <= N
X != Y
1 <= Z <= 50
*/
using namespace std;
struct KSD
{
int v,len,next;
}e[M];
int head[N],cnt;
int ru[2500],rv[2500],rl[2500];
inline void add(int u,int v,int len)
{
e[++cnt].v=v;
e[cnt].len=len;
e[cnt].next=head[u];
head[u]=cnt;
}
inline void Add(int u,int v,int len){add(u,v,len),add(v,u,0);}
int s,t,d[N];
queue<int>q;
bool bfs()
{
while(!q.empty())q.pop();
memset(d,0,sizeof d);
int i,u,v;
q.push(s),d[s]=1;
while(!q.empty())
{
u=q.front(),q.pop();
for(i=head[u];i;i=e[i].next)
{
if(!d[v=e[i].v]&&e[i].len)
{
d[v]=d[u]+1;
if(v==t)return 1;
q.push(v);
}
}
}
return 0;
}
int dinic(int x,int flow)
{
if(x==t)return flow;
int remain=flow,i,v,k;
for(i=head[x];i&&remain;i=e[i].next)
{
if(d[v=e[i].v]==d[x]+1&&e[i].len)
{
k=dinic(v,min(remain,e[i].len));
if(!k)d[v]=0;
e[i].len-=k,e[i^1].len+=k;
remain-=k;
}
}
return flow-remain;
}
int n,m,p,maxflow;
int id[105][P],num;
void build(int p)
{
int i,j,k;
for(i=1;i<=n;i++)Add(id[p-1][i],id[p][i]=++num,inf);
for(i=1;i<=m;i++)Add(id[p-1][ru[i]],id[p][rv[i]],rl[i]);
Add(id[p][n],t,inf);
}
int main()
{
freopen("test.in","r",stdin);
int i,j,k;
int a,b,c;
scanf("%d%d%d",&n,&m,&p),cnt=1;
if(n==1)
{
puts("0");
return 0;
}
s=0,t=num=1,Add(s,2,p);
for(i=1;i<=m;i++)scanf("%d%d%d",&ru[i],&rv[i],&rl[i]);
for(i=1;i<=n;i++)id[1][i]=++num;
for(i=2;i<=n+p;i++)
{
build(i);
while(bfs())maxflow+=dinic(s,inf);
if(maxflow==p)
{
printf("%d\n",i-1);
return 0;
}
}
puts("-1");
return 0;
}
郑重声明:本站内容如果来自互联网及其他传播媒体,其版权均属原媒体及文章作者所有。转载目的在于传递更多信息及用于网络分享,并不代表本站赞同其观点和对其真实性负责,也不构成任何其他建议。