LibreOJ #539. 「LibreOJ NOIP Round #1」旅游路线（倍增+二分）

　　哎一开始看错题了啊T T...最近状态一直不对...最近很多傻逼题都不会写了T T

　　考虑距离较大肯定不能塞进状态...钱数<=n^2能够承受， 油量再塞就不行了...显然可以预处理出点i到j走ci步的最长距离（一开始以为一条路耗油为路的长度T T）， 这个是经典题， 倍增求就好了...然后就可以转移了呀T T

　　最后二分就好了呀T T...

　　我怎么这么菜啊， 还写了好久T T...

　　一开始还写成n^4logn TLE了半天没查到错

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=, inf=1e9;
int n, m, C, T, l, r, up, x, y, z;
int f[maxn][maxn][maxn], g[maxn][maxn*maxn], mi[maxn], p[maxn], c[maxn], dis[maxn][maxn], tmp[maxn], dist[maxn];
inline void read(int &k)
{
    int f=; k=; char c=getchar();
    while(c<'' || c>'') c=='-'&&(f=-), c=getchar();
    while(c<='' && c>='') k=k*+c-'', c=getchar();
    k*=f;
}
inline int max(int a, int b){return a>b?a:b;}
inline void solve(int x, int step)
{
    for(int i=;i<=n;i++) tmp[i]=dist[i]=(i==x)?:-inf;
    for(int now=;now<=up;now++)
    if(step&(<<now))
    {
        for(int i=;i<=n;i++)tmp[i]=dist[i];
            for(int k=;k<=n;k++)
                for(int j=;j<=n;j++)
                    dist[j]=max(dist[j], tmp[k]+f[k][j][now]);
    }
    for(int i=;i<=n;i++) dis[x][i]=dist[i];
}
int main()
{
    read(n); read(m); read(C); read(T);
    for(int i=;i<=n;i++) read(p[i]), read(c[i]);
    up=floor(log2(C));
    for(int i=;i<=n;i++) for(int j=;j<=n;j++) for(int k=;k<=up;k++) f[i][j][k]=(i==j)?:-inf;
    for(int i=;i<=m;i++) read(x), read(y), read(z), f[x][y][]=z;
    for(int l=;l<=up;l++)
        for(int k=;k<=n;k++)
            for(int i=;i<=n;i++)
                for(int j=;j<=n;j++)
                    f[i][j][l]=max(f[i][j][l], f[i][k][l-]+f[k][j][l-]);
    for(int i=;i<=n;i++) for(int j=;j<=n;j++) for(int k=;k<=up;k++) f[i][j][k]=max(f[i][j][k], f[i][j][k-]);
    for(int i=;i<=n;i++) solve(i, c[i]);
    for(int i=;i<=n*n;i++)
        for(int j=;j<=n;j++)
        {
            if(i) g[j][i]=g[j][i-];
            if(i>=p[j]) for(int k=;k<=n;k++) g[j][i]=max(g[j][i], g[k][i-p[j]]+dis[j][k]);
        }
    while(T--)
    {
        read(x); read(y); read(z);
        if(g[x][y]<z){puts("-1"); continue;}
        int l=, r=y;
        while(l<r)
        {
            int mid=(l+r)>>;
            if(g[x][mid]>=z) r=mid;
            else l=mid+;
        }
        printf("%d\n", y-l);
    }
}