ECNU OJ 3353 塗黑板(線段樹離散化)

kewlgrl發表於2018-03-24

3353. 塗黑板

DESCRIPTIONSTATISTICSDISCUSSION

Time limit per test: 3.0 seconds

Memory limit: 256 megabytes

林吉吉是一位富有藝術氣息的紳士，讀論文之餘他最喜歡做的事情就是塗教室裡的黑板。
黑板可以看作一塊二維平面，林吉吉想要在黑板上塗黑 n 個矩形，
第 i 個矩形利用其左下角座標 (li,di) 和右上角座標 (ri,ui) 表示。

在塗黑了所有之後，林吉吉想要知道黑板上被塗黑的面積究竟有多少。

Input

第一行為資料組數 T

每組資料第一行為矩陣數目 n，接下來 n 行為 li,di,ri,ui

資料約束：
對於 40% 的資料保證 n≤1020≤li≤ri≤1020≤di≤ui≤102 
對於 70% 的資料保證 n≤1020≤li≤ri≤1090≤di≤ui≤109 
對於所有資料保證 T≤100n≤1030≤li≤ri≤1090≤di≤ui≤109 座標為整數

Output

每組資料輸出一行整數，黑板被塗黑的面積

Examples

input
2
2
1 1 10 10
5 5 15 15
2
1 1 5 5
10 10 15 15
output
156
41

偷的HDU 1542的程式碼，直接改了改輸入輸出，double變ll。

浮點數先要離散化;然後把矩形分成兩條邊,上邊和下邊,對橫軸建樹,然後從下到上掃描上去,

用cnt表示該區間下邊比上邊多幾個。


#include<bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define MAXN 100010
#include <cstdio>
#include <cstring>
#include <cctype>
#include <algorithm>
#define lson l , m , rt << 1
#define rson m + 1 , r , rt << 1 | 1
const int maxn = 2222;
int cnt[maxn << 2];
long long sum[maxn << 2];
long long X[maxn];
struct Seg
{
    long long h, l, r;
    int s;
    Seg() {}
    Seg(long long a,long long b,long long c,int d) : l(a), r(b), h(c), s(d) {}
    bool operator < (const Seg &cmp) const
    {
        return h < cmp.h;
    }
} ss[maxn];
void PushUp(int rt,int l,int r)
{
    if (cnt[rt]) sum[rt] = X[r+1] - X[l];
    else if (l == r) sum[rt] = 0;
    else sum[rt] = sum[rt<<1] + sum[rt<<1|1];
}
void update(int L,int R,int c,int l,int r,int rt)
{
    if (L <= l && r <= R)
    {
        cnt[rt] += c;
        PushUp(rt, l, r);
        return ;
    }
    int m = (l + r) >> 1;
    if (L <= m) update(L, R, c, lson);
    if (m < R) update(L, R, c, rson);
    PushUp(rt, l, r);
}
int Bin(long long key,int n,long long X[])
{
    int l = 0, r = n - 1;
    while (l <= r)
    {
        int m = (l + r) >> 1;
        if (X[m] == key) return m;
        if (X[m] < key) l = m + 1;
        else r = m - 1;
    }
    return -1;
}
int main()
{
    #ifdef ONLINE_JUDGE
#else
    freopen("F:/cb/read.txt","r",stdin);
//freopen("F:/cb/out.txt","w",stdout);
#endif
    /*ios::sync_with_stdio(false);
    cin.tie(0);*/
    int n, cas = 1;
    int ca;
    scanf("%d",&ca);
    while(ca--)
    {
        scanf("%d",&n);
        int m = 0;
        while (n --)
        {
            long long a, b, c, d;
            scanf("%lld%lld%lld%lld",&a,&b,&c,&d);
            X[m] = a;
            ss[m++] = Seg(a, c, b, 1);
            X[m] = c;
            ss[m++] = Seg(a, c, d, -1);
        }
        sort(X, X + m);
        sort(ss, ss + m);
        int k = 1;
        for (int i = 1 ; i < m ; i ++)
        {
            if (X[i] != X[i-1]) X[k++] = X[i];
        }
        memset(cnt, 0, sizeof(cnt));
        memset(sum, 0, sizeof(sum));
        long long ret = 0;
        for (int i = 0 ; i < m - 1 ; i ++)
        {
            int l = Bin(ss[i].l, k, X);
            int r = Bin(ss[i].r, k, X) - 1;
            if (l <= r) update(l, r, ss[i].s, 0, k - 1, 1);
            ret += sum[1] * (ss[i+1].h - ss[i].h);
        }
        printf("%lld\n",ret);
    }
    return 0;
}

相關文章