BZOJ 1560 [JSOI2009] 火星藏宝图

Description

Input

Output

Sample Input

Sample Output

HINT

Source

神奇的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;
}