Pagini recente » Cod sursa (job #2817568) | Cod sursa (job #1666005) | Cod sursa (job #2400053) | Cod sursa (job #3219150) | Cod sursa (job #3144953)
#include <bits/stdc++.h>
using namespace std;
ifstream fin("cowfood.in");
ofstream fout("cowfood.out");
void maxSelf(int &x, int y) {
if(y > x) {
x = y;
}
}
const int MOD = 3210121;
void addSelf(int &x, int y) {
x += y;
if(x >= MOD) {
x -= MOD;
}
if(x < 0) {
x += MOD;
}
}
int add(int x, int y) {
addSelf(x, y);
return x;
}
void multSelf(int &x, int y) {
x = (long long) x * y % MOD;
}
int mult(int x, int y) {
multSelf(x, y);
return x;
}
int fact(int n) {
static vector<int> data(1, 1);
for(int i = data.size(); i <= n; i++) {
data.emplace_back(mult(data[i - 1], i));
}
return data[n];
}
int inv(int n) {
static vector<int> data(2, 1);
for(int i = data.size(); i <= n; i++) {
data.emplace_back(mult(data[MOD % i], MOD - MOD / i));
}
return data[n];
}
int ifact(int n) {
static vector<int> data(2, 1);
for(int i = data.size(); i <= n; i++) {
data.emplace_back(mult(data[i - 1], inv(i)));
}
return data[n];
}
int comb(int n, int k) {
return mult(fact(n), mult(ifact(k), ifact(n - k)));
}
int main() {
int k, s, n;
fin >> k >> s >> n;
vector<vector<int>> v(n, vector<int>(k));
for(auto &row: v) {
for(auto &it: row) {
fin >> it;
}
}
int cnt = 0;
for(int mask = 0; mask < (1 << n); mask++) {
vector<int> r(k, 1);
for(int i = 0; i < n; i++) {
if((mask >> i) & 1) {
for(int j = 0; j < k; j++) {
maxSelf(r[j], v[i][j]);
}
}
}
int desired = s - accumulate(r.begin(), r.end(), 0);
if(desired < 0) {
continue;
}
if(__builtin_popcount(mask) & 1) {
addSelf(cnt, -comb(desired + k, k));
} else {
addSelf(cnt, comb(desired + k, k));
}
}
fout << cnt << '\n';
return 0;
}