Cf1114f xdtree bitmark

维护一个乘积和一个bitmark表示有那些质因子

注意 mp[i]|=(1ll<<j); 1后面加ll；

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>

using namespace std;

int n,q;

int a[400005];

int isp[305];
int primes[85];
long long afac[85];
long long mp[305];
int cnt;

long long mod=1e9+7;
long long powd(long long a,long long b)
{
    long long ans=1;
    while(b>0)
    {
        if(b%2==1)
            ans=ans*a%mod;
        a=a*a%mod;
        b/=2;
    }
    return ans;
}

int init()
{
    cnt=0;
    for(int i=2;i<=300;i++)
    {
        if(isp[i]==0)
            primes[cnt++]=i;
        for(int j=0;j<cnt;j++)
        {
            if(i*primes[j]>300)
                break;
            isp[i*primes[j]]=1;
            if(i%primes[j]==0)
                break;
        }
    }

    for(int i=0;i<cnt;i++)
    {
        afac[i]=powd(primes[i],mod-2);
    }
    for(int i=1;i<=300;i++)
    {
        for(int j=0;j<cnt;j++)
        {
            if(i%primes[j]==0)
                mp[i]|=(1ll<<j);
        }
    }
}

struct node
{
    int l,r;
    long long sum,mk;
    int f;
    long long fsum,fmk;
}tr[400005 * 5];

node up(node rt,node ls,node rs)
{
    node ntr=rt;
    ntr.sum=(ls.sum*rs.sum)%mod;
    ntr.mk=(ls.mk|rs.mk);

    return ntr;
}

void build(int i,int l,int r)
{
    tr[i].l=l;tr[i].r=r;
    tr[i].f=0;tr[i].fsum=1;tr[i].fmk=0;
    if(l==r)
    {
        tr[i].sum=a[l];tr[i].mk=mp[a[l]];
    }
    else
    {
        int mid=(l+r)/2;
        build(i*2,l,mid);build(i*2+1,mid+1,r);
        tr[i]=up(tr[i],tr[i*2],tr[i*2+1]);
    }
    return;
}

void down(int i)
{
    tr[i].f=0;tr[i*2].f=tr[i*2+1].f=1;

    tr[i*2].fsum=tr[i*2].fsum*tr[i].fsum%mod;tr[i*2].fmk|=tr[i].fmk;
    tr[i*2].sum=tr[i*2].sum*(powd(tr[i].fsum,(tr[i*2].r-tr[i*2].l+1)))%mod;tr[i*2].mk|=tr[i].fmk;

    tr[i*2+1].fsum=tr[i*2+1].fsum*tr[i].fsum%mod;tr[i*2+1].fmk|=tr[i].fmk;
    tr[i*2+1].sum=tr[i*2+1].sum*(powd(tr[i].fsum,(tr[i*2+1].r-tr[i*2+1].l+1)))%mod;tr[i*2+1].mk|=tr[i].fmk;

    tr[i].fmk=0;tr[i].fsum=1;
}

void update(int i,int ql,int qr,int x)
{
    int l=tr[i].l,r=tr[i].r;
    if(l==ql && r==qr)
    {
        tr[i].f=1;tr[i].fsum=tr[i].fsum*x%mod;tr[i].fmk|=mp[x];
        tr[i].sum=tr[i].sum*powd(x,r-l+1)%mod;tr[i].mk|=mp[x];
    } else
    {
        if(tr[i].f==1)
            down(i);

        int mid=(l+r)/2;
        if(qr<=mid)
            update(i*2,ql,qr,x);
        else if(ql>mid)
            update(i*2+1,ql,qr,x);
        else
        {
            update(i*2,ql,mid,x);
            update(i*2+1,mid+1,qr,x);
        }
        tr[i]=up(tr[i],tr[i*2],tr[i*2+1]);
    }
}

node query(int i,int ql,int qr)
{
    int l=tr[i].l,r=tr[i].r;
    if(l==ql && r==qr)
    {
        return tr[i];
    } else
    {
        if(tr[i].f==1)
            down(i);

        int mid=(l+r)/2;
        if(qr<=mid)
            return query(i*2,ql,qr);
        else if(ql>mid)
            return query(i*2+1,ql,qr);
        else
        {
            node ans;
            node a1=query(i*2,ql,mid);
            node a2=query(i*2+1,mid+1,qr);
            return up(ans,a1,a2);
        }
    }
}
int main() {
    //init();
    while(~scanf("%d%d",&n,&q))
    {
        init();

        for(int i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
        }

        build(1,1,n);
        while(q--)
        {
            char ch[15];
            scanf("%s",ch);
            if(ch[0]=='M')
            {
                int l,r,x;
                scanf("%d%d%d",&l,&r,&x);
                update(1,l,r,x);
            } else{
                int l,r;
                scanf("%d%d",&l,&r);
                node ans=query(1,l,r);
                long long anss=ans.sum;
                long long mk=ans.mk;
                for(int i=0;i<cnt;i++)
                {
                    if(((mk>>i)&1)!=0)
                    {
                       anss=anss*(1-afac[i])%mod;
                    }
                }
                anss=(anss+mod)%mod;
                printf("%lld\n",anss);
            }
        }
    }
    return 0;
}