51nod 1488 帕斯卡小三角_51nod1488-程序员宅基地

技术标签: 单调栈  二分  单调队列优化DP  

1488 帕斯卡小三角

· f(1,j)=a[j], 1≤j≤n.
· f(i,j)=min(f(i-1,j),f(i-1,j-1))+a[j], 2≤i≤n, i≤j≤n.
a是一个长度为n的数组。现在有若干个询问,输入x,y,求f(x,y)

Input
单组测试数据。
第一行有一个整数n (1 ≤ n ≤ 10^5),表示a的长度。
第二行有n个整数a[1], a[2],…,a[n](0 ≤ a[i] ≤ 10^4)。
接下来一行有一个整数m(1 ≤ m ≤ 10^5)表示询问次数。
接下来m行每行有两个整数xi, yi(1 ≤ xi ≤ yi ≤ n)。
Output
输出m行,对应f(xi,yi)。
Input示例
6
2 2 3 4 3 4
4
4 5
3 4
3 4
2 3
Output示例
12
9
9
5


【分析】

一道很神奇的题目。

首先我们需要开启脑洞,观察到对于一个 (x,y) 的答案,如果把它的转移过程在三角形上的路径标出来,那么这条路径一定是由第一行某个点 (1,s) ,竖直向下转移到 (yx+1,s) ,然后再沿一条斜率为-1的直线转移到 (x,y)

然后得到朴素算法:用s数组记录a数组的前缀和,扫一遍y坐标左边的点 (1,s)(s<=y) ,不断更新答案。

观察到如果 (1,k) 位置比 (1,j) 位置答案更优,且 k<j 。那么可以列出一个奇葩的不等式。

化简得到 (s[k]a[k]k)(s[j]a[j]j)a[k]a[j]<xy

这是一个可以斜率优化的东西。

所以我们可以维护一个上凸壳。将询问的以 y <script type="math/tex" id="MathJax-Element-10">y</script>为关键字排序,每遇到一组询问,在凸壳中二分出合法的左端点,然后再在合法区间内二分出最优答案。

早上和中午为维护上凸壳还是下凸壳的问题纠结了很久2333…笨


【代码】

//51nod 1488 帕斯卡小三角
#include<bits/stdc++.h>
#define D double
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mxn=100005;
int n,m,T,top;
struct query {
   int x,y,id;} q[mxn];
int a[mxn],s[mxn],sta[mxn],ans[mxn];
inline bool comp(query a,query b)
{
    return a.y<b.y;
}
inline int read()
{
    int x=0;char ch=getchar();
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch>='0' && ch<='9') x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
    return x;
}
inline D g(int k,int j)
{
    return ((D)(s[k]-a[k]*k)-(D)(s[j]-a[j]*j))/(D)(a[k]-a[j]);
}
inline int find(int x)
{
    int l=1,r=top;
    while(l<r)
    {
        int mid=l+r>>1;
        if(sta[mid]<x) l=mid+1;
        else r=mid;
    }
    return l;
}
int main()
{
    int i,j;
    n=read();
    fo(i,1,n) a[i]=read();
    fo(i,1,n) s[i]=s[i-1]+a[i];
    m=read();
    fo(i,1,m)
    {
        q[i].x=read();
        q[i].y=read();
        q[i].id=i;
    }
    sort(q+1,q+m+1,comp);
    for(i=1,j=1;i<=n;i++)
    {
        while(top && a[sta[top]]>=a[i]) top--;
        while(top>1 && g(sta[top],i)>=g(sta[top-1],i)) top--;
        sta[++top]=i;
        while(q[j].y==i && j<=m)
        {
            int l=find(q[j].y-q[j].x),r=top;
            while(l<r)
            {
                int mid=l+r>>1;
                if(g(sta[mid],sta[mid+1])<q[j].x-q[j].y) r=mid;
                else l=mid+1;
            }
            l=sta[l],r=q[j].y;
            ans[q[j].id]=s[r]-s[l]+(q[j].x-(r-l))*a[l];
            j++;
        }
    }
    fo(i,1,m) printf("%d\n",ans[i]);
    return 0;
}
版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/chai_jing/article/details/73648205

智能推荐

oracle 12c 集群安装后的检查_12c查看crs状态-程序员宅基地

文章浏览阅读1.6k次。安装配置gi、安装数据库软件、dbca建库见下:http://blog.csdn.net/kadwf123/article/details/784299611、检查集群节点及状态:[root@rac2 ~]# olsnodes -srac1 Activerac2 Activerac3 Activerac4 Active[root@rac2 ~]_12c查看crs状态

解决jupyter notebook无法找到虚拟环境的问题_jupyter没有pytorch环境-程序员宅基地

文章浏览阅读1.3w次,点赞45次,收藏99次。我个人用的是anaconda3的一个python集成环境,自带jupyter notebook,但在我打开jupyter notebook界面后,却找不到对应的虚拟环境,原来是jupyter notebook只是通用于下载anaconda时自带的环境,其他环境要想使用必须手动下载一些库:1.首先进入到自己创建的虚拟环境(pytorch是虚拟环境的名字)activate pytorch2.在该环境下下载这个库conda install ipykernelconda install nb__jupyter没有pytorch环境

国内安装scoop的保姆教程_scoop-cn-程序员宅基地

文章浏览阅读5.2k次,点赞19次,收藏28次。选择scoop纯属意外,也是无奈,因为电脑用户被锁了管理员权限,所有exe安装程序都无法安装,只可以用绿色软件,最后被我发现scoop,省去了到处下载XXX绿色版的烦恼,当然scoop里需要管理员权限的软件也跟我无缘了(譬如everything)。推荐添加dorado这个bucket镜像,里面很多中文软件,但是部分国外的软件下载地址在github,可能无法下载。以上两个是官方bucket的国内镜像,所有软件建议优先从这里下载。上面可以看到很多bucket以及软件数。如果官网登陆不了可以试一下以下方式。_scoop-cn

Element ui colorpicker在Vue中的使用_vue el-color-picker-程序员宅基地

文章浏览阅读4.5k次,点赞2次,收藏3次。首先要有一个color-picker组件 <el-color-picker v-model="headcolor"></el-color-picker>在data里面data() { return {headcolor: ’ #278add ’ //这里可以选择一个默认的颜色} }然后在你想要改变颜色的地方用v-bind绑定就好了,例如:这里的:sty..._vue el-color-picker

迅为iTOP-4412精英版之烧写内核移植后的镜像_exynos 4412 刷机-程序员宅基地

文章浏览阅读640次。基于芯片日益增长的问题,所以内核开发者们引入了新的方法,就是在内核中只保留函数,而数据则不包含,由用户(应用程序员)自己把数据按照规定的格式编写,并放在约定的地方,为了不占用过多的内存,还要求数据以根精简的方式编写。boot启动时,传参给内核,告诉内核设备树文件和kernel的位置,内核启动时根据地址去找到设备树文件,再利用专用的编译器去反编译dtb文件,将dtb还原成数据结构,以供驱动的函数去调用。firmware是三星的一个固件的设备信息,因为找不到固件,所以内核启动不成功。_exynos 4412 刷机

Linux系统配置jdk_linux配置jdk-程序员宅基地

文章浏览阅读2w次,点赞24次,收藏42次。Linux系统配置jdkLinux学习教程,Linux入门教程(超详细)_linux配置jdk

随便推点

matlab(4):特殊符号的输入_matlab微米怎么输入-程序员宅基地

文章浏览阅读3.3k次,点赞5次,收藏19次。xlabel('\delta');ylabel('AUC');具体符号的对照表参照下图:_matlab微米怎么输入

C语言程序设计-文件(打开与关闭、顺序、二进制读写)-程序员宅基地

文章浏览阅读119次。顺序读写指的是按照文件中数据的顺序进行读取或写入。对于文本文件,可以使用fgets、fputs、fscanf、fprintf等函数进行顺序读写。在C语言中,对文件的操作通常涉及文件的打开、读写以及关闭。文件的打开使用fopen函数,而关闭则使用fclose函数。在C语言中,可以使用fread和fwrite函数进行二进制读写。‍ Biaoge 于2024-03-09 23:51发布 阅读量:7 ️文章类型:【 C语言程序设计 】在C语言中,用于打开文件的函数是____,用于关闭文件的函数是____。

Touchdesigner自学笔记之三_touchdesigner怎么让一个模型跟着鼠标移动-程序员宅基地

文章浏览阅读3.4k次,点赞2次,收藏13次。跟随鼠标移动的粒子以grid(SOP)为partical(SOP)的资源模板,调整后连接【Geo组合+point spirit(MAT)】,在连接【feedback组合】适当调整。影响粒子动态的节点【metaball(SOP)+force(SOP)】添加mouse in(CHOP)鼠标位置到metaball的坐标,实现鼠标影响。..._touchdesigner怎么让一个模型跟着鼠标移动

【附源码】基于java的校园停车场管理系统的设计与实现61m0e9计算机毕设SSM_基于java技术的停车场管理系统实现与设计-程序员宅基地

文章浏览阅读178次。项目运行环境配置:Jdk1.8 + Tomcat7.0 + Mysql + HBuilderX(Webstorm也行)+ Eclispe(IntelliJ IDEA,Eclispe,MyEclispe,Sts都支持)。项目技术:Springboot + mybatis + Maven +mysql5.7或8.0+html+css+js等等组成,B/S模式 + Maven管理等等。环境需要1.运行环境:最好是java jdk 1.8,我们在这个平台上运行的。其他版本理论上也可以。_基于java技术的停车场管理系统实现与设计

Android系统播放器MediaPlayer源码分析_android多媒体播放源码分析 时序图-程序员宅基地

文章浏览阅读3.5k次。前言对于MediaPlayer播放器的源码分析内容相对来说比较多,会从Java-&amp;amp;gt;Jni-&amp;amp;gt;C/C++慢慢分析,后面会慢慢更新。另外,博客只作为自己学习记录的一种方式,对于其他的不过多的评论。MediaPlayerDemopublic class MainActivity extends AppCompatActivity implements SurfaceHolder.Cal..._android多媒体播放源码分析 时序图

java 数据结构与算法 ——快速排序法-程序员宅基地

文章浏览阅读2.4k次,点赞41次,收藏13次。java 数据结构与算法 ——快速排序法_快速排序法