P2605 [ZJOI2010]基站选址

P2605 [ZJOI2010]基站选址

 

题目描述

详见:P2605 [ZJOI2010]基站选址

Solution

首先不难想到一个  O(n^2k)  的DP。

f[i][j]  表示前i个村庄选择了j个基站的总费用。

f[i][j]=min(f[i][k]+cost(k,i))+c[i]

考虑如何优化这个转移。

对于村庄i,我们记录它覆盖范围内最靠前的村庄  fi[i]  和最靠后的村庄  la[i] 。

倘若在la[i]之后的村庄建立基站,则若上一个基站在fi[i]之前,那么就会多一个w[i]的费用。

增加费用的状态是形如f[1][j-1]...f[\;\;fi[i]\;\;][j-1] 的连续一段状态。

我们只需要让f滚动,维护区间加,询问区间最小值,线段树维护即可。

 

Ps:luogu题解中很多代码求  la[i]   时使用的是lower\_\;bound 再特判减1,但这样如果出现相等的值时,会找到相等的值中第一个值,与  la[i]  的实际意义不符,但由于我太蒟,不知道会不会对答案产生影响qwq。

Code

#include<bits/stdc++.h>
using namespace std;
const int MAXN=1e5+50;
const int INF=0x3f3f3f3f;
int d[MAXN],a[MAXN],s[MAXN],w[MAXN],fi[MAXN],la[MAXN],f[MAXN];
vector<int> e[MAXN];
inline int read()
{int x=0,f=1; char c=getchar();while (c<'0'||c>'9') { if(c=='-') f=-1; c=getchar(); }while (c>='0'&&c<='9') { x=(x<<3)+(x<<1)+(c^48); c=getchar(); }return x*f;
}
struct Segment_Tree
{struct segnode{int l,r,tag,mn; } tree[MAXN<<2];void up(int x){ tree[x].mn=min(tree[x<<1].mn,tree[x<<1|1].mn); }void down(int x){if (tree[x].tag!=0){int q=tree[x].tag;tree[x<<1].tag+=q,tree[x<<1].mn+=q;tree[x<<1|1].tag+=q,tree[x<<1|1].mn+=q;tree[x].tag=0;}}void build(int x,int l,int r){tree[x].tag=tree[x].mn=0;if ((tree[x].l=l)==(tree[x].r=r)) {tree[x].mn=f[l];return;}int mid=(l+r)>>1;build(x<<1,l,mid);build(x<<1|1,mid+1,r);up(x);}int query(int x,int l,int r){if (tree[x].l>=l&&tree[x].r<=r)  return tree[x].mn;down(x);int mid=(tree[x].l+tree[x].r)>>1;if (r<=mid) return query(x<<1,l,r);else if (l>mid) return query(x<<1|1,l,r);else return min(query(x<<1,l,mid),query(x<<1|1,mid+1,r));}void change(int x,int l,int r,int y){if (tree[x].l>=l&&tree[x].r<=r){tree[x].mn+=y;tree[x].tag+=y;return;}down(x);int mid=(tree[x].l+tree[x].r)>>1;if (r<=mid) change(x<<1,l,r,y);else if (l>mid) change(x<<1|1,l,r,y);else change(x<<1,l,mid,y),change(x<<1|1,mid+1,r,y);up(x);}	
} segment;
int main()
{int n=read(),m=read(); for (int i=2;i<=n;i++) d[i]=read();for (int i=1;i<=n;i++) a[i]=read();for (int i=1;i<=n;i++) s[i]=read();for (int i=1;i<=n;i++) w[i]=read();n++; d[n]=w[n]=INF;for (int i=1;i<=n;i++){fi[i]=lower_bound(d+1,d+n+1,d[i]-s[i])-d;la[i]=upper_bound(d+1,d+n+1,d[i]+s[i])-d-1;//cout<<i<<" "<<fi[i]<<" "<<la[i]<<endl;e[la[i]].push_back(i);}int cnt=0,ans=0;for (int i=1;i<=n;i++){f[i]=cnt+a[i];for (int k=0;k<e[i].size();k++) cnt+=w[e[i][k]];}ans=f[n];//for (int i=1;i<=n;i++) cout<<i<<" "<<f[i]<<endl;cout<<endl;for (int j=2;j<=m+1;j++){segment.build(1,1,n);for (int i=1;i<=n;i++){if (j>i) f[i]=a[i];else f[i]=segment.query(1,j-1,i-1)+a[i];for (int k=0;k<e[i].size();k++) if (fi[e[i][k]]>1) segment.change(1,1,fi[e[i][k]]-1,w[e[i][k]]);}ans=min(ans,f[n]);}printf("%d\n",ans);return 0;
}

 

 

 

 

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

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

相关文章

CF980D Perfect Groups

CF980D Perfect Groups 题意&#xff1a; 将一个串划分为多个子集&#xff08;不要求连续&#xff09;&#xff0c;要求同一子集内两任意元素的积为平方数 定义一个串的答案为所需的最少子集个数 一个长度为 n 的串有 n(n1)2\frac{n(n1)}{2}2n(n1)​个非空子串&#xff0c;…

P1600 天天爱跑步

P1600 天天爱跑步 题目描述 详见&#xff1a;P1600 天天爱跑步 Solution 树上差分LCA。 Code #include<bits/stdc.h> using namespace std; const int MAXN4e550; vector<int> e[MAXN]; struct wnode{int opt,x,c; }; vector<wnode> w[MAXN]; int dep[…

求斐波那契数列第n位的几种实现方式及性能对比

在每一种编程语言里&#xff0c;斐波那契数列的计算方式都是一个经典的话题。它可能有很多种计算方式&#xff0c;例如&#xff1a;递归、迭代、数学公式。哪种算法最容易理解&#xff0c;哪种算法是性能最好的呢&#xff1f;这里给大家分享一下我对它的研究和总结&#xff1a;…

Codeforces Round #729 (Div. 2)

Codeforces Round #729 (Div. 2) 题号题目知识点AOdd SetBPlus and MultiplyCStrange FunctionDPriority QueueE1Abnormal Permutation Pairs (easy version)E2Abnormal Permutation Pairs (hard version)

P2305 [NOI2014]购票

P2305 [NOI2014]购票 题目描述 详见&#xff1a;P2305 [NOI2014]购票 Solution 写出一个朴素的DP可以看出显然是一个斜率优化&#xff0c;且是在树上求解答案。 因此用点分治维护树上斜率优化DP。 设现在的重心为&#xff0c;根为&#xff0c;我们需要先递归处理所在的连…

ASP.NET Core 3.0 自动挡换手动挡:在 Middleware 中执行 Controller Action

最近由于发现奇怪的 System.Data.SqlClient 性能问题&#xff08;详见之前的博文&#xff09;&#xff0c;被迫提前了向 .NET Core 3.0 的升级工作&#xff08;3.0 Preview 5 中问题已被修复&#xff09;。郁闷的是&#xff0c;在刚开始对部分项目进行升级的时候就遇到了一个障…

CCPC秦皇岛gym102361A. Angle Beats

CCPC秦皇岛gym102361A. Angle Beats 题意&#xff1a; 给你n个点的坐标&#xff0c;现在有q次询问&#xff0c;每次询问给你一个坐标&#xff0c;问这个坐标可以与给定的n个点组成多少个不同的直角三角形 n<2000&#xff0c;q<2000 题解&#xff1a; 正解貌似是用极角…

P5327 [ZJOI2019]语言

P5327 [ZJOI2019]语言 题目描述 详见&#xff1a;P5327 [ZJOI2019]语言 简要题意&#xff1a;给定一棵树和一些链&#xff0c;问树上处于同一条链的不同点对数。 Solution 对于每一个点&#xff0c;考虑以它为端点的可行路径有哪些。 我们可以发现&#xff0c;可以到达的…

EasyOffice-.NetCore一行代码导入导出Excel,生成Word

Excel和Word操作在开发过程中经常需要使用&#xff0c;这类工作不涉及到核心业务&#xff0c;但又往往不可缺少。以往的开发方式在业务代码中直接引入NPOI、Aspose或者其他第三方库&#xff0c;工作繁琐&#xff0c;耗时多&#xff0c;扩展性差——比如基础库由NPOI修改为EPPlu…

CF1572B. Xor of 3

CF1572B. Xor of 3 题意&#xff1a; 给你个01序列&#xff0c;你有一种操作&#xff1a;每次选位置x&#xff0c;然后位置x&#xff0c;x1&#xff0c;x2的值变为三者的异或值。 现在要让所有的数都等于0&#xff0c;请输出存在的合法操作序列 题解&#xff1a; 首先如果有…

P4103 [HEOI2014]大工程

P4103 [HEOI2014]大工程 题目描述 详见&#xff1a;P4103 [HEOI2014]大工程 Solution 显然是虚树的板子题啊&#xff08;我也不造我为啥调了1h&#xff09; 直接建虚树&#xff0c;DP。 两两路径和很好求&#xff0c;只要考虑每条边的贡献&#xff0c;设为询问的关键点个…

【中部武汉】理想离家并不遥远

▼我叫Mike, 是一名“高端”程序猿。毕业于武汉某男生比例占据 70%的高校。我的家乡是湖北襄阳&#xff0c;一个历史悠久的文化名城。7年前&#xff0c;我还是一个刚毕业的毛头小子&#xff0c;在学校也算是半个风云人物&#xff0c;年轻气盛&#xff0c;心比天高&#xff0c;放…

P2495 [SDOI2011]消耗战

P2495 [SDOI2011]消耗战 题目描述 详见&#xff1a;P2495 [SDOI2011]消耗战 Solution 此题是虚树的经典版子题吧qwq。 是不是直接贴代码就行了啊&#xff08;反正不是luogu题解&#xff0c;没有人查&#xff09;。 大概就是先建出虚树&#xff08;一般给定一堆关键点的题…

使用Try.NET创建可交互.NET文档

原文地址&#xff1a;Create Interactive .NET Documentation with Try .NET[1] 原文作者&#xff1a;Maria 译文地址&#xff1a;https://www.cnblogs.com/lwqlun/p/10894497.html 译者&#xff1a;Lamond Lu背景当我们编写开发人员使用的文档时&#xff0c;我们需要捕捉他们的…

Link-Cut Tree

Link-Cut Tree 概述. LCT是一种支持动态维护树上路径信息的数据结构&#xff0c;其本质是实链剖分&#xff0c;通过其他数据结构维护实链的信息达到维护路径及一些子树信息的效果&#xff08;通常为splay&#xff09; 刚开始学的时候感觉很难&#xff0c;但做过几道题了解套…

2021-2022ACM赛季小总结

最近一个月比赛月&#xff0c;做题量和频率太高&#xff0c;写博客浪费时间&#xff0c;所以一直没写。。 当前总结&#xff1a; CCPC桂林铜&#xff0c;银牌题线段树差一行调出&#xff0c;可惜又无奈&#xff0c;缺少运气 ICPC沈阳铜&#xff0c;银牌题的tarjan没想出来&…

ASP.NET Core模块化前后端分离快速开发框架介绍

源码地址GitHub&#xff1a;https://github.com/iamoldli/NetModular演示地址地址&#xff1a;http://118.24.75.170:6220/账户&#xff1a;admin密码&#xff1a;admin前端框架演示地址(临时)地址&#xff1a;http://progqx5cu.bkt.clouddn.com/skins/index.html#/账户&#x…

[POJ2888] Magic Bracelet

[POJ2888] Magic Bracelet 题目描述 简要题意&#xff1a;给圆上个点染色&#xff0c;颜色有种&#xff0c;其中对颜色不能相邻&#xff0c;循环同构&#xff0c;多组数据&#xff0c;询问染色方案数。 Solution 大概就是一道挺显然的Burnside题&#xff08;一般染色&#x…

手写编译器

假期没事想搞搞编译器&#xff0c;但学校课程未教过编译原理&#xff0c;自己只好在b站看了stanford的CS143《编译原理》&#xff0c;并结合相关资料&#xff0c;打算寒假期间自己手写一个简易的编译器 ReactjsbootStrap开发自制编译语言Monkey的编译器&#xff1a;创建简易的页…

使用命令行导出 SQL Server 数据层应用程序

点击上方蓝字关注“汪宇杰博客”我们可以使用 SSMS 导出 SQL Server 数据库的数据层应用程序。我在本地机器上使用这种方法已经有好几年了。如果不知道什么是 DAC&#xff0c;您可以参考 Microsoft 文档&#xff1a;https://docs.microsoft.com/en-us/sql/relational-databases…