【bzoj2731】[HNOI2012]三角形覆盖问题

题目链接

Description

二维平面中，给定 N个等腰直角三角形（每个三角形的两条直角边分别
平行于坐标轴，斜边从左上到右下）。我们用三个非负整数( x, y, d)来描
述这样一个三角形，三角形三个顶点的坐标分别为$(x, y)$, $(x + d, y)$和$(x, y+d)$。要求计算这 N个三角形所覆盖的总面积。例如，下图有 3 个三角形，覆盖的总面积为 11.0。

Input

输入文件第一行为一个正整数N，表示三角形的个数。接下来的 N
行每行有用空格隔开的三个非负整数， x, y , d，描述一个三角
形的顶点坐标，分别为
$(x, y)$,$(x + d, y)$, $( x, y+d)$，
其中 x, y, d 满足$0≤ x, y, d≤1000000$。
对于50%的数据，$1≤ N≤500$；
100%的数据，$1≤N≤10000$。

Output

仅包含一行，为一个实数S，表示所有三角形所覆盖的总面积，输出恰
好保留一位小数。输入数据保证$S≤2^{31}$ 。

Sample Input

3 .

1 1 4

2 0 2

3 2 2

Sample Output

11.0

题解

膜拜vfk当年的暴力满分解。
http://vfleaking.blog.163.com/blog/static/17480763420123301447215/
用双向链表实现一下就好了。
输入的d居然有为0的。。。

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

const int N = 10000 + 10;

struct T{
    int x, y, d, l, r;
    bool operator < (const T &a) const{
        return y < a.y;
    }
}t[N];
int n, nxt[N], pre[N], mx;
long long ans;
int hd, tl, now, last;
int cov[1000010];

void init(){
    scanf("%d", &n);
    for(int i = 1; i <= n; i++){
        scanf("%d%d%d", &t[i].x, &t[i].y, &t[i].d);
        t[i].l = t[i].x; t[i].r = t[i].x + t[i].d - 1;
        mx = max(mx, t[i].y + t[i].d);
    }
    sort(t+1, t+n+1);
}

void del(int k){
    nxt[pre[k]] = nxt[k];
    pre[nxt[k]] = pre[k];
}

bool check(int u){
    if(!t[u].d)return false;
    for(int i = nxt[hd]; i != tl; i = nxt[i])
        if(t[i].l <= t[u].l && t[u].r <= t[i].r) return false;
    pre[nxt[hd]] = u;
    nxt[u] = nxt[hd];
    nxt[hd] = u;
    pre[u] = hd;
    return true;
}

void work(){
    hd = 0, tl = n + 1;
    nxt[hd] = tl; pre[tl] = hd;
    for(int i = t[1].y, j = 1; i <= mx; i++){
        now = last;
        for(int k = nxt[hd]; k != tl; k = nxt[k]){
            cov[t[k].r]--;
            if(cov[t[k].r] == 0) now--;
            t[k].r--;
            if(t[k].r < t[k].l) del(k);
        }
        ans += now + last;
        while(j <= n && t[j].y == i){
            if(check(j))
            for(int k = t[j].l; k <= t[j].r; k++){
                if(!cov[k]) now++;
                cov[k]++;
            }
            j++;
        }
        last = now;
    }
    printf("%.1lf\n", ans / 2.0);
}

int main(){
    init();
    work();
    return 0;
}