【坐标离散化】矩形面积求并

传送门:CodeVS3044


这道题如果数据范围小而且是整数的话实际很简单,只需附加到坐标系中,暴力出被覆盖的格子即可 但是由于数据是小数而且很大,直接暴力显然不支持,于是就把行列坐标分别排序,然后按大小一一对应,换句话说把第k小的行坐标直接用k表示,然后暴力,最后求面积时在把点坐标反带回去算边长


详见代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
int n,b[105][5];
double a[105][5],x[205],y[205],ans;
bool f[205][205];
void init(int n){
for (int i=1;i<=n;i++){
a[i][1]=a[i][2]=a[i][3]=a[i][4]=0;
b[i][1]=b[i][2]=b[i][3]=b[i][4]=0;
}
for (int i=1;i<=2*n;i++)
x[i]=y[i]=0;
for (int i=1;i<=2*n;i++)
for (int j=1;j<=2*n;j++)
f[i][j]=0;
ans=0;
}
int find(double s[205],double N){
for (int i=1;i<=2*n;i++)
if (s[i]==N) return i;
return -1;
}
void solve(){
init(n);
for (int i=1;i<=n;i++){
scanf("%lf%lf%lf%lf",&a[i][1],&a[i][2],&a[i][3],&a[i][4]);
x[i*2-1]=a[i][1];
y[i*2-1]=a[i][2];
x[i*2]=a[i][3];
y[i*2]=a[i][4];
}
sort(x+1,x+2*n+1);
sort(y+1,y+2*n+1);
for (int i=1;i<=n;i++){
b[i][1]=find(x,a[i][1]);
b[i][2]=find(y,a[i][2]);
b[i][3]=find(x,a[i][3]);
b[i][4]=find(y,a[i][4]);
}
for (int i=1;i<=n;i++)
for (int X=b[i][1];X<b[i][3];X++)
for (int Y=b[i][2];Y<b[i][4];Y++)
f[X][Y]=1;
for (int i=1;i<=n*2;i++)
for (int j=1;j<=n*2;j++)
if (f[i][j])
ans+=(x[i+1]-x[i])*(y[j+1]-y[j]);
printf("%.2f\n",ans);
}
int main(){
while (scanf("%d",&n) && n) solve();
return 0;
}

最后经柏神提醒这道题应该是可以线段树优化的 应该是一行一行求吧