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=110, 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=1; k=0; char c=getchar();
    while(c<'0' || c>'9') c=='-'&&(f=-1), c=getchar();
    while(c<='9' && c>='0') k=k*10+c-'0', 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=1;i<=n;i++) tmp[i]=dist[i]=(i==x)?0:-inf;
    for(int now=0;now<=up;now++)
    if(step&(1<<now))
    {
        for(int i=1;i<=n;i++)tmp[i]=dist[i];
            for(int k=1;k<=n;k++)
                for(int j=1;j<=n;j++)
                    dist[j]=max(dist[j], tmp[k]+f[k][j][now]);
    }
    for(int i=1;i<=n;i++) dis[x][i]=dist[i];
}
int main()
{
    read(n); read(m); read(C); read(T);
    for(int i=1;i<=n;i++) read(p[i]), read(c[i]);
    up=floor(log2(C));
    for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) for(int k=0;k<=up;k++) f[i][j][k]=(i==j)?0:-inf;
    for(int i=1;i<=m;i++) read(x), read(y), read(z), f[x][y][0]=z;
    for(int l=1;l<=up;l++)
        for(int k=1;k<=n;k++)
            for(int i=1;i<=n;i++)
                for(int j=1;j<=n;j++)
                    f[i][j][l]=max(f[i][j][l], f[i][k][l-1]+f[k][j][l-1]);
    for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) for(int k=1;k<=up;k++) f[i][j][k]=max(f[i][j][k], f[i][j][k-1]);
    for(int i=1;i<=n;i++) solve(i, c[i]);
    for(int i=0;i<=n*n;i++)
        for(int j=1;j<=n;j++)
        {
            if(i) g[j][i]=g[j][i-1];
            if(i>=p[j]) for(int k=1;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=0, r=y;
        while(l<r)
        {
            int mid=(l+r)>>1;
            if(g[x][mid]>=z) r=mid;
            else l=mid+1;
        }
        printf("%d
", y-l);
    }
}