[THUWC2017]在美妙的数学王国中畅游

[THUWC2017]在美妙的数学王国中畅游

e和sin信息不能直接合并

泰勒展开，大于21次太小，认为是0，保留前21次多项式即可

然后就把e，sin ，kx+b都变成多项式了，pushup合并

上LCT

// luogu-judger-enable-o2
// luogu-judger-enable-o2
#include<bits/stdc++.h>
#define reg register int
#define il inline
#define fi first
#define se second
#define mk(a,b) make_pair(a,b)
#define numb (ch^'0')
#define ld double
using namespace std;
typedef long long ll;
template<class T>il void rd(T &x){
    char ch;x=;bool fl=false;
    while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true);
    for(x=numb;isdigit(ch=getchar());x=x*+numb);
    (fl==true)&&(x=-x);
}
template<class T>il void output(T x){if(x/)output(x/);putchar(x%+'');}
template<class T>il void ot(T x){if(x<) putchar('-'),x=-x;output(x);putchar(' ');}
template<class T>il void prt(T a[],int st,int nd){for(reg i=st;i<=nd;++i) ot(a[i]);putchar('\n');}

namespace Miracle{
const int N=1e5+;
const int K=;
int n,m;
struct Poly{
    ld a[K];
    Poly(){memset(a,,sizeof a);}
    il ld &operator[](const int &x){return a[x];}
    il const ld &operator[](const int &x) const {return a[x];}
    Poly friend operator +(const Poly &a,const Poly &b){
        Poly c;
        for(reg i=;i<K;++i){
            c.a[i]=a.a[i]+b.a[i];
        }
        return c;
    }
    ld calc(ld x){
        ld tmp=,ret=;
        for(reg i=;i<K;++i){
            ret=ret+a[i]*tmp;
            tmp=tmp*x;
        }
        return ret;
    }
};
int C[][];
ld jie[K],ma[K],mb[K];
Poly exp(double a,double b){
    Poly ret;
    ma[]=;mb[]=;
    for(reg i=;i<K;++i){
        ma[i]=ma[i-]*a;
        mb[i]=mb[i-]*b;
    }
    for(reg i=;i<K;++i){
        for(reg j=;j<=i;++j){
            ret.a[j]+=C[i][j]*ma[j]*mb[i-j]*jie[i];
        }
    }
    return ret;
}
Poly sin(double a,double b){
    Poly ret;
    ma[]=;mb[]=;
    for(reg i=;i<K;++i){
        ma[i]=ma[i-]*a;
        mb[i]=mb[i-]*b;
    }
    for(reg i=;i<K;i+=){
        for(reg j=;j<=i;++j){
            if(((i-)/)&){
                ret.a[j]-=C[i][j]*ma[j]*mb[i-j]*jie[i];
            }else{
                ret.a[j]+=C[i][j]*ma[j]*mb[i-j]*jie[i];
            }
        }
    }
    return ret;
}
Poly init(int typ,double a,double b){
    if(typ==) return sin(a,b);
    if(typ==) return exp(a,b);
    Poly ret;
    ret.a[]=a;
    ret.a[]=b;
    return ret;
}
struct node{
    int ch[],fa;
    int rev;
    Poly f,sum;
}t[N];
#define ls t[x].ch[0]
#define rs t[x].ch[1]
void pushup(int x){
    if(!x) return;
    t[x].sum=t[ls].sum+t[rs].sum+t[x].f;
}
bool nrt(int x){
    return t[t[x].fa].ch[]==x||t[t[x].fa].ch[]==x;
}
void rev(int x){
    t[x].rev^=;
    swap(ls,rs);
}
void pushdown(int x){
    if(t[x].rev){
        rev(ls);rev(rs);
        t[x].rev=;
    }
}
void rotate(int x){
    int y=t[x].fa,d=t[y].ch[]==x;
    t[t[y].ch[d]=t[x].ch[!d]].fa=y;
    if(nrt(y)) t[t[x].fa=t[y].fa].ch[t[t[y].fa].ch[]==y]=x;
    else t[x].fa=t[y].fa;
    t[t[x].ch[!d]=y].fa=x;
    pushup(y);
}
int sta[N];
void splay(int x){
    int y=x,z=;
    sta[++z]=x;
    while(nrt(y)) y=t[y].fa,sta[++z]=y;
    while(z) pushdown(sta[z--]);
    while(nrt(x)){
        y=t[x].fa,z=t[y].fa;
        if(nrt(y)){
            rotate((t[y].ch[]==x)==(t[z].ch[]==y)?y:x);
        }
        rotate(x);
    }
    pushup(x);
}
void access(int x){
    for(reg y=;x;y=x,x=t[x].fa){
        splay(x);t[x].ch[]=y;pushup(x);
    }
}
void makert(int x){
    access(x);splay(x);rev(x);
}
void link(int x,int y){
    makert(x);t[x].fa=y;
}
void split(int x,int y){
    makert(x);access(y);splay(y);
}
void cut(int x,int y){
    split(x,y);
    t[y].ch[]=t[x].fa=;
    pushup(y);
}
int findrt(int x){
    access(x);splay(x);
    while(ls) x=ls;
    splay(x);
    return x;
}
ld query(int x,int y,ld v){
    split(x,y);
    return t[y].sum.calc(v);
}
char s[];
int main(){
    char haha[];
    rd(n);rd(m);scanf("%s",haha+);

    int typ;
    ld a,b;
    int x,y;
    C[][]=;
    for(reg i=;i<=K;++i){
        C[i][]=;
        for(reg j=;j<=K;++j){
            C[i][j]=C[i-][j]+C[i-][j-];
        }
    }
    jie[]=jie[]=;
    for(reg i=;i<=K;++i){
        jie[i]=jie[i-]/i;
    }
    for(reg i=;i<=n;++i){
        rd(typ);scanf("%lf%lf",&a,&b);
        t[i].f=init(typ,a,b);
        pushup(i);
    }
    while(m--){
        scanf("%s",s+);
        if(s[]=='a'){
            rd(x);rd(y);
            ++x;++y;
            link(x,y);
        }else if(s[]=='d'){
            rd(x);rd(y);
            ++x;++y;
            cut(x,y);
        }else if(s[]=='m'){
            rd(x);rd(typ);scanf("%lf%lf",&a,&b);
            ++x;
            access(x);splay(x);
            t[x].f=init(typ,a,b);
            pushup(x);
        }else {
//            cout<<" tr "<<endl;
            rd(x);rd(y);scanf("%lf",&a);
            ++x;++y;
//            cout<<x<<" "<<y<<" a "<<a<<endl;
            if(findrt(x)==findrt(y)) printf("%.10lf\n",query(x,y,a));
            else puts("unreachable");
        }
    }
    return ;
}

}
signed main(){
    Miracle::main();
    return ;
}

/*
   Author: *Miracle*
   Date: 2019/4/11 20:11:40
*/

多项式可以支持合并