#include using namespace std; typedef long long int ll; typedef pair P; typedef vector VI; typedef vector VVI; const ll MOD = 1000000007; const ll INF = 1e18; #define REP(i, n) for (int i = 0; i < n; i++) #define ALL(v) v.begin(), v.end() ll power(ll x, ll y) { if (y==0) return 1; else if (y==1) return x%MOD; else if (y%2==0) { ll pow=power(x,y/2); return (pow*pow)%MOD; } else { ll pow=power(x,y/2); return ((pow*pow)%MOD)*x%MOD; } } ll divid(ll x, ll y) { return ((x%MOD)*power(y,MOD-2))%MOD; } int main(){ ll n, k; cin >> n >> k; VI a(n); REP(i,n) cin >> a[i]; VI c(n+1,0); ll p=1; REP(i,n){ c[i]=p; p=(p*(k+i))%MOD; p=divid(p,i+1); } VI lc(n+1,0); REP(i,n) lc[i+1]=(lc[i]+c[i])%MOD; ll ans=0; REP(i,n){ ans+=lc[i+1]*lc[n-i]%MOD*a[i]%MOD; ans%=MOD; } cout << ans << endl; return 0; }