#include<cmath>

#include<cstdio>

#include<cstring>

#include<iostream>

#include<algorithm>

#define mod 10007

using namespace std;

int tag1[],tag2[],lump[],a[];

int n,sz;

void reset(int x)

{

    for(int i=(x-)*sz+;i<=min(x*sz,n);i++)

    {

        a[i]=(a[i]*tag2[x]+tag1[x])%mod;

    }

    tag1[x]=;

    tag2[x]=;

}

void add(int l,int r,int c)

{

    reset(lump[l]);

    for(int i=l;i<=min(lump[l]*sz,r);i++)

    {

        a[i]+=c;

        a[i]%=mod;

    }

    if(lump[l]!=lump[r])

    {

        reset(lump[r]);

        for(int i=(lump[r]-)*sz+;i<=r;i++)

        {

            a[i]+=c;

            a[i]%=mod;

        }

    }

    for(int i=lump[l]+;i<=lump[r]-;i++)

    {

        tag1[i]+=c;

        tag1[i]%=mod;

    }

}

void mul(int l,int r,int c)

{

    reset(lump[l]);

    for(int i=l;i<=min(lump[l]*sz,r);i++)

    {

        a[i]*=c;

        a[i]%=mod;

    }

    if(lump[l]!=lump[r])

    {

        reset(lump[r]);

        for(int i=(lump[r]-)*sz+;i<=r;i++)

        {

            a[i]*=c;

            a[i]%=mod;

        }

    }

    for(int i=lump[l]+;i<=lump[r]-;i++)

    {

        tag1[i]*=c;

        tag1[i]%=mod;

        tag2[i]*=c;

        tag2[i]%=mod;

    }

}

int main()

{

    int opt,l,r,c;

    scanf("%d",&n);

    sz=sqrt(n);

    for(int i=;i<=n;i++)

    {

        scanf("%d",&a[i]);

        lump[i]=(i-)/sz+;

    }

    for(int i=;i<=lump[n];i++)

    {

        tag2[i]=;

    }

    for(int i=;i<=n;i++)

    {

        scanf("%d%d%d%d",&opt,&l,&r,&c);

        if(!opt)

        {

            add(l,r,c);

        }

        else

        {

            if(opt==)

            {

                mul(l,r,c);

            }

            else

            {

                printf("%d\n",(a[r]*tag2[lump[r]]+tag1[lump[r]])%mod);

            }

        }

    }

}