解法
首先可以观察到所有矩形排成了类似楼梯的形状,所以如果对x排序,y也是有序的。
先按x从小到大排序
然后考虑dp:
f
[
i
]
表
示
最
后
一
个
矩
形
选
择
的
是
第
i
个
的
最
大
价
值
f[i]表示最后一个矩形选择的是第i个的最大价值
f[i]表示最后一个矩形选择的是第i个的最大价值
转移:
f
[
i
]
=
m
a
x
j
=
1
i
−
1
(
f
[
i
]
,
f
[
j
]
+
y
[
i
]
∗
(
x
[
i
]
−
x
[
j
]
)
−
a
[
i
]
)
f[i]=max_{j=1}^{i-1}(f[i],f[j]+y[i]*(x[i]-x[j])-a[i])
f[i]=maxj=1i−1(f[i],f[j]+y[i]∗(x[i]−x[j])−a[i])
考虑优化这个dp:
f
[
j
]
+
y
[
i
]
∗
x
[
i
]
−
y
[
i
]
∗
x
[
j
]
−
a
[
i
]
f[j]+y[i]*x[i]-y[i]*x[j]-a[i]
f[j]+y[i]∗x[i]−y[i]∗x[j]−a[i]这个转移式中:只有
y
[
i
]
∗
x
[
j
]
y[i]*x[j]
y[i]∗x[j]这一项与i,j都有关,然后又因为y[i]具有单调性,所以可以考虑斜率优化:
我们观察和j有关的项,可以认为
x
[
i
]
x[i]
x[i]是横坐标,
f
[
i
]
f[i]
f[i]是纵坐标,这样每个决策可以表示成(x[j],f[j])的点,然后我们的dp转移变成用斜率为y[i]的直线切之前的决策形成的凸包。具体的式子就是
f
[
j
]
=
x
[
j
]
∗
y
[
i
]
+
d
p
[
i
]
f[j]=x[j]*y[i]+dp[i]
f[j]=x[j]∗y[i]+dp[i]
f
[
i
]
=
d
p
[
i
]
+
y
[
i
]
∗
x
[
i
]
−
a
[
i
]
f[i]=dp[i]+y[i]*x[i]-a[i]
f[i]=dp[i]+y[i]∗x[i]−a[i]
不满足条件的情况是当前队头的凸包直线斜率大于
y
[
i
]
y[i]
y[i],因为
y
[
i
]
y[i]
y[i]递减,如果当前队头的凸包直线斜率大于
y
[
i
]
y[i]
y[i],那么这条直线以后都不会被选到。
队尾的处理是维护凸包,就是直接比斜率,我们这里维护的是上凸壳(因为
y
[
i
]
y[i]
y[i]是单调递减的)。
#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=1e6+5;
inline int read(){
char c=getchar();int t=0,f=1;
while((!isdigit(c))&&(c!=EOF)){if(c=='-')f=-1;c=getchar();}
while((isdigit(c))&&(c!=EOF)){t=(t<<3)+(t<<1)+(c^48);c=getchar();}
return t*f;
}
int n;
int dp[maxn];
struct node{
int x,y,a;
}p[maxn];
int h,t;
bool cmp(node a,node b){
return a.x<b.x;
}
inline double slope(int l,int r){
return 1.0*(dp[r]-dp[l])/(p[r].x-p[l].x);
}
int q[maxn];
signed main(){
n=read();
for(int i=1;i<=n;i++){
p[i].x=read(),p[i].y=read(),p[i].a=read();
}
sort(p+1,p+1+n,cmp);
h=t=1;
for(int i=1;i<=n;i++){
while(h+1<=t&&slope(q[h],q[h+1])>=p[i].y)h++;
int j=q[h];
dp[i]=max(dp[i-1],dp[j]+(p[i].x-p[j].x)*p[i].y-p[i].a);
while(h+1<=t&&slope(q[t-1],q[t])<=slope(q[t],i))t--;
q[++t]=i;
}
printf("%lld\n",dp[n]);
return 0;
}