链接:
题意:
给定一个竖直的平面和平面的宽度,在此范围内会有价值为v[i]的n个馅饼分别在t[i]时刻落在p[i]位置,假设每秒钟最多移动两个单位长度,初始位置任选,问接到馅饼的最大价值是多少。
思路:
首先明确是dp,然后考虑dp的数组怎么开。
其实是见过一个类似的题目,当时是给南阳师范学院校赛验题的时候见到的:
馅饼这道题显然是"轻羽飞扬"(hdu1176)的加强版,这道题将原题中0-10的位置update为W,移动速度改为<=2,并为物品加上了value。
当时的dp数组很好开,dp[i][j]表示第i秒,在j位置能获得的最大收益。
而在本题,先不说t的范围,仅仅对W开一维数组,空间上就不能承受(W <= 1e8)。
虽然如此,不如我们先把dp的式子写出来,然后考虑能不能优化。
首先根据时间排个序,dp[i]表示当第i个馅饼必须选时的最大收益。通过枚举1~i-1个馅饼,判断能否从第j个馅饼跑到第i个,如果可以,取max。这样做的时间复杂度是O(n^2),空间复杂度O(n),当n<=1e4时还可以尝试写一下,1e5肯定会T
那么要改变dp的思路,考虑,什么时候会导致状态转移?——当接到一个馅饼后跑到下一个馅饼时。
由于我们一秒可以移动两格,也就是每半秒可以移动一格,所以,将每个馅饼时间*2,与每秒最多移动一格等价。
假设接到第i个馅饼之后,能接到第j个馅饼,那么必然满足这个式子:
p[j].time - p[i].time >= |p[i].position - p[j].position|
拆开绝对值符号,分类讨论:
当p[i].pos > p[j].pos时,p[i].time + p[i].pos <= p[j].time + p[j].pos
不等号左侧减去一个较大值2*p[i].pos,右侧减去一个较小值2*p[j].pos
不等关系仍成立:p[i].time - p[i].pos <= p[j].time - p[j].pos
当p[i].pos < p[j].pos时,p[i].time - p[i].pos <= p[j].time - p[j].pos
不等号左侧加上一个较小值2*p[i].pos,右侧加上一个较大值2*p[j].pos
不等关系仍成立:p[i].time + p[i].pos <= p[j].time + p[j].pos
两式中的"="至多有一个成立(不存在两个完全相同的馅饼)
所以,问题转化为,令x = p[i].time+p[i].pos,y = p[i].time-p[i].pos,当I(xi,yi),J(xj,yj)为坐标系中两点且I在J的左下角(允许x值同或y值同)时,表示可以转移。
用树状数组维护,加上之前说的状态转移,并把y值离散化,根据大小关系分配id,此题可解。
代码:
#include <bits/stdc++.h>
#define lowB(x) (x&(-x))
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
static const int maxn = 100010;
static const int INF = 0x3f3f3f3f;
static const int mod = (int)1e9 + 7;
static const double eps = 1e-6;
static const double pi = acos(-1);
int a[maxn],ys[maxn],n;
struct node{
int time,pos,val,x,y;
bool operator <(const node& other)const{
if (x == other.x) return y>other.y;
return x < other.x;
}
}p[maxn];
inline int sum(int x){
int res = 0;
while(x){
res = max(res,a[x]);
x -= lowB(x);
}
return res;
}
inline void add(int x,int y){
while(x <= n+1){
a[x] = max(a[x],y);
x += lowB(x);
}
}
void redirect(){
#ifdef LOCAL
freopen("test.txt","r",stdin);
#endif
}
int main(){
redirect();
int W,tmp,ans = 0;
scanf("%d %d",&W,&n);
ys[0] = -W;
for(int i = 1;i <= n;i++){
scanf("%d %d %d",&p[i].time,&p[i].pos,&p[i].val);
p[i].time <<= 1;
p[i].x = p[i].time+p[i].pos;
p[i].y = p[i].time-p[i].pos;
ys[i] = p[i].y;
}
sort(p+1,p+1+n);
sort(ys+1,ys+1+n);
for(int i = 1;i <= n;i++){
int id = p[i].y = lower_bound(ys,ys+n+1,p[i].y)-ys+1;
tmp = p[i].val + sum(id);
ans = max(ans,tmp);
add(id,tmp);
}
printf("%d\n",ans);
return 0;
}