#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; constexpr int INF = 1001001001; // constexpr int mod = 1000000007; constexpr int mod = 998244353; template inline bool chmax(T& x, T y){ if(x < y){ x = y; return true; } return false; } template inline bool chmin(T& x, T y){ if(x > y){ x = y; return true; } return false; } struct ExtendedSieve{ vector table; const int n; ExtendedSieve(int n) : n(n) { table.assign(n + 1, 0); table[0] = table[1] = 1; for(int i = 2; i <= n; ++i){ if(table[i] > 0) continue; table[i] = i; for(int j = i + i; j <= n; j += i){ if(table[j] == 0) table[j] = i; } } } int operator[](const int& k){ return table[k]; } bool is_prime(int x){ if(x < 2) return false; return x == table[x]; } map prime_factorization(int x){ map res; while(x > 1){ ++res[table[x]]; x /= table[x]; } return res; } }; int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); int N, K; cin >> N >> K; vector A(N); vector> primes(N); ExtendedSieve sieve(100000); for(int i = 0; i < N; ++i){ cin >> A[i]; primes[i] = sieve.prime_factorization(A[i]); } set> st; for(int mask = 0; mask < 1 << N; ++mask){ if(__builtin_popcount(mask) < K) continue; int sum = 0; map prod; for(int i = 0; i < N; ++i){ if(mask >> i & 1){ sum += A[i]; for(auto& [p, c] : primes[i]){ prod[p] += c; } } } st.emplace(sieve.prime_factorization(sum)); st.emplace(prod); } cout << st.size() << '\n'; return 0; }