由于某人找了个单调栈的题解但是没研究透所以让我们来研究。。。。。。。。。。。。

首先先来考虑下面一种情况,假设第\(k\)次切割时,天数为\(d_k\),高度为\(b_k\),第\(k+1\)次切割时,天数为\(d_{k+1}\),高度为\(b_{k+1}\),那么我们定义一个切割速度,令\(v=\frac{b_{k+1}-b_k}{d_{k+1}-d_k}\),这个切割速度有什么含义呢,如果在\(d_k\)天时所有的草都是\(b_k\)高,那么生长速度\(>v\)的草都要割掉,注意这里的\(v\)极有可能是一个浮点数,而在写代码的时候最好还是避免浮点数的出现,所以要怎么转化?这个应该很显然,令\(v=\lfloor \frac{b_{k+1}-b_k}{d_{k+1}-d_k} \rfloor\)即可,因为判断的时候都是找的大于这个\(v\)的,而给出的生长速度都是整数,所以把这个\(v\)下取整后再判断大小也没有任何影响。怎么求出答案呢?很显然,就是\((d_{k+1}-d_k)*\)需要割掉的草的总生长速度\(+t\),\(t\)是什么?分类讨论一下,如果\(b_k>b_{k+1}\),那么除了长出来的草,还需要割掉\(b_k\)比\(b_{k+1}\)高的那一段,也就是\((b_k-b_{k+1})*\)需要割掉的草的总数量,如果\(b_k==b_{k+1}\),\(t=0\),如果如果\(b_k<b_{k+1}\),就不能长出来的全都割掉,因为只割到\(b_{k+1}\),所以需要把之前多算的那些减去,即\(-(b_{k+1}-b_k)*\)需要割掉的草的总数量,综上所述,\(t=(b_k-b_{k+1})*\)需要割掉的草的总数量,于是需要求的就是需要割掉的草的总生长速度和需要割掉的草的总数量,注意到\(m\)最大也就\(10^6\),所以开\(10^6\)个桶,第\(i\)个桶里边存生长速度为\(i\)的草的数量,然后就可以运用一下前缀和的思想求出上边的两个值。

好了基本思想有了,但是注意一个问题,并不是所有的需要切割的草在上一次切割时都恰好为\(b\),即上一次没有被切割但是下一次它长到了可以被切割的高度的草,把这些草。所以上边讨论了五百多字就白扯了吗,显然不是,还是有一部分草是满足上述的切割办法的,也就是两次都被切割的,所以只需要特殊考虑上一次没有被切割但是下一次它长到了可以被切割的高度的草,把这些草处理了就行。通过上一次的切割处理这个很不好处理,但是我们可以发现,如果我们定义第0天时割掉了所有高度大于0的草(对答案没有影响因为初始时都是0),那么假设当前时刻为\(d_k\),这些草在之前的某个时刻\(d_t(d_t \in[0,d_k)\ \ )\)一定会被割掉,于是\(d_t\)到\(d_k\)套用上述式子就行,聪明的你一定会发现这样会重复割草,没错,的确会,为了避免这种情况,假设上次割掉了生长速度超过\(last\)的草,那么这次只需要割掉生长速度超过切割速度且小于等于\(last\)的草即可,但是这样还是会出现重复割草的情况,因为可能前边已经割过的又被割了一次,所以需要记录前边每一次割草割过的最小切割速度,举个例子吧,假设记录的这个为\(to_i\),表示的意思就是在第\(i\)次割草的时候生长速度大于\(to_i\)的已经都被割掉了,所以如果此时\(k-i\)的切割速度大于\(to_i\)就需要直接\(break\)掉,避免了重复割草。

但是这样做的时间复杂度应该是\(O(m^2)\)的,需要优化一下。

注意到草的生长速度是有单调性的,在不去割草的情况下生长速度越大的高度一定越大废话,于是可以维护一个\(to\)值单调上升的单调栈。若栈顶那一次切割能切到的最小速度\(\leq\)当前能切到的最小速度,则前面切不到的依旧切不到,若栈顶那一次切割能切到的最小速度\(>\)当前能切到的最小速度,则先计算比栈顶那一次切割能切到的最小速度大的,并将栈顶\(pop\)掉,重复此过程至栈顶速度\(\leq\)当前能切到的最小速度。

但是这样会不会漏掉什么情况呢?答案是不会,因为每一次的割草都是割的连续的速度,只可能出现重复情况,而重复情况上边已经排除,所以这种做法是可以的。

由于本人蒟蒻,可能有些地方写的不严谨或者不清楚,又或者有的地方写错了。。。反正感觉不对的地方欢迎指出来。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int lqs=5e5+10,m=1e6;
int day,a[lqs<<1],n,stk[lqs],top;
ll s[lqs<<1],d[lqs],b[lqs],to[lqs];
int main(){
scanf("%d%d",&n,&day);
for(int w,i=1;i<=n;i++)
scanf("%d",&w),a[w]++;
for(int i=1;i<=day;i++)
scanf("%lld%lld",&d[i],&b[i]);
for(int i=1;i<=m;i++)
s[i]=s[i-1]+1ll*a[i]*i;
for(int i=1;i<=m;i++)
a[i]+=a[i-1];
top=1;
for(int i=1;i<=day;i++){
ll res=0;
int last=m,x=max(to[stk[top]],min((b[i]-b[stk[top]])/(d[i]-d[stk[top]]),1ll*m));
while(x<last){
res+=(d[i]-d[stk[top]])*(s[last]-s[x])+(b[stk[top]]-b[i])*(a[last]-a[x]);
last=x;
if(last>to[stk[top]])break;
if(top==0)break;top--;
x=max(to[stk[top]],min((b[i]-b[stk[top]])/(d[i]-d[stk[top]]),1ll*m));
}
to[i]=last;
if(to[i]<m)stk[++top]=i;
printf("%lld\n",res);
}
}

[PA2015]Siano 单调栈的更多相关文章

  1. 2018.07.23[PA2015]Siano(线段树)

    [PA2015]Siano 描述 Description 农夫Byteasar买了一片n亩的土地,他要在这上面种草. 他在每一亩土地上都种植了一种独一无二的草,其中,第i亩土地的草每天会长高a[i]厘 ...

  2. BZOJ1012: [JSOI2008]最大数maxnumber [线段树 | 单调栈+二分]

    1012: [JSOI2008]最大数maxnumber Time Limit: 3 Sec  Memory Limit: 162 MBSubmit: 8748  Solved: 3835[Submi ...

  3. BZOJ 4453: cys就是要拿英魂![后缀数组 ST表 单调栈类似物]

    4453: cys就是要拿英魂! Time Limit: 3 Sec  Memory Limit: 128 MBSubmit: 90  Solved: 46[Submit][Status][Discu ...

  4. BZOJ 3238: [Ahoi2013]差异 [后缀数组 单调栈]

    3238: [Ahoi2013]差异 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 2326  Solved: 1054[Submit][Status ...

  5. poj 2559 Largest Rectangle in a Histogram - 单调栈

    Largest Rectangle in a Histogram Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 19782 ...

  6. bzoj1510: [POI2006]Kra-The Disks(单调栈)

    这道题可以O(n)解决,用二分还更慢一点 维护一个单调栈,模拟掉盘子的过程就行了 #include<stdio.h> #include<string.h> #include&l ...

  7. BZOJ1057[ZJOI2007]棋盘制作 [单调栈]

    题目描述 国际象棋是世界上最古老的博弈游戏之一,和中国的围棋.象棋以及日本的将棋同享盛名.据说国际象棋起源于易经的思想,棋盘是一个8*8大小的黑白相间的方阵,对应八八六十四卦,黑白对应阴阳. 而我们的 ...

  8. 洛谷U4859matrix[单调栈]

    题目描述 给一个元素均为正整数的矩阵,上升矩阵的定义为矩阵中每行.每列都是严格递增的. 求给定矩阵中上升子矩阵的数量. 输入输出格式 输入格式: 第一行两个正整数n.m,表示矩阵的行数.列数. 接下来 ...

  9. POJ3250[USACO2006Nov]Bad Hair Day[单调栈]

    Bad Hair Day Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 17774   Accepted: 6000 Des ...

随机推荐

  1. 如何解决ubuntu 12.04重启后出现waiting for network configuration和网络标志消失问题

    如何解决ubuntu 12.04重启后出现waiting for network configuration和网络标志消失问题 作为菜鸟的我在学着设置网络后,重启电脑后显示 waiting forne ...

  2. Firewalld 的基本使用

    RHEL 7 系统中集成了多款防火墙管理工具,其中 firewalld(Dynamic Firewall Manager of Linux systems,Linux 系统的动态防火墙管理器)服务是默 ...

  3. zabbix内存百分比监控告警

    本文结合配置内存不足10%触发报警的需求,zabbix给我们提供的模板,里面都已经配置好了item和trigger.但是给我们的模板是当内存小于20M的时候才会触发报警,这样不能满足我们的需求,我们需 ...

  4. Flink 如何分流数据

    场景 分流方式 如何分流 使用Filter分流 使用Split分流 使用Side Output分流 场景 获取流数据的时候,通常需要根据所需把流拆分出其他多个流,根据不同的流再去作相应的处理. 举个例 ...

  5. 微信小程序实现连续扫码功能(uniapp)

    注:本文使用的是 uniapp 语法. 微信小程序提供了扫码API:wx.scanCode,但它只能扫一次码,想要实现连续扫码,需要借用 camera 组件.camera 组件不仅能拍照,还具有扫码功 ...

  6. 【Vulnhub】FristiLeaks v1.3

    靶机信息 下载连接 https://download.vulnhub.com/fristileaks/FristiLeaks_1.3.ova.torrent https://download.vuln ...

  7. 一起读《Java编程思想》(第四版)

    空余时间看<Java编程思想>(第四版)这本书,遇到不懂的知识点就记录在本博客内. 1.复用具体实现 Java代码复用的三种常见方式:继承.组合.代理. 1.继承:使用extends关键字 ...

  8. 几款一元单片机对比:CMS8S5880、STM8S003、N76E003

    大概17年开始,STM8S003的价格被贸易商炒货,变得很不稳定,一度上涨到2~3元,因为市场需求大增,小家电.无线充和一些简单功能的产品,本人就有在空气净化器.433M触摸开关.数据收发模块.红外控 ...

  9. 11.实战交付一套dubbo微服务到k8s集群(4)之使用Jenkins进行持续构建交付dubo服务的提供者

    1.登录到jenkins,新建一个项目 2.新建流水线 3.设置保留的天数及份数 4. 添加参数 # 参数 . name: git_repo type: string description: 项目在 ...

  10. 黑马程序员spring data jpa 2019年第一版本

    第一步首先创建一个maven工程,导入对于的pom依赖 <project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xs ...