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;
}

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章