BZOJ2636: crisis(可持久化线段树)

传送门：

解题思路：

题目描述是一大坑点，cancel后面是直接加ask或者redo的。

那么就可以愉快地可持久化了。

注意需要支持区间修改，那么就只需要在再次更新这个点的时候将标记储存在新的儿子中。

最后由下至上询问就好了。

代码：

 #include<cmath>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 typedef double dnt;
 const double eps=1e-;
 struct trnt{
     int ls;
     int rs;
     double dltx,dlty;
     double alpha;
     bool lunked;
 }tr[];
 struct pnt{
     double x,y;
     void insert(void)
     {
         scanf("%lf%lf",&x,&y);
         return ;
     }
 }p[];
 int siz;
 int cnt;
 int n,m;
 char cmd[];
 int root[];
 void move(int spc,double dx,double dy)
 {
     tr[spc].dltx+=dx;
     tr[spc].dlty+=dy;
     return ;
 }
 void patrol(int spc,double beta)
 {
     tr[spc].alpha+=beta;
     double nx,ny,x,y;
     x=tr[spc].dltx,y=tr[spc].dlty;
     nx=x*cos(beta)-y*sin(beta);
     ny=y*cos(beta)+x*sin(beta);
     tr[spc].dltx=nx;
     tr[spc].dlty=ny;
     return ;
 }
 void lunk(int spc)
 {
     tr[spc].alpha=tr[spc].dltx=tr[spc].dlty=;
     tr[spc].lunked=true;
     return ;
 }
 void pushdown(int spc,int lst)
 {
     tr[++siz]=tr[tr[lst].ls];
     tr[spc].ls=siz;
     tr[++siz]=tr[tr[lst].rs];
     tr[spc].rs=siz;
     if(tr[spc].lunked)
     {
         lunk(tr[spc].ls);
         lunk(tr[spc].rs);
     }
     patrol(tr[spc].ls,tr[spc].alpha);
     patrol(tr[spc].rs,tr[spc].alpha);
     move(tr[spc].ls,tr[spc].dltx,tr[spc].dlty);
     move(tr[spc].rs,tr[spc].dltx,tr[spc].dlty);
     tr[spc].lunked=false;
     tr[spc].alpha=tr[spc].dltx=tr[spc].dlty=;
     return ;
 }
 void build(int l,int r,int &spc)
 {
     spc=++siz;
     if(l==r)return ;
     int mid=(l+r)>>;
     build(l,mid,tr[spc].ls);
     build(mid+,r,tr[spc].rs);
     return ;
 }
 void Move(int l,int r,int ll,int rr,int spc,int lst,double dx,double dy)
 {
     if(!spc)return ;
     if(ll>r||l>rr)return ;
     if(ll<=l&&r<=rr)
     {
         move(spc,dx,dy);
         return ;
     }
     int mid=(l+r)>>;
     pushdown(spc,lst);
     Move(l,mid,ll,rr,tr[spc].ls,tr[lst].ls,dx,dy);
     Move(mid+,r,ll,rr,tr[spc].rs,tr[lst].rs,dx,dy);
     return ;
 }
 void Patrol(int l,int r,int ll,int rr,int spc,int lst,double beta)
 {
     if(!spc)return ;
     if(ll>r||l>rr)return ;
     if(ll<=l&&r<=rr)
     {
         patrol(spc,beta);
         return ;
     }
     int mid=(l+r)>>;
     pushdown(spc,lst);
     Patrol(l,mid,ll,rr,tr[spc].ls,tr[lst].ls,beta);
     Patrol(mid+,r,ll,rr,tr[spc].rs,tr[lst].rs,beta);
     return ;
 }
 void Lunk(int l,int r,int ll,int rr,int spc,int lst)
 {
     if(!spc)return ;
     if(ll>r||l>rr)return ;
     if(ll<=l&&r<=rr)
     {
         lunk(spc);
         return ;
     }
     int mid=(l+r)>>;
     pushdown(spc,lst);
     Lunk(l,mid,ll,rr,tr[spc].ls,tr[lst].ls);
     Lunk(mid+,r,ll,rr,tr[spc].rs,tr[lst].rs);
     return ;
 }
 void query(int l,int r,int spc,int pos,double &x,double &y)
 {
     if(!spc)return ;
     if(l>r)return ;
     int mid=(l+r)>>;
     if(pos<=mid)query(l,mid,tr[spc].ls,pos,x,y);
     else        query(mid+,r,tr[spc].rs,pos,x,y);
     if(tr[spc].lunked)x=y=;
     if(fabs(tr[spc].alpha)>eps)
     {
         double nx,ny;
         nx=x*cos(tr[spc].alpha)-y*sin(tr[spc].alpha);
         ny=y*cos(tr[spc].alpha)+x*sin(tr[spc].alpha);
         x=nx,y=ny;
     }
     x+=tr[spc].dltx;
     y+=tr[spc].dlty;
     return ;
 }
 int main()
 {
 //    freopen("zoo.in","r",stdin);
 //    freopen("zoo.out","w",stdout);
     scanf("%d",&n);
     for(int i=;i<=n;i++)p[i].insert();
     scanf("%d",&m);
     cnt=;
     build(,n,root[cnt]);
     while(m--)
     {
         scanf("%s",cmd+);
         if(cmd[]=='M')
         {
             int l,r;
             scanf("%d%d",&l,&r);
             if(l>r)std::swap(l,r);
             double dx,dy;
             scanf("%lf%lf",&dx,&dy);
             tr[++siz]=tr[root[cnt]];
             root[++cnt]=siz;
             Move(,n,l,r,root[cnt],root[cnt-],dx,dy);
         }
         if(cmd[]=='P')
         {
             int l,r;
             scanf("%d%d",&l,&r);
             double beta;
             scanf("%lf",&beta);
             if(l>r)std::swap(l,r);
             tr[++siz]=tr[root[cnt]];
             root[++cnt]=siz;
             Patrol(,n,l,r,root[cnt],root[cnt-],beta);
         }
         if(cmd[]=='L')
         {
             int l,r;
             scanf("%d%d",&l,&r);
             if(l>r)std::swap(l,r);
             tr[++siz]=tr[root[cnt]];
             root[++cnt]=siz;
             Lunk(,n,l,r,root[cnt],root[cnt-]);
         }
         if(cmd[]=='C')
         {
             int a;
             scanf("%d",&a);
             cnt-=a;
         }
         if(cmd[]=='R')
         {
             int a;
             scanf("%d",&a);
             cnt+=a;
         }
         if(cmd[]=='A')
         {
             int pos;
             scanf("%d",&pos);
             double x=p[pos].x,y=p[pos].y;
             query(,n,root[cnt],pos,x,y);
             printf("%.6lf %.6lf\n",x,y);
         }
     }
     return ;
 }