码迷,mamicode.com
首页 > 其他好文 > 详细

COJ969 WZJ的数据结构(负三十一)

时间:2015-08-03 12:57:45      阅读:207      评论:0      收藏:0      [点我收藏+]

标签:

 

WZJ的数据结构(负三十一)
难度级别:D; 运行时间限制:3000ms; 运行空间限制:262144KB; 代码长度限制:2000000B
试题描述
A国有两个主基站,供给全国的资源。定义一个主基站能覆盖的范围为:以该主基站为圆心,半径为r的圆(包括边界)。
如果一个子基站能被一个主基站覆盖,则它是激活的。
有N个事件,事件分两种:
1.新建一个坐标位于(x,y)的子基站。 
2.给出两个主基站的半径:r1、r2。询问处于非激活状态的子基站个数。
输入
第一行是四个正整数:x1、y1、x2、y2。表示两个主基站的坐标是(x1,y1)和(x2,y2)。
第二行是一个正整数N,表示有N个事件。
接下来的N行,每行三个正整数。
如果第一个数为1,则接下来两个数为x,y,表示新建一个坐标为(x,y)的新基站。
如果第一个数为2,则接下来两个数为r1,r2,表示询问当两个主基站的覆盖半径为r1和r2时,处于非激活状态的子基站数。
输出
对于每个询问输出答案。
输入示例
1 10 5 2
10
1 2 6
1 1 9
1 3 8
1 6 7
1 4 12
2 1 1
2 3 2
2 8 2
2 2 2
2 3 2
输出示例
4
3
0
4
3
其他说明
1<=x1,y1,x2,y2,x,y,r1,r2<=10^9
1<=N<=200000

妈妈我会写线段树分治了!

按论文的方法,操作对应一个时间区间,将这个区间打上标记,那么操作的标记数是O(nlogn)的,询问对应一个点,将点所在的区间打上标记,那么询问的标记数也是O(nlogn)的。最后对每一条线段离线做就可以了。

技术分享
#include<cstdio>
#include<cctype>
#include<queue>
#include<cstring>
#include<algorithm>
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define dwn(i,s,t) for(int i=s;i>=t;i--)
#define ren for(int i=first[x];i!=-1;i=next[i])
using namespace std;
inline int read() {
    int x=0,f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c==-) f=-1;
    for(;isdigit(c);c=getchar()) x=x*10+c-0;
    return x*f;
}
typedef long long ll;
const int maxn=200010;
ll x0,y0,x1,y1,tmp[maxn];
int n,N,first1[maxn*3],next1[maxn*15],id1[maxn*15],ToT1;
int first2[maxn*3],next2[maxn*15],id2[maxn*15],ToT2;
struct Query {
    int tp,id;
    ll x,y;
}Q[maxn];
struct Point {
    ll x,y;int id;
    bool operator < (const Point& a) const {
        if(x!=a.x) return x>a.x;
        if(y!=a.y) return y>a.y;
    }
}A[maxn];
void pre() {
    rep(i,1,n) tmp[i]=-Q[i].y;
    sort(tmp+1,tmp+n+1);
    rep(i,1,n) Q[i].y=lower_bound(tmp+1,tmp+n+1,-Q[i].y)-tmp;
}
int ans[maxn],sumv[maxn];
void add(int x,int v) {for(;x<=n;x+=x&-x) sumv[x]+=v;}
int sum(int x) {int res=0;for(;x;x-=x&-x) res+=sumv[x];return res;}
void query1(int o,int l,int r,int ql,int qr,int v) {
    if(ql<=l&&r<=qr) {
        id1[++ToT1]=v;next1[ToT1]=first1[o];first1[o]=ToT1;
    }
    else {
        int mid=l+r>>1,lc=o<<1,rc=lc|1;
        if(ql<=mid) query1(lc,l,mid,ql,qr,v);
        if(qr>mid) query1(rc,mid+1,r,ql,qr,v);
    }
} 
void query2(int o,int l,int r,int x,int v) {
    id2[++ToT2]=v;next2[ToT2]=first2[o];first2[o]=ToT2;
    if(l==r) return;
    int mid=l+r>>1,lc=o<<1,rc=lc|1;
    if(x<=mid) query2(lc,l,mid,x,v);
    else query2(rc,mid+1,r,x,v);
}
int main() {
    x0=read();y0=read();x1=read();y1=read();
    n=read();
    rep(i,1,n) {
        Q[Q[i].id=i].tp=read();
        ll x=read(),y=read();
        if(Q[i].tp==1) {
            Q[i].x=(x-x0)*(x-x0)+(y-y0)*(y-y0),Q[i].y=(x-x1)*(x-x1)+(y-y1)*(y-y1);
            query1(1,1,n,i,n,i);
        }
        else {
            Q[i].x=x*x,Q[i].y=y*y;
            query2(1,1,n,i,i);
        }
    }
    pre();
    rep(i,1,3*n) {
        int m=0;
        for(int j=first2[i];j;j=next2[j]) A[++m]=(Point){Q[id2[j]].x,Q[id2[j]].y,Q[id2[j]].id};
        for(int j=first1[i];j;j=next1[j]) A[++m]=(Point){Q[id1[j]].x,Q[id1[j]].y,0};
        sort(A+1,A+m+1);
        rep(j,1,m) if(!A[j].id) add(A[j].y,1); else ans[A[j].id]+=sum(A[j].y-1);
        rep(j,1,m) if(!A[j].id) add(A[j].y,-1);
    }
    rep(i,1,n) if(Q[i].tp==2) printf("%d\n",ans[i]);
    return 0;
}
View Code

CDQ分治常数比较小:

技术分享
#include<cstdio>
#include<cctype>
#include<queue>
#include<cstring>
#include<algorithm>
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define dwn(i,s,t) for(int i=s;i>=t;i--)
#define ren for(int i=first[x];i!=-1;i=next[i])
using namespace std;
inline int read() {
    int x=0,f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c==-) f=-1;
    for(;isdigit(c);c=getchar()) x=x*10+c-0;
    return x*f;
}
typedef long long ll;
const int maxn=200010;
ll x0,y0,x1,y1,tmp[maxn];
struct Query {
    int tp,id;
    ll x,y;
}Q[maxn];
struct Point {
    ll x,y;int id;
    bool operator < (const Point& a) const {
        if(x!=a.x) return x>a.x;
        if(y!=a.y) return y>a.y;
    }
}A[maxn];
int n,ans[maxn],sumv[maxn];
void add(int x,int v) {for(;x<=n;x+=x&-x) sumv[x]+=v;}
int sum(int x) {int res=0;for(;x;x-=x&-x) res+=sumv[x];return res;}
void solve(int L,int R) {
    if(L>=R) return;
    int mid=L+R>>1,m0=0,m=0;
    solve(L,mid);solve(mid+1,R);
    rep(i,L,mid) if(Q[i].tp==1) A[++m0]=(Point){Q[i].x,Q[i].y,0};
    if(!m0) return;m=m0;
    rep(i,mid+1,R) if(Q[i].tp==2) A[++m0]=(Point){Q[i].x,Q[i].y,Q[i].id};
    if(m0==m) return;m=m0;
    sort(A+1,A+m+1);
    rep(i,1,m) if(!A[i].id) add(A[i].y,1); else ans[A[i].id]+=sum(A[i].y-1);
    rep(i,1,m) if(!A[i].id) add(A[i].y,-1);
}
void pre() {
    rep(i,1,n) tmp[i]=-Q[i].y;
    sort(tmp+1,tmp+n+1);
    rep(i,1,n) Q[i].y=lower_bound(tmp+1,tmp+n+1,-Q[i].y)-tmp;
}
int main() {
    x0=read();y0=read();x1=read();y1=read();
    n=read();
    rep(i,1,n) {
        Q[Q[i].id=i].tp=read();
        ll x=read(),y=read();
        if(Q[i].tp==1) Q[i].x=(x-x0)*(x-x0)+(y-y0)*(y-y0),Q[i].y=(x-x1)*(x-x1)+(y-y1)*(y-y1);
        else Q[i].x=x*x,Q[i].y=y*y;
    }
    pre();solve(1,n);
    rep(i,1,n) if(Q[i].tp==2) printf("%d\n",ans[i]);
    return 0;
}
View Code

 

COJ969 WZJ的数据结构(负三十一)

标签:

原文地址:http://www.cnblogs.com/wzj-is-a-juruo/p/4699030.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!