HDU 1542 HDU 1225 (线段树扫描线）
关于扫描线，我是看的学长博客：
http://blog.csdn.net/lwt36/article/details/48908031
所以基本代码差不多。
HDU1542
题意：
最初始扫描线，想看模拟过程的可以去看上面那个博客，模拟一下就知道了。
其中有一点要注意，那就是线段树离散化后知后觉的区间，表示的是[l,r+1],
为什么要这样呢，你根据这题的样例模拟一下过程，你会发现，如果你采用的
是[l,r]的话，最右边的那个端点没有在更新是并没有被包括进来，而这个点
也是答案的一部分，这就造成了答案变小，所以一定要采用[l,r+1]这样的区间
存储方式，也就是线段树每一个节点便是的区间是[l,r+1]，否则错。代码：
//
// Created by CQU_CST_WuErli
// Copyright (c) 2016 CQU_CST_WuErli. All rights reserved.
//
// #include
#include
#include vector v;
vector<double> mp;
int n;
void pushup(int l,int r,int rt) {
if (num[rt]) sum[rt]=mp[r]-mp[l-1];
else if (l==r) sum[rt]=0;
else sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
void update(int L,int R,int val,int l,int r,int rt) {
if (L<=l && r<=R) {
num[rt]+=val;
pushup(l,r,rt);
return;
}
int mid=(l+r)>>1;
if (L<=mid) update(L,R,val,Lson);
if (R>mid) update(L,R,val,Rson);
pushup(l,r,rt);
}
int main(int argc, char const *argv[]) {
#ifdef LOCAL
freopen("C:\\Users\\john\\Desktop\\in.txt","r",stdin);
// freopen("C:\\Users\\john\\Desktop\\out.txt","w",stdout);
#endif
// BUG;
int kase=1;
while(SI(n)==1 && n) {
v.clear();mp.clear();
rep(i,1,n) {
double X1,Y1,X2,Y2;
cin >> X1 >> Y1 >> X2 >> Y2;
v.push_back(Seg(X1,X2,Y1,1));
v.push_back(Seg(X1,X2,Y2,-1));
mp.push_back(X1);mp.push_back(X2);
}
sort(v.begin(),v.end());
sort(mp.begin(),mp.end());
mp.resize(unique(mp.begin(),mp.end())-mp.begin());
// for (auto i: mp) cout << i << ' ';
// cout << endl;
CLR(num);CLR(sum);
double ans=0;
int m=mp.size();
rep(i,0,v.size()-2) {
// lookln(v[i].h);
int l=lower_bound(mp.begin(),mp.end(),v[i].l)-mp.begin();
int r=lower_bound(mp.begin(),mp.end(),v[i].r)-mp.begin();
l++,r++;
// cout << l << ' ' << r << endl;
update(l,r-1,v[i].val,1,m,1);
// lookln(sum[1]);
ans+=sum[1]*(v[i+1].h-v[i].h);
}
printf("Test case #%d\n",kase++);
printf("Total explored area: %.2f\n\n", ans);
}
return 0;
} HDU1225
题意：
和上一题有区别，要求的是覆盖两次以上的面积和，所以需要开两个数组，分别
表示覆盖至少一次和至少两次。在更新是要注意，因为没有pushdown这个操作
所以父亲节点的信息没有传下来，比如你分别更新了[1,5],[1,3],[4,5]，
这个时候，你在这三个区间都有覆盖一次的标记，但是父亲区间实际是被覆盖
两次的，所以更新是要注意。代码：
//
// Created by CQU_CST_WuErli
// Copyright (c) 2016 CQU_CST_WuErli. All rights reserved.
//
// #include
#include
#include vector v;
vector<double> mp;
void pushup(int l,int r,int rt) {
if (cnt[rt]>=2) two[rt]=one[rt]=mp[r]-mp[l-1];
else if (cnt[rt]==1) {
one[rt]=mp[r]-mp[l-1];
if (l==r) two[rt]=0;
else two[rt]=one[rt<<1]+one[rt<<1|1];
}
else {
if (l==r) one[rt]=two[rt]=0;
else {
one[rt]=one[rt<<1]+one[rt<<1|1];
two[rt]=two[rt<<1]+two[rt<<1|1];
}
}
}
void update(int L,int R,int val,int l,int r,int rt) {
if (L<=l && r<=R) {
cnt[rt]+=val;
pushup(l,r,rt);
return;
}
int mid=(l+r)>>1;
if (L<=mid) update(L,R,val,Lson);
if (R>mid) update(L,R,val,Rson);
pushup(l,r,rt);
}
int main(int argc, char const *argv[]) {
#ifdef LOCAL
freopen("C:\\Users\\john\\Desktop\\in.txt","r",stdin);
// freopen("C:\\Users\\john\\Desktop\\out.txt","w",stdout);
#endif
SI(T);
while (T--) {
// BUG;
double X1,X2,Y1,Y2;
SI(n);
mp.clear();
v.clear();
rep(i,1,n) {
scanf("%lf%lf%lf%lf",&X1,&Y1,&X2,&Y2);
v.push_back(Seg(X1,X2,Y1,1));
v.push_back(Seg(X1,X2,Y2,-1));
mp.push_back(X1);mp.push_back(X2);
}
sort(v.begin(),v.end());
sort(mp.begin(),mp.end());
mp.resize(unique(mp.begin(),mp.end())-mp.begin());
CLR(cnt);
CLR(one);CLR(two);
int m=mp.size();
double ans=0.0;
// BUG;
rep(i,0,v.size()-2) {
int l=lower_bound(mp.begin(),mp.end(),v[i].l)-mp.begin();
// BUG;
int r=lower_bound(mp.begin(),mp.end(),v[i].r)-mp.begin();
l++,r++;
if (l1,v[i].val,1,m,1);
// BUG;
ans+=two[1]*(v[i+1].h-v[i].h);
}
printf("%.2f\n",ans+eps);
}
return 0;
}