#include<bits/stdc++.h>
using namespace std;

#define int long long 
#define F first 
#define S second 
#define pb push_back

const int mol=1e9+7;

vector<int>node[5005];
bool vis[5005];
int dis[5005];
bool ok=1;
int t1=0,t2=0;

void dfs(int v){
    vis[v]=1;
    if(dis[v]==1) t2++;
    else t1++;
    for(int k:node[v]){
        if(!vis[k]) dis[k]=(dis[v]^1),dfs(k);
        else{
            if(dis[k]!=(dis[v]^1)){
                ok=0;
                return;
            }
        }
    }
}

int dp[2][5005];

void solve(){
    int n,m;
    cin>>n>>m;

    vector<int>a(n+1);
    for(int i=1;i<=n;i++) cin>>a[i];

    for(int i=1;i<=n;i++){
        for(int j=i+1;j<=n;j++){
            if(a[i]>a[j]){
                node[i].pb(j);
                node[j].pb(i);
            }
        }
    }

    dp[0][0]=1;
    int iu=1;
    for(int i=1;i<=n;i++){
        if(!vis[i]){
            t1=0,t2=0;
            dis[i]=1;
            dfs(i);
        }
        else continue;
        if(!ok){
            cout<<"0\n";
            return;
        }
        if(t1==0){
            for(int j=m;j>=1;j--) dp[iu^1][j]=dp[iu^1][j-1];
            dp[iu^1][0]=0;
            continue;
        }
        for(int j=m;j>=0;j--) dp[iu][j]=0;
        for(int j=m;j>=0;j--){
            if(j-t1>=0)
            dp[iu][j]=dp[iu^1][j-t1];
            if(j-t2>=0)
            dp[iu][j]+=dp[iu^1][j-t2];
            dp[iu][j]%=mol;
        }
        iu^=1;
    }

    iu^=1;
    int ans=0;
    for(int j=m;j>=0;j--){
        if(n-j>m) break;
        ans+=dp[iu][j];
        ans%=mol;
    }
    ans<<=1;
    ans%=mol;
    cout<<ans<<"\n";
}

signed main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    solve();
    return 0;
}