线段树扫描线模板

来源：互联网发布：阿里云售后自助服务编辑：程序博客网时间：2024/06/06 02:40
#include <cstdio>//hdu1542 求平面矩形面积交#include <algorithm>#include <vector>#define lson l,m,rt<<1#define rson m+1,r,rt<<1|1const int N=1222;using namespace std;int n,cnt;struct node{    double l,r,h;    int f;    node(double _l=0,double _r=0,double _h=0,int _f=0):l(_l),r(_r),h(_h),f(_f){};}a[N];bool cmp(node a,node b){return a.h<b.h;}double sum[N<<2];int flag[N<<2];vector<double> vt;int Getid(double x) {return lower_bound(vt.begin(),vt.end(),x)-vt.begin()+1;}void Pushup(int l,int r,int rt){    if(flag[rt]) sum[rt]=vt[r]-vt[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 x,int l,int r,int rt){    if(L<=l&&r<=R)    {        flag[rt]+=x;        Pushup(l,r,rt);        return;    }    int m=l+r>>1;    if(L<=m) Update(L,R,x,lson);    if(m<R) Update(L,R,x,rson);    Pushup(l,r,rt);}int main(){    int i,cas=0;    while(~scanf("%d",&n)&&n)    {        vt.clear();        cnt=0;        for(i=1;i<=n;i++)        {            double x1,x2,y1,y2;            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);            vt.push_back(x1);vt.push_back(x2);            a[++cnt]=node(x1,x2,y1,1);            a[++cnt]=node(x1,x2,y2,-1);        }        sort(vt.begin(),vt.end());        vt.erase(unique(vt.begin(),vt.end()),vt.end());        sort(a+1,a+1+cnt,cmp);        double ans=0;        for(i=1;i<=cnt;i++)        {            int l=Getid(a[i].l);            int r=Getid(a[i].r)-1;            Update(l,r,a[i].f,1,(int)vt.size(),1);            ans+=sum[1]*(a[i+1].h-a[i].h);        }        printf("Test case #%d\n",++cas);        printf("Total explored area: %.2f\n\n",ans);    }}/*********************************************/#include <cstdio>//hdu1828求矩形周长交#include <algorithm>#include <vector>#define lson l,m,rt<<1#define rson m+1,r,rt<<1|1const int N=100010;using namespace std;int n,cnt,sum[N<<2],flag[N<<2],num[N<<2],nl[N<<2],nr[N<<2];struct node{    int l,r,h,f;    node(int _l=0,int _r=0,int _h=0,int _f=0):l(_l),r(_r),h(_h),f(_f){}}a[N];bool Cmp(node a,node b){    if(a.h!=b.h) return a.h<b.h;    return a.f>b.f;}vector<int> vt;int Getid(int x) {return lower_bound(vt.begin(),vt.end(),x)-vt.begin()+1;}void pushup(int l,int r,int rt){    if(flag[rt])    {        sum[rt]=vt[r]-vt[l-1];        num[rt]=nl[rt]=nr[rt]=1;    }    else if(l==r) sum[rt]=num[rt]=nl[rt]=nr[rt]=0;    else    {        sum[rt]=sum[rt<<1]+sum[rt<<1|1];        nl[rt]=nl[rt<<1];        nr[rt]=nr[rt<<1|1];        num[rt]=num[rt<<1]+num[rt<<1|1]-nl[rt<<1|1]*nr[rt<<1];    }}void Update(int L,int R,int x,int l,int r,int rt){    if(L<=l&&r<=R)    {        flag[rt]+=x;        Pushup(l,r,rt);        return;    }    int m=l+r>>1;    if(L<=m) Update(L,R,x,lson);    if(m<R) Update(L,R,x,rson);    Pushup(l,r,rt);}int main(){    int i;    while(~scanf("%d",&n))    {        vt.clear();        cnt=0;        int x1,x2,y1,y2;        for(i=1;i<=n;i++)        {            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);            vt.push_back(x1),vt.push_back(x2);            a[++cnt]=node(x1,x2,y1,1);            a[++cnt]=node(x1,x2,y2,-1);        }        sort(vt.begin(),vt.end());        vt.erase(unique(vt.begin(),vt.end()),vt.end());        sort(a+1,a+1+cnt,Cmp);        int ans=0,last=0;        for(i=1;i<=cnt;i++)        {            int l=Getid(a[i].l);            int r=Getid(a[i].r)-1;            Update(l,r,a[i].f,1,(int)vt.size(),1);            ans+=abs(sum[1]-last)+num[1]*(a[i+1].h-a[i].h)*2;            last=sum[1];        }        printf("%d\n",ans);    }}
0 0