hdu5862Counting Intersections 扫描线+线段树

题目大意:给你n个平行于坐标轴的线段,求其交点数,n<=100000

题解:首先将线段离散化。 然后将所有线段按照横坐标的顺序,按照先插入再查找再删除的顺序进行操作。

对于横线 如果是左端点,则将其纵坐标加1,右端点则减1,对于竖线直接求和就行了。(线段树点要开多

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <map>
using namespace std;
#define maxn 1000100
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define getmid int mid=(l+r)>>1
int sum[maxn],save[maxn];
struct node
{
    int tp,x,y,v;
}po[maxn];
node get(int a,int b,int c,int d)
{
    node tmp;
    tmp.tp=a; tmp.x=b; tmp.y=c; tmp.v=d;
    return tmp;
}
bool cmp(node a,node b)
{
    if(a.x==b.x)return a.tp<b.tp;
    return a.x<b.x;
}
int ysiz,posiz,mpsiz;
map<int,int>mp;
int query(int l,int r,int rt,int st,int ed)
{
    if(st<=l&&r<=ed)
    {
        return sum[rt];
    }
    getmid;
    int res1=0,res2=0;
    if(st<=mid) res1=query(lson,st,ed);
    if(ed>mid) res2=query(rson,st,ed);
    return res1+res2;
}
void update(int l,int r,int rt,int pos,int v)
{
    if(l==r)
    {
        sum[rt]+=v;
        return ;
    }
    getmid;
    if(pos<=mid) update(lson,pos,v);
    else update(rson,pos,v);
    sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
int n,m;
int main()
{
    int cas;
    scanf("%d",&cas);
    while(cas--)
    {
        mp.clear();
        int x1,x2,y1,y2;
        scanf("%d",&n);
        ysiz=0,posiz=0,mpsiz=0;
        for(int i=1;i<=n;i++)
        {
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            if(x1==x2)
            {
                if(y1>y2) swap(y1,y2);
                posiz++; po[posiz]={1,x1,y1,y2};

                ysiz++;save[ysiz]=y1;
                ysiz++;save[ysiz]=y2;
            }
            else
            {
                if(x1>x2)swap(x1,x2);
                posiz++;po[posiz]={0,x1,y1,1};
                posiz++;po[posiz]={0,x2+1,y2,-1};
                
                ysiz++;save[ysiz]=y1;
            }
        }
        long long ans=0;
        sort(save+1,save+1+ysiz);
        for(int i=1;i<=ysiz;i++)
        {
            if(mp[save[i]]) continue;
            
            mpsiz++;mp[save[i]]=mpsiz;
        }
        sort(po+1,po+1+posiz,cmp);
        memset(sum,0,sizeof(sum));
        for(int i=1;i<=posiz;i++)
        {
            if(po[i].tp==0)
            {
                int tmp=mp[po[i].y];
                update(1,mpsiz,1,tmp,po[i].v);
            }
            else
            {
                int xl=mp[po[i].y];
                int xr=mp[po[i].v];
                ans+=query(1,mpsiz,1,xl,xr);
            }
        }
        printf("%lld\n",ans);
    }
}

一点)


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值