【POI2008】STR

最新推荐文章于 2024-01-25 09:23:53 发布
最新推荐文章于 2024-01-25 09:23:53 发布
阅读量649 收藏
点赞数
分类专栏： 数据结构 解题报告 poi 文章标签： poi 解题报告
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/samjia2000/article/details/47427989
版权
解题报告 同时被 3 个专栏收录
69 篇文章 0 订阅
订阅专栏
数据结构
22 篇文章 0 订阅
订阅专栏
poi
1 篇文章 0 订阅
订阅专栏

题目大意
给出一个平面内的n个点，有一系列询问形如: x0,y0,x1,y1 ，输出平面的点中更接近 (x0,y0) 的个数，更接近 (x1,y1) 的个数，与两点距离相等的点的个数。

主要思想
可以按照询问分6种情况讨论!
情况一
情况二
情况三
情况四
情况五
情况六
以上六种情况中， p1 表示第一个点， p2 表示第二个点，他们之间的关系有六种，前三种分别为，他们组成的长方形的长与纵轴平行，他们组成的长方形的长与横轴平行，他们组成了一个正方形，后三种类似。每个图被分为三个部分：
1. 实线与阴影部分
2. 实线与阴影部分外偏 p1 的部分
3. 实线与阴影部分外偏 p2 的部分
对于每一种情况，将分出的部分再分成几块来算，用某个数据结构（比如树状数组）来维护就好了。
下面附一下代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<numeric>
#include<queue>
#include<functional>
#include<set>
#include<map>

#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define MAXN 100010

typedef double db;

using namespace std;

int get(){
    char ch;
    while (ch=getchar(),(ch<'0'||ch>'9')&& ch!='-');
    char c=ch;
    int s=(c!='-')*c-48;
    while (ch=getchar(),ch>='0'&&ch<='9')s=s*10+ch-48;
    return c=='-'?-s:s;
}

int n,m,z,p;
struct point{
    int x,y;
    db v;
}a[MAXN],b[MAXN];
int hx[MAXN],kx,hy[MAXN],ky,tree[MAXN],c[MAXN];
struct question{
    db x1,x2,y1,y2;
    int ty,ans1,ans2,ans3;
    db v1,v2;
    bool bz;
}q[MAXN];

void gettype(int x1,int y1,int x2,int y2,int i){
    q[i].bz=0;
    if(y1<y2){
        x1^=x2^=x1^=x2;
        y1^=y2^=y1^=y2;
        q[i].bz=1;
    }
    q[i].x1=x1;q[i].y1=y1;q[i].x2=x2;q[i].y2=y2;
    if (x1<x2){
        if (y1-y2>x2-x1)q[i].ty=1;
        if (y1-y2<x2-x1)q[i].ty=2;
        if (y1-y2==x2-x1)q[i].ty=3;
    }
    else{
        if (y1-y2>x1-x2)q[i].ty=4;
        if (y1-y2<x1-x2)q[i].ty=5;
        if (y1-y2==x1-x2)q[i].ty=6;
    }
}

bool cmp1(point i,point j){
    if (i.y!=j.y)return i.y<j.y;
    return i.x<j.x;
}

bool cmp2(point i,point j){
    if (i.x!=j.x)return i.x<j.x;
    return i.y<j.y;
}

bool cmp3(int i,int j){
    if (q[i].v1!=q[j].v1)return q[i].v1<q[j].v1;
    return q[i].v2<q[j].v2;
}

bool cmp4(point i,point j){
    if (i.v!=j.v)return i.v<j.v;
    return i.x<j.x;
}

void add(int x){
    while (x<=z){
        tree[x]++;
        x+=x & -x;
    }
}

int getsum(int x){
    int ans(0);
    while (x){
        ans+=tree[x];
        x=x-(x&-x);
    }
    return ans;
}

int askx(db x){
    int l=1,r=kx,ans(0);
    while (l<=r){
        int mid=(l+r)/2;
        if (hx[mid]>x)r=mid-1;
        else{
            l=mid+1;
            ans=mid;
        }
    }
    return ans;
}

int asky(db x){
    int l=1,r=ky,ans(0);
    while (l<=r){
        int mid=(l+r)/2;
        if (hy[mid]>x)r=mid-1;
        else{
            l=mid+1;
            ans=mid;
        }
    }
    return ans;
}

void type1(){
    int h1=1,tot=0;
    fo(i,1,p)
    if (q[i].ty==1)c[++tot]=i;
    sort(a+1,a+1+z,cmp1);
    ky=0;
    fo(i,1,z)
    if (a[i].y>a[i-1].y)hy[++ky]=a[i].y;
    sort(a+1,a+1+z,cmp2);
    kx=0;
    fo(i,1,z)
    if (a[i].x>a[i-1].x)hx[++kx]=a[i].x;
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].x1;
        q[x].v2=q[x].y1-db(q[x].x2-q[x].x1+q[x].y1-q[x].y2)/2;
    }
    sort(c+1,c+1+tot,cmp3);
    fo(i,1,z){
        while (h1<=tot&&q[c[h1]].v1<=a[i].x){
            int x=c[h1++];
            int t1=getsum(asky(q[x].v2)),t2=getsum(asky(q[x].v2-1)),t3=getsum(ky);
            q[x].ans2+=t2;
            q[x].ans1+=t3-t1;
            if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
            else q[x].ans2+=t1-t2;
        }
        add(asky(a[i].y));
    }
    while(h1<=tot){
        int x=c[h1++];
        int t1=getsum(asky(q[x].v2)),t2=getsum(asky(q[x].v2-1)),t3=getsum(ky);
        q[x].ans2+=t2;
        q[x].ans1+=t3-t1;
        if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
        else q[x].ans2+=t1-t2;
    }

    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].x2;
        q[x].v2=q[x].y2+db(q[x].x2-q[x].x1+q[x].y1-q[x].y2)/2;
    }
    sort(c+1,c+1+tot,cmp3);
    fo(i,1,z)tree[i]=0;
    h1=tot;
    fd(i,z,1){
        while (h1&&q[c[h1]].v1>=a[i].x){
            int x=c[h1--];
            int t1=getsum(asky(q[x].v2)),t2=getsum(asky(q[x].v2-1)),t3=getsum(ky);
            q[x].ans2+=t2;
            q[x].ans1+=t3-t1;
            if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
            else q[x].ans2+=t1-t2;
        }
        add(asky(a[i].y));
    }
    while(h1){
        int x=c[h1--];
        int t1=getsum(asky(q[x].v2)),t2=getsum(asky(q[x].v2-1)),t3=getsum(ky);
        q[x].ans2+=t2;
        q[x].ans1+=t3-t1;
        if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
        else q[x].ans2+=t1-t2;
    }

    fo(i,1,z)a[i].v=a[i].y-a[i].x;
    sort(a+1,a+1+z,cmp4);
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].v2-q[x].v1;
        q[x].v2=0;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=1;
    fo(i,1,z)tree[i]=0;
    fo(i,1,z){
        while (h1<=tot&&q[c[h1]].v1<=a[i].v){
            int x=c[h1++];
            if (q[x].x2!=q[x].x1)
            q[x].ans2+=getsum(askx(q[x].x2))-getsum(askx(q[x].x1-1));
        }
        add(askx(a[i].x));
    }
    while (h1<=tot){
        int x=c[h1++];
        if (q[x].x2!=q[x].x1)
        q[x].ans2+=getsum(askx(q[x].x2))-getsum(askx(q[x].x1-1));
    }

    fo(i,1,z)tree[i]=0;
    h1=tot;
    fd(i,z,1){
        while (h1&&q[c[h1]].v1>=a[i].v){
            int x=c[h1--];
            q[x].ans1+=getsum(askx(q[x].x2))-getsum(askx(q[x].x1-1));
        }
        add(askx(a[i].x));
    }
    while (h1){
        int x=c[h1--];
        if (q[x].x2!=q[x].x1)
        q[x].ans1+=getsum(askx(q[x].x2))-getsum(askx(q[x].x1-1));
    }

    h1=1;
    int l;
    a[0].v=a[1].v-1;
    a[z+1].v=a[z].v+1;
    fo(i,1,z+1){
        if (a[i].v>a[i-1].v){
            while(h1<=tot&&q[c[h1]].v1<a[i-1].v)h1++;
            while(h1<=tot&&q[c[h1]].v1==a[i-1].v){
                int ll=l,rr=i-1,a1(-1),a2(-1),x=c[h1++];
                while (ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x<q[x].x1)ll=mid+1;
                    else{
                        rr=mid-1;
                        a1=mid;
                    }
                }
                ll=l;
                rr=i-1;
                while(ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x>q[x].x2)rr=mid-1;
                    else{
                        ll=mid+1;
                        a2=mid;
                    }
                }
                if (a1==-1||a2==-1)continue;
                q[x].ans3+=a2-a1+1;
            }
            l=i;
        }
    }
}

void type2(){
    int h1,tot(0);
    fo(i,1,p)
    if (q[i].ty==2)c[++tot]=i;

    sort(a+1,a+1+z,cmp1);
    fo(i,1,tot){
        int x=c[i];
        q[x].v2=q[x].x2-db(q[x].x2-q[x].x1+q[x].y1-q[x].y2)/2;
        q[x].v1=q[x].y2;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=1;
    fo(i,1,z)tree[i]=0;
    fo(i,1,z){
        while(h1<=tot&&q[c[h1]].v1<=a[i].y){
            int x=c[h1++];
            int t1=getsum(askx(q[x].v2)),t2=getsum(askx(q[x].v2-1)),t3=i-1;
            int l1=q[x].ans1,l2=q[x].ans2,l3=q[x].ans3;
            q[x].ans1+=t2;
            q[x].ans2+=t3-t1;
            if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
            else q[x].ans1+=t1-t2;
        }
        add(askx(a[i].x));
    }
    while(h1<=tot){
        int x=c[h1++];
        int t1=getsum(askx(q[x].v2)),t2=getsum(askx(q[x].v2-1)),t3=z;
            int l1=q[x].ans1,l2=q[x].ans2,l3=q[x].ans3;
        q[x].ans1+=t2;
        q[x].ans2+=t3-t1;
        if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
        else q[x].ans1+=t1-t2;
    }

    fo(i,1,z)tree[i]=0;
    fo(i,1,tot){
        int x=c[i];
        q[x].v2=q[x].x1+db(q[x].x2-q[x].x1+q[x].y1-q[x].y2)/2;
        q[x].v1=q[x].y1;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=tot;
    fd(i,z,1){
        while(h1&&q[c[h1]].v1>=a[i].y){
            int x=c[h1--];
            int t1=getsum(askx(q[x].v2)),t2=getsum(askx(q[x].v2-1)),t3=z-i;
            int l1=q[x].ans1,l2=q[x].ans2,l3=q[x].ans3;
            q[x].ans2+=t3-t1;
            q[x].ans1+=t2;
            if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
            else q[x].ans1+=t1-t2;
        }
        add(askx(a[i].x));
    } 
    while(h1){
        int x=c[h1--];
        int t1=getsum(askx(q[x].v2)),t2=getsum(askx(q[x].v2-1)),t3=z;
            int l1=q[x].ans1,l2=q[x].ans2,l3=q[x].ans3;
        q[x].ans2+=t3-t1;
        q[x].ans1+=t2;
        if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
        else q[x].ans1+=t1-t2;
    }

    fo(i,1,z)a[i].v=a[i].y-a[i].x;
    sort(a+1,a+1+z,cmp4);
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].v1-q[x].v2;
        q[x].v2=0;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=1;
    fo(i,1,z)tree[i]=0;
    fo(i,1,z){
        while (h1<=tot&&q[c[h1]].v1<=a[i].v){
            int x=c[h1++];
            q[x].ans2+=getsum(asky(q[x].y1))-getsum(asky(q[x].y2-1));
        }
        add(asky(a[i].y));
    }
    while (h1<=tot){
        int x=c[h1++];
        q[x].ans2+=getsum(asky(q[x].y1))-getsum(asky(q[x].y2-1));
    }

    fo(i,1,z)tree[i]=0;
    h1=tot;
    fd(i,z,1){
        while (h1&&q[c[h1]].v1>=a[i].v){
            int x=c[h1--];
            q[x].ans1+=getsum(asky(q[x].y1))-getsum(asky(q[x].y2-1));
        }
        add(asky(a[i].y));
    }
    while (h1){
        int x=c[h1--];
        q[x].ans1+=getsum(asky(q[x].y1))-getsum(asky(q[x].y2-1));
    }

    h1=1;
    int l=1;
    a[0].v=a[1].v-1;
    a[z+1].v=a[z].v+1;
    fo(i,1,z+1){
        if (a[i].v>a[i-1].v){
            while(h1<=tot&&q[c[h1]].v1<a[i-1].v)h1++;
            while(h1<=tot&&q[c[h1]].v1==a[i-1].v){
                int ll=l,rr=i-1,a1(-1),a2(-1),x=c[h1++];
                db len=(q[x].x2-q[x].x1+q[x].y1-q[x].y2)/2;
                while (ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x<q[x].x2-len)ll=mid+1;
                    else{
                        rr=mid-1;
                        a1=mid;
                    }
                }
                ll=l;
                rr=i-1;
                while(ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x>q[x].x1+len)rr=mid-1;
                    else{
                        ll=mid+1;
                        a2=mid;
                    }
                }
                if (a1==-1||a2==-1)continue;
                q[x].ans3+=a2-a1+1;
            }
            l=i;
        }
    }
}

void type3(){
    int h1,tot(0);
    fo(i,1,p)
    if (q[i].ty==3)c[++tot]=i;

    fo(i,1,z)tree[i]=0;
    sort(a+1,a+1+z,cmp2);
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].x1;
        q[x].v2=q[x].y2;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=1;
    fo(i,1,z){
        while(h1<=tot&&q[c[h1]].v1<=a[i].x){
            int x=c[h1++];
            int t1=getsum(asky(q[x].v2)),t2=getsum(ky);
            q[x].ans3+=t1;
            q[x].ans1+=t2-t1;
        }
        add(asky(a[i].y));
    }
    while(h1<=tot){
        int x=c[h1++];
        int t1=getsum(asky(q[x].v2)),t2=getsum(ky);
        q[x].ans3+=t1;
        q[x].ans1+=t2-t1;
    }

    fo(i,1,z)tree[i]=0;
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].x2;
        q[x].v2=q[x].y1;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=tot;
    fd(i,z,1){
        while(h1&&q[c[h1]].v1>=a[i].x){
            int x=c[h1--];
            int t1=getsum(ky),t2=getsum(asky(q[x].v2-1));
            q[x].ans2+=t2;
            q[x].ans3+=t1-t2;
        }
        add(asky(a[i].y));
    }
    while(h1){
        int x=c[h1--];
        int t1=getsum(ky),t2=getsum(asky(q[x].v2-1));
        q[x].ans2+=t2;
        q[x].ans3+=t1-t2;
    }

    fo(i,1,z)a[i].v=a[i].y-a[i].x;
    sort(a+1,a+1+z,cmp4);
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].v2-q[x].v1;
        q[x].v2=0;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=1;
    fo(i,1,z)tree[i]=0;
    fo(i,1,z){
        while (h1<=tot&&q[c[h1]].v1<=a[i].v){
            int x=c[h1++];
            q[x].ans2+=getsum(askx(q[x].x2))-getsum(askx(q[x].x1));
            q[x].ans3+=getsum(askx(q[x].x1))-getsum(askx(q[x].x1-1));
        }
        add(askx(a[i].x));
    }
    while (h1<=tot){
        int x=c[h1++];
        q[x].ans2+=getsum(askx(q[x].x2))-getsum(askx(q[x].x1));
            q[x].ans3+=getsum(askx(q[x].x1))-getsum(askx(q[x].x1-1));
    }
    fo(i,1,z)tree[i]=0;
    h1=tot;
    fd(i,z,1){
        while (h1&&q[c[h1]].v1>=a[i].v){
            int x=c[h1--];
            q[x].ans1+=getsum(askx(q[x].x2-1))-getsum(askx(q[x].x1-1));
            q[x].ans3+=getsum(askx(q[x].x2))-getsum(askx(q[x].x2-1));
        }
        add(askx(a[i].x));
    }
    while (h1){
        int x=c[h1--];
        q[x].ans1+=getsum(askx(q[x].x2-1))-getsum(askx(q[x].x1-1));
        q[x].ans3+=getsum(askx(q[x].x2))-getsum(askx(q[x].x2-1));
    }

    h1=1;
    int l;
    a[0].v=a[1].v-1;
    a[z+1].v=a[z].v+1;
    fo(i,1,z+1){
        if (a[i].v>a[i-1].v){
            while(h1<=tot&&q[c[h1]].v1<a[i-1].v)h1++;
            while(h1<=tot&&q[c[h1]].v1==a[i-1].v){
                int ll=l,rr=i-1,a1(-1),a2(-1),x=c[h1++];
                while (ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x<q[x].x1)ll=mid+1;
                    else{
                        rr=mid-1;
                        a1=mid;
                    }
                }
                ll=l;
                rr=i-1;
                while(ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x>q[x].x2)rr=mid-1;
                    else{
                        ll=mid+1;
                        a2=mid;
                    }
                }
                if (a1==-1||a2==-1)continue;
                q[x].ans3+=a2-a1+1;
            }
            l=i;
        }
    }
}

void type4(){
    int h1,tot(0);
    fo(i,1,p)
    if (q[i].ty==4)c[++tot]=i;

    sort(a+1,a+1+z,cmp2);
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].x2;
        q[x].v2=q[x].y2+db(q[x].x1-q[x].x2+q[x].y1-q[x].y2)/2;
    }
    fo(i,1,z)tree[i]=0;
    sort(c+1,c+1+tot,cmp3);
    h1=1;
    fo(i,1,z){
        while(h1<=tot&&q[c[h1]].v1<=a[i].x){
            int x=c[h1++];
            int t1=getsum(asky(q[x].v2)),t2=getsum(asky(q[x].v2-1)),t3=getsum(ky);
            q[x].ans1+=t3-t1;
            q[x].ans2+=t2;
            if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
            else q[x].ans2+=t1-t2;
        }
        add(asky(a[i].y));
    }
    while(h1<=tot){
        int x=c[h1++];
        int t1=getsum(asky(q[x].v2)),t2=getsum(asky(q[x].v2-1)),t3=getsum(ky);
        q[x].ans1+=t3-t1;
        q[x].ans2+=t2;
        if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
        else q[x].ans2+=t1-t2;
    }
    fo(i,1,z)tree[i]=0;
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].x1;
        q[x].v2=q[x].y1-db(q[x].x1-q[x].x2+q[x].y1-q[x].y2)/2;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=tot;
    fd(i,z,1){
        while (h1&&q[c[h1]].v1>=a[i].x){
            int x=c[h1--];
            int t1=getsum(asky(q[x].v2)),t2=getsum(asky(q[x].v2-1)),t3=getsum(ky);
            q[x].ans1+=t3-t1;
            q[x].ans2+=t2;
            if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
            else q[x].ans2+=t1-t2;
        }
        add(asky(a[i].y));
    }
    while (h1){
        int x=c[h1--];
        int t1=getsum(asky(q[x].v2)),t2=getsum(asky(q[x].v2-1)),t3=getsum(ky);
        q[x].ans1+=t3-t1;
        q[x].ans2+=t2;
        if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
        else q[x].ans2+=t1-t2;
    }

    fo(i,1,z)tree[i]=0;
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].v1+q[x].v2;
        q[x].v2=0;
    }
    sort(c+1,c+1+tot,cmp3);
    fo(i,1,z)a[i].v=a[i].x+a[i].y;
    sort(a+1,a+1+z,cmp4);
    h1=1;
    fo(i,1,z){
        while (h1<=tot&&q[c[h1]].v1<=a[i].v){
            int x=c[h1++];
            q[x].ans2+=getsum(askx(q[x].x1))-getsum(askx(q[x].x2-1));
        }
        add(askx(a[i].x));
    }
    while (h1<=tot){
        int x=c[h1++];
        q[x].ans2+=getsum(askx(q[x].x1))-getsum(askx(q[x].x2-1));
    }

    fo(i,1,z)tree[i]=0;
    h1=tot;
    fd(i,z,1){
        while (h1&&q[c[h1]].v1>=a[i].v){
            int x=c[h1--];
            q[x].ans1+=getsum(askx(q[x].x1))-getsum(askx(q[x].x2-1));
        }
        add(askx(a[i].x));
    }
    while (h1){
        int x=c[h1--];
        q[x].ans1+=getsum(askx(q[x].x1))-getsum(askx(q[x].x2-1));
    }

    h1=1;
    int l;
    a[0].v=a[1].v-1;
    a[z+1].v=a[z].v+1;
    fo(i,1,z+1){
        if (a[i].v>a[i-1].v){
            while(h1<=tot&&q[c[h1]].v1<a[i-1].v)h1++;
            while(h1<=tot&&q[c[h1]].v1==a[i-1].v){
                int ll=l,rr=i-1,a1(-1),a2(-1),x=c[h1++];
                while (ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x<q[x].x2)ll=mid+1;
                    else{
                        rr=mid-1;
                        a1=mid;
                    }
                }
                ll=l;
                rr=i-1;
                while(ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x>q[x].x1)rr=mid-1;
                    else{
                        ll=mid+1;
                        a2=mid;
                    }
                }
                if (a1==-1||a2==-1)continue;
                q[x].ans3+=a2-a1+1;
            }
            l=i;
        }
    }
}

void type5(){
    int h1,tot(0);
    fo(i,1,p)
    if (q[i].ty==5)c[++tot]=i;

    sort(a+1,a+1+z,cmp1);
    fo(i,1,tot){
        int x=c[i];
        q[x].v2=q[x].x2+db(q[x].x1-q[x].x2+q[x].y1-q[x].y2)/2;
        q[x].v1=q[x].y2;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=1;
    fo(i,1,z)tree[i]=0;
    fo(i,1,z){
        while(h1<=tot&&q[c[h1]].v1<=a[i].y){
            int x=c[h1++];
            int t1=getsum(askx(q[x].v2)),t2=getsum(askx(q[x].v2-1)),t3=i-1;
            q[x].ans2+=t2;
            q[x].ans1+=t3-t1;
            if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
            else q[x].ans2+=t1-t2;
        }
        add(askx(a[i].x));
    }
    while(h1<=tot){
        int x=c[h1++];
        int t1=getsum(askx(q[x].v2)),t2=getsum(askx(q[x].v2-1));
        q[x].ans2+=t2;
        q[x].ans1+=z-t1;
        if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
        else q[x].ans2+=t1-t2;
    }

    fo(i,1,z)tree[i]=0;
    fo(i,1,tot){
        int x=c[i];
        q[x].v2=q[x].x1-db(q[x].x1-q[x].x2+q[x].y1-q[x].y2)/2;
        q[x].v1=q[x].y1;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=tot;
    fd(i,z,1){
        while(h1&&q[c[h1]].v1>=a[i].y){
            int x=c[h1--];
            int t1=getsum(askx(q[x].v2)),t2=getsum(askx(q[x].v2-1)),t3=z-i;
            q[x].ans2+=t2;
            q[x].ans1+=t3-t1;
            if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
            else q[x].ans2+=t1-t2;
        }
        add(askx(a[i].x));
    } 
    while(h1){
        int x=c[h1--];
        int t1=getsum(askx(q[x].v2)),t2=getsum(askx(q[x].v2-1)),t3=z;
        q[x].ans2+=t2;
        q[x].ans1+=z-t1;
        if (int(q[x].v2)==q[x].v2)q[x].ans3+=t1-t2;
        else q[x].ans2+=t1-t2;
    }

    fo(i,1,z)tree[i]=0;
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].v1+q[x].v2;
        q[x].v2=0;
    }
    sort(c+1,c+1+tot,cmp3);
    fo(i,1,z)a[i].v=a[i].x+a[i].y;
    sort(a+1,a+1+z,cmp4);
    h1=1;
    fo(i,1,z){
        while (h1<=tot&&q[c[h1]].v1<=a[i].v){
            int x=c[h1++];
            q[x].ans2+=getsum(asky(q[x].y1))-getsum(asky(q[x].y2-1));
        }
        add(asky(a[i].y));
    }
    while (h1<=tot){
        int x=c[h1++];
        q[x].ans2+=getsum(asky(q[x].y1))-getsum(asky(q[x].y2-1));
    }

    fo(i,1,z)tree[i]=0;
    h1=tot;
    fd(i,z,1){
        while (h1&&q[c[h1]].v1>=a[i].v){
            int x=c[h1--];
            q[x].ans1+=getsum(asky(q[x].y1))-getsum(asky(q[x].y2-1));
        }
        add(asky(a[i].y));
    }
    while (h1){
        int x=c[h1--];
        q[x].ans1+=getsum(asky(q[x].y1))-getsum(asky(q[x].y2-1));
    }


    h1=1;
    int l=1;
    a[0].v=a[1].v-1;
    a[z+1].v=a[z].v+1;
    fo(i,1,z+1){
        if (a[i].v>a[i-1].v){
            while(h1<=tot&&q[c[h1]].v1<a[i-1].v)h1++;
            while(h1<=tot&&q[c[h1]].v1==a[i-1].v){
                int ll=l,rr=i-1,a1(-1),a2(-1),x=c[h1++];
                db len=(q[x].x1-q[x].x2+q[x].y1-q[x].y2)/2;
                while (ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x<q[x].x1-len)ll=mid+1;
                    else{
                        rr=mid-1;
                        a1=mid;
                    }
                }
                ll=l;
                rr=i-1;
                while(ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x>q[x].x2+len)rr=mid-1;
                    else{
                        ll=mid+1;
                        a2=mid;
                    }
                }
                if (a1==-1||a2==-1)continue;
                q[x].ans3+=a2-a1+1;
            }
            l=i;
        }
    }
}

void type6(){
    int h1,tot(0);
    fo(i,1,p)
    if (q[i].ty==6)c[++tot]=i;

    fo(i,1,z)tree[i]=0;
    sort(a+1,a+1+z,cmp2);
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].x2;
        q[x].v2=q[x].y1;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=1;
    fo(i,1,z){
        while(h1<=tot&&q[c[h1]].v1<=a[i].x){
            int x=c[h1++];
            int t1=getsum(asky(q[x].v2-1)),t3=i-1;
            q[x].ans2+=t1;
            q[x].ans3+=t3-t1;
        }
        add(asky(a[i].y));
    }
    while(h1<=tot){
        int x=c[h1++];
        int t1=getsum(asky(q[x].v2-1)),t3=z;
        q[x].ans2+=t1;
        q[x].ans3+=t3-t1;
    }

    fo(i,1,z)tree[i]=0;
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].x1;
        q[x].v2=q[x].y2;
    }
    sort(c+1,c+1+tot,cmp3);
    h1=tot;
    fd(i,z,1){
        while(h1&&q[c[h1]].v1>=a[i].x){
            int x=c[h1--];
            int t1=getsum(ky),t2=getsum(asky(q[x].v2));
            q[x].ans3+=t2;
            q[x].ans1+=t1-t2;
        }
        add(asky(a[i].y));
    }
    while(h1){
        int x=c[h1--];
        int t1=getsum(ky),t2=getsum(asky(q[x].v2));
        q[x].ans3+=t2;
        q[x].ans1+=t1-t2;
    }

    fo(i,1,z)tree[i]=0;
    fo(i,1,tot){
        int x=c[i];
        q[x].v1=q[x].v1+q[x].v2;
        q[x].v2=0;
    }
    sort(c+1,c+1+tot,cmp3);
    fo(i,1,z)a[i].v=a[i].x+a[i].y;
    sort(a+1,a+1+z,cmp4);
    h1=1;
    fo(i,1,z){
        while (h1<=tot&&q[c[h1]].v1<=a[i].v){
            int x=c[h1++];
            q[x].ans2+=getsum(askx(q[x].x1-1))-getsum(askx(q[x].x2-1));
            q[x].ans3+=getsum(askx(q[x].x1))-getsum(askx(q[x].x1-1));
        }
        add(askx(a[i].x));
    }
    while (h1<=tot){
        int x=c[h1++];
        q[x].ans2+=getsum(askx(q[x].x1-1))-getsum(askx(q[x].x2-1));
        q[x].ans3+=getsum(askx(q[x].x1))-getsum(askx(q[x].x1-1));
    }

    fo(i,1,z)tree[i]=0;
    h1=tot;
    fd(i,z,1){
        while (h1&&q[c[h1]].v1>=a[i].v){
            int x=c[h1--];
            q[x].ans1+=getsum(askx(q[x].x1))-getsum(askx(q[x].x2));
            q[x].ans3+=getsum(askx(q[x].x2))-getsum(askx(q[x].x2-1));
        }
        add(askx(a[i].x));
    }
    while (h1){
        int x=c[h1--];
        q[x].ans1+=getsum(askx(q[x].x1))-getsum(askx(q[x].x2));
        q[x].ans3+=getsum(askx(q[x].x2))-getsum(askx(q[x].x2-1));
    }

    h1=1;
    int l;
    a[0].v=a[1].v-1;
    a[z+1].v=a[z].v+1;
    fo(i,1,z+1){
        if (a[i].v>a[i-1].v){
            while(h1<=tot&&q[c[h1]].v1<a[i-1].v)h1++;
            while(h1<=tot&&q[c[h1]].v1==a[i-1].v){
                int ll=l,rr=i-1,a1(-1),a2(-1),x=c[h1++];
                while (ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x<q[x].x2)ll=mid+1;
                    else{
                        rr=mid-1;
                        a1=mid;
                    }
                }
                ll=l;
                rr=i-1;
                while(ll<=rr){
                    int mid=(ll+rr)/2;
                    if (a[mid].x>q[x].x1)rr=mid-1;
                    else{
                        ll=mid+1;
                        a2=mid;
                    }
                }
                if (a1==-1||a2==-1)continue;
                q[x].ans3+=a2-a1+1;
            }
            l=i;
        }
    }
}

void getans(){
    fo(i,1,p)
    if (q[i].bz)printf("%d %d %d\n",q[i].ans2,q[i].ans1,q[i].ans3);
    else printf("%d %d %d\n",q[i].ans1,q[i].ans2,q[i].ans3);
}

int main(){
    n=get();m=get();z=get();p=get();
    fo(i,1,z)a[i].x=b[i].x=get(),a[i].y=b[i].y=get();
    fo(i,1,p){
        int x1=get(),y1=get(),x2=get(),y2=get();
        gettype(x1,y1,x2,y2,i);
    }
    type1();
    type2();
    type3();
    type4();
    type5();
    type6();
    getans();
}
samjia2000
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
专栏目录
POI2008BLO 吐槽
huyuncong的专栏
08-14 765
1123: [POI2008]BLO Time Limit: 10 Sec  Memory Limit: 162 MB Submit: 359  Solved: 126 [Submit][Status][Discuss] Description Byteotia城市
POI操作Excel
11-05
Apache POI 是一个开源项目，专门用于处理 Microsoft Office 格式文件，如 Word、Excel 和 PowerPoint。在 Java 开发中，POI 提供了丰富的 API，使得开发者可以方便地读写 Excel 文档。这里我们将详细探讨如何使用 ...
[POI2008] STA-Station
skyflying266的博客
05-22 229
[POI2008] STA-Station
POI 2008 Blockade
最新发布
cgj0128_1的博客
01-25 440
作为抗议的一个措施，程序员们计划封锁 Byteotia 的一个城镇，阻止人进入，离开甚至通过这个城镇。如图中点1，当删去与其相连的边后，它与其它所有点都不能互相到达，即(1,2),(1,3)......(2,1),(3,1),(4,1),(5,1)无法互通，其数目为。写一个程序，从标准输入读取 Byteotia 的路网，输出对于每个确定的城镇，如果程序员们封锁了这个城镇，将有多少次拜访不能进行？Byteotia 城市有 n 个城镇，m 条双向道路。每条道路连接两个不同的城镇，没有重复的道路，所有城镇连通。
[POI2008] STA-Station (学学换根dp)
xianqiuyigedao的博客
03-31 413
题目 题意: 给定一个 n 个点的树，请求出一个结点，使得以这个结点为根时，所有结点的深度之和最大。 一个结点的深度之定义为该节点到根的简单路径上边的数量。 思路: 换根dp. 1.先以任意一个点为根dfs1，求出f[root]。 2.考虑换根。看由u换成其子节点j之后的影响，进行dfs2。 从而在O(n)的时间内完成换根，得到。 f[j] = f[u] - s[j] + n - s[j] 时间复杂度: O(n) 代码: // Problem: P3478 [POI2008] STA-Station //
POI2008 Blockade
weixin_30274627的博客
08-04 84
为获得更好的阅读体验请访问我的 Blog Luogu 如果节点 \(i\) 不是割点，那么只有 \(2*(n-1)\) 个点对不连通。 如果节点 \(i\) 是割点，设在搜索树上节点 \(i\) 的子节点里有 \(k\) 个节点 \(p_1,p_2,\dots,p_k\) 满足割点判断法则 \(Low[p_t] \geq Dfn[x]\) 那么这时候将节点 \(i\) 其连着的边删除后整个图分成 ...
Excel POI读取封装(文件+示范代码)
03-02
Excel POI读取封装(文件+示范代码) package org.excel.service; import java.io.FileInputStream; import java.io.FileNotFoundException; import java.io.FileWriter; import java.io.IOException; import java....
poi解析EXCEL2010 java代码.doc
06-11
10. **循环遍历**: 代码通过双重循环遍历所有行和列，将数据存储在一个二维字符串数组`str`中。然后，这些数据可以进一步处理并映射到`VHrUser`对象。 11. **Date格式转换**: 代码中使用了`DateFormat`和`...
TensorFlow实现模型断点训练,checkpoint模型载入方式
12-17
深度学习中，模型训练一般都需要很长的时间，由于很多原因，导致模型中断... # saver.restore(sess, './log/' + "model_savemodel.cpkt-" + str(20000)) sess.run(tf.global_variables_initializer()) ckpt = tf.tr
导出excel错误处理
11-25
执行当前 Web 请求期间，出现未处理的异常。请检查堆栈跟踪信息，以了解有关该错误以及代码中导致错误的出处的详细信息。 ... ASP.NET 未被授权访问所请求的资源。请考虑授予 ASP.NET 请求标识访问此资源的权限。...
【POI2008】BLO
weixin_30693183的博客
05-04 171
这道题是在基本的tarjan求割点的算法上进一步加深，加上了计数问题。 基本思路是显然的，我们必定要跑一遍tarjan求割点，之后我们就要分类讨论。 若这个点不是割点，那么把他的边去掉之后这个点就与原来的联通块分开了，一共有(n-1)对，因为(x,y),(y,x)算两对，所以答案为n-1<<1 若这个点是割点，那么把他的边去掉之后这张图就变成了若干部分：这个点自己、在搜索树中这个...
【POI 2008, BLO】割点
weixin_51979465的博客
08-10 109
对于每一个点，如果这个点是割点，那么删了他之后对这个图来说是没有影响的，所以说这个点的贡献就是2*(n-1)，如果这个点是割点的话，那么对于把它删去的所有连通块都乘一下就行，注意求那个割边的过程，就是如果与这个点相连的点没有返租边的话那就是割点了，下面看代码吧。有n个点m条边无向连通图的图，保证没有重边和自环。对于每个点，输出将这个点的所有边删除之后，有多少点对不能互相连通。这里的点对是有顺序的，也就是(u,v)和(v,u)需要被统计两次。第一行包含两个数n,m(1≤n≤105,1≤m≤5×105)。...
BZOJ1129: [POI2008]Per
L_0_Forever_LF的专栏
12-08 827
模数m不是质数很麻烦qwq 先把他分解成m=∏Mi=1pikim=∏i=1Mpikim=\prod_{i=1}^Mpi^{ki} 用每个pikipikipi^{ki}做模数计算最后EXCRT合并 那么现在模数Mod=pikiMod=pikiMod=pi^{ki} 计算s的排名，按位枚举i，计算1~i-1位与s相同，第i位&lt; s的序列数，最后+1 令p[i]表示位置i在序列里排第p[...
P3469 [POI2008]BLO-Blockade
NOIP提高组日志
07-08 184
https://www.luogu.org/problemnew/show/P3469 #include <iostream> #include <cstdio> #include <cstring> #include <stack> #include <algorithm> using namespace std; const int...
BZOJ1125: [POI2008]Poc
L_0_Forever_LF的专栏
12-08 508
我是真的没有任何码力可言qwq splay维护串的hash值 每次更改hash值，在splay上给所有那个值的点打个答案标记 特殊处理交换同一个串code：#include<set> #include<map> #include<deque> #include<queue> #include<stack> #include<cmath> #include<ctime> #include<bit
[POI2008]BLO
a5163273的博客
05-14 107
Description Byteotia城市有n个 towns m条双向roads. 每条 road 连接 两个不同的 towns ,没有重复的road. 所有towns连通。 Input 输入n<=100000 m<=500000及m条边 Output 输出n个数，代表如果把第i个点去掉，将有多少对点不能互通。 Sample Input 5 5 1 2 2 3 1 3 3 ...
2008 POI各省类型 统计分析 2008poi 2008信息点
QQ_2801061513的博客
02-23 296
2008-餐饮服务 281906 2008-道路附属设施 10214 2008-地名地址信息 4561834 2008-风景名胜 37629 2008-公共设施 27113 2008-公司企业 850646 2008-购物服务 614680 2008-交通设施服务 32574 2008-金融保险服务 178688 2008-科教文化服务 243571 2008-汽车服务 210662 2008-商务住宅 131547 2008-生活服务 239899 20
洛谷3467 [POI2008]PLA-Postering 单调栈
olahiuj的博客
04-15 515
题目描述 All the buildings in the east district of Byteburg were built in accordance with the old arbitecture: they stand next to each other with no spacing inbetween. Together they form a very long cha
python爬取高德poi数据
11-11
以下是使用Python爬取高德POI数据的步骤： 1. 导入必要的库，如requests、pandas等。 2. 构造请求URL，包括高德开放... fnew.write(str(r[0])+','+str(r[1])+'|'+str(r[2])+','+str(r[3])+'\t'+str(count)+'\n') ```
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值