Pagini recente » Cod sursa (job #1494227) | Cod sursa (job #1800374) | Cod sursa (job #748758) | Cod sursa (job #3223290) | Cod sursa (job #2809342)
#include <fstream>
#include <iostream>
#include <vector>
using namespace std;
ifstream fin( "cowfood.in" );
ofstream fout( "cowfood.out" );
const int MAXS = 10005;
const int MAXK = 35;
const int MAXN = 25;
const int MOD = 3210121;
int cnt[MAXN][MAXK];
int comb[MAXS + MAXK][MAXK];
void precalc() {
for ( int i = 1; i < MAXS + MAXK; ++i ) {
comb[i][1] = i;
for ( int j = 2; j < MAXK; ++j ) {
comb[i][j] = (comb[i - 1][j] + comb[i - 1][j - 1]) % MOD;
}
}
}
int sum[MAXN];
//tle
int sb( int n, int k ) {
return comb[n + k - 1][k - 1];
}
int main() {
int k, s, n, bad = 0;
fin >> k >> s >> n;
for ( int i = 0; i < n; ++i ) {
for ( int j = 0; j < k; ++j ) {
fin >> cnt[i][j];
}
}
precalc();
sum[0] = 1;
for ( int i = 1; i <= s; ++i ) {
sum[i] = (sum[i-1] + sb(i, k)) % MOD;
}
for ( int mask = 1; mask < (1 << n); ++mask ) {
vector<int> a(k);
int h = 0;
for ( int i = 0; i < n; ++i ) {
if ( (mask >> i) & 1 ) {
for ( int j = 0; j < k; ++j ) {
a[j] = max(a[j], cnt[i][j]);
}
++h;
}
}
int sp = 0;
for ( auto val : a ) {
sp += val;
}
if ( s >= sp ) {
if ( h & 1 ) {
bad = (bad + sum[s - sp]) % MOD;
} else {
bad = (bad - sum[s - sp]) % MOD;
}
}
}
bad = (bad + 1 + s * k) % MOD;
fout << (sum[s] - bad + MOD) % MOD;
fin.close();
fout.close();
return 0;
}