poj1177 Picture

最新推荐文章于 2018-09-21 22:54:15 发布

BinJourney

最新推荐文章于 2018-09-21 22:54:15 发布

阅读量548

点赞数

分类专栏： ACM 文章标签：计算几何扫描线坐标离散化线段树

本文链接：https://blog.csdn.net/lj94093/article/details/43876217

版权

ACM 专栏收录该内容

137 篇文章

订阅专栏

本文介绍了一种使用扫描线算法结合线段树的数据结构来解决几何问题的方法。通过具体的代码实现，展示了如何处理线段的交叠和计算被线段覆盖区域的长度。适用于需要高效处理大量线段覆盖问题的场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

基本上照抄的，学习一下扫描线的线段树

#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<math.h>
using namespace std;
typedef long long ll;
struct Line{
    int x,y1,y2,flag;
};
/*********************************

*   l-左子树节点 r-右子树节点
*   num-子树中有多少条线段
*   len-当前结点覆盖y方向的总长度
*   cover-是否被覆盖 也就是入边+1，出边-1
*   lb-左端点是否被覆盖 rb-右端点是否被覆盖
**********************************/
//每个节点表示的是该节点l对应的y值和r对应的y值之间的部分
struct Node {
    int l,r,num,len,cover;
    bool lb,rb;
};
Node node[20010];
Line l[10010];
int y[10010];
void add_line(int x1,int x2,int y1,int y2,int &cnt){
    //入边
    l[cnt].x=x1;
    l[cnt].y1=y1;
    l[cnt].y2=y2;
    l[cnt].flag=1;
    y[cnt++]=y1;
    //出边
    l[cnt].x=x2;
    l[cnt].y1=y1;
    l[cnt].y2=y2;
    l[cnt].flag=-1;
    y[cnt++]=y2;
}
//入边在前
int cmp(Line a,Line b){
    if (a.x==b.x) return a.flag>b.flag;
    return a.x<b.x;
}
void Build(int t,int l,int r){
    node[t].l=l;
    node[t].r=r;
    node[t].num=0;
    if (l==r-1) return;
    int mid=(r+l)>>1;
    Build(t*2,l,mid);
    Build(t*2+1,mid,r);
}
void Update_len(int t){
    //>0即该节点对应扫描线对应的两个y值被完全是覆盖了
    if (node[t].cover>0) {
        node[t].num=node[t].lb=node[t].rb=1;
        node[t].len=y[node[t].r]-y[node[t].l];
    }
    //完全没被覆盖
    else if(node[t].l==node[t].r-1) node[t].rb=node[t].lb=node[t].len=node[t].num=0;
    //部分覆盖
    else {
        node[t].rb=node[2*t+1].rb;
        node[t].lb=node[2*t].lb;
        node[t].len=node[2*t].len+node[2*t+1].len;
        node[t].num=node[2*t].num+node[2*t+1].num-node[2*t+1].lb*node[2*t].rb;
    }
}
void Update(int t,Line p){
    if (y[node[t].l]>=p.y1 &&y[node[t].r]<=p.y2 ){
        node[t].cover+=p.flag;
        Update_len(t);
        return;
    }
    if (node[t].l==node[t].r-1) return;
    int mid=(node[t].l+node[t].r)>>1;
    if (p.y1<=y[mid]) Update(2*t,p);//中线下方
    if (p.y2>y[mid]) Update(2*t+1,p);//中线上方
    Update_len(t);
}
int solve(int n,int cnt,Line *l){
    memset(node,0,sizeof(node));
    Build(1,0,cnt-1);
    int ans=0,lines=0,last=0;
    for(int i=0;i<n;i++){
        Update(1,l[i]);
        if (i>0) ans+=2*lines*(l[i].x-l[i-1].x);//平等于x轴的长度
        ans+=abs(node[1].len-last);
        last=node[1].len;
        lines=node[1].num;
    }
    return ans;
}
int main() {
    #ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
    #endif // ONLINE_JDUGE
    int n,x1,x2,y1,y2,cnt=0;
    scanf("%d",&n);
    for(int i=0;i<n;i++){
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        add_line(x1,x2,y1,y2,cnt);
    }

    sort(y,y+cnt);
    sort(l,l+cnt,cmp);
    int t=cnt;
    t=unique(y,y+cnt)-y;
    int ans=solve(cnt,t,l);
    printf("%d\n",ans);
    return 0;
}