POJ 1151 Atlantis(扫描线)

题目原链接:http://poj.org/problem?id=1151

题目中文翻译:

POJ 1151 Atlantis

Time Limit: 1000MS		Memory Limit: 10000K
Total Submissions: 25769		Accepted: 9477

Description

有几个古希腊文本包含传说中的亚特兰蒂斯岛的描述。 其中一些文本甚至包括岛屿部分地图。 但不幸的是，这些地图描述了亚特兰蒂斯的不同区域。 您的朋友Bill必须知道地图的总面积。 你（不明智地）自告奋勇写了一个计算这个数量的程序。

Input

输入包含几个测试用例。 每个测试用例都以一行包含一个整数n（1 <= n <= 100）开始，指示可用的地图。以下n行描述了每个地图。 这些行中的每一行包含四个数字x1; y1; x2; y2（0 <= x1 <x2 <= 100000; 0 <= y1 <y2 <= 100000），不一定是整数。 值（x1; y1）和（x2; y2）是地图左上角和右下角的坐标。

输入文件以包含单个0的行作为终止。不处理它。

Output

对于每个测试用例，您的程序应输出一个部分。 每个部分的第一行必须是“Test case #k”，其中k是测试用例的编号（从1开始）。 第二个必须是“Total explored area：a”，其中a是总探索区域（即此测试用例中所有矩形的覆盖区域），精确到小数点后两位数。

在每个测试用例后输出一个空行。

Sample Input

2

10 10 20 20

15 15 25 25.5

0

Sample Output

Test case #1

Total explored area: 180.00

解题思路:

本人太菜,无法描述,请看大佬详解:AKIOI

AC代码:

(由此大佬博客借鉴而来：code)

 #include<cstdio>
 #include<iostream>
 #include<algorithm>

 using namespace std;

 struct kkk{//线段树
     int l,r;//线段树的左右整点
     int c;//c用来记录重叠情况
     double cnt,lf,rf;//cnt用来计算实在的长度，lf,rf分别是对应的左右真实的浮点数端点
 }s[];

 struct k2{
     double x,y1,y2;
     int f;
 }l[];
 //把一段段平行于y轴的线段表示成数组 ，
 //x是线段的x坐标，y1,y2线段对应的下端点和上端点的坐标
 //一个矩形 ，左边的那条边f为1，右边的为-1，
 //用来记录重叠情况，可以根据这个来计算，kkk节点中的c 

 double y[];//记录y坐标的数组

 bool cmp(k2 a,k2 b) {
     return a.x < b.x;
 }

 void build(int t,int l,int r) {
     s[t].l = l;s[t].r = r;
     s[t].cnt = s[t].c = ;
     s[t].lf = y[l];
     s[t].rf = y[r];
     if(l +  == r) return ;
     int mid = (l + r) >> ;
     build(t << ,l,mid);
     build(t <<  | ,mid,r);
 }

 void calen(int t) {//计算长度
     if(s[t].c > ) {// t对应对节点有线段覆盖
         s[t].cnt = s[t].rf - s[t].lf;
         return ;
     }
     // 现在是t对应的线段没有完整的被覆盖，但是他的孩子节点可能部分被覆盖
     if(s[t].l +  == s[t].r) s[t].cnt = ;//线段树叶子结点，代表一个点，特判，长度为0
     else s[t].cnt = s[t<<].cnt + s[t<<|].cnt;//否则，用孩子的和来表示
 }

 void update(int t,k2 e) {//加入线段e，后更新线段树
     if(e.y1 == s[t].lf && e.y2 == s[t].rf) {//如果正好找到区间
         s[t].c += e.f;
         calen(t);
         return ;
     }
     if(e.y2 <= s[t<<].rf) update(t<<,e);//下传左儿子
     else if(e.y1 >= s[t<<|].lf) update(t<<|,e);//下传右儿子
     else {//左右儿子都下传
         k2 tmp = e;
         tmp.y2 = s[t<<].rf;
         update(t<<,tmp);
         tmp = e;
         tmp.y1 = s[t<<|].lf;
         update(t<<|,tmp);
     }
     calen(t);
 }

 int main() {
     int i,n,t,aa = ;
     double x1,y1,x2,y2;
     while(scanf("%d",&n),n) {
         aa++;
         t = ;
         for(int i = ;i <= n; i++) {
             scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
             l[t].x = x1;
             l[t].y1 = y1;
             l[t].y2 = y2;
             l[t].f = ;
             y[t++] = y1;
             l[t].x = x2;
             l[t].y1 = y1;
             l[t].y2 = y2;
             l[t].f = -;
             y[t++] = y2;
         }
         sort(l+,l+t,cmp);
         sort(y+,y+t);
         build(,,t-);//建树
         update(,l[]);//下传lazy标记
         double res = ;
         for(int i = ;i < t; i++) {
             res += s[].cnt * (l[i].x - l[i-].x);
             update(,l[i]);
         }
         printf("Test case #%d\nTotal explored area: %.2f\n\n",aa,res);
     }
 }