1. 程式人生 > >LibreOJ #539. 「LibreOJ NOIP Round #1」旅遊路線

LibreOJ #539. 「LibreOJ NOIP Round #1」旅遊路線

技術分享 play tin 旅遊 opened getch ostream int -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\n", y-l); } }
View Code

LibreOJ #539. 「LibreOJ NOIP Round #1」旅遊路線