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;
}