#include <bits/stdc++.h>
#define int long long
#define pii pair<int, int>
using namespace std;
const int mod = 1e9 + 7;
int f_pow(int x, int y) {
    if (y == 1) return x;
    if (y == 0) return 1;
    int t = f_pow(x, y / 2) % mod;
    if (y % 2)
        return t * t % mod * x % mod;
    else
        return t * t % mod;
}
struct node{
    int val, mi, left;
};
bool cmp1(const node x,const node y){
    return x.val < y.val;
}
bool cmp2(const node x,const node y){
    return x.left < y.left;
}
signed main() {
    // ios::sync_with_stdio(false);
    // cin.tie(0);
    // cout.tie(0);
    int n, k, tot = 0;
    cin >> n >> k;
    vector<node> a(n);
    for (int i = 0; i < n; i++) {
        cin >> a[i].val;
        a[i].mi = log2(a[i].val);
        tot += a[i].mi;
    }
    vector<int> mult(n);  // 储存这里即将乘2的多少次方
    int ava = (tot + k) / n;
    sort(a.begin(), a.end(),cmp1);
    int m;
    for (m = n - 1; m > 0; m--) {
        ava = (tot + k) / (m + 1);
        tot -= a[m].mi;
        if (a[m].mi < ava) break;
    }
    for (int i = 0; i <= m; i++) {
        mult[i] = ava - a[i].mi;
        k -= mult[i];
    }
    sort(a.begin(), a.begin()+m+1, cmp2);
    for (int i = 0; i <= m; i++) {
        if (k == 0) break;
        mult[i]++;
        k--;
    }
    // puts("mult:");for (int i = 0; i <= m; i++) cout << mult[i] << ' ';
    int ans = 0;
    // printf("\nnums:\n");
    for (int i = 0; i < n; i++) {
        a[i].val = a[i].val * f_pow(2, mult[i]) % mod;
        // cout << a[i].val << ' ';
        ans = (ans + a[i].val) % mod;
    }
    cout << ans % mod;
    return 0;
}