poj 1389

题目概述

给定若干个矩形左下角和右上角顶点的坐标，求其面积并
矩形的边一定平行或垂直与坐标轴

时限

1000ms/3000ms

输入

每行四个非负整数，代表矩形两个顶点的坐标，每组数据到四个-1结束，输入到连续两行四个-1结束

限制

1<=矩形数<=1000;0<=坐标值<=50000

输出

每行一个数，为所求面积并

样例输入

0 0 4 4
1 1 5 2
1 1 2 5
-1 -1 -1 -1
0 0 2 2
1 1 3 3
2 2 4 4
-1 -1 -1 -1
-1 -1 -1 -1　

样例输出

18
10

讨论

计算几何，扫描线+离散化+线段树，这个题和刚写出题解的poj 1151除了输入输出格式外几乎一致，额不打算再赘述一遍，注释也会少一些
到现在仍然有疑问，线段树到底需要多大的空间？开8倍都不够，10倍才过，平时的题3倍就够了
另外这回采用的是poj 1151中第二版的写法，追求精简

题解状态

316K，16MS，C++，1663B

题解代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define INF 0x3f3f3f3f
#define MAXN 1002
#define memset0(a) memset(a,0,sizeof(a))
#define EPS 1e-8

struct Ln//水平线的结构 扫描线也是水平的
{
    int l, r, y, f;//水平线的左右端点横坐标 纵坐标 以及上下底边的标记
    bool operator<(Ln &b)const
    {
        return y < b.y;
    }
}lns[MAXN * 2];
int N;//矩形总数
int hashx[MAXN * 2];//将横坐标离散化的数组
int L[MAXN * 10], R[MAXN * 10], c[MAXN * 10], len[MAXN * 10];//这四个是线段树用的 节点表示的范围（离散下标） 节点覆盖值 节点有效宽度
void build(int i, int l, int r)//建树 i是内部参数 后面两个是节点表示的范围
{
    L[i] = l, R[i] = r, c[i] = 0;
    if (l + 1 != r) {
        build(i * 2, l, (l + r) / 2);
        build(i * 2 + 1, (l + r) / 2, r);
    }
}
void modify(int i, Ln &a)//更新操作
{
    if (a.l == hashx[L[i]] && a.r == hashx[R[i]])
        c[i] += a.f;
    else if (a.l >= hashx[L[i * 2 + 1]])
        modify(i * 2 + 1, a);
    else if (a.r <= hashx[R[i * 2]])
        modify(i * 2, a);
    else {
        Ln b = a;
        b.r = hashx[R[i * 2]];
        modify(i * 2, b);
        b = a;
        b.l = hashx[L[i * 2 + 1]];
        modify(i * 2 + 1, b);
    }
    if (c[i])//计算有效长度
        len[i] = hashx[R[i]] - hashx[L[i]];
    else
        len[i] = len[i * 2] + len[i * 2 + 1];
}
int fun()
{
    int x1, y1, x2, y2;
    while (scanf("%d%d%d%d", &x1, &y1, &x2, &y2) && x1 >= 0) {//input//数据范围都是非负的 有一个负的就停止读入
        N++;//先+1后赋值以配合线段树下标从1开始的特性
        lns[N].l = x1, lns[N].r = x2, lns[N].y = y1, lns[N].f = 1, hashx[N] = x1;
        N++;
        lns[N].l = x1, lns[N].r = x2, lns[N].y = y2, lns[N].f = -1, hashx[N] = x2;
    }
    sort(lns + 1, lns + N + 1);
    sort(hashx + 1, hashx + N + 1);
    build(1, 1, N);
    double area = 0;
    for (int p = 1; p <= N; p++) {
        modify(1, lns[p]);
        area += len[1] * (lns[p + 1].y - lns[p].y);
    }
    return area;
}
int main(void)
{
    //freopen("vs_cin.txt", "r", stdin);
    //freopen("vs_cout.txt", "w", stdout);

    char c;
    while (~scanf("%c", &c) && c != '-') {//数据范围非负 因而只要碰到负号就停
        ungetc(c, stdin);//如果不是负号 那就是个数字 还得退回输入流
        printf("%d\n", fun());
        getchar();//吞掉行尾换行符 方便识别下一行行首字符
        N = 0;//N不能忘记归零 这是本题唯一一个需要清零的变量
    }
}

EOF