description

• 题面较长，这里给出题目链接

solution

• 考虑预处理出\(f[i][j]\)表示在第\(i\)个点加满油后，从第\(i\)个点出发，至多消耗\(j\)元钱走过的最大路程，那么对于每一个询问就可以二分答案\(O(logq)\)查询了

• 可以得出转移方程\(f[i][k]=\max(f[j][k-p[j]]+g[i][j],f[i][k])\)，其中\(g[i][j]\)表示从在\(i\)点加满油后从\(i\)走到\(j\)能走过的最大路程

• \(g\)可以使用倍增\(floyd\)预处理出\(h[k][i][j]\)表示从\(i\)走到\(j\),至多走\(2^k\)步的最长路程后求出。

• 倍增\(floyd\)代码：

  for(int k=1;k<=16;++k){ memcpy(h[k],h[k-1],sizeof(h[k])); for(int w=1;w<=n;++w) for(int i=1;i<=n;++i) for(int j=1;j<=n;++j) if(h[k-1][i][w]>inf&&h[k-1][w][j]) h[k][i][j]=max(h[k][i][j],h[k-1][i][w]+h[k-1][w][j]);
  }

• 于是我们就可以\(O(n^3logC)\)完成预处理，总时间复杂度\(O(n^3logC+Tlogq)\)

code

#include<bits/stdc++.h>
using namespace std;
const int N=110;
int n,m,C,T,p[N],c[N],b[N],inf;
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
    return x*f;
}
int f[N][N*N];//在i加了油后出发，用j块钱能走的最大路程
int g[N][N];//在i处加了油后从i到j的最长路
int h[20][N][N];//从i走到j,至多走1<<k条路的最长路
int main(){
    n=read();m=read();C=read();T=read();
    for(int i=1;i<=n;++i) p[i]=read(),c[i]=read();
    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 u=read(),v=read(),w=read();
        h[0][u][v]=max(h[0][u][v],w);
    }
    for(int k=1;k<=16;++k){
        memcpy(h[k],h[k-1],sizeof(h[k]));
        for(int w=1;w<=n;++w)
            for(int i=1;i<=n;++i)
                for(int j=1;j<=n;++j)
                    if(h[k-1][i][w]>inf&&h[k-1][w][j]) h[k][i][j]=max(h[k][i][j],h[k-1][i][w]+h[k-1][w][j]);
    }
    for(int i=1;i<=n;++i){
        memset(g[i],192,sizeof(g[i]));g[i][i]=0;
        int s=min(C,c[i]);
        for(int k=0;k<=16;++k){
            if((s>>k)&1){
                memset(b,inf,sizeof(b));
                for(int j=1;j<=n;++j){
                    if(g[i][j]>inf){
                        for(int w=1;w<=n;++w)
                            if(h[k][j][w]>inf) b[w]=max(b[w],g[i][j]+h[k][j][w]);
                    }
                }
                memcpy(g[i],b,sizeof(b));
            }
        }
    }
    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 u=read(),q=read(),d=read();
        int l=p[u],r=q;
        while(l<r){
            int mid=(l+r)>>1;
            if(f[u][mid-p[u]]>=d) r=mid;
            else l=mid+1;
        }
        if(f[u][l-p[u]]<d) puts("-1");
        else printf("%d\n",q-l);
    }
    return 0;
}