LOJ#539.「LibreOJ NOIP Round #1」旅游路线[DP 倍增Floyd]

题目链接：传送门

题目大意:详见原题，比较复杂

题目思路：首先考虑有多少个因素，有出发点，钱的数量，油耗，走的距离。因为这些数据的范围较大，不可能全部放入数组进行dp推导。

因此我们考虑f[i][j]表示从 i出发，已经在i加了油，带了j 块钱，能走的最大距离。 

考虑如何转移，注意到油量比较大，所有不可能把它记到状态里。那我们就枚举下一次在哪里加油： 

 

 f[i][k]=max(f[i][k],f[j][k-p[j]]+g[i][j]);

接下来，第二部分我们需要计算g[i][j].

其中 g[i][j] 表示从 i走到 k，至多走 min{C,ci} 条边的最长路。 
现在我们只需解决怎么预处理出 gg 数组。 
即，第三部分，我们需要倍增 floydfloyd 。

预处理出 h[k][i][j] 表示从i 走到 j ,至多走 2^k 步的最长路，然后就能 O(n3logC) 求了。

#include<bits/stdc++.h>
using namespace std;
const int maxn =180;
typedef long long ll;
int n,m,C,T;
int p[maxn],c[maxn];//费用和油量
int s,q,d;
int h[18][maxn][maxn];
int INF;
int a[maxn],b[maxn];
int g[maxn][maxn];
int f[maxn][maxn*maxn]; //fi,j 表示从 i 出发，已经在 i 加了油，带了 j 块钱，能走多远。
void dp(int s)//预处理g数组，g[i][j] 表示从 i 走到 j，至多走 min{C,ci} 条边的最长路。
{
    memset(a,192,sizeof a);
    a[s]=0;
    for(int k=0;k<=16;k++)
    {
        if((min(C,c[s])>>k)&1)
        {
            memset(b,192,sizeof b);
            for(int i=1;i<=n;i++)
            {
                if(a[i]>INF)
                for(int j=1;j<=n;j++)
                {
                    if(h[k][i][j]>INF)
                    {
                        b[j]=max(b[j],a[i]+h[k][i][j]);
                    }
                }
            }
            memcpy(a,b,sizeof b);
        }
    }
    for(int i=1;i<=n;i++) g[s][i]=a[i];
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&C,&T);
    for(int i=1;i<=n;i++)scanf("%d%d",&p[i],&c[i]);
    memset(h,192,sizeof h);
    INF=h[0][0][0];
    for(int i=1;i<=n;i++)h[0][i][i]=0;
    for(int i=1;i<=m;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        h[0][x][y]=max(h[0][x][y],z);
    }
    for(int k=1;k<=16;k++)//倍增Floyd的预处理
    {
        memcpy(h[k],h[k-1],sizeof h[k]);
        for(int l=1;l<=n;l++)
        {
            for(int i=1;i<=n;i++)
            {
                for(int j=1;j<=n;j++)
                {
                    if(h[k-1][i][l]>INF&&h[k-1][l][j])
                    {
                        h[k][i][j]=max(h[k][i][j],h[k-1][i][l]+h[k-1][l][j]);
                    }
                }
            }
        }
    }
    for(int i=1;i<=n;i++)
    {
        dp(i);
    }
    for(int k=0;k<=n*n;k++)
    {
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                f[i][k]=max(f[i][k],g[i][j]);
            }
        }
    }
    for(int k=0;k<=n*n;k++)
    {
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                if(k>=p[j])
                {
                    f[i][k]=max(f[i][k],f[j][k-p[j]]+g[i][j]);
                }
            }
        }
    }
    while(T--)
    {
        int s,q,d;
        scanf("%d%d%d",&s,&q,&d);
        int ans=-1;
        int l=p[s],r=q;
        while(l<=r)
        {
            int mid=(l+r)>>1;
            if(f[s][mid-p[s]]>=d)
            {
                r=mid-1;
                ans=q-mid;
            }
            else
            {
                l=mid+1;
            }
        }
        printf("%d\n",max(-1,ans));
    }
    return 0;
}