BZOJ1570Blue Mary的旅行

最新推荐文章于 2019-03-12 10:55:52 发布

DCrusher

最新推荐文章于 2019-03-12 10:55:52 发布

阅读量520

点赞数

分类专栏：最大流

本文链接：https://blog.csdn.net/DCrusher/article/details/50899107

版权

最大流专栏收录该内容

19 篇文章 0 订阅

订阅专栏

1570: [JSOI2008]Blue Mary的旅行
Time Limit: 15 Sec Memory Limit: 162 MB
Submit: 288 Solved: 150
Description
在一段时间之后，网络公司终于有了一定的知名度，也开始收到一些订单，其中最大的一宗来自B市。Blue Mary决定亲自去签下这份订单。为了节省旅行经费，他的某个金融顾问建议只购买U航空公司的机票。U航空公司的所有航班每天都只有一班，并且都是上午出发当天下午到达的，所以他们每人每天只能坐一班飞机。经过调查，他们得到了U航空公司经营的所有航班的详细信息，这包括每一航班的出发地，目的地以及最多能买到的某一天出发的票数。(注意: 对于一个确定的航班，无论是哪一天，他们最多能买到的那一天出发的票数都是相同的。) Blue Mary注意到他们一定可以只乘坐U航空公司的航班就从A市到达B市，但是，由于每一航班能买到的票的数量的限制，他们所有人可能不能在同一天到达B市。所以现在Blue Mary需要你的帮助，设计一个旅行方案使得最后到达B市的人的到达时间最早。
Input
第一行包含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的机票。) 输入保证从一个城市到另一个城市的单向航班最多只有一个。
Output
仅有一行，包含一个正整数，表示最后到达B市的人的最早到达时间。假设他们第一次乘飞机的那一天是第一天。
Sample Input
3 3 5
1 2 1
2 3 5
3 1 4
Sample Output
6
HINT
约定：
2 <= N <= 50
1 <= M <= 2450
1 <= T <= 50
1 <= X，Y <= N
X != Y
1 <= Z <= 50
最大流分层建图。。
虽然一眼最大流，但是yy了一下建图并没有rush出正解。。
于是看了Po姐姐的建图思路。。
题目大意：给定一张有向图，每条边每天最多经过有限次，一个人每天只能经过一条边，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) 都不是很大可以跑出来
第一次（相当于默写出）最大流。。
附上本蒟蒻的代码：

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define T 5009
#define inf 0x3f3f3f3f
int n,m,t,head,tail,dis[1000001],h[1000001],q[1000001],cnt=1,ans,sum;
struct kx
{
    int to,next,v;
}edge[1000001];
struct data
{
    int x,y,z;
}e[100001];

int read()
{
    int w=0,c=1; char ch=getchar();
    while (ch<'0' || ch>'9')
      {
        if (ch=='-') c=-1;
        ch=getchar();
      }
    while (ch>='0' && ch<='9')
      w=w*10+ch-'0',ch=getchar();
    return w*c;
}

void add(int u,int v,int w)
{
    cnt++,edge[cnt].next=h[u],h[u]=cnt,edge[cnt].to=v,edge[cnt].v=w;
    cnt++,edge[cnt].next=h[v],h[v]=cnt,edge[cnt].to=u,edge[cnt].v=0;
}

bool bfs()
{
    int j,p;
    memset(dis,-1,sizeof(dis));
    head=0,tail=1;
    dis[0]=q[0]=0;
    while (head<tail)
      {
        head++,j=q[head],p=h[j];
        while (p)
          {
            if (dis[edge[p].to]<0 && edge[p].v>0)
              {
                dis[edge[p].to]=dis[j]+1;
                tail++,q[tail]=edge[p].to;
              }
            p=edge[p].next;
          }
      }
    if (dis[T]>0)
      return true;
    else
      return false;
}

int dfs(int x,int f)
{
    int w,used=0,i=h[x];
    if (x==T)
      return f;
    while (i)
      {
        if (edge[i].v && dis[edge[i].to]==dis[x]+1)
          {
            w=f-used,w=dfs(edge[i].to,min(w,edge[i].v));
            edge[i].v-=w,edge[i^1].v+=w;
            used+=w;
            if (used==f)
              return f;
          }
        i=edge[i].next;
      }
    if (!used)
      dis[x]=-1;
    return used;
}

int main()
{
    int i,j;
    n=read(),m=read(),t=read();
    for (i=1;i<=m;i++)
      e[i].x=read(),e[i].y=read(),e[i].z=read();
    add(0,1,t);
    for (i=1;i<=n+t;i++)
      {
        for (j=1;j<=m;j++)
          add(i*n-n+e[j].x,i*n+e[j].y,e[j].z);
        for (j=1;j<=n;j++)
          add(i*n-n+j,i*n+j,inf);
        add(i*n+n,T,inf);
        while (bfs())
          while (sum=dfs(0,inf))
            ans+=sum;
        if (ans==t)
          {
            printf("%d",i);
            break;
          }
      }
    return 0;
}