hdu_5634

直接用线段树就可以做，

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <set>
#include <queue>
#include <map>
#include <cmath>
#include <vector>
using namespace std;
#define N 500010
#define pi acos(-1.0)
#define inf 0x3f3f3f3f
#define pb(x) push_back((x))
typedef __int64 lll;
typedef unsigned long long ull;
int phi[10000010];
int a[1200010];
lll sum[1200010];
void phi_table(int n){//欧拉打表
    phi[1]=1;
    for(int i=2;i<=n;i++){
        if(!phi[i]) 
          for(int j=i;j<=n;j+=i){
              if(!phi[j]) phi[j]=j;
              phi[j]=phi[j]/i*(i-1);
          }
    }
}

void build(int root,int l,int r){//建造线段树
    if(l==r){
        scanf("%d",&a[root]);
        sum[root]=a[root];
        return ; 
    }
    int mid=(l+r)>>1;
    build(root<<1,l,mid);
    build(root<<1|1,mid+1,r);
    sum[root]=sum[root<<1]+sum[root<<1|1];
    if(a[root<<1]==a[root<<1|1]) a[root]=a[root<<1];
}

void pushdown(int x, int l, int mid, int r){//区间拆分向下传递操作
    a[x << 1] = a[x << 1 | 1] = a[x];  
    sum[x << 1] = (lll)a[x] * (mid - l + 1);  
    sum[x << 1 | 1] = (lll)a[x] * (r - mid);  
    a[x] = 0;  
}
void gan1(int root,int l,int r,int ll,int rr){//第一个操作
    if(l>=ll&&r<=rr){
        if(a[root]){
            a[root]=phi[a[root]];
            sum[root]=(lll)a[root]*(r-l+1);
            return;
        }
        int mid=(r+l)>>1;
        gan1(root<<1,l,mid,ll,rr);
        gan1(root<<1|1,mid+1,r,ll,rr);
        sum[root]=sum[root<<1]+sum[root<<1|1];
        if(a[root<<1]==a[root<<1|1]) a[root]=a[root<<1];
        else a[root]=0;
        return;
    }
    int mid=(l+r)>>1;
    if(a[root]) pushdown(root,l,mid,r);//如果这个区间的值已经统一，那么在进行更改的时候，如果要向下操作，就要将这个统一的区间拆分成两个统一的小区间，以免给以后造成麻烦
    if(ll<=mid) gan1(root<<1,l,mid,ll,rr);
    if(rr>mid) gan1(root<<1|1,mid+1,r,ll,rr);
    sum[root]=sum[root<<1]+sum[root<<1|1];
    if(a[root<<1]==a[root<<1|1]) a[root]=a[root<<1];
    else a[root]=0;
}

void gan2(int root,int l,int r,int ll,int rr,int x){//第二个操作

    if(l>=ll&&r<=rr){
        a[root]=x;
        sum[root]=(lll)a[root]*(r-l+1);
        return ;
    }
    int mid=(l+r)>>1;
    if(a[root]) pushdown(root,l,mid,r);//如果这个区间的值已经统一，那么在进行更改的时候，如果要向下操作，就要将这个统一的区间拆分成两个统一的小区间，以免给以后造成麻烦
    if(ll<=mid) gan2(root<<1,l,mid,ll,rr,x);
    if(rr>mid) gan2(root<<1|1,mid+1,r,ll,rr,x);
    sum[root]=sum[root<<1]+sum[root<<1|1];
    if(a[root<<1]==a[root<<1|1]) a[root]=a[root<<1];
    else a[root]=0;
}

lll out(int root,int l,int r,int ll,int rr){
    if(l>=ll&&r<=rr) return sum[root];
    lll ans=0;
    int mid=(l+r)>>1;
    if(a[root]) return sum[root]/(r-l+1)*(min(r,rr)-max(l,ll)+1);//如果这个区间的所有值都一样，那么直接就可以得出结果。
    if(ll<=mid) ans+=out(root<<1,l,mid,ll,rr);
    if(rr>mid) ans+=out(root<<1|1,mid+1,r,ll,rr);
    return ans;
}

int main(){
    int t,l,r,x,n,m;
    phi_table(1e7);//欧拉打表
    scanf("%d",&t);
    while(t--){
        memset(sum,0,sizeof(sum));
        memset(a,0,sizeof(a));
        scanf("%d%d",&n,&m);
        build(1,1,n);
        int tar;
        while(m--){
            scanf("%d",&tar);
            if(tar==1){
                scanf("%d%d",&l,&r);
                gan1(1,1,n,l,r);
            }
            else if(tar==2){
                scanf("%d%d%d",&l,&r,&x);
                gan2(1,1,n,l,r,x);
            }
            else{
                scanf("%d%d",&l,&r);
                printf("%I64d\n",out(1,1,n,l,r));
            }
        }
    }
    return 0;
}