厦门大学“网宿杯“17届程序设计竞赛决赛 E芜湖起飞(三分+最短路)

题目描述：
安徽芜湖有n个机场，一共有m条线路在空管部门报备。
每条线路单向连接两个机场，并且需要的通行时间每天都可能不一样。具体来说，设目前是第x天，那么第i条线路所需要的通行时间为kix+bi
一年一共有H天，也就是说，x取[0,H]中的整数。现在大司马想从1号机场在一天内换乘任意多次航班前往n号机场，他总是选择用时最短的方式，现在他想知道哪一天需要花最长的时间。
思路：和普通的最短路不同，这个题目的路径是变化的，因为bi和ki是确定的，所以路径长度只随着天数的变化而变化。题目要求的是0-h天里花费时间最长的那一天。因为路径变化并没有规律，首先想到的就是暴力。但是h最大可达到1e9所以可以排除了。后来考虑到bi和ai都是固定的，想了一下最短的情况应该就是第0天或者第h天(瞎想的没有什么理论依据)。那么答案应该在0-h中间，0-h的最短路可以看作是一个向上凸的二次函数曲线，可以用三分的方法求出最高点。三分设左边界为l，右边界为r，ml=（l+r)>>1，mr=（r+ml）>>1;如果ml天的最短路大于mr天的最短路，那么答案可能的区间可以缩小到[l,mr]，反之区间缩小为[ml,r]。mr-ml<10后，直接暴力把答案精确的求出来就行了

#include<bits/stdc++.h>
using namespace std;
typedef long long int ll;
const int maxm=114614;
const int maxn=114614;
const ll inf=0x3f3f3f3f3f3f3f;
int head[maxn],vis[maxn];
int cnt,n,m;
ll dis[maxn],ans;
struct node
{
    int to,next,l;
}e[maxm],ee[maxm];
struct edge
{
    int u,v,a,b;
}sj[maxn];
void add(int x,int y,int l)
{
    e[cnt].to=y;
    e[cnt].next=head[x];
    e[cnt].l=l;
    head[x]=cnt++;
}
void spfa(int s)
{
    queue<int> q;
    for(int i=0;i<=n;i++)
    {
        vis[i]=0;
        dis[i]=inf;
    }
    dis[s]=0;
    vis[s]=1;
    q.push(s);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=head[u];i!=-1;i=e[i].next)
        {
            int v=e[i].to;
            if(dis[v]>dis[u]+e[i].l)
            {
                dis[v]=dis[u]+e[i].l;
                if(!vis[v])
                {
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
}
ll solve(int x)
{
    cnt=0;
    memset(head,-1,sizeof(head));
    for(int i=1;i<=m;i++)
    {
        add(sj[i].u,sj[i].v,sj[i].a*x+sj[i].b);
    }
    spfa(1);
    ans=max(ans,dis[n]);
    return dis[n];
}
int main()
{
    int h,u,v,k,b;
    scanf("%d%d%d",&n,&m,&h);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d%d",&sj[i].u,&sj[i].v,&sj[i].a,&sj[i].b);
    }
    int l=0,r=h;
    while(r-l>10)
    {
        int lm=(l+r)>>1;
        int rm=(lm+r)>>1;
        if(solve(lm)>solve(rm))r=rm;
        else l=lm;

    }
    for(int i=l;i<=r;i++)solve(i);
    printf("%lld\n",ans);
    return 0;
}