洛谷 P3063 [USACO12DEC]牛奶的路由Milk Routing

　　　　洛谷 P3063 [USACO12DEC]牛奶的路由Milk Routing

题目描述

Farmer John's farm has an outdated network of M pipes (1 <= M <= 500) for pumping milk from the barn to his milk storage tank. He wants to remove and update most of these over the next year, but he wants to leave exactly one path worth of pipes intact, so that he can still pump milk from the barn to the storage tank.

The pipe network is described by N junction points (1 <= N <= 500), each of which can serve as the endpoint of a set of pipes. Junction point 1 is the barn, and junction point N is the storage tank. Each of the M bi-directional pipes runs between a pair of junction points, and has an associated latency (the amount of time it takes milk to reach one end of the pipe from the other) and capacity (the amount of milk per unit time that can be pumped through the pipe in steady state). Multiple pipes can connect between the same pair of junction points.

For a path of pipes connecting from the barn to the tank, the latency of the path is the sum of the latencies of the pipes along the path, and the capacity of the path is the minimum of the capacities of the pipes along the path (since this is the "bottleneck" constraining the overall rate at which milk can be pumped through the path). If FJ wants to send a total of X units of milk through a path of pipes with latency L and capacity C, the time this takes is therefore L + X/C.

Given the structure of FJ's pipe network, please help him select a single path from the barn to the storage tank that will allow him to pump X units of milk in a minimum amount of total time.

农民约翰的农场有一套老旧的管网，管网由M条管道(1<=M<=500)构成，用于将牛奶从谷仓运到储奶罐。 他想在明年移除和更新大部分管道，但他想原封不动地保留一条完整的路径，这样他仍然可以把牛奶从谷仓输送到储罐。

管网由N个节点(1<=N<=500)组成，每个点都可以作为一组管道的端点。结点1是谷仓，结点N是储罐。M条双向管道中的每一条都连接一对节点，并且都有一个延迟值(牛奶达到管的另一端的用时)和容量值(单位时间内可以稳定通过管道的牛奶量)。多条管道可以连接同一对节点。

对于一条连接谷仓与储罐的路径，路径的延迟等于沿途所有管道的延迟之和，路径的容量等于沿途管道最小的容量(因为这是制约牛奶运送的“瓶颈”)。如果约翰通过一条延迟为L、容量为C的管道运送X个单位的牛奶，需要的时间为L+X/C。

给出约翰的管网结构，请帮助他选择一条路径，使得他从谷仓到储罐运送X个单位牛奶的总时间最少。

输入输出格式

输入格式：

* Line 1: Three space-separated integers: N M X (1 <= X <= 1,000,000).

* Lines 2..1+M: Each line describes a pipe using 4 integers: I J L C. I and J (1 <= I,J <= N) are the junction points at both ends of the pipe. L and C (1 <= L,C <= 1,000,000) give the latency and capacity of the pipe.

第1行：三个空格分隔的整数:N M X(1<=X<=1000000)。

第2行到第M+1行：每一行描述一条管道，有4个整数：I J L C。I和J(1<=I,J<=N)是这条管道连接的两个点。L和C(1<=L，C<=1000000)是这条管道的延迟和容量。

输出格式：

* Line 1: The minimum amount of time it will take FJ to send milk along a single path, rounded down to the nearest integer.

第1行:约翰沿着一条路径送牛奶花费的最少的时间，向下取整到最近的整数。

输入输出样例

输入样例#1： 复制

3 3 15 
1 2 10 3 
3 2 10 2 
1 3 14 1 

输出样例#1： 复制

27 

说明

FJ wants to send 15 units of milk through his pipe network. Pipe #1 connects junction point 1 (the barn) to junction point 2, and has a latency of 10 and a capacity of 3. Pipes #2 and #3 are similarly defined.

The path 1->3 takes 14 + 15/1 = 29 units of time. The path 1->2->3 takes 20 + 15/2 = 27.5 units of time, and is therefore optimal.

约翰想要通过管网运送15个单位的牛奶。管道1连接节点1（谷仓）和节点2，延迟为10，容量为3。管道2和管道3也以相似的方式来定义。

路径1->3花费14+15/1=29个单位的时间。路径1->2->3花费20+15/2=27.5个单位的时间，用时最少。

 

思路：spfa，但要注意题目中的限制条件（已在上文标出）

难度：普及/提高-

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#define ll long long
#define rp (i-1)*n+j
#define cp (i-1)*n+j+n*n
#define inf 50000000
#define re register
using namespace std;
struct po {
    int nxt, to, w, dis;
}edge[20001];
int n, m, k, num, s, t;
int head[20001], dis[20001], c[20001], b[20001], minn=909909090;
inline void add_edge(int from, int to, int w, int dis) {
    edge[++num].nxt = head[from];
    edge[num].to = to;
    edge[num].w = w;
    edge[num].dis = dis;
    head[from] = num;
}
inline void add(int from, int to, int w, int dis) {
    add_edge(from, to, w, dis);
    add_edge(to, from, w, dis);
}
inline void spfa(int x) {
    queue<int> q;
    memset(dis, 50, sizeof(dis));
    memset(b, 0, sizeof(b));
    q.push(s); dis[s] = 0; b[s] = 1;
    while(!q.empty()) {
        int u = q.front(); q.pop(); b[u] = 0;
        for(re int i = head[u]; i; i = edge[i].nxt) {
            int v = edge[i].to;
            if(dis[v] > dis[u]+edge[i].dis && edge[i].w >= x) {
                dis[v] = dis[u]+edge[i].dis;
                if(!b[v]) { q.push(v); b[v]=1; }
            }
        }
    }
}
int main() {
    cin >> n >> m >> k;
    for(re int i = 1; i <= m; i++) {
        int x, y, d, l;
        cin >> x >> y >> d >> l;
        add(x, y, l, d);
        c[i] = l;
    }
    sort(c+1, c+m+1);
    s = 1, t = n;
    for(re int i = 1; i <= m; i++) {
        spfa(c[i]);
        minn = min(minn,dis[t]+k/c[i]);
    }
    cout<<minn;
    return 0;
}