#include<iostream>
#include<cstring>
#include<vector>
using namespace std;
using LL = long long;

const int mod = 1e9 + 7;

void add(int &a, int b){
    a += b;
    if (a >= mod) a -= mod;
}

int mul(int a, int b){
    return 1LL * a * b % mod;
}

int qpow(int a, int b, int mod){
    int res = 1;
    while (b){
        if (b & 1) res = 1LL * res * a % mod;
        a = 1LL * a * a % mod;
        b >>= 1;
    }
    return res;
}

const int N = 44, M = 25, T = N - M - 1;
int dp[N + 1][1 << T], ndp[N + 1][1 << T];

int main(){

#ifdef LOCAL
    freopen("data.in", "r", stdin);
    freopen("data.out", "w", stdout);
#endif

    cin.tie(0);
    cout.tie(0);
    ios::sync_with_stdio(0);

    LL a[50]{};

    auto insert = [&](LL x){
        for(int i = N - 1; i >= 0; i--){
            if (x >> i & 1){
                if (a[i]) x ^= a[i];
                else{
                    for(int j = N - 1; j > i; j--){
                        if (a[j] >> i & 1){
                            a[j] ^= x;
                        }
                    }
                    a[i] = x;
                    return;
                }
            }
        }
    };

    int n, k;
    cin >> n >> k;
    for(int i = 0; i < n; i++){
        LL x;
        cin >> x;
        insert(x);
    }
    vector<LL> p;
    for(int i = N - 1; i >= 0; i--){
        if (a[i] > 0){
            p.push_back(a[i]);
        }
    }

    int pows[50]{};
    for(int i = 1; i <= N; i++) pows[i] = qpow(i, k, mod);

    int ans = 0;
    if (p.size() <= M){
        const int m = p.size();
        vector<LL> val(1 << m);
        for(int j = 0; j < 1 << m; j++){
            if (j != 0){
                int t = __lg(j);
                val[j] = val[j ^ (1 << t)] ^ p[t];
            }
            add(ans, pows[__builtin_popcountll(val[j])]);
        }
        ans = mul(ans, qpow(2, n - p.size(), mod));
    }
    else{
        LL state = 0;
        for(auto x : p) state |= 1LL << __lg(x);
        vector<int> bit;
        for(int i = 0; i < N; i++){
            if (!(state >> i & 1)){
                bit.push_back(i);
            }
        }
        const int m = bit.size();
        dp[0][0] = 1;
        for(auto x : p){
            for(int i = 0; i <= p.size(); i++){
                memset(ndp[i], 0, sizeof ndp[i]);
            }
            int val = 0;
            for(int i = 0; i < bit.size(); i++){
                if (x >> bit[i] & 1){
                    val |= 1 << i;
                }
            }
            for(int i = 0; i <= p.size(); i++){
                for(int j = 0; j < 1 << m; j++){
                    add(ndp[i][j], dp[i][j]);
                    add(ndp[i + 1][j ^ val], dp[i][j]);
                }
            }
            swap(dp, ndp);
        }
        for(int i = 0; i <= p.size(); i++){
            for(int j = 0; j < 1 << m; j++){
                add(ans, mul(dp[i][j], pows[i + __builtin_popcount(j)]));
            }
        }
        ans = mul(ans, qpow(2, n - p.size(), mod));
    }
    cout << ans << '\n';

}