Codeforces Round #343 (Div. 2) D. Babaei and Birthday Cake（线段树+离散化优化DP）

题目链接：点击打开链接

题意：给出n个圆柱体的地面半径和高， 要求只能有一个直接放在桌子上， 其他的要放在他上面， 第i个能放在第j个上面的条件是：当且仅当第i个的体积大于第j个且j < i 。 求能叠起来的最大体积。

思路：一看就是一个DP， 而且状态很容易表示， d[i]表示到第i个为止能得到的最大总体积。   转移到 max(d[j]) + a[i], （j < i && a[i] > a[j]）。  但是n非常大， 显然要优化， 因为第二层循环所做的事情就是在i之前找一个满足a[j] > a[i]的最大dp[j]。  那么为了同时满足两个条件， 可以把其中一个条件（每个的体积）当成线段树下标， 变相维护了一个条件， 又因为小标必须是整数， 且体积太大， 我们想到了离散化。  那么这题就很简单了， 离散化每个的体积， 用体积作为线段树下标， 维护一个区间最大值， 这个值就是之前已经加进去的dp[j]。 

细节参见代码：

#include<cstdio>#include<cstring>#include<algorithm>#include<iostream>#include<string>#include<vector>#include<stack>#include<bitset>#include<cstdlib>#include<cmath>#include<set>#include<list>#include<deque>#include<map>#include<queue>#define Max(a,b) ((a)>(b)?(a):(b))#define Min(a,b) ((a)<(b)?(a):(b))using namespace std;typedef long long ll;typedef long double ld;const ld eps = 1e-9, PI = 3.1415926535897932384626433832795;const int mod = 1000000000 + 7;const int INF = int(1e9);const ll INF64 = ll(1e18);const int maxn = 100000 + 10;int T,n,m;ll d[maxn],b[maxn];ll maxv[maxn<<2];struct node {    ll r, h, v;}a[maxn];void PushUp(int o) {    maxv[o] = max(maxv[o<<1], maxv[o<<1|1]);}void build(int l, int r, int o) {    int m = (l + r) >> 1;    maxv[o] = 0;    if(l == r) return ;    build(l, m, o<<1);    build(m+1, r, o<<1|1);}void update(int L, int R, int v, int l, int r, int o) {    int m = (l + r) >> 1;    if(L <= l && r <= R) {        if(maxv[o] < d[v]) {            maxv[o] = d[v];        }        return ;    }    if(L <= m) update(L, R, v, l, m, o<<1);    if(m < R) update(L, R, v, m+1, r, o<<1|1);    PushUp(o);}ll query(int L, int R, int l, int r, int o) {    int m = (l + r) >> 1;    if(L <= l && r <= R) {        return maxv[o];    }    ll ans = 0;    if(L <= m) ans = max(ans, query(L, R, l, m, o<<1));    if(m < R) ans = max(ans, query(L, R, m+1, r, o<<1|1));    PushUp(o);    return ans;}int main() {    scanf("%d",&n);    for(int i=1;i<=n;i++) {        scanf("%I64d%I64d",&a[i].r,&a[i].h);        a[i].v = a[i].r * a[i].r * a[i].h;        b[i] = a[i].v;    }    sort(b+1, b+1+n);    int m = unique(b+1, b+1+n) - b - 1;    build(1, m, 1);    ll cur = 0;    for(int i=1;i<=n;i++) {        int L = lower_bound(b+1, b+1+m, a[i].v) - b;        ll v;        if(L > 1) v = query(1, L-1, 1, m, 1);        else v = 0;        d[i] = (v + a[i].v);        update(L, L, i, 1, m, 1);        cur = max(cur, d[i]);    }    double ans = PI * cur;    printf("%.8f\n",ans);    return 0;}