luogu P3613 睡觉困难综合征

题目地址
做法和原题差不多原题地址
显然每一个节点开一个大小为K的数组储存答案肯定是不行的
但是 k k <=64,所以可以开一个unsignedlonglong储存答案
但是合并答案有点恶心
每个节点储存两个值 A A B A A 表示INF经过操作后的值, B B 表示0经过操作后的值
对于一个节点
首先 LA L A LB L B 为左边的答案, RA R A RB R B 为右边的答案
以A为例
对于每一位i
如果 LA L A 的第 i i 位为0那么A的第 i i 位的值为RA的第 i i
如果LA的第 i i 位为1那么A的第 i i 位的值为RB的第 i i

A=((notLA)andRA)or(LAandRB)

同理

B=((notLB)andRA)or(LBandRB) B = ( ( n o t L B ) a n d R A ) o r ( L B a n d R B )

那么这题就可以树链剖分或者LCT( 大常数树)
个人比较推荐树链剖分,因为LCT不仅难调而且还不好写并且常数还巨大
因为LCT有翻转操作所以还要维护 RevA R e v A RevB R e v B 表示反向操作的答案
总之能不写LCT就不写
ps:但是我还是写的LCT

#include <cstdio>
#include <cctype>
#include <cmath>
#include <algorithm>
#define File(S) freopen(S".in","r",stdin),freopen(S".out","w",stdout)
using namespace std;
namespace Data_stream{
    const int IN_LEN=1<<17;
    const int OUT_LEN=1<<17;
    char buf[IN_LEN],obuf[OUT_LEN],*pcur=obuf,*HEAD=buf,*TAIL=buf,ch;
    int BUF,PS=0,Buf[30];
    bool sign;
    inline char gec(){
        if(HEAD==TAIL){
            TAIL=(HEAD=buf)+fread(buf,1,IN_LEN,stdin);
            if(HEAD==TAIL)return EOF;
        }
        return *HEAD++;
    }
    inline void puc(char ch){
        if(pcur==obuf+OUT_LEN){
            fwrite(obuf,1,OUT_LEN,stdout);
            pcur=obuf;
        }
        *pcur++=ch;
    }
    template <typename T>
    inline int read(T&a){
        for(sign=0,ch=gec();!isdigit(ch)&&ch!=EOF;ch=gec())
            if(ch=='-')sign=1;
        if(ch==EOF)return EOF;
        for(a=0;isdigit(ch);ch=gec())a=((a+(a<<2))<<1)+(ch^'0');
        if(sign)a=-a;
        return 0;
    }
    template <typename T>
    inline void print(T x){
        if(x<0)puc('-'),x=-x;
        for(;x;x=BUF)
            BUF=x/10,Buf[++PS]=x-((BUF+(BUF<<2))<<1);
        while(PS)puc('0'^Buf[PS--]);
    }
    inline void flush(){
        fwrite(obuf,1,pcur-obuf,stdout);
    }
}
using namespace Data_stream;
#define Hf(x) (x->f?((x->f->s[0]==x)||(x->f->s[1]==x)):0)
#define Get(type) (type=(~type&s[1]->Fz)|(type&s[1]->Fe))
#define Get_(type) (type=(~type&s[0]->Rfz)|(type&s[0]->Rfe))
typedef unsigned long long Ul;
const int N=2e5+10;
const Ul INF=0xffffffffffffffff;//Max unsigned long long 
inline void Calc(Ul&v,const Ul&val,const int&opt){(!opt)?(v&=val):((opt&1)?(v|=val):(v^=val));}
struct Tree{
    Tree(int opt=0,Ul v=0):opt(opt),v(v){Fz=Fe=Rfz=Rfe=0;rev=0;s[0]=s[1]=f=0;} 
    int opt;
    Ul Fz,Fe,Rfz,Rfe,v;
    bool rev;
    Tree *s[2],*f;  
    inline void Push(){
        if(!rev)return;
        rev=0;swap(s[0],s[1]);
        swap(Fz,Rfz);swap(Fe,Rfe);
        if(s[0])s[0]->rev^=1;
        if(s[1])s[1]->rev^=1;
    }
    inline void Update(){
        if(s[0])Fz=s[0]->Fz,Fe=s[0]->Fe;
        else Fz=0,Fe=INF;
        Calc(Fz,v,opt);Calc(Fe,v,opt);
        if(s[1])Get(Fz),Get(Fe),Rfz=s[1]->Rfz,Rfe=s[1]->Rfe;
        else Rfz=0,Rfe=INF;
        Calc(Rfz,v,opt);Calc(Rfe,v,opt);
        if(s[0])Get_(Rfz),Get_(Rfe);
    }
    void*operator new(size_t);
}t[N],*cur=t,*vec[N],*stack[N];
inline void*Tree::operator new(size_t){return cur++;} 
int n,m,k,top;
inline void Rot(Tree *p){
    Tree *f=p->f;int v=(p!=f->s[0]);
    if(f->s[v]=p->s[v^1])f->s[v]->f=f;
    if(Hf(f))f->f->s[f!=f->f->s[0]]=p;
    p->f=f->f;p->s[v^1]=f;f->f=p;
    f->Update();p->Update();
} 
inline void Splay(Tree *p){
    Tree *last;
    stack[++top]=p->s[0];stack[++top]=p->s[1];
    for(last=p;Hf(p);p=p->f)stack[++top]=p->f->s[p==p->f->s[0]],stack[++top]=p;
    stack[++top]=p;
    for(;top;top--)if(stack[top])stack[top]->Push();
    for(p=last;Hf(p);Rot(p)){
        Tree *f=p->f;
        if(Hf(f)) Rot(((p!=f->s[0])^(f!=f->f->s[0]))?p:f);
    }
}
inline void Access(Tree *x){for(Tree *last=0;x;last=x,x=x->f){Splay(x);x->s[1]=last;x->Update();}} 
inline void MakeRoot(Tree *x){Access(x);Splay(x);x->rev^=1;}
inline void link(Tree*x,Tree*y){MakeRoot(x);x->f=y;}
inline void cut(Tree*x,Tree*y){MakeRoot(x);Access(y);Splay(x);x->s[1]=y->f=0;x->Update();}
int main()
{
    int x,y,type;
    Ul z;
    read(n);read(m);read(k);
    for(int i=1;i<=n;i++)read(t[i].opt),t[i].opt--,read(t[i].v),t[i].Update();
    for(int i=1;i<n;i++)read(x),read(y),link(t+x,t+y); 
    while(m--){
        read(type);read(x);read(y);read(z);
        if(type&1){
            MakeRoot(t+x);Access(t+y);Splay(t+y);
            Ul Fz=t[y].Fz,Fe=t[y].Fe,Ans=0;
            for(int i=k-1;~i;i--){
                Ul P=(Ul)1<<i; 
                if((Fe&P)&&!(Fz&P)&&(z>=P))Ans|=P,z-=P;
                else Ans|=Fz&P;
            }
            printf("%llu\n",Ans);
        }
        else{
            Splay(t+x);t[x].opt=y-1;t[x].v=z;
            t[x].Update();
        }
    }
    return 0;
}
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值