CF1114F Please, another Queries on Array?（线段树，数论，欧拉函数，状态压缩）

这题我在考场上也是想出了正解的……但是没调出来。

题目链接：CF原网

题目大意：给一个长度为 $n$ 的序列 $a$，$q$ 个操作：区间乘 $x$，求区间乘积的欧拉函数模 $10^9+7$ 的值。

$1\le n\le 4\times 10^5,1\le q\le 2\times 10^5,1\le a_i,x\le 300$。时限 5.5s，空限 256MB。

---

明显线段树。

有一个想法是维护区间积的欧拉函数，但是这样时间复杂度和代码复杂度都很高……

我的做法是维护区间积。而欧拉函数，就是看看区间中包含什么质因子，然后除一下乘一下好了。

区间积就不用说了。

包含什么质因子？难道要开bool数组吗？时间复杂度很高……

经过后台黑科技操作发现 $300$ 以内的质数只有 $62$ 个。明摆着状压的节奏！

好的，这题做完了。细节的东西在代码中都有。

对于我的代码实现来说：（以下令 $k=62$）

建树 $O(kn)$。

合并节点 $O(1)$。

下推标记 $O(\log n)$。

区间乘 $O(\log^2 n+k)$。

查询欧拉函数 $O(\log n+k)$。

总时间复杂度应该是 $O((n+q)k+q\log^2n)$。其实跑得不慢，我跑得最慢的点是1934ms。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=,mod=;
#define lson o<<1,l,mid
#define rson o<<1|1,mid+1,r
#define FOR(i,a,b) for(int i=(a);i<=(b);i++)
#define ROF(i,a,b) for(int i=(a);i>=(b);i--)
#define MEM(x,v) memset(x,v,sizeof(x))
inline int read(){
    char ch=getchar();int x=,f=;
    while(ch<'' || ch>'') f|=ch=='-',ch=getchar();
    while(ch>='' && ch<='') x=x*+ch-'',ch=getchar();
    return f?-x:x;
}
int n,q,pri[],pl,a[maxn],inv[],f[],tag1[maxn*];    //tag1表示区间要乘多少
ll tag2[maxn*];    //tag2表示区间会多出哪些质因子（也是压缩过的）
//为什么要两个标记呢？不能直接对tag1分解质因子吗？
//因为tag1乘几遍就会被取模，这样看起来质因子就变了。所以额外加一个tag2表示真的质因子集合。
bool vis[];
void init(){
    FOR(i,,){
        if(!vis[i]) pri[++pl]=i;
        for(int j=;j<=pl && i*pri[j]<=;j++){
            vis[i*pri[j]]=true;
            if(i%pri[j]==) break;
        }
    }
    inv[]=;
    FOR(i,,) inv[i]=mod-1ll*(mod/i)*inv[mod%i]%mod;
    FOR(i,,pl) f[i]=1ll*inv[pri[i]]*(pri[i]-)%mod;
    //f[i]表示除以p[i]，再乘上p[i]-1，便于计算欧拉函数
}
inline int qpow(int a,int b){
    int ans=;
    for(;b;b>>=,a=1ll*a*a%mod) if(b&) ans=1ll*ans*a%mod;
    return ans;
}
struct node{
    int pro;ll has;
}nd[maxn*];    //一个线段树节点，pro是区间积，has是区间包含哪些质因子（压缩过的）
void pushup(node &o,node l,node r){    //合并
    o.has=l.has|r.has;    //直接取或
    o.pro=1ll*l.pro*r.pro%mod;
}
void setmult(int o,int l,int r,int x,ll y){    //对第o个节点（管辖[l,r]）区间乘x，质因子多了y
    tag1[o]=1ll*tag1[o]*x%mod;
    tag2[o]|=y;
    nd[o].pro=1ll*nd[o].pro*qpow(x,r-l+)%mod;    //记得乘r-l+1次方
    nd[o].has|=y;
}
void pushdown(int o,int l,int r){    //下传标记
    if(!tag2[o]) return;
    int mid=(l+r)>>;
    setmult(lson,tag1[o],tag2[o]);
    setmult(rson,tag1[o],tag2[o]);
    tag1[o]=;tag2[o]=;    //记得tag1闲置时是1
}
void build(int o,int l,int r){
    tag1[o]=;tag2[o]=;
    if(l==r){
        nd[o].pro=a[l];
        FOR(i,,pl)    //记录质因子集合
            if(a[l]%pri[i]==) nd[o].has|=1ll<<(i-);
        return;
    }
    int mid=(l+r)>>;
    build(lson);build(rson);
    pushup(nd[o],nd[o<<],nd[o<<|]);
}
void mult(int o,int l,int r,int ql,int qr,int x,ll y){    //外面调用时先把质因子集合弄好，会省时间
    if(l>=ql && r<=qr){
        setmult(o,l,r,x,y);    //直接设上
        return;
    }
    pushdown(o,l,r);
    int mid=(l+r)>>;
    if(mid>=ql) mult(lson,ql,qr,x,y);
    if(mid<qr) mult(rson,ql,qr,x,y);
    pushup(nd[o],nd[o<<],nd[o<<|]);
}
node query(int o,int l,int r,int ql,int qr){
    if(l>=ql && r<=qr) return nd[o];
    pushdown(o,l,r);
    int mid=(l+r)>>;
    if(mid<ql) return query(rson,ql,qr);
    if(mid>=qr) return query(lson,ql,qr);
    node ans;
    pushup(ans,query(lson,ql,qr),query(rson,ql,qr));    //合并两边
    return ans;
}
int main(){
    init();
    n=read();q=read();
    FOR(i,,n) a[i]=read();
    build(,,n);
    FOR(i,,q){
        char op[];
        scanf("%s",op);
        int l=read(),r=read();
        if(op[]=='M'){    //乘操作
            int x=read();ll y=;
            FOR(i,,pl) if(x%pri[i]==) y|=1ll<<(i-);    //先处理质因子集合
            mult(,,n,l,r,x,y);
        }
        else{    //求欧拉函数操作
            node ans=query(,,n,l,r);
            int s=ans.pro;    //区间积
            FOR(i,,pl) if(ans.has&(1ll<<(i-))) s=1ll*s*f[i]%mod;    //区间含有第i个质数，那就要除以p[i]，再乘上p[i]-1
            printf("%d\n",s);
        }
    }
}