#include <bits/stdc++.h>

using namespace std;

using i64=long long;
using i128=__int128_t;
using u32=unsigned int;
using u64=unsigned long long;
using db=double;
using pii=pair<int,int>;
using vi=vector<int>;

template<typename T>
inline T read(){
    T x=0,f=0;char ch=getchar();
    while(!isdigit(ch)) f|=(ch=='-'),ch=getchar();
    while(isdigit(ch)) x=x*10+(ch-'0'),ch=getchar();
    return f?-x:x;
}
template<typename T>
inline void write(T x){
    if(x<0) {putchar('-'),x=-x;}
    static int bit[50];
    int tp=0;
    do bit[++tp]=x%10,x/=10; while(x);
    while(tp) putchar(bit[tp--]+'0');
}

#define rdi read<int>
#define rdi64 read<i64>
#define fi first
#define se second
#define pb push_back

const int N=1e5+10,MX=1e9;

int n,q,a[N];
vi e[N];
i128 sum2;

inline int lg2(int x) {return !x?-1:__lg(x);}

int f[N][18],dep[N];
int dfn[N],low[N],ti;
void dfs(int x,int fa){
    dfn[x]=++ti,dep[x]=dep[fa]+1;
    f[x][0]=fa;
    for(int i=1;i<=lg2(dep[x]);i++) f[x][i]=f[f[x][i-1]][i-1];
    for(auto y:e[x]){
        if(y==fa) continue;
        dfs(y,x);
    }
    low[x]=ti;
}

int lca(int x,int y){
    if(dep[x]<dep[y]) swap(x,y);
    for(int i=lg2(dep[x]-dep[y]);i>=0;i--)
        if(dep[f[x][i]]>=dep[y]) x=f[x][i];
    if(x==y) return x;
    for(int i=lg2(dep[x]);i>=0;i--)
        if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
    return f[x][0];
}

struct Data{
    int s0;i64 s1;i128 s2;
    Data() {s0=s1=s2=0;}
    Data(int s0,i64 s1,i128 s2):s0(s0),s1(s1),s2(s2){}
    Data(int v){
        if(v>=0) s0=1,s1=v,s2=(i64)v*v;
        else s0=-1,s1=v,s2=(i64)-v*v;
    }
    Data operator + (const Data &rhs) const{
        return {s0+rhs.s0,s1+rhs.s1,s2+rhs.s2};
    }
    Data operator - (const Data &rhs) const{
        return {s0-rhs.s0,s1-rhs.s1,s2-rhs.s2};
    }
    Data operator * (int b) const{
        return {s0*b,s1*b,s2*b};
    }
};

int lst;

struct Op{
    int op,x,y;
    i64 z;
}p[N];

struct SGT{
#define lson (t[nw].ls)
#define rson (t[nw].rs)
#define mid ((l+r)>>1)
    struct Node{int ls,rs;Data d;}t[N*64];
    int tot;
    void add(int &nw,int pr,int l,int r,int x,int v){
        t[nw=++tot]=t[pr],t[nw].d=t[nw].d+v;
        if(l==r) return;
        if(x<=mid) add(lson,t[pr].ls,l,mid,x,v);
        else add(rson,t[pr].rs,mid+1,r,x,v);
    }
    i128 query(int rt1,int rt2,int rt3,int rt4,int l,int r,
            i64 lim,int x,int y,int lc,int tid,
            Data d={0,0,0},i128 sum=0){
        auto qinfo = [&](int typ){
            auto insub = [&](int x,int y){
                return dfn[x]<=dfn[y]&&dfn[y]<=low[x];
            };
            auto on_path = [&](int u){
                return insub(lc,u)&&(insub(u,x)||insub(u,y));
            };

            Data ret;int l1,r1;
            if(typ==1) ret=t[rt1].d+t[rt2].d-t[rt3].d-t[rt4].d,l1=l,r1=r;
            else if(typ==2) ret=t[t[rt1].ls].d+t[t[rt2].ls].d-t[t[rt3].ls].d-t[t[rt4].ls].d,l1=l,r1=mid;
            else if(typ==3) ret=t[t[rt1].rs].d+t[t[rt2].rs].d-t[t[rt3].rs].d-t[t[rt4].rs].d,l1=mid+1,r1=r;
            for(int i=lst+1;i<=tid;i++)
                if(p[i].op==1&&on_path(p[i].x)){
                    if(l1<=p[i].y&&p[i].y<=r1) ret=ret-p[i].y;
                    if(l1<=p[i].z&&p[i].z<=r1) ret=ret+p[i].z;
                }
            return ret;
        };
        if(l==r){
            d=d+qinfo(1);
            i64 rest=lim-(d.s1-(i64)l*d.s0);
            i128 res=(i128)l*l*d.s0+sum;
            if(l) res-=(i128)(2*l-1)*rest;
            return res+(sum2-d.s2-sum);
        }
        Data tmp=qinfo(3),nxt=tmp+d;
        if(nxt.s1-(i64)mid*nxt.s0<=lim) return query(t[rt1].ls,t[rt2].ls,t[rt3].ls,t[rt4].ls,l,mid,lim,x,y,lc,tid,nxt,sum);
        else return query(t[rt1].rs,t[rt2].rs,t[rt3].rs,t[rt4].rs,mid+1,r,lim,x,y,lc,tid,d,sum+qinfo(2).s2);
    }
    void clear() {tot=0;}
#undef lson
#undef rson
#undef mid
}t;

int rt[N];
void rebuild(int x,int fa){
    t.add(rt[x],rt[fa],0,MX,a[x],a[x]);
    for(auto y:e[x]){
        if(y==fa) continue;
        rebuild(y,x);
    }
}

void init(){
    n=rdi(),q=rdi();
    for(int i=1;i<n;i++){
        int x=rdi(),y=rdi();
        e[x].pb(y),e[y].pb(x);
    }
    for(int i=1;i<=n;i++) a[i]=rdi();

    dfs(1,0);
    for(int i=1;i<=n;i++)
        sum2+=(i64)a[i]*a[i];

    static int a1[N];
    copy(a+1,a+n+1,a1+1);
    for(int i=1;i<=q;i++){
        p[i].op=rdi();
        if(p[i].op==1){
            p[i].x=rdi();
            p[i].y=a[p[i].x],p[i].z=rdi();
            a[p[i].x]=p[i].z;
        }
        else{
            p[i].x=rdi(),p[i].y=rdi(),p[i].z=rdi64();
        }
    }
    copy(a1+1,a1+n+1,a+1);
}

i128 query(int x,int y,i64 k,int tid){
    int lc=lca(x,y);
    return t.query(rt[x],rt[y],rt[lc],rt[f[lc][0]],0,MX,k,x,y,lc,tid);
}

void solve(){
    init();
    int bsiz=sqrt(n)*3.75;
    for(int l=1,r;l<=q;l=r+1){
        r=min(q,l+bsiz-1);
        t.clear(),rebuild(1,0),lst=l-1;
        for(int i=l;i<=r;i++){
            if(p[i].op==1) sum2=sum2-(i64)p[i].y*p[i].y+(i64)p[i].z*p[i].z,a[p[i].x]=p[i].z;
            else write(query(p[i].x,p[i].y,p[i].z,i)),putchar('\n');
        }
    }
}

int main(){
#ifdef LOCAL
    freopen("1.in","r",stdin);
    freopen("1.out","w",stdout);
#endif
    solve();
    return 0;
}