洛谷P1884 [USACO12FEB]Overplanting S （矩形切割）

一种矩形切割的做法：

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 typedef long long LL;
 4 const int maxn=1005;
 5 struct node{//矩形的结构体
 6     LL x1,y1,x2,y2;
 7 }a[maxn];
 8 int n,cnt=0;
 9
10 void Cut(LL id,LL x1,LL y1,LL x2,LL y2){//矩形切割
11     LL k1,k2,k3,k4;
12     k1=max(a[id].x1,x1);
13     k2=min(a[id].x2,x2);
14     k3=min(a[id].y1,y1);
15     k4=max(a[id].y2,y2);
16     if(a[id].x1<k1) a[++cnt]=(node){a[id].x1,a[id].y1,k1,a[id].y2};
17     if(a[id].x2>k2) a[++cnt]=(node){k2,a[id].y1,a[id].x2,a[id].y2};
18     if(a[id].y1>k3) a[++cnt]=(node){k1,a[id].y1,k2,k3};
19     if(a[id].y2<k4) a[++cnt]=(node){k1,k4,k2,a[id].y2};
20 }
21
22 int read(){
23     int x=0,f=1;
24     char c=getchar();
25     while(c<'0'||c>'9'){if(c=='-') f=-1;c=getchar();}
26     while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
27     return x*f;
28 }
29
30 int main(){
31     n=read();
32     LL x1,y1,x2,y2;
33     x1=read();y1=read();x2=read();y2=read();
34     a[++cnt]=(node){x1,y1,x2,y2};//先加入一个矩形
35     for(int i=2;i<=n;i++){
36         x1=read();y1=read();x2=read();y2=read();
37         for(int j=1;j<=cnt;j++){
38             if(a[j].x1>=x2||a[j].x2<=x1||a[j].y1<=y2||a[j].y2>=y1)continue;
39             Cut(j,x1,y1,x2,y2);//若相交，用当前矩形去切割
40             a[j--]=a[cnt--];//删除原矩形，用最后一个矩形来覆盖它
41         }
42         a[++cnt]=(node){x1,y1,x2,y2};//加入新矩形
43     }
44     LL ans=0;
45     for(int i=1;i<=cnt;i++)
46         ans+=(a[i].x2-a[i].x1)*(a[i].y1-a[i].y2);
47     cout<<ans;
48 }

好像更普遍的解法是扫描线+线段树，这里先留一个坑，以后再来补上这种做法。