线段树+欧拉函数——cf1114F

调了半天，写线段树老是写炸

/*
两个操作
1.区间乘法
2.区间乘积询问欧拉函数 

欧拉函数计算公式
    phi(mul(ai))=mul(ai) * (p1-1)/p1 * (p2-1)/p2 * .. * (pk-1)/pk
因为只有300以内的质数（62个）用一个long long来状态压缩
因此线段树结点维护住区间的质数状态集合S，区间的乘积

操作1 [l,r] x：把x质因数分解，然后更新S，然后再更新乘积,
操作2 [l,r]：询问到区间的状态集合S，区间的乘积，再求逆元进行除法

先把62个质数的逆元求出来 

线段树结点维护区间乘积，区间质数集合S，然后两个lazy标记
*/
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define mod 1000000007
#define maxn 400005
ll n,q;
ll p[]={,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
    ,,,,,,,,,,,,,,,,,,,,,,,,};
ll invp[];
ll Pow(ll a,ll b){
    ll res=;
    while(b){
        if(b%)res=res*a%mod;
        b>>=;a=a*a%mod;
    }
    return res;
}

ll a[maxn];

#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
//mul是区间乘积，l_mul是区间乘的lazy
ll msk[maxn<<],mul[maxn<<],l_msk[maxn<<],l_m[maxn<<];
inline ll calc(ll x){
    ll res=;
    for(ll i=;i<;i++)
        if(x%p[i]==) res|=((ll)<<i);
    return res;
}
void pushup(int rt){
    mul[rt]=mul[rt<<] * mul[rt<<|] % mod;
    msk[rt]=msk[rt<<] | msk[rt<<|];
}
void pushdown(int l,int r,int rt){
    if(l_msk[rt]!=){
        l_msk[rt<<] |= l_msk[rt]; msk[rt<<] |= l_msk[rt];
        l_msk[rt<<|] |= l_msk[rt]; msk[rt<<|] |= l_msk[rt];
        l_msk[rt]=;
    }
    if(l_m[rt]!=){//整个区间都要乘
        int m=l+r>>;
        l_m[rt<<]=l_m[rt<<] * l_m[rt]%mod;
        mul[rt<<]=Pow(l_m[rt],m-l+) * mul[rt<<]%mod;
        l_m[rt<<|]=l_m[rt<<|] * l_m[rt]%mod;
        mul[rt<<|]=Pow(l_m[rt],r-m) * mul[rt<<|]%mod;
        l_m[rt]=;
    }
}
void build(int l,int r,int rt){
    l_msk[rt]=;l_m[rt]=;
    if(l==r){
        mul[rt]=a[l];
        msk[rt]=calc(a[l]);
        return ;
    }
    int m=l+r>>;
    build(lson);build(rson);
    pushup(rt);
}
void up_mul(int L,int R,ll x,int l,int r,int rt){
    if(L<=l && R>=r){
        l_m[rt]=l_m[rt]*x%mod;
        mul[rt]=mul[rt]*Pow(x,r-l+)%mod;
        return;
    }
    pushdown(l,r,rt);
    int m=l+r>>;
    if(L<=m)up_mul(L,R,x,lson);
    if(R>m)up_mul(L,R,x,rson);
    pushup(rt);
}
void up_S(int L,int R,ll S,int l,int r,int rt){
    if(L<=l && R>=r){
        l_msk[rt]|=S;
        msk[rt]|=S;
        return;
    }
    pushdown(l,r,rt);
    int m=l+r>>;
    if(L<=m) up_S(L,R,S,lson);
    if(R>m) up_S(L,R,S,rson);
    pushup(rt);
}
ll q_mul(int L,int R,int l,int r,int rt){
    if(L<=l && R>=r)return mul[rt];
    pushdown(l,r,rt);
    int m=l+r>>;
    ll res=;
    if(L<=m)
        res=res*q_mul(L,R,lson)%mod;
    if(R>m)
        res=res*q_mul(L,R,rson)%mod;
    return res;
}
ll q_S(int L,int R,int l,int r,int rt){
    if(L<=l && R>=r)return msk[rt];
    pushdown(l,r,rt);
    int m=l+r>>;
    ll S=;
    if(L<=m)
        S|=q_S(L,R,lson);
    if(R>m)
        S|=q_S(L,R,rson);
    return S;
}
int main(){
    for(int i=;i<;i++)
        invp[i]=Pow(p[i],mod-);
    cin>>n>>q;
    for(int i=;i<=n;i++)scanf("%lld",&a[i]);

    build(,n,);

    char op[];ll l,r,x;
    while(q--){
        scanf("%s",op);
        if(op[]=='M'){
            scanf("%lld%lld%lld",&l,&r,&x);
            up_mul(l,r,x,,n,);
            ll S=calc(x);
            up_S(l,r,S,,n,);
        }
        else {
            scanf("%lld%lld",&l,&r);
            ll S=q_S(l,r,,n,);
            ll mul=q_mul(l,r,,n,);
            for(ll i=;i<;i++)
                if(S & ((ll)<<i))
                    mul=(mul*invp[i]%mod*(p[i]-)%mod);
            cout<<mul<<endl;
        }
    }
    return ;
}