标签:blog io for ar 问题 div amp 时间
(BZOJ挂了,还没在BZOJ测,先是在wikioi测过了,,)
囧。在军训时立志要学lct!!!这是一道lct的裸题,只有access操作(10行都没有啊亲。。。缩行大法的话,我就不说了。。)(link操作相当于水过),其实lct很简单。。想想都有点小激动。。。。。。
lct用splay维护的话,一下就写好了。。但是我在写lct的时候,发现了一些我原来splay的老问题,我原来也知道了的,就是将null的ch给赋值了,因为在rot操作里没有特判,所以导致了null的孩子被赋值了,导致我的lct调试了2天T_T
这不科学啊!!(但是想想了还好了,,因为lct在splay的过程中,是有父亲而父亲不一定有这个孩子的,因为lct维护的是森林。。所以在splay里有加特判,即父亲是否有这个孩子
还有个地方卡了我的调试,就是access操作,本来我是没有加ch[1]->f=null的,之前调试的时候以为是这里出问题,但是证明我很蠢。。。。这里一定不能这样搞啊!因为这样就将孩子和父亲断了联系。。。囧。
本题有个特性,即已经有序了的,我们可以从后往前看这个数组,想象成森林,某根的孩子就是在原数组的左边(深度比他深,即越后的元素深度越浅,因为lct在找preferred path的时候splay是按深度构图的,那么孩子就是在父亲的左边,而羊是从左往右走的,所以在更新的时候,不能断了左边的联系T_T,就是我之前的那个错误。。。。将孩子的父亲断了。。T_T。)
link-cut tree具体看论文吧,很容易看懂的,真心话。就是用一颗splay维护某个点的路径,即access操作,我相信多做题就行了的
。。找到错误真开心~调试出来了真开心~
有时间在学习完整套lct然后写篇专题。。囧。。最近看的数学真心感觉有点难啊,数论怎么那么hentai啊?公式推来推去就是为了将O(n)变成O(1)么。这东西能在考场上想到吗囧。好吧,我是蒟蒻所以想不到。
#include <cstdio> #define read(x) x=getint() using namespace std; inline int getint() { char c; int ret=0; for(c=getchar(); c<‘0‘ || c>‘9‘; c=getchar()); for(; c>=‘0‘ && c<=‘9‘; c=getchar()) ret=ret*10+c-‘0‘; return ret; } const int N=200005; struct node* null; struct node { node *f, *ch[2]; int s; void pushup() { s=1+ch[0]->s+ch[1]->s; } bool d() { return f->ch[1]==this; } void setc(node* c, bool d) { ch[d]=c; c->f=this; } bool check() { return f==null || (f->ch[0]!=this && f->ch[1]!=this); } //check操作是lct特地的,因为为了省空间,我们不开n颗splay,只用节点关联就行了,这样就无法避免父亲没有这个孩子的情况。所以还要特判 }*nd[N]; void rot(node* r) { node* f=r->f; bool d=r->d(); if(f->check()) r->f=f->f; else f->f->setc(r, f->d()); //这里一定要这样写,不然会让null的孩子改变,这样在splay的循环就会死循环T_T f->setc(r->ch[!d], d); r->setc(f, !d); f->pushup(); } inline void splay(node* r) { while(!r->check()) if(r->f->check()) rot(r); else r->d()==r->f->d()?(rot(r->f), rot(r)):(rot(r), rot(r)); r->pushup(); } inline void access(node* f) { for(node* c=null; f!=null; f=f->f) { splay(f); f->setc(c, 1); f->pushup(); c=f; } } inline void link(node* c, node* f) { access(c); splay(c); c->ch[0]->f=null; c->ch[0]=null; c->f=f; c->pushup(); } inline void init() { null=new node; null->s=0; null->f=null->ch[0]=null->ch[1]=null; } int main() { init(); int n, t; read(n); for(int i=0; i<n; ++i) { nd[i]=new node; nd[i]->s=1; nd[i]->ch[0]=nd[i]->ch[1]=nd[i]->f=null; } for(int i=0; i<n; ++i) { read(t); if(i+t<n) nd[i]->f=nd[i+t]; } int m, a, b, c; read(m); for(int i=0; i<m; ++i) { read(a); read(b); if(a==1) { access(nd[b]); splay(nd[b]); printf("%d\n", nd[b]->s); } else { read(c); if(b+c<n) link(nd[b], nd[b+c]); else link(nd[b], null); } } return 0; }
【BZOJ】2002: [Hnoi2010]Bounce 弹飞绵羊(lct),布布扣,bubuko.com
【BZOJ】2002: [Hnoi2010]Bounce 弹飞绵羊(lct)
标签:blog io for ar 问题 div amp 时间
原文地址:http://www.cnblogs.com/iwtwiioi/p/3912395.html