#include<bits/stdc++.h>
using namespace std;
/*
 * Author: Yang
 * Date: 2024-12-03
 */
//#pragma GCC optimize(2)
using ll=long long;
using uint=unsigned int;
using ull=unsigned ll;
using pii=pair<int,int>;
#define all(v) begin(v),end(v)
#define fst first
#define scd second
#define rep(i,s,t) for(int i=(s);i<=(t);++i)
#define irep(i,t,s) for(int i=(t);i>=(s);--i)
template<typename T>
bool chmax(T &x,T y){return x<y?(x=y,1):0;}
template<typename T>
bool chmin(T &x,T y){return y<x?(x=y,1):0;}
const int mod=1e9+7;
int qp(int a,ll b,int p=mod){
    int r=1;
    while(b){
        if(b&1) r=1ll*r*a%p;
        a=1ll*a*a%p;b>>=1;
    }
    return r;
};
int solve_2_to_n(vector<ll> a, int k){
    constexpr int m=45;
    vector<int> pwk(m+1,0);
    rep(i,1,m){
        pwk[i]=qp(i,k);
    }
    array<int,2> n;
    n[0]=a.size()/2;
    n[1]=a.size()-a.size()/2;
    vector<ll> S[2];
    rep(t,0,1){
        S[t].resize(1<<n[t],0);
        rep(x,1,(1<<n[t])-1){
            S[t][x]=S[t][x^(x&-x)]^a[__lg(x&-x)+t*n[0]];
        }
    }
    vector<int> pc(1<<(m/2),0);
    rep(x,1,(1<<(m/2))-1){
        pc[x]=pc[x^(x&-x)]+1;
    }
    int ans=0;
    rep(x,0,(1<<n[0])-1){
        rep(y,0,(1<<n[1])-1){
            ll v=S[0][x]^S[1][y];
            int c=pc[v>>(m/2)]+pc[v&((1ll<<(m/2))-1)];
            (ans+=pwk[c])%=mod;
        }
    }
    return ans;
}
template<int m0=45>
int solve(vector<ll> a,int k){
    int m=m0;
    vector<int> pwk(m+1,0);
    rep(i,1,m0){
        pwk[i]=qp(i,k);
    }
    auto rm_digit=[&](int i){
        m--;
        ll const msk=(1ll<<i)-1;
        for(auto &x:a){
            x=(x&msk)|(x>>1)&~msk;
        }
    };
    int n=a.size();
    vector<int> piv(a.size()),pivs;
    int p=0;
    rep(i,0,n-1){
        while(p<m0 && !any_of(i+all(a),[&](ll x){return x>>p&1;})) p++;
        if(!a[i]>>p&1){
            rep(j,i+1,n-1){
                if(a[j]>>p&1){
                    a[i]^=a[j];
                    break;
                }
            }
        }
        rep(j,i+1,n-1){
            if(a[j]>>p&1) a[j]^=a[i];
        }
        piv[i]=p;
        if(p<m0) pivs.push_back(p);
    }
    irep(i,(int)pivs.size()-1,0) rm_digit(pivs[i]);
    auto trans=[&](valarray<int> const &v, ll x)->valarray<int> {
        valarray<int> r(0,1<<m);
        rep(i,0,(1<<m)-1) r[i^x]=v[i];
        return r;
    };
    valarray<valarray<int>> f(valarray<int>(0,1<<m),n+1);
    f[0][0]=1;
    // cout<<m<<'\n';
    rep(i,0,(int)a.size()-1){
        // cout<<i<<' '<<a[i]<<','<<piv[i]<<'\n';
        irep(c,n-1,0){
            bool dc=piv[i]<m0;
            // cout<<' '<<c<<"->"<<c+dc<<'\n';
            f[c+dc]+=trans(f[c],a[i]);
            f[c+dc]%=mod;
        }
    }
    int ans=0;
    rep(c,0,n){
        rep(x,0,(1<<m)-1){
            // if(f[c][x]) cout<<c<<' '<<x<<" "<<f[c][x]<<' '<<pwk[__builtin_popcount(x)+c]<<'\n';
            (ans+=1ll*f[c][x]*pwk[__builtin_popcount(x)+c]%mod)%=mod;
        }
    }
    return ans;
}
signed main(){
    //ios::sync_with_stdio(false);
    int n,k;
    cin>>n>>k;
    vector<ll> a;
    rep(i,1,n){
        ll x;cin>>x;
        a.push_back(x);
    }
    cout<<(n<=25?solve_2_to_n(a,k):solve(a,k))<<"\n";
}