BZOJ 4811([Ynoi2017]由乃的OJ-树链剖分)

最新推荐文章于 2019-02-21 20:15:13 发布

nike0good

最新推荐文章于 2019-02-21 20:15:13 发布

阅读量931

点赞数

分类专栏：树链剖分

本文链接：https://blog.csdn.net/nike0good/article/details/70984151

版权

树链剖分专栏收录该内容

7 篇文章 0 订阅

订阅专栏

Description

由乃正在做她的OJ。现在她在处理OJ上的用户排名问题。OJ上注册了n个用户，编号为1～”，一开始他们按照编号
排名。由乃会按照心情对这些用户做以下四种操作，修改用户的排名和编号：然而由乃心情非常不好，因为Deus天
天问她题。。。因为Deus天天问由乃OI题，所以由乃去学习了一下OI，由于由乃智商挺高，所以OI学的特别熟练她
在RBOI2016中以第一名的成绩进入省队，参加了NOI2016获得了金牌保送
Deus：这个题怎么做呀？
yuno：这个不是NOI2014的水题吗。。。
Deus：那如果出到树上，多组链询问，带修改呢？
yuno：诶。。。？？？
Deus：这题叫做睡觉困难综合征哟~
虽然由乃OI很好，但是她基本上不会DS，线段树都只会口胡，比如她NOI2016的分数就是100+100+100+0+100+100。
。。NOIP2017的分数是100+0+100+100+0+100所以她还是只能找你帮她做了。。。
给你一个有n个点的树，每个点的包括一个位运算opt和一个权值x，位运算有&,l,^三种，分别用1,2,3表示。
每次询问包含三个数x,y,z,初始选定一个数v。然后v依次经过从x到y的所有节点，每经过一个点i，v就变成v opti
xi，所以他想问你，最后到y时，希望得到的值尽可能大，求最大值？给定的初始值v必须是在[0,z]之间。每次修
改包含三个数x,y,z，意思是把x点的操作修改为y，数值改为z

Input

第一行三个数n，m，k。k的意义是每个点上的数，以及询问中的数值z都 <2^k。之后n行
每行两个数x,y表示该点的位运算编号以及数值
之后n - 1行，每行两个数x,y表示x和y之间有边相连
之后m行，每行四个数，Q,x,y,z表示这次操作为Q(1位询问，2为修改)，x，y，z意义如题所述
0 <= n , m <= 100000 , k <= 64
Output

对于每个操作1，输出到最后可以造成的最大刺激度v
Sample Input

5 5 3

1 7

2 6

3 7

3 6

3 1

1 2

2 3

3 4

1 5

1 1 4 7

1 1 3 5

2 1 1 3

2 3 3 3

1 1 3 2
Sample Output

我们存2个值，之前是0、1的情况下，经过一堆操作后分别会变为几
然后就可以压位，然后转到树链刨分上。
注意向上和向下的合并要分别计算

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<vector>
#include<iomanip> 
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 (o<<1)
#define Rson ((o<<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 (2139062143)
#define F (100000007)
#define MAXN (100000+10)
#define MAXM (200000+10)
#define pb push_back
#define mp make_pair
#pragma comment(linker, "/STACK:1024000000,1024000000") 
typedef long long ll;
typedef unsigned long long ull;
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+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
inline ull read()
{
    ull x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
struct node{
    ull x0,x1;
    node(){}
    node(ull _a,ull _b):x0(_a),x1(_b){}
    friend node operator+(node a,node b) {
        node c;
        c.x0=(a.x0&b.x1)|((~a.x0)&b.x0);  
        c.x1=(a.x1&b.x1)|((~a.x1)&b.x0);  
        return c;
    }
    node(int op,ull x) {
        switch(op) {
            case 1:{
                x0=0&x;
                x1=(~0)&x;
                break;
            }
            case 2:{
                x0=0|x;
                x1=(~0)|x;
                break;
            }
            case 3:{
                x0=0^x;
                x1=(~0)^x;
                break;
            }
        }
    }
};

int op[MAXN]; 
ull a[MAXN];

node vl[MAXN<<2],vr[MAXN<<2];
void pushup(int o) {
    vl[o]=vl[Lson]+vl[Rson];
    vr[o]=vr[Rson]+vr[Lson];
}
void build(int l,int r,int o) {
    if (l==r) {
        vl[o]=vr[o]=node(op[l],a[l]);   return ;
    }
    int m=(l+r)>>1;
    build(l,m,Lson),build(m+1,r,Rson);
    pushup(o);
}
void update(int l,int r,int o,int p) {
    if (l==r) {
        vl[o]=vr[o]=node(op[l],a[l]);   return ;
    }
    int m=(l+r)>>1;
    if (p<=m) update(l,m,Lson,p); 
    else update(m+1,r,Rson,p);
    pushup(o);
}
node queryL(int l,int r,int o,int L,int R) {
    if(L<=l && r<=R ) return vl[o];
    int m=(l+r)>>1;
    if(R<=m) return queryL(l,m,Lson,L,R);
    else if(L>m) return queryL(m+1,r,Rson,L,R); 
    else return queryL(l,m,Lson,L,R)+queryL(m+1,r,Rson,L,R);
}
node queryR(int l,int r,int o,int L,int R) {
    if(L<=l && r<=R ) return vr[o];
    int m=(l+r)>>1;
    if(R<=m) return queryR(l,m,Lson,L,R);
    else if(L>m) return queryR(m+1,r,Rson,L,R); 
    else return queryR(m+1,r,Rson,L,R)+queryR(l,m,Lson,L,R);
}


int n,m;
struct Tree{
    void mem(){MEM(Pre) siz=1;}
    int edge[MAXM],Next[MAXM],Pre[MAXN],siz;
    void addedge(int u,int v)
    {
        edge[++siz]=v;
        Next[siz]=Pre[u];
        Pre[u]=siz;
    }
    void addedge2(int u,int v){addedge(u,v);addedge(v,u);}
    bool vis[MAXN]; 
    int cnt,id[MAXN];
    int son[MAXN],dep[MAXN],sz[MAXN],top[MAXN],pre[MAXN],q[MAXN];
    void build()
    {
        MEM(vis) cnt=0; MEM(id)
        MEM(son) MEM(dep) MEM(sz) MEM(top) MEM(pre) MEM(q) 
        int r=1;
        vis[dep[1]=q[1]=1]=1;
        For(i,r)
        {
            int u=q[i];
            Forp(u)
            {
                int v=edge[p];
                if (vis[v]) continue; else vis[v]=1;
                dep[ q[++r]=v ]=dep[u]+1;
                pre[v]=u;
            }
        }
        ForD(i,r) {
            sz[pre[q[i]]] += ++sz[q[i]];
            if (sz[son[pre[q[i]]]]<sz[q[i]] ) son[pre[q[i]]] = q[i];
        }
        For(i,r) {
            if (!top[q[i]])
                for(int x=q[i];x;x=son[x]) {
                    top[x]=q[i];
                    id[x]=++cnt; 
                }
        }

    }
    int lca(int a,int b)
    {
        while(1) {
            if (top[a]==top[b]) return dep[a]<=dep[b] ? a:b;
            if (dep[top[a]]<dep[top[b]]) swap(a,b);
            a=pre[top[a]];
        }
    }

    node Ask(int a,int b)
    {
        node vl(3,0),vr(3,0);
        while (top[a]^top[b]) {
            if (dep[top[a]]>dep[top[b]]) {
                vl=vl+queryR(1,n,1,id[top[a]],id[a]);
                a=pre[top[a]];
            } else {
                vr=queryL(1,n,1,id[top[b]],id[b])+vr;
                b=pre[top[b]];
            }
        }   
        if (dep[a]>dep[b]) {
            return vl+queryR(1,n,1,id[b],id[a]) + vr;
        } else {
            return vl+queryL(1,n,1,id[a],id[b])+vr;
        }
    }

}S;
int op2[MAXN];
ull a2[MAXN];
int main()
{
//  freopen("bzoj4811.data","r",stdin);
//  freopen("bzoj4811.out","w",stdout);

    n=read(),m=read();int k=read();
    S.mem();
    For(i,n) op2[i]=read(),a2[i]=read();
    For(i,n-1)  {
        S.addedge2(read(),read());
    }
    S.build();
    For(i,n) {
        op[S.id[i]]=op2[i];
        a[S.id[i]]=a2[i];
    }
    build(1,n,1);

    For(i,m) {
        ull op=read(),x=read(),y=read(),z=read();
        if (op==1) {
            node t=S.Ask(x,y);
            ull ans=0,now=0;
            RepD(i,k-1) {
                if(t.x0&(1ull<<i)){ 
                    ans+=(1ull<<i);
                }else if ((t.x1&(1ull<<i))&&now+(1ull<<i)<=z) {
                    ans+=(1ull<<i);
                    now+=(1ull<<i);
                }
            }
            printf("%llu\n",ans);
        } else {
            ::op[S.id[x]]=y;
            ::a[S.id[x]]=z;
            update(1,n,1,S.id[x]);
        }
    }
    return 0;
}

nike0good

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
BZOJ 4811([Ynoi2017]由乃的OJ-树链剖分)

Description由乃正在做她的OJ。现在她在处理OJ上的用户排名问题。OJ上注册了n个用户，编号为1～”，一开始他们按照编号排名。由乃会按照心情对这些用户做以下四种操作，修改用户的排名和编号：然而由乃心情非常不好，因为Deus天天问她题。。。因为Deus天天问由乃OI题，所以由乃去学习了一下OI，由于由乃智商挺高，所以OI学的特别熟练她在RBOI2016中以第一名的成绩进入省队，参
复制链接

扫一扫

专栏目录