#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

int N, Q;
vector<ll> sum_, mx, mn, lz;
vector<char> hasLz;

void pull(int nd){ sum_[nd]=sum_[2*nd]+sum_[2*nd+1]; mx[nd]=max(mx[2*nd],mx[2*nd+1]); mn[nd]=min(mn[2*nd],mn[2*nd+1]); }

void applyA(int nd,int l,int r,ll v){ sum_[nd]=v*(ll)(r-l+1); mx[nd]=mn[nd]=v; lz[nd]=v; hasLz[nd]=1; }

void push(int nd,int l,int r){ if(hasLz[nd]){ int m=(l+r)/2; applyA(2*nd,l,m,lz[nd]); applyA(2*nd+1,m+1,r,lz[nd]); hasLz[nd]=0; } }

void build(int nd,int l,int r,vector<ll>&a){ if(l==r){sum_[nd]=mx[nd]=mn[nd]=a[l]; return;} int m=(l+r)/2; build(2*nd,l,m,a); build(2*nd+1,m+1,r,a); pull(nd); }

void updA(int nd,int l,int r,int ql,int qr,ll v){ if(qr<l||r<ql) return; if(ql<=l&&r<=qr){applyA(nd,l,r,v); return;} push(nd,l,r); int m=(l+r)/2; updA(2*nd,l,m,ql,qr,v); updA(2*nd+1,m+1,r,ql,qr,v); pull(nd); }

ll isq(ll x){ ll v=(ll)sqrtl((long double)x); while((v+1)*(v+1)<=x) v++; while(v>0&&v*v>x) v--; return v; }

void updS(int nd,int l,int r,int ql,int qr){ if(qr<l||r<ql||mx[nd]<=1) return; if(ql<=l&&r<=qr&&mn[nd]==mx[nd]){ applyA(nd,l,r,isq(mx[nd])); return; } push(nd,l,r); int m=(l+r)/2; updS(2*nd,l,m,ql,qr); updS(2*nd+1,m+1,r,ql,qr); pull(nd); }

ll qry(int nd,int l,int r,int ql,int qr){ if(qr<l||r<ql) return 0; if(ql<=l&&r<=qr) return sum_[nd]; push(nd,l,r); int m=(l+r)/2; return qry(2*nd,l,m,ql,qr)+qry(2*nd+1,m+1,r,ql,qr); }

int main(){
  ios_base::sync_with_stdio(false); cin.tie(0);
  cin>>N>>Q;
  vector<ll> A(N); for(auto&x:A) cin>>x;
  sum_.assign(4*N,0); mx.assign(4*N,0); mn.assign(4*N,0); lz.assign(4*N,0); hasLz.assign(4*N,0);
  build(1,0,N-1,A);
  while(Q--){
    int t; cin>>t;
    if(t==0){ int l,r; cin>>l>>r; if(l<r) cout<<qry(1,0,N-1,l,r-1)<<"\n"; else cout<<0<<"\n"; }
    else if(t==1){ int l,r; ll x; cin>>l>>r>>x; if(l<r) updA(1,0,N-1,l,r-1,x); }
    else { int l,r; cin>>l>>r; if(l<r) updS(1,0,N-1,l,r-1); }
  }
}