如果只有第3个操作，那么这就是个sd题，随便lct搞搞就过去了
然后就是一个神仙东西 taylor公式
我不会，看gsy博客https://www.cnblogs.com/zhoushuyu/p/8148732.html
#include<bits/stdc++.h>
#define il inline
#define vd void
typedef long long ll;
#define maxn 100010
#define k 16
il int gi(){
int x=0,f=1;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-')f=-1;
ch=getchar();
}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return x*f;
}
int F[maxn];double A[maxn],B[maxn];
int ch[maxn][2],fa[maxn];
double f[maxn][k];
bool rev[maxn];
il bool isrt(int x){return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x;}
il vd Rev(int x){if(x)rev[x]^=1,std::swap(ch[x][0],ch[x][1]);}
il vd down(int x){
if(!isrt(x))down(fa[x]);
if(rev[x])Rev(ch[x][0]),Rev(ch[x][1]),rev[x]=0;
}
il vd upd(int x){
for(int i=0;i<16;++i)f[x][i]=f[ch[x][0]][i]+f[ch[x][1]][i];
if(F[x]==1){
double val=1,Sin=sin(B[x]),Cos=cos(B[x]);
for(int i=0;i<k;i+=4){
f[x][i]+=Sin*val,val*=A[x];
f[x][i+1]+=Cos*val,val*=A[x];
f[x][i+2]-=Sin*val,val*=A[x];
f[x][i+3]-=Cos*val,val*=A[x];
}
}else if(F[x]==2){
double val=exp(B[x]);
for(int i=0;i<k;++i)f[x][i]+=val,val*=A[x];
}else f[x][0]+=B[x],f[x][1]+=A[x];
}
il vd rotate(int x){
int y=fa[x],z=fa[y],o=ch[y][1]==x;
fa[x]=z;if(!isrt(y))ch[z][y==ch[z][1]]=x;
ch[y][o]=ch[x][!o];fa[ch[x][!o]]=y;
ch[x][!o]=y,fa[y]=x;
upd(y);
}
il vd splay(int x){
down(x);
for(int y=fa[x],z=fa[y];!isrt(x);rotate(x),y=fa[x],z=fa[y])
if(!isrt(y))rotate(((ch[y][1]==x)==(ch[z][1]==y))?y:x);
upd(x);
}
il vd access(int x){for(int y=0;x;x=fa[y=x])splay(x),ch[x][1]=y,upd(x);}
il vd makert(int x){access(x),splay(x),Rev(x);}
il vd link(int x,int y){makert(x),fa[x]=y;}
il vd cut(int x,int y){makert(x),access(y),splay(y),ch[y][0]=fa[x]=0,upd(y);}
il int find(int x){access(x);splay(x);while(ch[x][0])x=ch[x][0];return x;}
double p[k];
int main(){
#ifndef ONLINE_JUDGE
freopen("4546.in","r",stdin);
freopen("4546.out","w",stdout);
#endif
p[0]=1;for(int i=1;i<k;++i)p[i]=p[i-1]*i;
int n=gi(),m=gi();scanf("%*s");
for(int i=1;i<=n;++i)F[i]=gi(),scanf("%lf%lf",&A[i],&B[i]),upd(i);
char opt[10];
int u,v;
while(m--){
scanf("%s",opt);
if(opt[0]=='a')u=gi(),v=gi(),link(u+1,v+1);
else if(opt[0]=='d')u=gi(),v=gi(),cut(u+1,v+1);
else if(opt[0]=='m')u=gi()+1,makert(u),F[u]=gi(),scanf("%lf%lf",&A[u],&B[u]),upd(u);
else{
u=gi()+1,v=gi()+1;
double X,ans=0,pX;scanf("%lf",&X);pX=1;
if(find(u)!=find(v)){puts("unreachable");continue;}
makert(u),access(v),splay(v);
for(int i=0;i<k;++i)ans+=f[v][i]*pX/p[i],pX*=X;
printf("%.8e\n",ans);
}
}
return 0;
}
虽然我不会，但我在BZOJ上的时间吊打了gsy