10pts求条，符合所有进食后入

#include<bits/stdc++.h>
using namespace std;
#define int long long
#define mod 1000000007
int c1[100005],c2[100005],n;
int q1(int x)
{
    int ans=0;
    for(int i=x;i;i-=(i&-i))
        ans=(ans+c1[i])%mod;
    return ans;
}
int q2(int x)
{
    int ans=0;
    for(int i=x;i;i-=(i&-i))
        ans=(ans+c2[i]);
    return ans;
}
void setp(int x,int k)
{
    int oldval=(q1(x)-q1(x-1)+mod)%mod;
    for(int i=x;i<=n;i+=(i&-i))
    {
        c2[i]=(c2[i]-oldval*oldval%mod+mod)%mod;
        c2[i]=(c2[i]+k*k)%mod;
        c1[i]=(c1[i]-oldval+mod)%mod;
        c1[i]=(c1[i]+k)%mod;
    }
}
constexpr int inv(int x)
{
    int a=x,b=mod-2,ret=1;
    while(b)
    {
        if(b&1)
            ret=(ret*a)%mod;
        a=(a*a)%mod;
        b>>=1;
    }
    return ret;
}
signed main()
{
	int m;
    cin>>n>>m;
    for(int i=1;i<=n;i++)
    {
        int t;
        cin>>t;
        setp(i,t);
    }
    for(int i=1;i<=m;i++)
    {
        int op,x,y;
        cin>>op>>x>>y;
        if(op==1)
            setp(x,y);
        else
        {
            int invn=inv(y-x+1);
            int sum=(q1(y)-q1(x-1)+mod)%mod;
            int sum2=(q2(y)-q2(x-1)+mod)%mod;
            int avg=invn*(sum+mod)%mod;
            cout<<(invn*(sum2-(2*avg*sum)%mod+(y-x+1)*avg*avg%mod+mod)%mod+mod)%mod<<'\n';
        }
    }
}