P7294-[USACO21JAN]Minimum Cost Paths P【单调栈】

正题

题目链接:https://www.luogu.com.cn/problem/P7294


题目大意

n×mn\times mn×m的网格,当你在(x,y)(x,y)(x,y)时你有两种选择

  1. 花费x2x^2x2的代价向右移动
  2. 花费cyc_ycy的代价向下移动

qqq次询问(1,1)(1,1)(1,1)走到(x,y)(x,y)(x,y)的最小代价。

1≤n≤109,1≤m,q≤2×1051\leq n\leq 10^9,1\leq m,q\leq 2\times 10^51n109,1m,q2×105


解题思路

假设我们开始都是直接走最上面向右的路(也就是代价都是121^212)。

然后考虑在某个位置(x,y)(x,y)(x,y)要走到(n,m)(n,m)(n,m)时向下走会产生的贡献为cy+(m−y)×(2x+1)c_y+(m-y)\times (2x+1)cy+(my)×(2x+1)(后面要抬m−ym-ymy个横着走,然后从x2x^2x2(x+1)2(x+1)^2(x+1)2要到加2x+12x+12x+1)。

然后拆一下就是cy+m−y+2xm−2xyc_y+m-y+2xm-2xycy+my+2xm2xy。发现斜率−2y-2y2y是按照yyy递增递减的,而且我们要求选出的若干个cyc_ycyyyy一定要递增,但是这样的话我们选出来的一定是递增的所以我们只需要考虑对于每个xxx选择一个yyy使得最小化cy+(m−y)×(2x+1)c_y+(m-y)\times (2x+1)cy+(my)×(2x+1)

按照询问排序,一个一个加入新的cyc_ycy,按照斜率维护一个上凸壳,然后在凸壳上面二分就好了。

时间复杂度O(m+qlog⁡m)O(m+q\log m)O(m+qlogm)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=2e5+10;
struct node{ll n,m,id;
}q[N];
ll n,m,t,top,b[N],k[N],s[N],sum[N],l[N],r[N],ans[N];
ll calc(ll i,ll j)//i<j
{return (b[j]-b[i]+k[i]-k[j]-1)/(k[i]-k[j]);}
ll cap(ll i,ll j){ll x=calc(i,j);if(x<=l[i])return 1;return 0;
}
ll getf(ll x,ll l,ll r)
{return (r+l)*(r-l+1)/2*k[x]+b[x]*(r-l+1);}
bool cmp(node x,node y)
{return x.m<y.m;}
signed main()
{
//	freopen("data.in","r",stdin);
//	freopen("data.out","w",stdout);scanf("%lld%lld",&n,&m);for(ll i=1;i<=m;i++){scanf("%lld",&b[i]);b[i]=b[i]-i;k[i]=2*(-i);}scanf("%lld",&t);for(ll i=1;i<=t;i++){scanf("%lld%lld",&q[i].n,&q[i].m);q[i].id=i;}sort(q+1,q+1+t,cmp);for(ll i=1,z=1;i<=m;i++){r[i]=n;while(top>0&&cap(s[top],i))top--;if(top)l[i]=max(calc(s[top],i),1ll);else l[i]=1;if(l[i]<=r[i]){r[s[top]]=l[i]-1;s[++top]=i;sum[top-1]=((top>1)?sum[top-2]:0)+getf(s[top-1],l[s[top-1]],r[s[top-1]]);sum[top]=sum[top-1]+getf(i,l[i],r[i]);}	while(z<=t&&q[z].m<=i){ll qn=q[z].n-1,L=1,R=top;while(L<=R){ll mid=(L+R)>>1;if(r[s[mid]]<qn)L=mid+1;else R=mid-1;}ans[q[z].id]=sum[L-1]+getf(s[L],l[s[L]],qn)+(qn+1)*qn*i+i*qn+i-1;z++;}}for(ll i=1;i<=t;i++)printf("%lld\n",ans[i]);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318468.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

.net core i上 K8S(五).netcore程序的hostip模式

正文上一章讲了pod的管理&#xff0c;今天再分享一个pod的访问方式1.Pod的HostIP模式Pod的HostIP模式&#xff0c;可以通过宿主机访问pod内的服务&#xff0c;创建yaml文件如下apiVersion: v1 kind: Pod metadata: name: netcore-podlabels:app: netcorepod spec:containers: …

状态压缩:枚举子集(最优组队)(ybtoj)(动态规划)

解析 很裸的状压dp 但是直接暴力的话状态2n,枚举2n 乘在一起会T诶 怎么办呢&#xff1f; 使用下面这个循环&#xff0c;就可以保证只会枚举当前状态s的子集 for(int i(s-1)&s;i;i(i-1)&s){........ }证明 举举例子就挺明显了 为什么不重不漏呢&#xff1f; 首先i肯…

P7295-[USACO21JAN]Paint by Letters P【平面图欧拉公式】

正题 题目链接:https://www.luogu.com.cn/problem/P7295 题目大意 给出n∗mn*mn∗m的网格&#xff0c;每个格子上有字母&#xff0c;相同字母的四联通相邻格子为连通&#xff0c;每次询问一个子矩阵求连通块个数。 1≤n,m,q≤10001\leq n,m,q\leq 10001≤n,m,q≤1000 解题思路…

【活动(深圳)】告别2018之12.22 大湾区.NET Meet 大会 ,同时有网络直播

今年的 Connect(); 主题更加聚焦开发者工具生产力、开源&#xff0c;以及无服务器&#xff08;Serverless&#xff09;云服务。Visual Studio 2019 AI 智能加持的 IntelliCode、实时代码协作共享 Live Share&#xff0c;.NET Core 3.0的预览版本附带了大量更新&#xff0c;旨在…

带花树算法

对于一般的二分图匹配我们肯定会想到匈牙利算法&#xff0c;但是如果图中出现奇环怎么办&#xff1f;此时匈牙利算法就不可以了&#xff0c;就需要另一个算法&#xff1a;带花树算法 主要就是为了解决奇环的问题 我们匹配时会发现&#xff0c;如果存在奇环&#xff0c;传统的匈…

最短路径(状压dp)(ybtoj)

解析 “bug总有de完的一天” 头疼 暴力写的话本题显然复杂度是假的 有一个很好的思路优化时间复杂度 先用dp[k][i]表示**从第k个关键点到任意i点的最短路 跑k遍 SPFA或迪杰斯特拉 即可 然后转移时可以只考虑关键点 使状态转移数大大降低 细节 头真疼 边界条件&#xff1a; …

P3352-[ZJOI2016]线段树【dp】

正题 题目链接:https://www.luogu.com.cn/problem/P3352 题目大意 nnn个数字的一个序列&#xff0c;每次随机选择一个区间让这个区间所有数等于这个区间的最大值&#xff0c;重复qqq次&#xff0c;对每个位置求所有情况下这个位置的值的和。 1≤n,q≤4001\leq n,q\leq 4001≤n…

P2403 [SDOI2010]所驼门王的宝藏(强连通分量)(拓扑排序)

文章目录题目描述解析代码洛谷传送门题目描述 解析 看题目要求很容易想到强连通分量缩点加拓扑dp 但是问题在于存图 第一感就是和暴力和每个点连边 但那样无论点数和边数都很爆炸 随后我们发现这个图非常稀疏 所以我们可以只连有宝藏的点 然而这样边数会被一行横门这样的数据…

Xamarin.Forms之UserDialogs 重制版本

在 forms 里面&#xff0c;目前使用比较多的弹出组件是 Acr.UserDialogs &#xff0c;但是这个组件有些小问题&#xff0c;比如 loading .hide 会同时把 toast 给一起关掉&#xff0c;android 下的 toast 希望是 安卓原生的toast 样子&#xff0c;而不是 底部弹出一个横条&…

H.Minimum-cost Flow

H.Minimum-cost Flow 题目&#xff1a; 其实就是给出每条边的单位费用&#xff0c;q次查询&#xff0c;每次查询改变所有边的容量&#xff08;所有边容量一样&#xff09;&#xff0c;问最后流出1流量的最小花费是多少&#xff1f; 题解&#xff1a; 暴力做法肯定是每次询问…

P6177-Count on a tree II/[模板]树分块

正题 题目链接:https://www.luogu.com.cn/problem/P6177 题目大意 nnn个点的一棵树mmm次询问树上颜色。 强制在线 1≤n≤4105,1≤m≤105,0≤vali<2311\leq n\leq 4\times 10^5,1\leq m\leq 10^5,0\leq val_i<2^{31}1≤n≤4105,1≤m≤105,0≤vali​<231 解题思路 把所…

eShopOnContainers 看微服务③:Identity Service

引言通常&#xff0c;服务所公开的资源和 API 必须仅限受信任的特定用户和客户端访问。那进行 API 级别信任决策的第一步就是身份认证——确定用户身份是否可靠。在微服务场景中&#xff0c;身份认证通常统一处理。一般有两种实现形式&#xff1a;基于API 网关中心化认证&#…

P2168 [NOI2015] 荷马史诗(哈夫曼编码树)

传送门 文章目录题目描述解析细节代码题目描述 解析 其实就是构造一棵树 另所有带权点都为叶子节点 其代价为权值与深度的乘积 求最小代价及最小代价下的最小深度 可以看成一开始有n棵小树 每次把k棵合并在一起 最后合成一棵大树就行了 每次合并的代价是k棵树的权值和 看到这…

B-Suffix Array

B-Suffix Array 题意&#xff1a; 一个字符串只含有a和b&#xff0c;先给出b数组的构造方式&#xff1a; 对于每个位置i来说&#xff1a; 如果存在一个位置j&#xff0c;使得j<i,且s[j] s[i],则b[i]i-j否则b[i]0 现在对字符串每个后缀都构造B数组&#xff0c;并按照字典…

Wannafly挑战赛10F-小H和遗迹【Trie,树状数组】

正题 题目链接:https://ac.nowcoder.com/acm/contest/72/F 题目大意 nnn个字符串&#xff0c;包括小写字母和#\##。其中#\##可以替换为任意字符串。求有多少对字符串可能相同。 保证每个字符串至少有一个#\##。 2≤n≤500000,∑i1n∣si∣≤1062\leq n\leq 500000,\sum_{i1}^n…

线段树合并、分裂

基本概念&#xff1a; 如果需要维护许多个大小为 \(10^5\) 级别的多重集&#xff0c;可以看做给每一个多重集建立一棵线段树。线段树的合并、分裂就是多重集的累加、分开。 这里使用动态开点的方式存储线段树树。 如果一个节点为空&#xff0c;那么它的编号为 \(0\) 。 变量释义…

(神奇的)虚树(初步了解)

参考文献&#xff1a; 博文1 博文2 博文3 引入 在一类树上动态规划问题中,题目给出的询问往往包含树上的很多各节点,并保证总的点数规模小于某个值. 如果我们直接在整颗树上进行dp的话,时间复杂度与询问的次数有关,这显然是不可接受的,如果我们可以找到一种动态规划的方法,使…

P2048 [NOI2010] 超级钢琴(RMQ 贪心)

文章目录题目描述解析代码传送门题目描述 解析 首先&#xff0c;如果只有一个和弦&#xff0c;那么问题显然简单了 用前缀和结合ST表随便做做即可 然而 这次要求前k大的 怎么办呢&#xff1f; 参照之前有一道序列合并的做法 我们想到&#xff0c;可以先建一个优先队列&#xf…

微服务架构基础之Service Mesh

ServiceMesh(服务网格) 概念在社区里头非常火&#xff0c;有人提出 2018 年是 ServiceMesh 年&#xff0c;还有人提出 ServiceMesh 是下一代的微服务架构基础。那么到底什么是 ServiceMesh&#xff1f;它的诞生是为了解决什么问题&#xff1f;企业是否适合引入 ServiceMesh&…

[SDOI2011]消耗战

[SDOI2011]消耗战 题意&#xff1a; 给出n个点的一棵带有边权的树,以及q个询问.每次询问给出k个点,询问这使得这k个点与1点不连通所需切断的边的边权和最小是多少. 题解&#xff1a; 树型dp虚树 dp[x]:切断x及其子树上询问点的最小代价 预处理出minv[pos]代表从11到pos路径…