P3063 【[USACO12DEC]牛奶的路由Milk Routing】【Spfa算法】

题目链接

【Spfa算法】

  这道题最初的时候，我想到的做法是最优比例生成树，后来想一想，找不到该如何求设那个答案mid，于是看了讨论，基本都是用最短路来求这道题，于是写了Spfa算法。

  这里除了Spfa算法，最重要的是路径的代替，我求到达每个点的最短路，然后把两点间的路径用L+X/min(目前节点的C，该点的对应原C)来表示，就是所求ans。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 1e9+7
using namespace std;
typedef long long ll;
const int maxN=505;
int N,M,X;      //节点1～N，M条边，X为总共需要运送的量
struct node
{
    int nex, l, c;
    node(int a=0, int b=0, int d=0):nex(a), l(b), c(d) {}
};
vector<node> vt[maxN];
struct point
{
    int l,c,time;
    point(int a=0, int b=0, int c=0):l(a), c(b), time(c) {}
}dis[maxN];
queue<int> Q;
int calc(int L, int C) { return L+X/C; }
void Spfa(int u)
{
    bool inq[maxN]; memset(inq, false, sizeof(inq));        //是否在队列中
    Q.push(u);
    inq[u]=true;
    while(!Q.empty())
    {
        u=Q.front();
        Q.pop();
        inq[u]=false;
        int len=(int)vt[u].size();
        for(int i=0; i<len; i++)
        {
            int v=vt[u][i].nex;
            int k=calc(dis[u].l+vt[u][i].l, min(dis[u].c, vt[u][i].c));
            if(dis[v].time>k)
            {
                dis[v].time=k;
                dis[v].l=dis[u].l+vt[u][i].l;
                dis[v].c=min(dis[u].c, vt[u][i].c);
                if(!inq[v])
                {
                    Q.push(v);
                    inq[v]=true;
                }
            }
        }
    }
    printf("%d\n",dis[N].time);
}
int main()
{
    while(scanf("%d%d%d",&N,&M,&X)!=EOF)
    {
        for(int i=1; i<=N; i++) { vt[i].clear(); dis[i]=( i==1?point(0, INF, 0):point(INF, INF, INF) ); }
        while(!Q.empty()) Q.pop();
        for(int i=1; i<=M; i++)
        {
            int e1,e2,e3,e4;
            scanf("%d%d%d%d",&e1,&e2,&e3,&e4);
            vt[e1].push_back(node(e2, e3, e4));
            vt[e2].push_back(node(e1, e3, e4));
        }
        Spfa(1);
    }
    return 0;
}

 

 

我的洛谷博客