POJ 1151 矩形切割

【题目】http://poj.org/problem?id=1151

【题意】求矩形面积并。

【解题方法】可以线段树扫描线，当然这里我是来试一试矩形切割的。

【AC 代码】在POJ和HDU均测试通过。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 105;
struct node{
    double x1,y1,x2,y2,sum;
}T[N];
int n;
void Cover(double x1,double y1,double x2,double y2,int k,int c)
{
    while(k<n&&(x1>=T[k].x2||x2<=T[k].x1||y1>=T[k].y2||y2<=T[k].y1)) k++;
    if(k>=n){
        T[c].sum+=(x2-x1)*(y2-y1);
        return ;
    }
    if(x1<T[k].x1){
        Cover(x1,y1,T[k].x1,y2,k+1,c);
        x1=T[k].x1;
    }
    if(x2>T[k].x2){
        Cover(T[k].x2,y1,x2,y2,k+1,c);
        x2=T[k].x2;
    }
    if(y1<T[k].y1){
        Cover(x1,y1,x2,T[k].y1,k+1,c);
        y1=T[k].y1;
    }
    if(y2>T[k].y2){
        Cover(x1,T[k].y2,x2,y2,k+1,c);
        y2=T[k].y2;
    }
}
int main()
{
    int ks=1;
    while(scanf("%d",&n)!=EOF&&n)
    {
        for(int i=0; i<n; i++){
            scanf("%lf%lf%lf%lf",&T[i].x1,&T[i].y1,&T[i].x2,&T[i].y2);
            T[i].sum=0;
        }
        for(int i=n-1; i>=0; i--){
            Cover(T[i].x1,T[i].y1,T[i].x2,T[i].y2,i+1,i);
        }
        double ans=0.0;
        for(int i=0; i<n; i++){
            ans+=T[i].sum;
        }
        printf("Test case #%d\n",ks++);
        printf("Total explored area: %.2f\n\n",ans);
    }
    return 0;
}