第一次写LCT,各种模板加入。。。以后都只遇到有新意的题目再更新了

这道题就是LCT,但是,难在一个回退的操作。这时,可以通过改变执行顺序,先把要回退后再做的操作先执行了,再回退到之前的执行。这时,建立一棵操作树,使用DFS可以回退的特点,就能完成这一功能。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <set>
#pragma comment(linker, "/STACK:1024000000,1024000000")
#define LL long long
using namespace std; const int MAX=300050;
const LL MOD=1e9+7;
int fa[MAX],ch[MAX][2],son[MAX],nxt[MAX],rev[MAX];
int n,q;
int opt_x[MAX],opt_u[MAX],opt_v[MAX],opt_d[MAX];
LL val[MAX],mul[MAX];
LL ans[MAX];
set<long long>st;
int stack[MAX]; void update(int x){
int l=ch[x][0],r=ch[x][1];
mul[x]=val[x];
if(l) mul[x]=mul[x]*mul[l]%MOD;
if(r) mul[x]=mul[x]*mul[r]%MOD;
} void pushdown(int x){
int l=ch[x][0],r=ch[x][1];
if(rev[x]){
rev[x]^=1;rev[l]^=1;rev[r]^=1;
swap(ch[x][0],ch[x][1]);
}
} bool isroot(int x){
return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x;
} void rotate(int x){
int y=fa[x],z=fa[y],l,r;
if(ch[y][0]==x)l=0;else l=1;r=l^1;
if(!isroot(y))
{
if(ch[z][0]==y)ch[z][0]=x;
else ch[z][1]=x;
}
fa[x]=z;fa[y]=x;fa[ch[x][r]]=y;
ch[y][l]=ch[x][r];ch[x][r]=y;
update(y);update(x);
} void splay(int u){
int top=0; stack[++top]=u;
for(int i=u;!isroot(i);i=fa[i]) stack[++top]=fa[i];
while(top) pushdown(stack[top--]);
while(!isroot(u)){
int y=fa[u],z=fa[y];
if(!isroot(y)){
if(ch[y][0]==u^ch[z][0]==y)rotate(u);
else rotate(y);
}
rotate(u);
}
} bool connect(int u,int v){
while(fa[u]) u=fa[u];
while(fa[v]) v=fa[v];
return u==v;
} void access(int u){
for(int t=0;u;t=u,u=fa[u]){
splay(u),ch[u][1]=t,update(u);
}
} void makeroot(int x){
access(x),splay(x),rev[x]^=1;
} bool link(int u,int v){
if(u>v) swap(u,v);
if(u==v) return false;
if(connect(u,v)) return false;
makeroot(u),fa[u]=v;
st.insert(MAX*u+v);
return true;
} bool cut(int u,int v){
if(u>v) swap(u,v);
if(u==v) return false;
if(st.find(MAX*u+v)==st.end()) return false;
makeroot(u),access(v),splay(v);
fa[u]=0,ch[v][0]=0;
st.erase(MAX*u+v);
return true;
} void dfs(int x){
bool done=false;
LL tmp;
int u=opt_u[x],v=opt_v[x];
if(opt_x[x]==1) done=link(u,v);
if(opt_x[x]==2) done=cut(u,v);
if(opt_x[x]==4){
if(u==v) ans[opt_d[x]]=val[u];
else{
if(connect(u,v)){
makeroot(u);
access(v),splay(v);
ans[opt_d[x]]=mul[v];
}
else{
ans[opt_d[x]]=0;
}
}
}
if(opt_x[x]==5){
done=true;
splay(u);
tmp=val[u];
val[u]=v*(v-1)/2%MOD;
update(u);
}
for(int i=son[x];i;i=nxt[i]) dfs(i);
if(done){
if(opt_x[x]==1) done=cut(u,v);
if(opt_x[x]==2) done=link(u,v);
if(opt_x[x]==5){
splay(u);
val[u]=tmp;
update(u);
}
}
} int main(){
int T;
LL k;
scanf("%d",&T);
while(T--){
scanf("%d%d",&n,&q);
for(int i=1;i<=n;i++){
scanf("%lld",&k);
val[i]=mul[i]=k*(k-1)/2%MOD;
}
for(int i=0;i<=n+10;i++) rev[i]=0,fa[i]=0,ch[i][0]=ch[i][1]=0;
for(int i=0;i<=q+10;i++) nxt[i]=son[i]=0;
st.clear();
int tot=0;
for(int i=1;i<=q;i++){
scanf("%d",&opt_x[i]);
if(opt_x[i]==3){
scanf("%d",&opt_u[i]);
nxt[i]=son[opt_u[i]],son[opt_u[i]]=i;
}
else{
scanf("%d%d",&opt_u[i],&opt_v[i]);
nxt[i]=son[i-1],son[i-1]=i;
if(opt_x[i]==4){
opt_d[i]=++tot;
}
}
}
dfs(0);
for(int i=1;i<=tot;i++){
printf("%lld\n",ans[i]);
}
}
return 0;
}

  

HDU 5467的更多相关文章

  1. HDU 1284 钱币兑换问题 母函数、DP

    题目链接:HDU 1284 钱币兑换问题 钱币兑换问题 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (J ...

  2. HDOJ 2111. Saving HDU 贪心 结构体排序

    Saving HDU Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total ...

  3. 【HDU 3037】Saving Beans Lucas定理模板

    http://acm.hdu.edu.cn/showproblem.php?pid=3037 Lucas定理模板. 现在才写,noip滚粗前兆QAQ #include<cstdio> #i ...

  4. hdu 4859 海岸线 Bestcoder Round 1

    http://acm.hdu.edu.cn/showproblem.php?pid=4859 题目大意: 在一个矩形周围都是海,这个矩形中有陆地,深海和浅海.浅海是可以填成陆地的. 求最多有多少条方格 ...

  5. HDU 4569 Special equations(取模)

    Special equations Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u S ...

  6. HDU 4006The kth great number(K大数 +小顶堆)

    The kth great number Time Limit:1000MS     Memory Limit:65768KB     64bit IO Format:%I64d & %I64 ...

  7. HDU 1796How many integers can you find(容斥原理)

    How many integers can you find Time Limit:5000MS     Memory Limit:32768KB     64bit IO Format:%I64d ...

  8. hdu 4481 Time travel(高斯求期望)(转)

    (转)http://blog.csdn.net/u013081425/article/details/39240021 http://acm.hdu.edu.cn/showproblem.php?pi ...

  9. HDU 3791二叉搜索树解题(解题报告)

    1.题目地址: http://acm.hdu.edu.cn/showproblem.php?pid=3791 2.参考解题 http://blog.csdn.net/u013447865/articl ...

随机推荐

  1. 微信小程序获取地理位置

    小程序只支持获取当前位置的经纬度,并不能直接获取到地理名称,需要通过第三方来逆地址解析,这里我选择的是腾讯位置服务 在使用前需要去申请key,这里是地址:https://lbs.qq.com/cons ...

  2. 关于swoole 和golang 的压力测试结果

    一.环境介绍 linux centos7 php7.1.18 go1.12.1 2核4G内存 二.代码 swoole代码 <?php $http = new swoole_http_server ...

  3. $P3931 SAC E一道难题 Tree$

    problem #include <bits/stdc++.h> #define rep(i,j,n) for(register int i=j;i<=n;i++) #define ...

  4. WP8开发常用解决方案收集

    我其实不怎么做wp的东西.但是偶尔还是会用到, 但是wp8开发的资料确实难找.特开此贴,记录一些常见的解决方案 1.水平滑动动画(比如app首次使用说明就可以用这个做) http://www.cnbl ...

  5. scala 变量定义,基本操作符

    scala是jvm语言,它将面向对象与函数风格结合的很好,它可以访问任何java类库并很好的结合使用. 它可以使程序更简单,同时可利用并发的威力. scala基本语法: package com.tes ...

  6. NHibernate学习(零)-本次学习遇到的错误汇总

    问题一: "System.TypeInitializationException"类型的未经处理的异常在 KimismeDemo.exe 中发生 其他信息: "NHibe ...

  7. Select2插件ajax方式加载数据并刷新页面数据回显

    今天在优化项目当中,有个要在下拉框中搜索数据的需求:最后选择使用selec2进行开发: 官网:http://select2.github.io/ 演示: 准备工作: 文件需要引入select2.ful ...

  8. Product Device Lot

    Product是指产品: 这个Product可以在不同的设备类型上生产, 同一类型的设备也可能硬件有差异,所以会有相对应的Device(Recipe): 同一Product(或同一Device)由于数 ...

  9. TCP协议滑动窗口(一)——控制大批量数据传输速率

    窗口大小:TCP头中一个16位的域,表示当前可用接受缓冲区大小.在每个TCP对等段连接初始化时,告诉对方自己的窗口大小(不一定是满额,假如满额65201字节,可能暂时通告5840字节).若客户端接受数 ...

  10. html5——背景

    背景大小 background-size: % %;//宽 高 background-size: 100px 100px;//宽 高 background-size: cover;//覆盖,图片会完全 ...