#include<bits/stdc++.h>
using namespace std;
#define ll long long
const ll mod=1000000007;

int a[4004];
pair<int,int> aa[8000000];
int main(){
    // int n; cin>>n;
    // for(int i=0;i<n;++i)cin>>a[i];
    int n=4000; srand(time(0));
    for(int i=0;i<n;++i)a[i]=rand();
    int ptr=0;
    for(int k=0;k<n;++k){
        set<int> st;
        for(int ii=0;ii<k;++ii)st.insert(a[ii]);
        for(int i=0;i<n-k;++i){
            st.insert(a[i+k]);
            auto it=st.end(); --it;
            aa[ptr]={*it,ptr}; ++ptr;
            st.erase(a[i]);
        }
    }
    ll ans=0;
    sort(aa,aa+ptr,greater<pair<int,int>>());
    set<int> rng; rng.insert({-1,ptr});
    for(int i=0;i<ptr;++i){
        // cout<<aa[i].first<<" "<<aa[i].second<<endl;
        auto it2=rng.lower_bound(aa[i].second);
        auto it1=it2; --it1;
        ll dis=(aa[i].second-*it1)*1ll*(*it2-aa[i].second);
        ans=(ans+aa[i].first*dis%mod)%mod;
        rng.insert(aa[i].second);
    }
    cout<<ans<<endl;
}