zoukankan      html  css  js  c++  java
  • D. Babaei and Birthday Cake--- Codeforces Round #343 (Div. 2)

    D. Babaei and Birthday Cake
    time limit per test
    2 seconds
    memory limit per test
    256 megabytes
    input
    standard input
    output
    standard output

    As you know, every birthday party has a cake! This time, Babaei is going to prepare the very special birthday party's cake.

    Simple cake is a cylinder of some radius and height. The volume of the simple cake is equal to the volume of corresponding cylinder. Babaei has n simple cakes and he is going to make a special cake placing some cylinders on each other.

    However, there are some additional culinary restrictions. The cakes are numbered in such a way that the cake number i can be placed only on the table or on some cake number j where j < i. Moreover, in order to impress friends Babaei will put the cake i on top of the cake j only if the volume of the cake i is strictly greater than the volume of the cake j.

    Babaei wants to prepare a birthday cake that has a maximum possible total volume. Help him find this value.

    Input

    The first line of the input contains a single integer n (1 ≤ n ≤ 100 000) — the number of simple cakes Babaei has.

    Each of the following n lines contains two integers ri and hi (1 ≤ ri, hi ≤ 10 000), giving the radius and height of the i-th cake.

    Output

    Print the maximum volume of the cake that Babaei can make. Your answer will be considered correct if its absolute or relative error does not exceed 10 - 6.

    Namely: let's assume that your answer is a, and the answer of the jury is b. The checker program will consider your answer correct, if .

    Examples
    input
    2
    100 30
    40 10
    output
    942477.796077000
    input
    4
    1 1
    9 7
    1 4
    10 7
    output
    3983.539484752
    Note

    In first sample, the optimal way is to choose the cake number 1.

    In second sample, the way to get the maximum volume is to use cakes with indices 1, 2 and 4.

    线段树维护比当前值a[i]小的最大dp[j]的值。j<i

    我们可以用体积作为线段树的横坐标(就那个意思).

    也可以直接上。(就是官方题解)

    /* ***********************************************
    Author        :guanjun
    Created Time  :2016/2/24 9:56:13
    File Name     :cf343d.cpp
    ************************************************ */
    #include <iostream>
    #include <cstring>
    #include <cstdlib>
    #include <stdio.h>
    #include <algorithm>
    #include <vector>
    #include <queue>
    #include <set>
    #include <map>
    #include <string>
    #include <math.h>
    #include <stdlib.h>
    #include <iomanip>
    #include <list>
    #include <deque>
    #include <stack>
    #define ull unsigned long long
    #define ll long long
    #define mod 90001
    #define INF 0x3f3f3f3f
    #define maxn 101000
    #define cle(a) memset(a,0,sizeof(a))
    const ull inf = 1LL << 61;
    const double eps=1e-5;
    using namespace std;
    priority_queue<int,vector<int>,greater<int> >pq;
    struct Node{
    int x,y;
    };
    struct cmp{
        bool operator()(Node a,Node b){
            if(a.x==b.x) return a.y> b.y;
            return a.x>b.x;
        }
    };
    
    bool cmp(int a,int b){
        return a>b;
    }
    struct node{
        ll Max;
        int l,r;
        int getmid(){return (l+r)/2;}
    }nod[maxn*4];
    void push_up(int i){
        nod[i].Max=max(nod[i<<1].Max,nod[i<<1|1].Max);
    }
    void build(int i,int l,int r){
        nod[i].l=l,nod[i].r=r;
        if(l==r){
            nod[i].Max=0; return ;
        }
        int mid=(l+r)/2;
        build(i<<1,l,mid);
        build(i<<1|1,mid+1,r);
    }
    void update(int i,int k,ll w){
        if(nod[i].l==k&&nod[i].r==k){
            nod[i].Max=max(nod[i].Max,w);return ;
        }
        int mid=nod[i].getmid();
        if(k<=mid)update(i<<1,k,w);
        else update(i<<1|1,k,w);
        push_up(i);
    }
    double query(int i,int l,int r){
        if(r<l)return 0;
        if(nod[i].l==l&&nod[i].r==r){
            return nod[i].Max;
        }
        int mid=nod[i].getmid();
        if(r<=mid)return query(i<<1,l,r);
        else if(l>mid)return query(i<<1|1,l,r);
        else return max(query(i<<1,l,mid),query(i<<1|1,mid+1,r));
    }
    map<ll,int>mp;
    ll a[maxn],b[maxn];
    int main()
    {
        #ifndef ONLINE_JUDGE
        freopen("in.txt","r",stdin);
        #endif
        //freopen("out.txt","w",stdout);
        int n;
        while(cin>>n){
            ll x,y;
            mp.clear();
            for(int i=1;i<=n;i++){
                cin>>x>>y;
                b[i]=a[i]=x*x*y;
            }
            sort(b+1,b+1+n);
            int num=0;
            for(int i=1;i<=n;i++){
                if(!mp[b[i]])mp[b[i]]=++num;
            }
            build(1,1,n);
            //puts("YEs");
            for(int i=1;i<=n;i++){
                //cout<<mp[a[i]]<<endl;
                ll dp=query(1,1,mp[a[i]]-1)+a[i];
                update(1,mp[a[i]],dp);
            }
            printf("%.8lf
    ",3.1415926*nod[1].Max);
        }
        return 0;
    }
    View Code
    /* ***********************************************
    Author        :guanjun
    Created Time  :2016/2/23 21:49:30
    File Name     :cf343d.cpp
    ************************************************ */
    #include <iostream>
    #include <cstring>
    #include <cstdlib>
    #include <stdio.h>
    #include <algorithm>
    #include <vector>
    #include <queue>
    #include <set>
    #include <map>
    #include <string>
    #include <math.h>
    #include <stdlib.h>
    #include <iomanip>
    #include <list>
    #include <deque>
    #include <stack>
    #define ull unsigned long long
    #define ll long long
    #define mod 90001
    #define INF 0x3f3f3f3f
    #define maxn 101100
    #define cle(a) memset(a,0,sizeof(a))
    #define X first
    #define Y second
    const ull inf = 1LL << 61;
    const double eps=1e-5;
    using namespace std;
    priority_queue<int,vector<int>,greater<int> >pq;
    pair<ll,int>a[maxn];
    struct node{
        int l,r;
        ll Max;
    }nod[maxn*4];
    
    void push_up(int i){
        nod[i].Max=max(nod[i<<1].Max,nod[i<<1|1].Max);
    }
    void build(int i,int l,int r){
        nod[i].l=l;
        nod[i].r=r;
        if(l==r){
            nod[i].Max=0;
            return ;
        }
        int mid=(l+r)/2;
        build(i<<1,l,mid);
        build(i<<1|1,mid+1,r);
        push_up(i);
    }
    void update(int i,int k,ll w){
        if(nod[i].l==k&&nod[i].r==k){
            nod[i].Max=max(nod[i].Max,w);
            return ;
        }
        int mid=(nod[i].l+nod[i].r)/2;
        if(k<=mid)update(i<<1,k,w);
        else update(i<<1|1,k,w);
        push_up(i);
    }
    double query(int i,int l,int r){
        if(nod[i].l==l&&nod[i].r==r){
            return nod[i].Max;
        }
        int mid=(nod[i].l+nod[i].r)/2;
        if(r<=mid)return query(i<<1,l,r);
        else if(l>mid)return query(i<<1|1,l,r);
        else return max(query(i<<1,l,mid),query(i<<1|1,mid+1,r));
    }
    ll dp[maxn];
    int main()
    {
        #ifndef ONLINE_JUDGE
        freopen("in.txt","r",stdin);
        #endif
        //freopen("out.txt","w",stdout);
        int n;
        while(cin>>n){
            ll x,y; 
            for(int i=1;i<=n;i++){
                scanf("%I64d %I64d",&x,&y);
                a[i].X=x*x*y;
                a[i].Y=-i;
            }
            build(1,1,n);
            sort(a+1,a+1+n);
            for(int i=1;i<=n;i++){
                int id=-a[i].Y;
                //cout<<id<<endl;
                dp[id]=query(1,1,id)+a[i].X;
                
                update(1,id,dp[id]);
            }
            printf("%.8lf
    ",nod[1].Max*3.1415926);
        }
        return 0;
    }
    View Code

    注意id 赋值为-i。这样才能保证相同的体积出现时我们取最后那一个。这里不用-i也行,也可以自己写一个比较函数。

  • 相关阅读:
    leetcode 14. 最长公共前缀
    leetcode13. 罗马数字转整数
    leetcode 21.合并两个有序链表
    leetcode14-最长公共前缀
    leetcode13-罗马数字转整数
    leetcode11- 盛最多水的容器
    leetcode9-回文数
    leetcode7- 整数反转
    leetcode 1-两数之和
    文件的上传
  • 原文地址:https://www.cnblogs.com/pk28/p/5212093.html
Copyright © 2011-2022 走看看