Pagini recente » Cod sursa (job #3241387) | Cod sursa (job #1291836) | Cod sursa (job #328294) | Cod sursa (job #2195938) | Cod sursa (job #1014039)
#include <fstream>
#include <iostream> //
#include <vector>
#include <algorithm>
#include <bitset>
using namespace std;
int M, N, R, C, MAX, total;
vector<vector<int> > mat;
vector<int> sum_row, sum_col;
bitset<8000> sol;
ifstream in("elimin.in");
ofstream out("elimin.out");
void input() {
in>>M>>N>>R>>C;
sum_row.resize(M);
sum_col.resize(N);
for (int i = 0; i < M; ++i) {
vector<int> row(N, 0);
mat.push_back(row);
int val, sum = 0;
for (int j = 0; j < N; ++j) {
in>>val;
mat[i][j] = val;
sum += val;
sum_col[j] += val;
}
total += sum;
sum_row[i] = sum;
}
}
void columnSort(int total) {
int i, j;
for (j = 0; j < N; ++j) {
int sum = 0;
for (i = 0; i < M; ++i) {
if (sol[i] == 0) {
sum += mat[i][j];
}
}
sum_col[j] = sum;
}
sort(sum_col.begin(), sum_col.end());
for (i = 0; i < C; ++i) {
total -= sum_col[i];
}
if (total > MAX) {
MAX = total;
}
}
void rowSort(int total) {
int i, j;
for (i = 0; i < M; ++i) {
int sum = 0;
for (j = 0; j < N; ++j) {
if (sol[j] == 0) {
sum += mat[i][j];
}
}
sum_row[i] = sum;
}
sort(sum_row.begin(), sum_row.end());
for (i = 0; i < R; ++i) {
total -= sum_row[i];
}
if (total > MAX) {
MAX = total;
}
}
void backPeLinii(const int &pos, const int &lvl, const int &total) {
if (lvl == R) {
columnSort(total);
return;
}
for (int i = pos + 1; (i < M) && (M - i >= R - lvl); ++i) {
if (total - sum_row[i] > MAX) {
sol[i] = 1;
backPeLinii(i, lvl + 1, total - sum_row[i]);
sol[i] = 0;
}
}
}
void backPeColoane(const int &pos, const int &lvl, const int &total) {
if (lvl == C) {
rowSort(total);
return;
}
for (int i = pos + 1; (i < N) && (N - i >= C - lvl); ++i) {
if (total - sum_col[i] > MAX) {
sol[i] = 1;
backPeColoane(i, lvl + 1, total - sum_col[i]);
sol[i] = 0;
}
}
}
void output() {
out<<MAX;
}
void solve() {
if (M < N) {
backPeLinii(-1, 0, total);
} else {
backPeColoane(-1, 0, total);
}
}
int main() {
input();
solve();
output();
return 0;
}