BZOJ 1560 [JSOI2009] 火星藏宝图

Description

Input

Output

Sample Input

4 10
1 1 20
10 10 10
3 5 60
5 3 30

Sample Output

-4

HINT

Source

JSOI2009Day2

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

神奇的DP～

类似于斜率优化的做法，先将所有点按照先x后y的顺序排序，再按顺序DP。用pos[i]表示第i列目前的最低位置，f[i]表示第i列目前的最大值。

不需要开long long；把(1,1)加入到数组中排序会WA，不知道是为什么～

#include<cstdio>#include<iostream>#include<algorithm>using namespace std;int n,m,pos[1001],f[1001],now;struct node{int x,y,val;}a[200001];bool operator < (node u,node v){return u.x==v.x ? u.y<v.y:u.x<v.x;}int read(){int x=0,f=1;char ch=getchar();while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}while(ch>='0' && ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f;}int main(){n=read();m=read();a[1]=(node){m,m,0};n++;for(int i=2;i<=n;i++) a[i].x=read(),a[i].y=read(),a[i].val=read();sort(a+1,a+n+1);pos[1]=1;for(int i=1;i<=n;i++){now=-0x7fffffff;for(int j=1;j<=a[i].y;j++)  if(pos[j]) now=max(now,f[j]-(a[i].y-j)*(a[i].y-j)-(a[i].x-pos[j])*(a[i].x-pos[j]));pos[a[i].y]=a[i].x;f[a[i].y]=now+a[i].val;}printf("%d\n",f[m]);return 0;}