#include using namespace std; #define REP(i,a,b) for(i=a;i'9')break;*x=(*x)*10+k-'0';}if(m)(*x)=-(*x);} void reader(ll *x){int k,m=0;*x=0;for(;;){mygc(k);if(k=='-'){m=1;break;}if('0'<=k&&k<='9'){*x=k-'0';break;}}for(;;){mygc(k);if(k<'0'||k>'9')break;*x=(*x)*10+k-'0';}if(m)(*x)=-(*x);} void reader(double *x){scanf("%lf",x);} int reader(char c[]){int i,s=0;for(;;){mygc(i);if(i!=' '&&i!='\n'&&i!='\r'&&i!='\t'&&i!=EOF) break;}c[s++]=i;for(;;){mygc(i);if(i==' '||i=='\n'||i=='\r'||i=='\t'||i==EOF) break;c[s++]=i;}c[s]='\0';return s;} template void reader(T *x, S *y){reader(x);reader(y);} template void reader(T *x, S *y, U *z){reader(x);reader(y);reader(z);} template void reader(T *x, S *y, U *z, V *w){reader(x);reader(y);reader(z);reader(w);} void writer(int x, char c){int s=0,m=0;char f[10];if(x<0)m=1,x=-x;while(x)f[s++]=x%10,x/=10;if(!s)f[s++]=0;if(m)mypc('-');while(s--)mypc(f[s]+'0');mypc(c);} void writer(ll x, char c){int s=0,m=0;char f[20];if(x<0)m=1,x=-x;while(x)f[s++]=x%10,x/=10;if(!s)f[s++]=0;if(m)mypc('-');while(s--)mypc(f[s]+'0');mypc(c);} void writer(double x, char c){printf("%.15f",x);mypc(c);} void writer(const char c[]){int i;for(i=0;c[i]!='\0';i++)mypc(c[i]);} void writer(const char x[], char c){int i;for(i=0;x[i]!='\0';i++)mypc(x[i]);mypc(c);} template void writerLn(T x){writer(x,'\n');} template void writerLn(T x, S y){writer(x,' ');writer(y,'\n');} template void writerLn(T x, S y, U z){writer(x,' ');writer(y,' ');writer(z,'\n');} template void writerArr(T x[], int n){int i;if(!n){mypc('\n');return;}rep(i,n-1)writer(x[i],' ');writer(x[n-1],'\n');} char memarr[17000000]; void *mem = memarr; #define MD 1000000007 int get_inv(ll a, int md){ll t=a,s=md,u=1,v=0,e;while(s){e=t/s;t-=e*s;u-=e*v;swap(t,s);swap(u,v);}if(u<0)u+=md;return u;} struct mint{ static unsigned md, W, R, Rinv, mdninv, RR; unsigned val; mint(){}mint(int a){val=mulR(a);}mint(unsigned a){val=mulR(a);}mint(ll a){val=mulR(a);}mint(ull a){val=mulR(a);} unsigned setmod(unsigned m){int i;unsigned t;W=32;md=m;R=(1ULL<>W);if(t>=md)t-=md;return t;}unsigned reduce(ull T){unsigned m=(unsigned)T*mdninv;unsigned t=(unsigned)((T+(ull)m*md)>>W);if(t>=md)t-=md;return t;} unsigned get(){return reduce(val);} mint&operator+=(mint a){val+=a.val;if(val>=md)val-=md;return*this;}mint&operator-=(mint a){if(val>=1;a*=a;}return r;} bool operator==(int a){return get()==a;} bool operator!=(int a){return get()!=a;} }; unsigned mint::md, mint::W, mint::R, mint::Rinv, mint::mdninv, mint::RR; mint operator+(int a, mint b){return mint(a)+=b;}mint operator-(int a, mint b){return mint(a)-=b;}mint operator*(int a, mint b){return mint(a)*=b;}mint operator/(int a, mint b){return mint(a)/=b;} mint operator+(ll a, mint b){return mint(a)+=b;}mint operator-(ll a, mint b){return mint(a)-=b;}mint operator*(ll a, mint b){return mint(a)*=b;}mint operator/(ll a, mint b){return mint(a)/=b;} mint mval[10000], minv[10000]; void mint_init(int md=MD, mint val[]=mval, int vals=10000, mint inv[]=minv, int invs=10000){int i;val[0].setmod(md);val[0].val=0;REP(i,1,vals){val[i].val=val[i-1].val+mint::R;if(val[i].val >=md)val[i].val-=md;}inv[1].val=1;REP(i,2,invs){inv[i].val=md-((ll)(md/i)*inv[md%i].val%md);}REP(i,1,invs)inv[i].val=(ull)inv[i].val*mint::R%md;} template struct fenwick{ int size, memory; T *data; void malloc(int mem){memory=mem;data=(T*)std::malloc(sizeof(T)*mem);} void *malloc(int mem, void *workMemory){memory=mem;data=(T*)workMemory;return (void*)(data+mem);} void free(void){memory=0;free(data);} void init(int N){size=N;memset(data,0,sizeof(T)*N);} void add(int k, T val){while(k=0)res+=data[k],k=(k&(k+1))-1;return res;} T range(int a, int b){if(b==-1)b=size-1;return get(b)-get(a-1);} int kth(T k){int i=0,j=size,c;T v;while(i=0;k--) if(d[k]d[k]) break; j=d[i], d[i]=d[k], d[k]=j; intReverse(d+k+1, size-k-1); return i; } int N; ll K; int P[100010]; ll fact[100010]; ll sa[100010], sum[100010]; int main(){ int i, j, k, m; fenwick t; mint res = 0; ll tmp; mint_init(); reader(&N,&K); rep(i,N) reader(P+i), P[i]--; if(0){ int i, j, k; int PP[100]; rep(i,N) PP[i] = P[i]; int res = 0; rep(k,K){ rep(i,N) REP(j,i+1,N) if(PP[i] > PP[j]) res++; next_permutation(PP, PP+N); } writerLn(res); } fact[0] = 1; REP(i,1,N+10){ if(fact[i-1] > 3000000000000000000LL / i) fact[i] = 3000000000000000000LL; else fact[i] = fact[i-1] * i; } if(K >= fact[N]){ tmp = K / fact[N]; K -= fact[N] * tmp; res += mint((ull)tmp) * (fact[N-1] * N / 2); } t.malloc(N); t.init(N); rep(i,N) t.add(i,1); rep(i,N){ sa[i] = t.get(P[i]-1); t.add(P[i], -1); } sum[0] = 0; rep(i,N) sum[i+1] = sum[i] + sa[i]; while(K){ for(i=N-1;i;i--){ m = N-i; if(P[i-1] > P[i]) break; if(fact[m+1] > K) break; } K -= fact[m]; res += mint((ull)fact[m]) * mint((ull)sum[N-m]) + mint((ull)(fact[m-1] * m / 2)); reverse(P+N-m, P+N); i = intNextPermutation(P, N); if(i < N-50){ t.init(N); rep(i,N) t.add(i,1); rep(i,N){ sa[i] = t.get(P[i]-1); t.add(P[i], -1); } sum[0] = 0; rep(i,N) sum[i+1] = sum[i] + sa[i]; } else { i = max(0, i-10); REP(j,i,N){ sa[j] = 0; REP(k,j+1,N) if(P[j] > P[k]) sa[j]++; sum[j+1] = sum[j] + sa[j]; } } } writerLn((int)res); return 0; }