【题意】给定正边权有向图,车油量上限C,每个点可以花费pi加油至min(C,ci),走一条边油-1,T次询问s点出发带钱q,旅行路程至少为d的最多剩余钱数。

n<=100,m<=1000,C<=10^5,q<=n^2

【算法】动态规划

【题解】官方题解

虽然不是DAG,但是由于q很小的特点,将q加入状态就满足DP的无后效性了。

令f[i][q]表示当前在i点并在i点加油,加油前钱数为q的最大路程。(q<pi时,f[i][q]=0)

假设下一加油点为j,转移方程:f[i][q]=max{f[j][q-pi]+w(i,j,ci)},其中w(i,j,c)表示i点到j点,油量为c的最大路程。(ci=min(ci,C))

现在主要问题是预处理w(i,j,ci),发现每走一条边c只会减少1,符合倍增变化规则一致的特点,考虑图上倍增。

虽然不是DAG,但是c加入状态就满足DP的无后效性了。

g(i,j,k)表示i点到j点,油量为2^k的最大路程,显然g(i,j,k)=max{g(i,x,k-1)+g(x,j,k-1)},x是中转点。

对于w(i,j,ci),将ci拆分二进制,每次枚举x作为中转点后直接取倍增数组g计算答案。

【倍增的思想是很经典的,将需要的ci拆分二进制后将1的位用倍增数组堆起来。但是应用到图上每次就都需要遍历全图作为可能的中转点,最后找到最优答案。】

最后得到了f[i][q],对每个询问在f[s]上二分到第一个大于等于d的f[s][q],q就是答案。

复杂度O(n^4+n^3*log k+T*log n^2),瓶颈复杂度O(n^4)。(n^4过100……)

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cctype>
using namespace std;
int read(){
char c;int s=,t=;
while(!isdigit(c=getchar()))if(c=='-')t=-;
do{s=s*+c-'';}while(isdigit(c=getchar()));
return s*t;
}
const int maxn=,maxm=,maxk=,inf=0x3f3f3f3f;
int n,m,C,T,tot;
int g[maxn][maxn][maxk],p[maxn],c[maxn],A[maxn],B[maxn],w[maxn][maxn],f[maxn][maxn*maxn];
void cmax(int &a,int b){if(b>a)a=b;}
int main(){
n=read();m=read();C=read();T=read();
for(int i=;i<=n;i++)for(int j=;j<=n;j++)for(int k=;k<=;k++)g[i][j][k]=-inf;
for(int i=;i<=n;i++)p[i]=read(),c[i]=min(C,read()),g[i][i][]=;
for(int i=;i<=m;i++){
int u=read(),v=read(),w=read();
g[u][v][]=w;
}
for(int k=;k<=;k++)
for(int i=;i<=n;i++)
for(int j=;j<=n;j++)
for(int x=;x<=n;x++)cmax(g[i][j][k],g[i][x][k-]+g[x][j][k-]);
for(int i=;i<=n;i++){
for(int j=;j<=n;j++)A[j]=-inf;A[i]=;
for(int k=;k<=;k++)if((c[i]>>k)&){
for(int j=;j<=n;j++){
B[j]=-inf;
for(int x=;x<=n;x++)cmax(B[j],A[x]+g[x][j][k]);
}
for(int j=;j<=n;j++)A[j]=B[j];
}
for(int j=;j<=n;j++)w[i][j]=A[j];
}
for(int q=;q<=n*n;q++)
for(int x=;x<=n;x++)if(q>=p[x])
for(int y=;y<=n;y++)cmax(f[x][q],f[y][q-p[x]]+w[x][y]);
while(T--){
int s=read(),q=read(),d=read();
int pl=lower_bound(f[s],f[s]+n*n+,d)-f[s];
if(pl>q)printf("-1\n");else printf("%d\n",q-pl);
}
return ;
}

【LibreOJ】#539. 「LibreOJ NOIP Round #1」旅游路线的更多相关文章

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

    哎一开始看错题了啊T T...最近状态一直不对...最近很多傻逼题都不会写了T T 考虑距离较大肯定不能塞进状态...钱数<=n^2能够承受, 油量再塞就不行了...显然可以预处理出点i到j走c ...

  2. LOJ#539. 「LibreOJ NOIP Round #1」旅游路线

    n<=100,m<=1000的图,在此图上用油箱容量C<=1e5的车来旅行,旅行时,走一条边会耗一单伟油,在点i时,若油量<ci,则可以把油以pi的价格补到ci,pi<= ...

  3. LOJ #539. 「LibreOJ NOIP Round #1」旅游路线 倍增floyd + 思维

    考试的时候是这么想的: 求出每一个点花掉 $i$ 的花费向其他点尽可能走的最长距离,然后二分这个花费,找到第一个大于 $d$ 的就输出$.$然而,我这个记忆化搜索 $TLE$ 的很惨$.$这里讲一下正 ...

  4. 「LOJ 539」「LibreOJ NOIP Round #1」旅游路线

    description 题面较长,这里给出题目链接 solution 考虑预处理出\(f[i][j]\)表示在第\(i\)个点加满油后,从第\(i\)个点出发,至多消耗\(j\)元钱走过的最大路程,那 ...

  5. LibreOj #539. 「LibreOJ NOIP Round #1」旅游路线

    题目链接 做完这道题,我深知当一个问题复杂度过高的时候,把一些可以分离的操作都分散开,可以大幅度降低复杂度..... 发现无论有多少钱,每到一个点后扩展到的距离被限制在 \(min(C, c[i])\ ...

  6. 「LibreOJ NOIP Round #1」旅游路线

    Description T 城是一个旅游城市,具有 nnn 个景点和 mmm 条道路,所有景点编号为 1,2,...,n1,2,...,n1,2,...,n.每条道路连接这 nnn 个景区中的某两个景 ...

  7. LibreOJ #517. 「LibreOJ β Round #2」计算几何瞎暴力

    二次联通门 : LibreOJ #517. 「LibreOJ β Round #2」计算几何瞎暴力 /* LibreOJ #517. 「LibreOJ β Round #2」计算几何瞎暴力 叫做计算几 ...

  8. LibreOJ #528. 「LibreOJ β Round #4」求和

    二次联通门 : LibreOJ #528. 「LibreOJ β Round #4」求和 /* LibreOJ #528. 「LibreOJ β Round #4」求和 题目要求的是有多少对数满足他们 ...

  9. LibreOJ #527. 「LibreOJ β Round #4」框架

    二次联通门 : LibreOJ #527. 「LibreOJ β Round #4」框架 /* LibreOJ #527. 「LibreOJ β Round #4」框架 %% xxy dalao 对于 ...

随机推荐

  1. Java中的线程状态转换和线程控制常用方法

    Java 中的线程状态转换: [注]:不是 start 之后就立刻开始执行, 只是就绪了(CPU 可能正在运行其他的线程). [注]:只有被 CPU 调度之后,线程才开始执行, 当 CPU 分配给你的 ...

  2. 【SQLAlchemy】SQLAlchemy修改查询字段列名

    SQLAlchemy问题记录 company price quantity Microsoft Google Google Google 要实现脚本 select price, sum(quantit ...

  3. Qt快速入门学习笔记(画图篇)

    1.Qt中提供了强大的2D绘图系统,可以使用相同的API在屏幕和绘图设备上进行绘制,它主要基于QPainter.QPaintDevice和QPaintEngine这三个类.其中QPainter用来执行 ...

  4. C# 大文件的复制方法

    如何复制读取大文件,也许困惑了很多人很长时间,这个不知道怎么搞,的确让人头疼欲裂,知道了你就才发现原来那么简单,话不多说,直入正题```` static void Main(string[] args ...

  5. [CF1037H] Security

    题目链接 codeforces. 洛谷. Solution 按照套路,可以\(SAM\)上线段树合并求出\(endpos\)集合,然后随便贪心一下就好了. #include<bits/stdc+ ...

  6. UOJ228:基础数据结构练习题——题解

    http://uoj.ac/problem/228 参考:https://www.cnblogs.com/ljh2000-jump/p/6357583.html 考虑当整个区间的最大值开方==最小值开 ...

  7. BZOJ 1016 Windy 数 | 数位DP

    题目: http://www.lydsy.com/JudgeOnline/problem.php?id=1026 题解: f[i][j][1/0]表示枚举到第i位,这位开头是j,当前的数大于(1)或小 ...

  8. HDU 2710

    Max Factor Time Limit: 2000/1000 MS(Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total ...

  9. ES6 Set,WeakSet,Map,WeakMap

    1. Set Set是一个集合,里面的值都是唯一的,没有重复的.Set中可以是任何数据类型,并且添加数据时会进行严格比较,重复数据无法加入. 2. WeakSet 弱引用Set.只能存储对象,不能存储 ...

  10. setTimeout()的应用

    错误写法:setTimeout(window.close(),5000); 正确写法:setTimeout(window.close,5000); 或者 setTimeout(function(){ ...