覆盖的面积
Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 5595 Accepted Submission(s): 2810
Problem Description
给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积.
![](https://i-blog.csdnimg.cn/blog_migrate/d1adc6253c233d9968f03fe1d872f86d.jpeg)
Input
输入数据的第一行是一个正整数T(1<=T<=100),代表测试数据的数量.每个测试数据的第一行是一个正整数N(1<=N<=1000),代表矩形的数量,然后是N行数据,每一行包含四个浮点数,代表平面上的一个矩形的左上角坐标和右下角坐标,矩形的上下边和X轴平行,左右边和Y轴平行.坐标的范围从0到100000.
注意:本题的输入数据较多,推荐使用scanf读入数据.
注意:本题的输入数据较多,推荐使用scanf读入数据.
Output
对于每组测试数据,请计算出被这些矩形覆盖过至少两次的区域的面积.结果保留两位小数.
Sample Input
2
5
1 1 4 2
1 3 3 7
2 1.5 5 4.5
3.5 1.25 7.5 4
6 3 10 7
3
0 0 1 1
1 0 2 1
2 0 3 1
Sample Output
7.63
0.00
Author
Ignatius.L & weigang Lee
题意:
给出矩形左下和右上角的坐标,求矩形交的面积。
代码:
//cnt[rt]>=2的一定是覆盖过两次以上的,直接计算。cnt[rt]==1时,如果rt的左右儿子 //被覆盖过大于等于1次,那么再被他们父亲覆盖上就是覆盖过两次了,值是左右儿子被 //覆盖过一次的长度的和。cnt[rt]==0照常更新。sum2记录覆盖过两次以上的边的长度。 //sum1记录覆盖过一次的边的长度。 #include<iostream> #include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int maxn=2005; double sum1[maxn*4],sum2[maxn*4],mp[maxn]; int cnt[maxn*4]; struct node{ double l,r,h; int d; node(){} node(double a,double b,double c,int d):l(a),r(b),h(c),d(d){} bool operator < (const node &p)const{ if(h==p.h) return d>p.d; return h<p.h; } }nodes[maxn*4]; int Bsearch(double a,int b,double *c) { int l=0,r=b-1,mid; while(l<=r){ mid=(l+r)>>1; if(c[mid]==a) return mid; else if(c[mid]>a) r=mid-1; else l=mid+1; } return -1; } void Pushup(int l,int r,int rt) { if(cnt[rt]){ sum1[rt]=mp[r+1]-mp[l]; if(cnt[rt]==1) sum2[rt]=sum1[rt<<1]+sum1[rt<<1|1]; else sum2[rt]=mp[r+1]-mp[l]; } else if(l==r) sum1[rt]=sum2[rt]=0; else{ sum1[rt]=sum1[rt<<1]+sum1[rt<<1|1]; sum2[rt]=sum2[rt<<1]+sum2[rt<<1|1]; } } void Update(int ql,int qr,int v,int l,int r,int rt) { if(ql<=l&&qr>=r){ cnt[rt]+=v; Pushup(l,r,rt); return; } //if(l==r) return; int m=(l+r)>>1; if(ql<=m) Update(ql,qr,v,l,m,rt<<1); if(qr>m) Update(ql,qr,v,m+1,r,rt<<1|1); Pushup(l,r,rt); } int main() { int t,n; scanf("%d",&t); while(t--){ scanf("%d",&n); int m=0,nu=0; double x1,y1,x2,y2; memset(sum1,0,sizeof(sum1)); memset(sum2,0,sizeof(sum2)); memset(cnt,0,sizeof(cnt)); for(int i=0;i<n;i++){ scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2); nodes[m++]=node(x1,x2,y1,1); nodes[m++]=node(x1,x2,y2,-1); mp[nu++]=x1;mp[nu++]=x2; } sort(mp,mp+nu); nu=unique(mp,mp+nu)-mp; sort(nodes,nodes+m); double ans=0; for(int i=0;i<m-1;i++){ int lef=Bsearch(nodes[i].l,nu,mp); int rig=Bsearch(nodes[i].r,nu,mp)-1; if(lef<=rig) Update(lef,rig,nodes[i].d,0,nu-1,1); ans+=sum2[1]*(nodes[i+1].h-nodes[i].h); } printf("%.2lf\n",ans); } return 0; }