#include<iostream>

 #include<stdio.h>

 #include<algorithm>

 #include<iomanip>

 #include<cmath>

 #include<cstring>

 #include<vector>

 #define ll __int64

 #define pi acos(-1.0)

 #define MAX 100003

 #define M 10007

 #define lson i<<1

 #define rson i<<1|1

 using namespace std;

 struct tree

 {

     int l,r,sum[],add,mul;

     void init(int _l,int _r){

         l=_l;r=_r;mul=;add=;

     }

     void Add(int a){

         sum[]=(sum[]+a*a%M*a%M*len()%M+*a*sum[]%M+*a*a%M*sum[]%M)%M;

         sum[]=(sum[]+a*a%M*len()%M+*a*sum[]%M)%M;

         sum[]=(sum[]+a*len()%M)%M;

         add=(add+a)%M;

     }

     void Mul(int a){

         sum[]=sum[]*a%M*a%M*a%M;

         sum[]=sum[]*a%M*a%M;

         sum[]=sum[]*a%M;

         add=add*a%M;

         mul=mul*a%M;

     }

     int len(){

         return r-l+;

     }

     int mid(){

         return (l+r)>>;

     }

 }T[MAX*];

 void down(int i)

 {

     T[lson].Mul(T[i].mul);

     T[rson].Mul(T[i].mul);

     T[lson].Add(T[i].add);

     T[rson].Add(T[i].add);

     T[i].add=;

     T[i].mul=;

 }

 void updateson(int i)

 {

     for(int j=;j<;j++)

         T[i].sum[j]=(T[lson].sum[j]+T[rson].sum[j])%M;

 }

 void built(int i,int l,int r)

 {

     T[i].init(l,r);

     memset(T[i].sum,,sizeof(T[i].sum));

     if(l==r) return ;

     int m=T[i].mid();

     built(lson,l,m);

     built(rson,m+,r);

 }

 void update(int i,int l,int r,int mul,int add)

 {

     if(T[i].l==l&&T[i].r==r){

         T[i].Mul(mul);

         T[i].Add(add);

         return ;

     }

     down(i);

     int m=T[i].mid();

     if(r<=m) update(lson,l,r,mul,add);

     else if(l>m) update(rson,l,r,mul,add);

     else{

         update(lson,l,m,mul,add);

         update(rson,m+,r,mul,add);

     }

     updateson(i);

 }

 int query(int i,int l,int r,int p)

 {

     if(T[i].l==l&&T[i].r==r) return T[i].sum[p-];

     down(i);

     int m=T[i].mid();

     if(r<=m) return query(lson,l,r,p);

     else if(l>m) return query(rson,l,r,p);

     else return (query(lson,l,m,p)+query(rson,m+,r,p))%M;

 }

 int main(){

     int m,n,op,x,y,p;

     while(scanf("%d%d",&n,&m)&&(m+n)){

         built(,,n);

         while(m--){

             scanf("%d%d%d%d",&op,&x,&y,&p);

             if(op==) update(,x,y,,p);

             else if(op==) update(,x,y,p,);

             else if(op==) update(,x,y,,p);

             else printf("%d\n",query(,x,y,p)%M);

         }

     }

     return ;

 }