BZOJ2131 免费的馅饼 二维排序+离散化+树状数组dp

链接：

免费的馅饼

题意：

给定一个竖直的平面和平面的宽度，在此范围内会有价值为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;
}