当前位置: 代码迷 >> 综合 >> 51nod 1488 帕斯卡小三角
  详细解决方案

51nod 1488 帕斯卡小三角

热度:100   发布时间:2023-10-29 06:20:01.0

前言

文化课弱鸡选手又来做题啦。。
(⊙o⊙)…尽量保证一周一或两题?(大雾)
可能过多几个星期文化课补上来就可以多做点题了吧。。

题意

· 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)

题解

随意观察一下可以发现,这个东西其实就是你在第一行某一个找一个点x
然后你每一次可以往下走或者往右下走,问你走到这个点的代价最小
然后再胡乱观察一下,可以发现肯定是从起点一直往下走,然后剩下一段往右下走。。
然后问题就变成了要找一个s使得
c是前缀和
c[y]?c[s]+a[s]?(y?x+s)c[y]?c[s]+a[s]?(y?x+s)最小\
然后把式子整理一下
c[y]+a[s]?(y?x)+(a[s]?s?c[s])c[y]+a[s]?(y?x)+(a[s]?s?c[s])
很容易发现,如果一个a比前面的a要小,那么前面肯定是没有用的
这个就可以斜率优化了
维护一个下凸的凸包
每一次询问二分找最值就可以了
其实s的选择是有一个区间限制的
就是[max(0,y?x),y][max(0,y?x),y]
我们按y排序,每一次加入新的就可以了
然后听说不管左边界也可以AC。。然后我就偷工减料了。。
其实要管的话,就把凸包的那个单调栈换成set就可以了
时间复杂度还是一个loglog的,就是常数大一点点
一开始没有读优过不去啊
CODE:

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
using namespace std;
const int N=100005;
int a[N],c[N];
int ans[N];
struct qq {int x,y,id; }s[N];
int n,Q;
bool cmp (qq x,qq y)    {
   return x.y<y.y;}
int q[N];
int st,ed,now;
double Y (int x)    {
   return a[x]*x-c[x];}
double X (int x)    {
   return -a[x];}
double get_k (int x,int y)
{return (double)(Y(x)-Y(y))/(double)(X(x)-X(y));
}
int calc (int xx,int x,int y)//这个点行不行 
{return c[y]-c[xx]+a[xx]*(x-y+xx);
}
void solve ()
{st=1;ed=0;now=1;for (int u=1;u<=n;u++){while (st<=ed&&a[q[ed]]>=a[u]) ed--;while (st<ed&&get_k(q[ed-1],q[ed])<get_k(q[ed],u))  ed--;   q[++ed]=u;while (now<=Q&&s[now].y==u)//可以用到这个了 {int l=st,r=ed,lalal;while (l<=r){int mid=(l+r)>>1;if (mid==ed||calc(q[mid],s[now].x,s[now].y)<calc(q[mid+1],s[now].x,s[now].y)) {lalal=mid;r=mid-1;}else l=mid+1;}ans[s[now].id]=calc(q[lalal],s[now].x,s[now].y);now++;}}
}
int read()
{char ch=getchar();int x=0;while (ch<'0'||ch>'9') ch=getchar();while (ch>='0'&&ch<='9') {
   x=x*10+ch-'0';ch=getchar();}return x;
}
int main()
{n=read();for (int u=1;u<=n;u++)  {a[u]=read();c[u]=c[u-1]+a[u];}Q=read();for (int u=1;u<=Q;u++){s[u].x=read();s[u].y=read();s[u].id=u;}sort(s+1,s+1+Q,cmp);solve();for (int u=1;u<=Q;u++)  printf("%d\n",ans[u]);return 0;
}