1. 程式人生 > >629D Babaei and Birthday Cake (線段樹+DP裸題目)

629D Babaei and Birthday Cake (線段樹+DP裸題目)

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define read(x,y) scanf("%d%d",&x,&y)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
const int  maxn =1e5+5;
const int mod=1e9+7;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
題目大意:就是給個蛋糕序列,
大的可以放到小的上面。。T_T。
求最大權值。

太經典了。。。。最大權值遞增序列。
不過以防萬一我還是離散化了一下,
因為體積值可能爆12次方。

把體積值對映成大小關係就行了,
在往樹狀數組裡面丟DP值時再對映成其體積就行了,
特別裸,對於遍歷到的i位置在樹狀陣列中找到最大的(下標小於當前位)DP值進行更新,
更新完後再丟進去就行了。


*/
#define pi acos(-1)
ll r[maxn],h[maxn],a[maxn],b[maxn];
ll dp[maxn];
///樹狀陣列
ll bit[maxn<<1];
int lowbit(int x){return x&(-x);}
void add(int x,ll d)
{
    for(;x<maxn;bit[x]=max(bit[x],d),x+=lowbit(x));
}
ll sum(int x)
{
    ll ret=0;
    for(;x>0;ret=max(ret,bit[x]),x-=lowbit(x));
    return ret;
}
ll n,ans;

int main()
{
    memset(bit,0,sizeof(bit));
    scanf("%lld",&n);
    for(int i=0;i<n;i++)
    {
        scanf("%lld%lld",&r[i],&h[i]);
        b[i]=r[i]*r[i]*h[i];
    }
    sort(b,b+n);///離散化
    for(int i=0;i<n;i++)
    {
        ll tp=r[i]*r[i]*h[i];
        a[i]=lower_bound(b,b+n,tp)-b+1;///離散化大小關係
        dp[i]=sum(a[i]-1)+b[a[i]-1];
        ans=max(ans,dp[i]);
        add(a[i],dp[i]);
    }
    printf("%.9f\n",1.0*ans*pi);
    return 0;
}