zoukankan      html  css  js  c++  java
  • HDU-3015 Disharmony Trees [数状数组]

    Problem Description
    One day Sophia finds a very big square. There are n trees in the square. They are all so tall. Sophia is very interesting in them.

    She finds that trees maybe disharmony and the Disharmony Value between two trees is associated with two value called FAR and SHORT.

    The FAR is defined as the following:If we rank all these trees according to their X Coordinates in ascending order.The tree with smallest X Coordinate is ranked 1th.The trees with the same X Coordinates are ranked the same. For example,if there are 5 tree with X Coordinates 3,3,1,3,4. Then their ranks may be 2,2,1,2,5. The FAR of two trees with X Coordinate ranks D1 and D2 is defined as F = abs(D1-D2).

    The SHORT is defined similar to the FAR. If we rank all these trees according to their heights in ascending order,the tree with shortest height is ranked 1th.The trees with the same heights are ranked the same. For example, if there are 5 tree with heights 4,1,9,7,4. Then their ranks may be 2,1,5,4,2. The SHORT of two trees with height ranks H1 and H2 is defined as S=min(H1,H2).

    Two tree’s Disharmony Value is defined as F*S. So from the definition above we can see that, if two trees’s FAR is larger , the Disharmony Value is bigger. And the Disharmony value is also associated with the shorter one of the two trees. 

    Now give you every tree’s X Coordinate and their height , Please tell Sophia the sum of every two trees’s Disharmony value among all trees.
     
    Input
    There are several test cases in the input

    For each test case, the first line contain one integer N (2 <= N <= 100,000) N represents the number of trees.

    Then following N lines, each line contain two integers : X, H (0 < X,H <=1,000,000,000 ), indicating the tree is located in Coordinates X and its height is H.
     
    Output
    For each test case output the sum of every two trees’s Disharmony value among all trees. The answer is within signed 64-bit integer.
     
    Sample Input
    2 10 100 20 200 4 10 100 50 500 20 200 20 100
     
    Sample Output
    1 13

    题意:求∑( abs(Xi-Xj)*min(Hi,Hj) )

    解:离散化X和H,对H升序排序,由于计算对H取的是最小值,因此只要考虑每项对后续的贡献;

    对于每一项i,∑( abs(Xi-Xj) )可以表达为 (前项个数*Xi - 前项总和) + (后项总和 - 后项个数*Xi );因此用两个树状数组分别维护个数和总和,遍历一遍求和,对每项计算完贡献后在数状数组中删去。

    #include <iostream>
    #include <cstdio>
    #include <cstdlib>
    #include <cmath>
    #include <cstring>
    #include <algorithm>
    #include <string>
    #include <vector>
    #include <queue>
    #include <stack>
    #include <set>
    #include <map>
    #define INF 0x3f3f3f3f
    #define ll long long
    #define ull unsigned long long
    #define lowbit(x) (x&(-x))
    #define eps 0.00000001
    #define PI acos(-1)
    #define pn printf("
    ");
    using namespace std;
    
    const int maxn = 1e5 + 5;
    int n;
    struct node{
        ll x, h;
    }arr[maxn];
    
    ll cnt[maxn], sum[maxn];
    
    void update(ll pos,ll val, ll *c)
    {
        while(pos <= maxn)
        {
            c[pos] += val;
            pos += lowbit(pos);
        }
    }
    
    ll query(ll pos, ll*c)
    {
        ll ret = 0;
        while(pos > 0)
        {
            ret += c[pos];
            pos -= lowbit(pos);
        }
        return ret;
    }
    
    int main()
    {
        while(~scanf("%d", &n))
        {
            memset(cnt, 0, sizeof cnt);
            memset(sum ,0, sizeof sum);
            for(int i=0;i<n;i++)
                scanf("%lld%lld",&arr[i].x, &arr[i].h);
            sort(arr, arr+n, [](node a, node b){
                return a.x < b.x;
            });
            for(int i=0;i<n;i++)
            {
                int j = i+1;
                ll mk = arr[i].x;
                arr[i].x = i+1;
                while(j < n && arr[j].x == mk)
                {
                    arr[j].x = i+1;
                    j ++;
                }
                i = j - 1;
            }
            
            sort(arr, arr+n, [](node a, node b){
                return a.h < b.h;
            });
            for(int i=0;i<n;i++)
            {
                int j = i+1;
                ll mk = arr[i].h;
                arr[i].h = i+1;
                while(j < n && arr[j].h == mk)
                {
                    arr[j].h = i+1;
                    j ++;
                }
                i = j - 1;
            }
            
            for(int i=0;i<n;i++)
            {
                update(arr[i].x, 1, cnt);
                update(arr[i].x, arr[i].x, sum);
            }
            
    
            ll ans = 0;
            for(int i=0;i<n;i++)
            {
                ll pre_cnt = query(arr[i].x - 1, cnt);
                ll pst_cnt = n-i - query(arr[i].x, cnt);
                ll pre_sum = query(arr[i].x - 1, sum);
                ll pst_sum = query(n, sum) - query(arr[i].x, sum);
         
                update(arr[i].x, -1, cnt);
                update(arr[i].x, -arr[i].x, sum);
                
                ans += ( pre_cnt * arr[i].x - pre_sum + pst_sum - pst_cnt * arr[i].x) * arr[i].h;
            }
            printf("%lld
    ", ans);
        }
    }
  • 相关阅读:
    JavaWeb---总结(七)HttpServletResponse对象(一)
    JavaWeb---总结(八)HttpServletResponse对象(二)
    JavaWeb---总结(九)通过Servlet生成验证码图片
    hdu 1102(最小生成树)
    poj 2002(好题 链式hash+已知正方形两点求另外两点)
    hdu 2461(AC) & poj 3695(TLE)(离散化+矩形并)
    poj 1389(离散化+计算几何)
    poj 1151(离散化+矩形面积并)
    hdu 3264(枚举+二分+圆的公共面积)
    poj 2546(两圆公共面积)
  • 原文地址:https://www.cnblogs.com/HazelNut/p/10300447.html
Copyright © 2011-2022 走看看