Cod sursa(job #1014086)

Utilizator razvan9310FMI - Razvan Damachi razvan9310 Data 22 octombrie 2013 01:18:40
Problema Elimin Scor 90
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.71 kb
#include <fstream>
#include <vector>
#include <algorithm>
#include <bitset>
using namespace std;

int MAX;
vector<vector<int> > mat;
vector<int> sum_row, sum_col;
ifstream in("elimin.in");
ofstream out("elimin.out");

void input(int &M, int &N, int &R, int &C, int &total) {
  in>>M>>N>>R>>C;
  sum_row.resize(M);
  sum_col.resize(N);
  vector<int> row(N, 0);
  for (int i = 0; i < M; ++i) {
    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 &M, int &N, int &C, int total, bitset<7295> &sol) {
  int i, j;
  for (i = 0; i < N; ++i) {
    sum_col[i] = 0;
  }
  for (i = 0; i < M; ++i) {
    if (!sol[i]) {
      for (j = 0; j < N; ++j) {
        sum_col[j] += mat[i][j];
      }
    }
  }

  sort(sum_col.begin(), sum_col.end());
  for (i = 0; i < C && total > MAX; ++i) {
    total -= sum_col[i];
  }
  if (total > MAX) {
    MAX = total;
  }
}

void rowSort(int &M, int &N, int &R, int total, bitset<7295> &sol) {
  int i, j;
  for (i = 0; i < M; ++i) {
    int sum = 0;
    for (j = 0; j < N; ++j) {
      if (!sol[j]) {
        sum += mat[i][j];
      }
    }
    sum_row[i] = sum;
  }
  sort(sum_row.begin(), sum_row.end());
  for (i = 0; i < R && total > MAX; ++i) {
    total -= sum_row[i];
  }
  if (total > MAX) {
    MAX = total;
  }
}

void linii(int &M, int &N, int &R, int &C, int &total, bitset<7295> &sol) {
  if (!R) {
    columnSort(M, N, C, total, sol);
    return;
  }
  int size = 1 << M;
  for (int i = 0; i < size; ++i) {
    int j, local = total, count = 0;
    for (j = 0; j < M && count < R; ++j) {
      if ((i & (1 << j)) && (local - sum_row[j] > MAX)) {
        sol[j] = 1, local -= sum_row[j];
        ++count;
      } else {
        sol[j] = 0;
      }
    }
    if ((count == R) && ((i >> j) == 0)) {
      columnSort(M, N, C, local, sol);
    }
  }
}

void coloane(int &M, int &N, int &R, int &C, int &total, bitset<7295> &sol) {
  if (!C) {
    rowSort(M, N, R, total, sol);
    return;
  }
  int size = 1 << N;
  for (int i = 0; i < size; ++i) {
    int j, local = total, count = 0;
    for (j = 0; j < N && count < C; ++j) {
      if ((i & (1 << j)) && (local - sum_col[j] > MAX)) {
        sol[j] = 1, local -= sum_col[j];
        ++count;
      } else {
        sol[j] = 0;
      }
    }
    if ((count == C) && ((i >> j) == 0)) {
      rowSort(M, N, R, local, sol);
    }
  }
}

int main() {
  int M, N, R, C, total;
  bitset<7295> sol;
  input(M, N, R, C, total);
  if (M <= N) {
    linii(M, N, R, C, total, sol);//backPeLinii(-1, 0, total);
  } else {
    coloane(M, N, R, C, total, sol);//backPeColoane(-1, 0, total);
  }
  out<<MAX;
  return 0;
}