【zkw线段树套自顶向下splay】cf44g

最新推荐文章于 2019-10-12 13:22:37 发布
最新推荐文章于 2019-10-12 13:22:37 发布
阅读量1.0k 收藏
点赞数
分类专栏： 数据结构 codeforces
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/huyuncong/article/details/9934667
版权

考场上队友写这道题tle了，貌似是数组寻址花的时间过久加上splay常数过大

因此赛后re一遍，很久没写数据结构，调了老半天，最后9000+ms过了，虽然写得有点冗而且还是没改成指针，但是zkw线段树以及自顶向下splay的常数确实比较小

后来优化了一下删除，快了1000ms+

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
const int oo=1073741819;
using namespace std;
struct node{
    int y,id,l,r,m_id;
}c[7000000];
struct point{
    int x,y,w;
}b[200000];
struct rect{
    int lx,rx,ly,ry,z;
}a[200000];
int n,m,ss,ans[200000],st,ls[200000],rs[200000],root[5000000];
int w[500000],u[500000],p[500000],m1,tot;
int max(int x,int y)
{
    if (b[x].w<b[y].w) return x;
    return y;
}
void updata(int x)
{
    c[x].m_id=max(c[x].id,c[c[x].l].m_id);
    c[x].m_id=max(c[x].m_id,c[c[x].r].m_id);
}
int right(int x,int y)
{
    c[y].l=c[x].r,c[x].r=y,updata(y);
    return x;
}
int left(int x,int y)
{
    c[y].r=c[x].l,c[x].l=y,updata(y);
    return x;
}
void rconnect(int x)
{
    if (!rs[0]) {rs[rs[0]=1]=x;return ;}
    c[rs[rs[0]]].l=x,rs[++rs[0]]=x;
}
void lconnect(int x)
{
    if (!ls[0]) {ls[ls[0]=1]=x;return ;}
    c[ls[ls[0]]].r=x,ls[++ls[0]]=x;
}
int cmp(int x,int y,int id)
{
    if (y!=c[x].y) return (y<c[x].y) ? -1 : 1;
    if (id!=c[x].id) return (id<c[x].id) ? -1 : 1;
    return 0;
}
void splay(int &mid,int y,int id)
{
    int pd,tmp;
    ls[0]=rs[0]=0;
    for (;mid && ((pd=cmp(mid,y,id))!=0);) {
        if (pd<0) {
            if (c[mid].l && cmp(c[mid].l,y,id)<0) mid=right(c[mid].l,mid);
            rconnect(mid),tmp=c[mid].l,c[mid].l=0,mid=tmp;
        }
        else if (pd>0) {
            if (c[mid].r && cmp(c[mid].r,y,id)>0) mid=left(c[mid].r,mid);
            lconnect(mid),tmp=c[mid].r,c[mid].r=0,mid=tmp; 
        }
    }
    if (!mid) {
        mid=++ss;
        c[ss].y=y,c[ss].id=id;
    }
    if (c[mid].l) lconnect(c[mid].l);
    if (c[mid].r) rconnect(c[mid].r);
    for (int i=ls[0];i>=1;i--) updata(ls[i]);
    for (int i=rs[0];i>=1;i--) updata(rs[i]);
    if (ls[0]) c[mid].l=ls[1];
    if (rs[0]) c[mid].r=rs[1];
    updata(mid);
}
void origin()
{
    for (m1=1;m1<=tot+2;m1<<=1) ;
    ss=0;
}
void ins(int x,int y,int id)
{
    for (x+=m1;x;x>>=1) 
        splay(root[x],y,id);
}
void del(int x,int y,int id)
{
    for (x+=m1;x;x>>=1) {
//        cout<<ss<<endl;
//            cout<<ss<<endl;
//            cout<<root[x]<<' '<<c[root[x]].y<<' '<<c[root[x]].id<<' '<<y<<' '<<id<<endl;
//        c[root[x]].id=0;
        splay(root[x],y,id);
        if (!c[root[x]].r) root[x]=c[root[x]].l;
        else if (!c[root[x]].l) root[x]=c[root[x]].r;
        else {
            int k=c[root[x]].r;
            for (;c[k].l;k=c[k].l) ;
            splay(c[root[x]].r,c[k].y,c[k].id);
            c[c[root[x]].r].l=c[root[x]].l;
            root[x]=c[root[x]].r;
            updata(root[x]);
        }
    }
}
void doit(int x,int ly,int ry)
{
    splay(root[x],ry,m+1);
    splay(c[root[x]].l,ly,0);
    int rt=root[x];
//    cout<<rt<<' '<<c[rt].l<<' '<<c[rt].r<<' '<<c[c[rt].l].r<<endl;
    st=max(st,c[c[c[rt].l].r].m_id);
//    cout<<st<<endl;
}
void ask(int lx,int rx,int ly,int ry)
{
    lx+=m1-1,rx+=m1+1;
//    cout<<lx<<' '<<rx<<endl;
    for (;!((lx^rx)==1);lx>>=1,rx>>=1) {
        if ((lx&1)==0) doit(lx+1,ly,ry);
        if ((rx&1)==1) doit(rx-1,ly,ry);
    }
}
bool cmp2(int i,int j)
{
    return w[i]<w[j];
}
bool cmp1(int i,int j)
{
    return a[i].z<a[j].z;
}
int main()
{
//    freopen("g.in","r",stdin);
//    freopen("g.out","w",stdout);
    scanf("%d",&n);
    tot=0;
    for (int i=1;i<=n;i++) {
        scanf("%d%d%d%d%d",&a[i].lx,&a[i].rx,&a[i].ly,&a[i].ry,&a[i].z);
        w[++tot]=a[i].lx,p[tot]=i;
        w[++tot]=a[i].rx,p[tot]=i+n;
    }
    scanf("%d",&m);
    b[0].w=b[m+1].w=oo;
    for (int i=1;i<=m;i++) {
        scanf("%d%d",&b[i].x,&b[i].y);b[i].w=i;
        w[++tot]=b[i].x,p[tot]=i+n+n;
    }
//    cout<<b[194].x<<' '<<b[194].y<<endl;
    for (int i=1;i<=tot;i++) u[i]=i;
    sort(u+1,u+tot+1,cmp2);
    w[u[0]=0]=-1;
    for (int i=1,k=0;i<=tot;i++) {
        if (w[u[i]]!=w[u[i-1]]) k++;
        if (p[u[i]]<=n) a[p[u[i]]].lx=k;
        else if (p[u[i]]<=n+n) a[p[u[i]]-n].rx=k;
        else b[p[u[i]]-n-n].x=k;
    }
    origin();
    for (int i=1;i<=m;i++) ins(b[i].x,b[i].y,i);
    for (int i=1;i<=n;i++) u[i]=i;
    sort(u+1,u+n+1,cmp1); 
//    memset(p,0,sizeof(p));
//    cout<<ss<<endl;
//    for (int i=1;i<=n;i++) cout<<a[i].lx<<' '<<a[i].rx<<endl;
    for (int i=1;i<=n;i++) {
        st=0;
        ask(a[u[i]].lx,a[u[i]].rx,a[u[i]].ly,a[u[i]].ry);
//        cout<<st<<endl;
//        cout<<endl;
//        if (p[st]) cout<<oo<<' '<<st<<' '<<b[st].x<<' '<<b[st].y<<' '<<b[st].w<<endl;
//        if (st==194) cout<<endl;
        if (st==0) continue;
        p[st]=1;
        ans[st]=u[i];
//        cout<<b[st].y<<' '<<st<<endl;
        del(b[st].x,b[st].y,st);
    }
//    cout<<ss<<' '<<m1<<endl;
    for (int i=1;i<=m;i++) printf("%d\n",ans[i]);
    return 0;
}

快了1000ms+

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
const int oo=1073741819;
using namespace std;
struct node{
    int y,id,l,r,m_id;
}c[7000000];
struct point{
    int x,y,w;
}b[200000];
struct rect{
    int lx,rx,ly,ry,z;
}a[200000];
int n,m,ss,ans[200000],st,ls[200000],rs[200000],root[5000000];
int w[500000],u[500000],p[500000],m1,tot;
int max(int x,int y)
{
    if (b[x].w<b[y].w) return x;
    return y;
}
void updata(int x)
{
    c[x].m_id=max(c[x].id,c[c[x].l].m_id);
    c[x].m_id=max(c[x].m_id,c[c[x].r].m_id);
}
int right(int x,int y)
{
    c[y].l=c[x].r,c[x].r=y,updata(y);
    return x;
}
int left(int x,int y)
{
    c[y].r=c[x].l,c[x].l=y,updata(y);
    return x;
}
void rconnect(int x)
{
    if (!rs[0]) {rs[rs[0]=1]=x;return ;}
    c[rs[rs[0]]].l=x,rs[++rs[0]]=x;
}
void lconnect(int x)
{
    if (!ls[0]) {ls[ls[0]=1]=x;return ;}
    c[ls[ls[0]]].r=x,ls[++ls[0]]=x;
}
int cmp(int x,int y,int id)
{
    if (y!=c[x].y) return (y<c[x].y) ? -1 : 1;
    if (id!=c[x].id) return (id<c[x].id) ? -1 : 1;
    return 0;
}
void splay(int &mid,int y,int id)
{
    int pd,tmp;
    ls[0]=rs[0]=0;
    for (;mid && ((pd=cmp(mid,y,id))!=0);) {
        if (pd<0) {
            if (c[mid].l && cmp(c[mid].l,y,id)<0) mid=right(c[mid].l,mid);
            rconnect(mid),tmp=c[mid].l,c[mid].l=0,mid=tmp;
        }
        else if (pd>0) {
            if (c[mid].r && cmp(c[mid].r,y,id)>0) mid=left(c[mid].r,mid);
            lconnect(mid),tmp=c[mid].r,c[mid].r=0,mid=tmp; 
        }
    }
    if (!mid) {
        mid=++ss;
        c[ss].y=y,c[ss].id=id;
    }
    if (c[mid].l) lconnect(c[mid].l);
    if (c[mid].r) rconnect(c[mid].r);
    for (int i=ls[0];i>=1;i--) updata(ls[i]);
    for (int i=rs[0];i>=1;i--) updata(rs[i]);
    if (ls[0]) c[mid].l=ls[1];
    if (rs[0]) c[mid].r=rs[1];
    updata(mid);
}
void origin()
{
    for (m1=1;m1<=tot+2;m1<<=1) ;
    ss=0;
}
void ins(int x,int y,int id)
{
    for (x+=m1;x;x>>=1) 
        splay(root[x],y,id);
}
void del(int x,int y,int id)
{
    for (x+=m1;x;x>>=1) {
//        cout<<ss<<endl;
//            cout<<ss<<endl;
//            cout<<root[x]<<' '<<c[root[x]].y<<' '<<c[root[x]].id<<' '<<y<<' '<<id<<endl;
        splay(root[x],y,id);
        b[id].w=oo;
        updata(root[x]);
/*        splay(root[x],y,id);
        if (!c[root[x]].r) root[x]=c[root[x]].l;
        else if (!c[root[x]].l) root[x]=c[root[x]].r;
        else {
            int k=c[root[x]].r;
            for (;c[k].l;k=c[k].l) ;
            splay(c[root[x]].r,c[k].y,c[k].id);
            c[c[root[x]].r].l=c[root[x]].l;
            root[x]=c[root[x]].r;
            updata(root[x]);
        }*/
    }
}
void doit(int x,int ly,int ry)
{
    splay(root[x],ry,m+1);
    splay(c[root[x]].l,ly,0);
    int rt=root[x];
//    cout<<rt<<' '<<c[rt].l<<' '<<c[rt].r<<' '<<c[c[rt].l].r<<endl;
    st=max(st,c[c[c[rt].l].r].m_id);
//    cout<<st<<endl;
}
void ask(int lx,int rx,int ly,int ry)
{
    lx+=m1-1,rx+=m1+1;
//    cout<<lx<<' '<<rx<<endl;
    for (;!((lx^rx)==1);lx>>=1,rx>>=1) {
        if ((lx&1)==0) doit(lx+1,ly,ry);
        if ((rx&1)==1) doit(rx-1,ly,ry);
    }
}
bool cmp2(int i,int j)
{
    return w[i]<w[j];
}
bool cmp1(int i,int j)
{
    return a[i].z<a[j].z;
}
int main()
{
    freopen("g.in","r",stdin);
    freopen("g.out","w",stdout);
    scanf("%d",&n);
    tot=0;
    for (int i=1;i<=n;i++) {
        scanf("%d%d%d%d%d",&a[i].lx,&a[i].rx,&a[i].ly,&a[i].ry,&a[i].z);
        w[++tot]=a[i].lx,p[tot]=i;
        w[++tot]=a[i].rx,p[tot]=i+n;
    }
    scanf("%d",&m);
    b[0].w=b[m+1].w=oo;
    for (int i=1;i<=m;i++) {
        scanf("%d%d",&b[i].x,&b[i].y);b[i].w=i;
        w[++tot]=b[i].x,p[tot]=i+n+n;
    }
//    cout<<b[194].x<<' '<<b[194].y<<endl;
    for (int i=1;i<=tot;i++) u[i]=i;
    sort(u+1,u+tot+1,cmp2);
    w[u[0]=0]=-1;
    for (int i=1,k=0;i<=tot;i++) {
        if (w[u[i]]!=w[u[i-1]]) k++;
        if (p[u[i]]<=n) a[p[u[i]]].lx=k;
        else if (p[u[i]]<=n+n) a[p[u[i]]-n].rx=k;
        else b[p[u[i]]-n-n].x=k;
    }
    origin();
    for (int i=1;i<=m;i++) ins(b[i].x,b[i].y,i);
    for (int i=1;i<=n;i++) u[i]=i;
    sort(u+1,u+n+1,cmp1); 
//    memset(p,0,sizeof(p));
//    cout<<ss<<endl;
//    for (int i=1;i<=n;i++) cout<<a[i].lx<<' '<<a[i].rx<<endl;
    for (int i=1;i<=n;i++) {
        st=0;
        ask(a[u[i]].lx,a[u[i]].rx,a[u[i]].ly,a[u[i]].ry);
//        cout<<st<<endl;
//        cout<<endl;
//        if (p[st]) cout<<oo<<' '<<st<<' '<<b[st].x<<' '<<b[st].y<<' '<<b[st].w<<endl;
//        if (st==194) cout<<endl;
        if (st==0) continue;
        p[st]=1;
        ans[st]=u[i];
//        cout<<b[st].y<<' '<<st<<endl;
        del(b[st].x,b[st].y,st);
    }
//    cout<<ss<<' '<<m1<<endl;
    for (int i=1;i<=m;i++) printf("%d\n",ans[i]);
    return 0;
}



huyuncong
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
ZKW线段树源码-pascal
03-10
ZKW线段树的实现源码-pascal版，可保证其正确性并在logN时间复杂度类进行区间操作
ZKW线段树模板类
08-06
zkw线段树模板类，可动态统计区间最大值。代码略作修改即可动态统计区间和
ZKW线段树
06-10
线段树
【zkw线段树讲稿】统计的力量-线段树.pdf
04-14
【zkw线段树讲稿】统计的力量-线段树.pdf
统计的力量(zkw线段树)
06-09
传说中的ZKW线段树 比一般线段树快 这个是ZKW讲课的课件 线段树是二叉搜索树的一种。与一般的二叉搜索树不同的是，线段树保存的是所有元素可能取的值，而不是每个元素。通常情况下，线段树只用叶节点表示每个值。而...
【字符串新武器】后缀自动机
热门推荐
huyuncong的专栏
05-19 1万+
发链：http://neroysq.blogcn.com/articles/%E5%90%8E%E7%BC%80%E8%87%AA%E5%8A%A8%E6%9C%BA%E5%88%9D%E6%8E%A2.html http://blog.sina.com.cn/s/blog_7812e98601012cim.html 详细构造见上述链接，此处介绍性质与理解 后缀自动机具有两大性质，考虑转
【可持久化数据结构】函数式线段树
huyuncong的专栏
05-19 7049
clj的论文很不错 总体思想就是只赋值不修改，同时充分运用历史版本，也正因为这个特性，所以可以完成在线询问历史版本的功能 这东西不同于后缀自动机基本基于原有知识就可以有直观的理解，往往平常想题的时候想到某个思路认为无法实现就直接枪毙，但这些东西恰好用函数式编程就迎刃而解，如果用以前的思维方式很有可能直接滤过，囧的就是明明拿着钥匙却偏偏认为那道门打不开，所以对于此类题目关键就是要贴合函数式编程的
李其乐 免费的馅饼+郭华阳 新一代打鼹鼠
huyuncong的专栏
04-30 4109
今天做了一道郭华阳出的模拟题，发现竟然和冬令营的答辩试题如出一辙，只改了一点条件，冬令营时就被什么树状数组+平衡树，坐标移位弄得头晕脑胀，一点也没听懂，结果冤家路窄，这次又只拿了40分，T T重新翻ppt，其实只是优化了转移而已，学了解析几何后，对坐标轴理解更加深刻，这两道题也比较容易理解了。有m只鼹鼠出现在一个1 ~ n的长条上。对于第k只鼹鼠，我们用t[k], w[k], x[k]来描述。即在t[k]时刻，会有一只分值为w[k]的鼹鼠k出现在位置x[k]上。现在你有一个锤头，每一时刻它都停留在一个位置上
三维偏序 树套树&cdq分治
huyuncong的专栏
10-18 3830
一般的最长上升子序列是二维偏序，所以排序后用线段树维护可以达到nlogn程度，但是如果是三元组呢？ 因为多了一个偏序关系，所以维护起来相当麻烦。 roosephu说了一句经典的:一维快排，二维归并，三维树状数组，就是一维排序x后，再而在快排保证前一段x小于后一段x基础上归并排
动态维护凸包
huyuncong的专栏
02-03 3660
依次插入n个点，询问凸包面积。 利用granham思想，维护凸包，水平序要维护两个凸壳，比较麻烦，所以用极角序，每次查入用平衡树维护极角离他最近的点，此时到双向链表上开始模拟granham双向删点，直到无法删为止，因为每个点至多删一次，所以是nlogn的复杂度。面积维护用叉积即时维护即可。 依次插入n个半平面，询问半平面交面积。 与上题类似，维护法向量的极角，双向删半平面（当相邻一半平面两交
zkw线段树 运用
huyuncong的专栏
05-23 3523
poj 3468 经典的区间修改和区间和查询，需要灵活运用差分和前缀和。 设A为原数组，A'为差分数组（A'i：=Ai-A[I-1]）。 那么，我们如果要求Ai，便可从A'1+...A'I得出，由此看来A=SA'（前缀和）， 于是，我们想要快速知道Ai~j的和方法之一便是SAj-SAi，由上式可知，SA=SSA'（前缀和的前缀和） SSA'i=SA’1+...SA'i=i*A‘1+(I-
【dominator tree】 Lengauer-Tarjan algorithm
huyuncong的专栏
11-14 3067
题意：给定源点，求出源点到其他各点的关键点  Lengauer-Tarjan algorithm 按理说这也是个经典算法,跟lca的tarjan和强连通的tarjan都有极其相似之处，但是貌似并没有推广 感觉出题比较好出，先求个dominator tree，然后再在上面各种搞，虽然多半会出成一个拼接题... 首先有几个链接 一些概念 http://en.wikipedia.org/
【回文自动机】ural2040
huyuncong的专栏
11-16 2694
题意：依次添加字符询问增加的
zkw线段树 标记下放
huyuncong的专栏
10-23 2300
zkw线段树的标记上传好像只能做区间修改的rmq，但标记下放应用面就与朴素线段树差不多了，但是漂亮很多。 在进行修改或询问时，我们先将当前区间的标记全部下放，从左右开区间端点从顶向下跑一边下放就是，同时，对于修改后，再从左右区间把祖先更新，zkw的非递归性质发挥的淋漓尽致。其实很好编，下面两道题都是1A、 好吧，我又重温了zkw的博客，这种方法不太纯正，但是，暂时还不理解那种标记永久化，所以还
【后缀自动机】自动机<->后缀树<->后缀数组
huyuncong的专栏
03-27 2194
一直都说自动机建后缀树，一直没真正构出后缀树过... 其实建后缀树很简单，父亲边已经有了，关键是边代表的子串怎么求。从叶子一层一层向上，对于节点i，我们已经知道了它在原串的位置和逆序后缀长度，他的父亲的逆序后缀长度也知道，父亲又是i的逆序后缀的前缀，这就可以直接在原串定位了。 另外，通过父亲边还可以求出节点i可接受的最短子串，设i的最长子串长度为maxlen[i]，最短子串长度为minlen[
【ctsc2012】solution
huyuncong的专栏
05-19 2053
梭哈游戏 题意：给定一系列比较规则，并给出A，B当前手牌，讨论牌抽完后，A能赢B的概率 比较规则比较复杂，可以先判牌形，再将花色大小压成一个数，排序后我们就只要统计，在某副手牌之前的手牌中，有多少个与其无重复手牌，判重工作可以由容斥原理和hash表完成，对于每副手牌需要2^5判断复杂度为o（手牌数*2^5），手牌数最多只有c（51，5）5s时限还算充裕，对于hash，我们考虑用53进制压位，因
【fhqtreap】航空管制
huyuncong的专栏
04-05 1929
飞机因为航空管制晚点了一个多小时，
【动态树】弹飞绵羊
huyuncong的专栏
03-09 1752
简洁暴力动态树 类似于树链剖分，只不过重边不以子树大小唯一确定，而是随着询问适时改变，从而摆脱静态的限制，可以动态添边删边。 边分为实边与虚边，实边连的称为路径，用splay维护（以保证可以动态加删边），从一节点往根走时将路径上的虚边改为实边，通过splay加速上升，而两点lca为最后一个虚变实的节点。 具体实现时，每个节点要存父亲，而根节点的父亲为此链的最浅处连出去的虚边父亲，虚变实时将新
zkw线段树修正 标记上升
huyuncong的专栏
10-12 1743
恢复一下被吞的blog。。 今天才发现统计的力量很坑爹，同时发现了我写线段树的漏洞。 首先说我自己的问题，线段树最后一层的第一个节点不能用，但我以前用了而且一直没出问题，这次涉及区间修改才出现问题。 然后是统计的力量，里面关于区间修改，区间查询最大值的代码有bug，不但修改上传标记与查询回收标记一个是求最小值的，一个是最大值的，而且标记上传也不完全，查询时不能改开区间等等问题，幸亏有奥特曼的...
zkw线段树的python实现
最新发布
05-10
以下是 zkw 线段树的 Python 实现： ```python class ZkwSegmentTree: def __init__(self, n): self.sz = 1 while self.sz self.sz self.tree = [0] * (2 * self.sz) def update(self, i, val): i += ...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值