第一次写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. golang——常用内建函数

    (1)func len(v Type) int 返回长度,取决于具体类型:字符串返回字节数:channel返回缓存元素的个数: (2)func cap(v Type) int 返回容量,取决于具体类型 ...

  2. .Net Core开源小工具mssql2mysql,从mssql生成mysql脚本

    Microsoft SQL Server to MySQL 这个工具用于从MSSQL生成MySQL脚本,生成的脚本包含表结构和数据 安装 这是一个.Net Core的具具,所以需要先安装.net co ...

  3. 实现一个类似于收趣APP的功能

    近日想做一个类似于收趣APP软件的一个功能,将头条.微信等其他App的文章能够通过分享微信好友的方式分享到自己的平台软件中. 分享的方式有三种: 1.通过微信好友的方式,将文章分享给收趣. 2.复制文 ...

  4. 研磨JavaScript系列(一):回归简单

    想要理解JavaScript,你得首先放下对象和类的概念,回到数据和代码的本原.编程世界只有数据和代码两种基本元素,而这两种元素又有着纠缠不清的关系.JavaScript就是把数据和代码都简化到最原始 ...

  5. 【RTTI】java Class类详解

    RTTI (Run-Time Type Information)运行时类信息 Java的Class类是java反射机制的基础,通过Class类我们可以获得关于一个类的相关信息,下面我们来了解一下有关j ...

  6. CSS——img标签消除3px

    1.dispaly:block 2.float:left 这两种都可以消除3px

  7. JS——dom

    节点的获取 <script> var div = document.getElementById("box");//返回指定标签 var div = document. ...

  8. php入门学习相关函数

      1.join(): 定义和用法 join() 函数返回由数组元素组合成的字符串. join() 函数是 implode() 函数的别名. 注释:join() 函数接受两种参数顺序.但是由于历史原因 ...

  9. 集合Set、List、Map的遍历方法

    package com.shellway.javase; import java.util.ArrayList; import java.util.Collection; import java.ut ...

  10. Linq处理decimal字段汇总Sum()为NULL

    xxxxxxxx.Sum(f => f.jifen).GetValueOrDefault(0)