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<bits/stdc++.h>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cctype>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <map>
#include <queue>
#include <stack>
#include <set>
#include <algorithm>
#include <sstream>
#define CLR(x) memset(x,0,sizeof(x))
#define OFF(x) memset(x,-1,sizeof(x))
#define MEM(x,a) memset((x),(a),sizeof(x))
#define For_UVa if (kase!=1) cout << endl
#define BUG cout << "I am here" << endl
#define lookln(x) cout << #x << "=" << x << endl
#define SI(a) scanf("%d",&a)
#define SII(a,b) scanf("%d%d",&a,&b)
#define SIII(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define rep(flag,start,end) for(int flag=start;flag<=end;flag++)
#define Rep(flag,start,end) for(int flag=start;flag>=end;flag--)
#define Lson l,mid,rt<<1
#define Rson mid+1,r,rt<<1|1
#define Root 1,n,1
#define BigInteger bign
const int MAX_L=2005;// For BigInteger
const int INF_INT=0x3f3f3f3f;
const long long INF_LL=0x7fffffff;
const int MOD=1e9+7;
const double eps=1e-9;
const double pi=acos(-1);
typedef long long  ll;
using namespace std;

const int N=220;
int num[N<<2];
double sum[N<<2];
struct Seg {
    double l,r,h;
    int val;
    Seg(){}
    Seg(double l,double r,double h,int val):l(l),r(r),h(h),val(val){}
    bool operator < (const Seg& rhs) const {
        return h<rhs.h;
    }
};
vector<Seg> 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<bits/stdc++.h>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cctype>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <map>
#include <queue>
#include <stack>
#include <set>
#include <algorithm>
#include <sstream>
#define CLR(x) memset(x,0,sizeof(x))
#define OFF(x) memset(x,-1,sizeof(x))
#define MEM(x,a) memset((x),(a),sizeof(x))
#define For_UVa if (kase!=1) cout << endl
#define BUG cout << "I am here" << endl
#define lookln(x) cout << #x << "=" << x << endl
#define SI(a) scanf("%d",&a)
#define SII(a,b) scanf("%d%d",&a,&b)
#define SIII(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define rep(flag,start,end) for(int flag=start;flag<=end;flag++)
#define Rep(flag,start,end) for(int flag=start;flag>=end;flag--)
#define Lson l,mid,rt<<1
#define Rson mid+1,r,rt<<1|1
#define Root 1,n,1
#define BigInteger bign
const int MAX_L=2005;// For BigInteger
const int INF_INT=0x3f3f3f3f;
const long long INF_LL=0x7fffffff;
const int MOD=1e9+7;
const double eps=1e-9;
const double pi=acos(-1);
typedef long long  ll;
using namespace std;

const int N=2e4;

int T;
int n;
int cnt[N<<2];
double one[N<<2],two[N<<2];

struct Seg {
    double l,r,h;
    int val;
    Seg(){}
    Seg(double l,double r,double h,int val):l(l),r(r),h(h),val(val){}
    bool operator < (const Seg& rhs) const {
        return h<rhs.h;
    }
};
vector<Seg> 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 (l<r) update(l,r-1,v[i].val,1,m,1);
            // BUG;
            ans+=two[1]*(v[i+1].h-v[i].h);
        }
        printf("%.2f\n",ans+eps);

    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值