#include #define VARNAME(x) #x #define show(x) cerr << #x << " = " << x << endl using namespace std; using ll = long long; using ld = long double; template vector Vec(int n, T v) { return vector(n, v); } template auto Vec(int n, Args... args) { auto val = Vec(args...); return vector(n, move(val)); } template ostream& operator<<(ostream& os, const vector& v) { os << "sz:" << v.size() << "\n["; for (const auto& p : v) { os << p << ","; } os << "]\n"; return os; } template ostream& operator<<(ostream& os, const pair& p) { os << "(" << p.first << "," << p.second << ")"; return os; } constexpr ll MOD = (ll)1e9 + 9LL; constexpr ld PI = static_cast(3.1415926535898); template constexpr T INF = numeric_limits::max() / 10; ll power(const ll p, const ll n) { if (n == 0) { return 1; } if (n % 2 == 1) { return (power(p, n - 1) * p) % MOD; } else { const ll pp = power(p, n / 2); return (pp * pp) % MOD; } } ll power2(const ll n) { return power(2, n); } int main() { cin.tie(0); ios::sync_with_stdio(false); int N, K; cin >> N >> K; vector> A(N); for (int i = 0; i < N; i++) { ll a; cin >> a; A[i] = bitset<60>(a); } ll mean = 0; ll var = 0; for (int i = 0; i < K; i++) { int one = 0; for (int j = 0; j < N; j++) { if (A[j][i]) { one++; } } const ll fac = power2(K - 1 + i); const ll m_2 = N; const ll v_4 = (N - 2 * one) * (N - 2 * one); mean += (m_2 * fac) % MOD; var += (v_4 * ((fac * fac) % MOD)) % MOD; var %= MOD; mean %= MOD; } cout << mean << endl; cout << var << endl; return 0; }