设$c[i]=g[i]+\frac{i(i+1)}{2}-a[i]\times i-a[i]$,$d[i]=a[i]-i$

$f[i]$表示以$i$为结尾最多保留多少个建筑,则

$f[i]=\max(f[j])+1$,$j<i且d[j]\leq d[i]$

$g[i]$表示以$i$为结尾的最小代价,则

$g[i]=\min(i\times d[j]+c[j])+b[i]+a[i]+\frac{i(i-1)}{2}$,$j<i,f[j]+1=f[i]且d[j]\leq d[i]$

按$f$一组一组来处理,每组按序号排序,然后分治处理。

分治的时候,按$d$排序,用栈维护凸壳,查询时在凸壳上二分。

时间复杂度$O(n\log^2n)$。

#include<cstdio>
#include<algorithm>
#define N 100010
using namespace std;
typedef long long ll;
int n,m,i,j,k,a[N],b[N],d[N],e[N],bit[N],f[N],ans0,G[N],nxt[N];
int cl,cr,L[N],R[N],top,q[N];
ll c[N],g[N],ans1=1LL<<62;
struct E{int x,t;E(){}E(int _x,int _t){x=_x,t=_t;}}h[N];
inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
inline void umax(int&a,int b){if(a<b)a=b;}
inline void umin(ll&a,ll b){if(a>b)a=b;}
inline int lower(int x){
int l=1,r=n+1,mid,t;
while(l<=r)if(e[mid=(l+r)>>1]<=x)l=(t=mid)+1;else r=mid-1;
return t;
}
inline void ins(int x,int y){for(;x<=n+1;x+=x&-x)umax(bit[x],y);}
inline int ask(int x){int t=-N;for(;x;x-=x&-x)umax(t,bit[x]);return t;}
inline void add(int x,int y){nxt[y]=G[x];G[x]=y;}
inline bool cmp(int x,int y){return d[x]==d[y]?c[x]>c[y]:d[x]<d[y];}
inline double pos(int x,int y){return 1.0*(c[y]-c[x])/(d[x]-d[y]);}
inline void insert(int x){
if(top&&d[x]==d[q[top]])top--;
while(top>1&&pos(q[top],x)>pos(q[top],q[top-1]))top--;
q[++top]=x;
}
inline void query(int x){
if(!top)return;
int l=1,r=top-1,mid,t=top;
while(l<=r){
mid=(l+r)>>1;
if(x>pos(q[mid],q[mid+1]))r=(t=mid)-1;else l=mid+1;
}
umin(g[x],1LL*x*d[q[t]]+c[q[t]]);
}
void solve(int l,int r){
if(l==r)return;
int mid=(l+r)>>1,i,j;
solve(l,mid),solve(mid+1,r);
for(cl=0,i=l;i<=mid;i++)if(!h[i].t)L[++cl]=h[i].x;
for(cr=0,i=r;i>mid;i--)if(h[i].t)R[++cr]=h[i].x;
if(!cl||!cr)return;
sort(L+1,L+cl+1,cmp),sort(R+1,R+cr+1,cmp);
for(i=j=1,top=0;i<=cr;i++){
while(j<=cl&&d[L[j]]<=d[R[i]])insert(L[j++]);
query(R[i]);
}
}
int main(){
read(n);
for(i=1;i<=n;i++)read(a[i]);
for(i=1;i<=n;i++)read(b[i]),e[i+1]=d[i]=a[i]-i;
sort(e+1,e+n+2);
for(i=1;i<=n+1;i++)bit[i]=-N;
ins(lower(0),0);
for(i=1;i<=n;i++)g[i]=ans1,j=lower(d[i]),ins(j,f[i]=ask(j)+1);
ans0=ask(n+1);
for(i=0;i<=n;i++)G[i]=-1;
for(i=n;~i;i--)if(f[i]>=0)add(f[i],i);
for(i=0;i<ans0;i++){
j=G[i],k=G[i+1],m=0;
while(~j||~k){
if(j<0)h[++m]=E(k,1),k=nxt[k];
else if(k<0||j<k)h[++m]=E(j,0),j=nxt[j];
else h[++m]=E(k,1),k=nxt[k];
}
solve(1,m);
for(k=G[i+1];~k;k=nxt[k]){
g[k]+=b[k]+a[k]+1LL*k*(k-1)/2;
c[k]=g[k]+1LL*k*(k+1)/2-1LL*a[k]*k-a[k];
}
}
for(i=0;i<=n;i++)if(f[i]==ans0)umin(ans1,g[i]+1LL*(n-i)*a[i]+1LL*(n-i)*(n-i+1)/2);
return printf("%d %lld",ans0,ans1),0;
}

  

BZOJ2149 : 拆迁队的更多相关文章

  1. bzoj2149拆迁队 斜率优化dp+分治

    2149: 拆迁队 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 397  Solved: 177[Submit][Status][Discuss] ...

  2. Bzoj2149拆迁队:cdq分治 凸包

    国际惯例的题面:我们考虑大力DP.首先重新定义代价为:1e13*选择数量-(总高度+总补偿).这样我们只需要一个long long就能维护.然后重新定义高度为heighti - i,这样我们能选择高度 ...

  3. 【BZOJ2149】拆迁队(斜率优化DP+CDQ分治)

    题目: 一个斜率优化+CDQ好题 BZOJ2149 分析: 先吐槽一下题意:保留房子反而要给赔偿金是什么鬼哦-- 第一问是一个经典问题.直接求原序列的最长上升子序列是错误的.比如\(\{1,2,2,3 ...

  4. ●BZOJ 2149 拆迁队

    题链: http://www.lydsy.com/JudgeOnline/problem.php?id=2149 题解: 斜率优化DP,栈维护凸包,LIS,分治(我也不晓得是不是CDQ分治...) 一 ...

  5. 核心业务系统数据库平台迁移: Oracle -> MySQL

    为了对核心技术拥有更多的自主控制能力,为了解决数据库的线性扩展问题,为了尽量减少对商业软件的依赖,为了摆脱对高端硬件的依赖,为了… 基于以上多种原因,2年前,我们计划将公司某核心应用平台进行大手术:数 ...

  6. ZJOI2019一轮停课刷题记录

    Preface 菜鸡HL终于狗来了他的省选停课,这次的时间很长,暂定停到一试结束,不过有机会二试的话还是可以搞到4月了 这段时间的学习就变得量大而且杂了,一般以刷薄弱的知识点和补一些新的奇怪技巧为主. ...

  7. RuPengWang项目

    项目 Day1------------------------- 说明:建外键约束.ashx+Razor RupengWang创建三个类库Model DAL BLL后台:RupengWang.Admi ...

  8. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  9. 玩家福音:10款最佳Linux免费游戏

    “我能在Linux平台上游戏吗?”这类疑问正困扰游戏玩家,那么答案就是“快去Linux平台吧!”.开源组织一直以来坚持不懈为Linux操作系统开发不同类型的游戏,在Linux平台下的游戏完全不亚于其他 ...

随机推荐

  1. hdu1798(几何面积计算)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1798 题意:给出两个圆的圆心坐标与半径,求他们相交部分的大小 思路:有三种情况: 1. 两圆相离,ar ...

  2. NYOJ926(概率)

    题目链接:http://acm.nyist.net/JudgeOnline/problem.php?pid=926 设最终A获胜的概率为P,则B获胜的概率为1-P: 因此我们只需要考虑A获胜的概率即可 ...

  3. 3.2 STL中的函数对象类模板

    *: STL中有一些函数对象类模板,如下所示: 1)例如要求两个double类型的x 和y 的积,可以: multiplies<double>()(x,y); 该表达式的值就是x*y的值. ...

  4. GoLang搞一个基本的HTTP服务

    慢慢和python的对应一下看看. package main import ( "fmt" "net/http" "strings" &qu ...

  5. jQuery函数attr()和prop()的区别

    在jQuery中,attr()函数和prop()函数都用于设置或获取指定的属性,它们的参数和用法也几乎完全相同. 但不得不说的是,这两个函数的用处却并不相同.下面我们来详细介绍这两个函数之间的区别. ...

  6. [LeetCode] Isomorphic Strings

    Isomorphic Strings Total Accepted: 30898 Total Submissions: 120944 Difficulty: Easy Given two string ...

  7. PHP不同域名cookie共享(单点登录实现原理)

    PHP使用P3P完成COOKIE跨域操作实际实用中,类似的需求有,比如说我们有两个域名,我们想实现在一个域名登录后,能自动完成另一个域名的登录,也就是单点登录(SSO)功能.为了测试的方便,先编辑ho ...

  8. SQL常用查询语句及函数

    1.日期匹配_获取时间差 select datediff(dd,getdate(),'12/25/2006')  --计算从今天到12/25/2006还有多少个月 2.不能通过IP连接数据库 在数据库 ...

  9. 自己动手写RTP服务器——关于RTP协议

    转自:http://blog.csdn.net/baby313/article/details/7353605 本文会带领着你一步步动手实现一个简单的RTP传输服务器,旨在了解RTP流媒体传输协议以及 ...

  10. Effective C++ 之 Item 5:了解C++默默编写并调用哪些函数

    Effective C++ chapter 2. 构造 / 析构 / 赋值运算 (Constructors, Destructors, and Assignment Operators) Item 5 ...