BZOJ 2683/4066(简单题-kd-tree重构)

最新推荐文章于 2018-10-03 13:59:00 发布
最新推荐文章于 2018-10-03 13:59:00 发布
阅读量1.8k 收藏
点赞数
分类专栏： KD-Tree 文章标签： bzoj
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/nike0good/article/details/48949063
版权

初始(N*N)矩阵全部为0，维护2个操作：
1.矩阵加
2.矩阵求和
1<=N<=500000,操作数不超过200000个，内存限制20M

2683离线
4066强制在线，数据加强

对于离线，我们可以离散化，然后树套树，cdq，怎么开心怎么来
但是在线的情况，我们不能离散化，于是用kd-tree树上求和
Bzoj 2683

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=next[p])  
#define Lson (x<<1)
#define Rson ((x<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (1000000000)
#define F (100000007)
#define MAXN (200000+10)
typedef long long ll;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int n;

int cmp_d=0;
class node
{
public:
    int x[2];
    int l,r,minv[2],maxv[2];
    int w,sumv;
    node(){}
    node(int a,int b,int _w){MEM(x) l=r=0; w=sumv=_w; x[0]=a,x[1]=b; Rep(i,2) minv[i]=maxv[i]=x[i];}
    int& operator[](int i){return x[i]; } 
};

int cmp(node a,node b){return a[cmp_d]<b[cmp_d];    }

class KD_Tree
{
public: 
    node a[MAXN*3];
    KD_Tree()
    {   
    }

    void mem()
    {
    }

    void update(node& o)
    {
        o.sumv=o.w;
        if (o.l)
        {
            node p=a[o.l];
            Rep(i,2) o.minv[i]=min(o.minv[i],p.minv[i]);
            Rep(i,2) o.maxv[i]=max(o.maxv[i],p.maxv[i]);    
            o.sumv+=p.sumv;
        }
        if (o.r)
        {
            node p=a[o.r];
            Rep(i,2) o.minv[i]=min(o.minv[i],p.minv[i]);
            Rep(i,2) o.maxv[i]=max(o.maxv[i],p.maxv[i]);    
            o.sumv+=p.sumv;
        }

    }

    int build(int L,int R,int nowd)
    {
        int m=(L+R)>>1;

        cmp_d=nowd;
        nth_element(a+L+1,a+m+1,a+R+1,cmp);

        if (L^m) a[m].l=build(L,m-1,nowd^1);
        if (R^m) a[m].r=build(m+1,R,nowd^1);

        update(a[m]);

        return m;

    } 

    int root;
    void _build(int L,int R,int nowd) //1-n的节点 至少为1 
    {
        root=build(L,R,nowd);
    }

    void insert(int o,int k,int nowd)
    {
        int p=a[o].x[nowd];
        int p2=a[k].x[nowd];

        if (p2<=p) 
        {
            if (a[o].l) 
                insert(a[o].l,k,nowd^1);
            else a[o].l=k;
        }
        else
        {
            if (a[o].r)
                insert(a[o].r,k,nowd^1);        
            else a[o].r=k;

        }   

        update(a[o]);   

    }
    void _insert(int k,int nowd)
    {
        int p=root;
        insert(root,k,nowd);        
    }


    int _x1,_y1,_x2,_y2;
    int _ans;

    void ask(int o)
    {
        if (o==0) return;

        if (_x1<=a[o].minv[0] && a[o].maxv[0]<=_x2 && _y1<=a[o].minv[1] && a[o].maxv[1]<=_y2 ) {
            _ans+=a[o].sumv;return;
        }       
        if (_x1<=a[o].x[0] && a[o].x[0]<=_x2 && _y1<=a[o].x[1] && a[o].x[1]<=_y2 ) {
            _ans+=a[o].w;
        }       

        if (a[o].l) {
            int p=a[o].l;
            if (a[p].minv[0]<=_x2 && _x1<=a[p].maxv[0] && a[p].minv[1]<=_y2 && _y1<=a[p].maxv[1] ) 
                ask(p);
        }
        if (a[o].r) {
            int p=a[o].r;
            if (a[p].minv[0]<=_x2 && _x1<=a[p].maxv[0] && a[p].minv[1]<=_y2 && _y1<=a[p].maxv[1] ) 
                ask(p);
        }


    }

    int _ask(int x1,int y1,int x2,int y2)
    {
        _x1=x1;_y1=y1;_x2=x2;_y2=y2;
        ;_ans=0;
        ask(root);
        return _ans; 
    }

}S;

int main()
{


    int N;
    cin>>N;

    n=0;
    S.a[++n]=node(0,0,0);
    S._build(1,n,0);

    int p;
    int ans=0;
    while (scanf("%d",&p)==1 && p^3)
    {
        if (p==1) {
            int x,y,A;
            scanf("%d%d%d",&x,&y,&A);
            S.a[++n]=node(x,y,A);
            S._insert(n,0);
        } else {
            int x1,y1,x2,y2;
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);

            ans=S._ask(x1,y1,x2,y2);
            printf("%d\n",ans); 
        }

    } 
    return 0;
}

但是这样不能保证kd-tree的平衡，过不了4066，
于是我们对kd-tree重构

Bzoj 4066

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=next[p])  
#define Lson (x<<1)
#define Rson ((x<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (1000000000)
#define F (100000007)
#define MAXN (200000+10)
#define fac 0.65
typedef long long ll;
int n;

int cmp_d=0;
class node
{
public:
    int x[2];
    int l,r,minv[2],maxv[2];
    int w,sumv,siz;
    node(){}
    node(int a,int b,int _w){l=r=0; siz=1; w=sumv=_w; x[0]=a,x[1]=b; Rep(i,2) minv[i]=maxv[i]=x[i];}
    int& operator[](int i){return x[i]; } 
};

int cmp(node a,node b){return a[cmp_d]<b[cmp_d];    }


int cmp2(int i,int j); 

int p;
char c;
int read()
{
    while (c=getchar(),!isdigit(c));
    p=c-'0';
    while (isdigit(c=getchar())) p=p*10+c-'0'; return p;
}
class KD_Tree
{
public: 
    node a[MAXN];

    void update(node& o)
    {
        o.sumv=o.w;
        o.minv[0]=o.maxv[0]=o.x[0];o.minv[1]=o.maxv[1]=o.x[1]; 
        o.siz=1;
        if (o.l)
        {
            node p=a[o.l];
            Rep(i,2) o.minv[i]=min(o.minv[i],p.minv[i]);
            Rep(i,2) o.maxv[i]=max(o.maxv[i],p.maxv[i]);    
            o.sumv+=p.sumv;
            o.siz+=p.siz;
        }
        if (o.r)
        {
            node p=a[o.r];
            Rep(i,2) o.minv[i]=min(o.minv[i],p.minv[i]);
            Rep(i,2) o.maxv[i]=max(o.maxv[i],p.maxv[i]);    
            o.sumv+=p.sumv;
            o.siz+=p.siz;
        }

    }

    int build(int L,int R,int nowd,node *a)
    {
        int m=(L+R)>>1;

        cmp_d=nowd;
        nth_element(a+L+1,a+m+1,a+R+1,cmp);

        if (L^m) a[m].l=build(L,m-1,nowd^1,a);
        if (R^m) a[m].r=build(m+1,R,nowd^1,a);

        update(a[m]);

        return m;

    } 

    int po[MAXN],pt;
    void dfs(int x)
    {
        po[++pt]=x;
        if (a[x].l) dfs(a[x].l);
        if (a[x].r) dfs(a[x].r);
    }

    int rebuild(int L,int R,int nowd)
    {
        int m=(L+R)>>1;

        cmp_d=nowd;
        nth_element(po+L+1,po+m+1,po+R+1,cmp2);
        int now=po[m];
        a[now].l=a[now].r=0;
        if (L^m) a[now].l=rebuild(L,m-1,nowd^1);
        if (R^m) a[now].r=rebuild(m+1,R,nowd^1);

        update(a[now]);

        return now;

    } 

    int root;
    void _build(int L,int R,int nowd) //1-n的节点 至少为1 
    {
        root=build(L,R,nowd,a);
    }

    int insert(int o,int k,int nowd)
    {
        if (!o) return k;
        int p=a[o].x[nowd];
        int p2=a[k].x[nowd];
        int nx=0;
        if (p2<=p) 
        {
            a[o].l=insert(a[o].l,k,nowd^1);
            nx=a[o].l;
        }
        else
        {
            a[o].r=insert(a[o].r,k,nowd^1);
            nx=a[o].r;
        }   
        update(a[o]);   

        if (a[nx].siz>(double)a[o].siz*fac)
        {
            pt=0,dfs(o);
            o=rebuild(1,pt,nowd);   

        }
        return o;
    }
    void _insert(int k,int nowd)
    {
        int p=root;
        root = insert(root,k,nowd);     
    }


    int _x1,_y1,_x2,_y2;
    int _ans;

    void ask(int o)
    {
        if (o==0) return;

        if (_x1<=a[o].minv[0] && a[o].maxv[0]<=_x2 && _y1<=a[o].minv[1] && a[o].maxv[1]<=_y2 ) {
            _ans+=a[o].sumv;return;
        }       
        if (_x1<=a[o].x[0] && a[o].x[0]<=_x2 && _y1<=a[o].x[1] && a[o].x[1]<=_y2 ) {
            _ans+=a[o].w;
        }       

        if (a[o].l) {
            int p=a[o].l;
            if (a[p].minv[0]<=_x2 && _x1<=a[p].maxv[0] && a[p].minv[1]<=_y2 && _y1<=a[p].maxv[1] ) 
                ask(p);
        }
        if (a[o].r) {
            int p=a[o].r;
            if (a[p].minv[0]<=_x2 && _x1<=a[p].maxv[0] && a[p].minv[1]<=_y2 && _y1<=a[p].maxv[1] ) 
                ask(p);
        }


    }

    int _ask(int x1,int y1,int x2,int y2)
    {
        _x1=x1;_y1=y1;_x2=x2;_y2=y2;
        ;_ans=0;
        ask(root);
        return _ans; 
    }

}S;

int cmp2(int i,int j){return S.a[i].x[cmp_d]<S.a[j].x[cmp_d];   }

int main()
{
//  freopen("bzoj4066_data.in","r",stdin);


    int N=read();

    n=0;
    S.a[++n]=node(N/2,N/2,0);
    S._build(1,n,0);

    int p;
    int ans=0;
    int x,y,A;
    int x1,y1,x2,y2;

    while (scanf("%d",&p)==1 && p^3)
    {
//      cout<<"t"<<endl;
        if (p==1) {
            x=read(),y=read(),A=read();
            x^=ans;y^=ans;A^=ans;
            S.a[++n]=node(x,y,A);
            S._insert(n,0);
        } else {
            x1=read(),y1=read(),x2=read(),y2=read();
            x1^=ans,y1^=ans,x2^=ans,y2^=ans;
            ans=S._ask(x1,y1,x2,y2);
            printf("%d\n",ans); 
        }

    } 
    return 0;
}
nike0good
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
【BZOJ2683】简单题
cz_xuyixuan的博客
04-25 355
【题目链接】点击打开链接【三倍经验链接】【BZOJ1176】【Balkan2007】Mokia【BZOJ4066】简单题【思路要点】KDTree+替罪羊树式重构即可。时间复杂度\(O(N\sqrt{N})\)。【代码】#include&lt;bits/stdc++.h&gt; using namespace std; #define MAXN 200005 #define ALPHA 0.65 t...
BZOJ 2989 数列 —— kd-tree + 旋转坐标系
Scar_Halo
04-23 242
题目链接：点我啊╭(╯^╰)╮ 题目大意： 解题思路：     将 （i,a[i]）（i,a[i]）（i,a[i]） 视为一个点，查询即为曼哈顿距离 ≤k\le k≤k 的点数     考虑用 kd−treekd-treekd−tree 维护，但查询项是一个以 （i,a[i]）（i,a[i]）（i,a[i]） 为中心...
kd-tree讲解 & bzoj 2648 & 2716 & 3053 题解
阿蒋的专栏
06-24 6778
【KD-TREE介绍】在SYC1999大神的“蛊惑”下，我开始接触这种算法。 首先，大概的概念可以去百度百科。具体实现，我是看RZZ的代码长大的。 我们可以想象在平面上有N个点。首先，按横坐标排序找到最中间的那个点。然后水平划一条线，把平面分成左右两个部分。再递归调用左右两块。注意，在第二次（偶数次）调用的时候，是找到纵坐标中最中间的点，并垂直画一条线。 这样效率看上去很好。维护的时候有点像
BZOJ2626 JZPFAR及kd-tree入门
Rayment的博客
06-10 479
Problem BZOJ权限题 洛谷 Solution hdu4347和这道题也很类似。 我们用kd-tree来做，用一个堆来维护可能成为k远的几个点，不断更新答案。优先搜可能更远的子树，再根据估价函数判断需不需要搜另一个子树。注意重载小于号和估价函数不要写错了。注意当现在的节点少于k时仍然不能进入空子树，否则可能会影响答案。 Code #include &amp;amp;lt;algori...
BZOJ 4066 - KDtree + 重构
过路人丶的博客
10-03 278
题目链接:点击这里 解题思路: 这题被强制在线了如果是离线(BZOJ 2683)可以CDQ分治 + 树状数组做. 所以只能用KDtree了. 若每次都将新节点插入原DKT中,最终会导致树会非常的不平衡,而失去了意义,导致超时. 所以我们要设置一个阈值,当节点数到达阈值时就对KDT进行重构保证稳定. #include&lt;cstdio&gt; #include&lt;cs...
BZOJ4066 简单题
neither_nor
07-13 988
强制在线还卡内存，于是上KDT，插入的时候替罪羊式重建，真是tm玄学复杂度 #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define MAXN 2
Node.js-ZOJCH是BZOJ题库的离线版
08-09
【标题】"Node.js-ZOJCH是BZOJ题库的离线版"揭示了这个项目的核心功能，即提供一个基于Node.js平台的离线版本的BZOJ（BestCoder Judging System）题库。BZOJ是一个在线编程竞赛平台，而ZOJCH则允许用户在没有网络...
【(伪)Top tree】【BZOJ3153 Sone1 题解】-附件资源
03-05
【(伪)Top tree】【BZOJ3153 Sone1 题解】-附件资源
wangjunrui666#bzoj-problem##2508. 简单题1
最新发布
07-25
1. 平面上加入一条直线 2. 删除一条已加入的直线 3. 求一个点到平面上所有直线距离平方和最小，你需要输出这个最小值
[资料共享] BZOJ所有题 OJ 在线评测 题面 PDF版 非常良心-附件资源
03-05
[资料共享] BZOJ所有题 OJ 在线评测 题面 PDF版 非常良心-附件资源
BZOJ2683 简单题 [KD-tree]
lemonoil的博客
10-10 496
2683: 简单题Time Limit: 50 Sec  Memory Limit: 128 MBSubmit: 1635  Solved: 661[Submit][Status][Discuss]Description你有一个N*N的棋盘，每个格子内有一个整数，初始时的时候全部为0，现在需要维护两种操作：   命令
BZOJ4066 - 简单题
VisJiao的博客
01-12 393
k-d树上求和
[模版] 带部分重构的KD_Tree
visit_world的博客
04-24 1326
RT ， 从黄学长那里初步学习到写KD树的姿势后自己摸索着写了指针＋替罪羊版的，效率还可以….能够忍受。BZOJ 4066 简单题#include <cstdio> #include <cstring> #include <cstdlib> #include <algorithm>#define rep(i, x, y) for (int i = (x), _ = (y); i <= _; ++i)
【kd-tree】BZOJ4520 CQOI2016K远点对
蒟蒻的奋斗史
04-13 2104
传送门平生第一次在bzoj上排到了rank1，载图留恋 之所以写这篇博客是因为CQOI2016的SB出题人居然让什么三分和卡壳这样A了（其实三分感觉不是很好卡…） 以下是Claris老司机对卡壳错误的证明…（图片来至UOJ用户群…） 其实这是一个模板题… 第k远点对肯定在某一个点的前k远距离中… 然后就用kd-tree把某一个点的前k远距离都找出来… 然后再加一些弱弱滴剪枝就可以
bzoj 4066: 简单题 带部分重构的KDtree
beginend
12-10 347
题意你有一个N*N的棋盘，每个格子内有一个整数，初始时的时候全部为0，现在需要维护两种操作： 1 x y A将格子x,y里的数字加上A 2 x1 y1 x2 y2输出x1 y1 x2 y2这个矩形内的数字和 3 终止程序 n<=500000，操作数<=200000，强制在线，空间限制为20M。分析如果没有空间限制和强制在线的话，就可以直接离线cdq分治或二维线段树。现在要求强制在线，就只能用
【BZOJ4520】K远点对（KD-Tree)
小蒟蒻yyb的博客
05-23 495
题面 BZOJ 洛谷 题解 考虑暴力。 维护一个大小为KKK的小根堆，然后每次把两个点之间的距离插进去，然后弹出堆顶 这样子可以用KD−TreeKD−TreeKD-Tree优化： 如果当前平面内可以产生的最远距离小于堆顶，则证明这个平面内的点无法产生贡献 就不需要计算了 #include&lt;iostream&gt; #include&lt;cstdio&gt; #in...
bzoj4066 简单题
Elijahqi
01-12 221
(http://www.elijahqi.win/2018/01/12/bzoj4066-%E7%AE%80%E5%8D%95%E9%A2%98) Description 你有一个N*N的棋盘，每个格子内有一个整数，初始时的时候全部为0，现在需要维护两种操作： 命令 参数限制 内容 1 x y A 1 2 x1 y1 x2 y2 1 1 输出x1 y1 x2 y2
BZOJ2683: 简单题 kdtree
リーインカーネイション
12-27 608
题目大意：n*n矩阵，单点加，矩形查和 1<=N<=500000,操作数不超过200000个，内存限制20M。 正解是cdq分治，强行拿kdtree卡了过去。 1.n是矩形大小，不是操作数量，不要被坑了 2.kdtree是空树直接返回0，否则就崩了 3.nth_element的写法是(begin,nth,end,comp)，话说我之前那个崂山白花蛇草水错成那样都能过。。。 4.kdtre
bzoj 4358: permu （莫队+栈||KD-tree||莫队+线段树）
clover_hxy的博客
03-31 1548
题目描述传送门题目大意：给出一个长度为n的排列P(P1,P2,…Pn)，以及m个询问。每次询问某个区间[l,r]中，最长的值域 连续段长度。题解1：莫队+线段树用权值线段树维护区间的连续最大长度，左端连续最长，右端连续最长。 应该是比较好想，好写，好调的做法。 时间复杂度O(nn√logn)O(n\sqrt n log n) ,实际测速中貌似会TLE，于是利用毕生所学进行卡常，最终卡时A掉了。
BZOJ JOI 2013-2014春季训练合宿题目解析1
"BZOJ JOI 2013~2014 春季training合宿 系列题解1" 这篇题解涉及到的是一个编程竞赛中的题目，主要运用了哈希算法来解决。哈希算法在计算机科学中是一种高效的数据结构和算法，它能够快速查找、插入和删除元素，其...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值