#include <fstream>
#include <vector>
#include <bitset>
#include <algorithm>
using namespace std;

vector<vector<int> > mat;
vector<int> sum_row, sum_col;
int M, N, R, C, MAX, total;
bitset<8000> final_sol, partial_sol;

void input() {
  ifstream in("elimin.in");
  in>>M>>N>>R>>C;
  int i, j, val;
  for (i = 0; i < M; ++i) {
    vector<int> row(N, 0);
    mat.push_back(row);
    int sum = 0;
    for (j = 0; j < N; ++j) {
      in>>val;
      mat[i][j] = val;
      sum += val;
    }
    sum_row.push_back(sum);
    total += sum;
  }
  in.close();
}

void column_sort(int total) {
  int i, j;
  for (j = 0; j < N; ++j) {
    int sum = 0;
    for (i = 0; i < M; ++i) {
      if (final_sol[i] == 0) {
        sum += mat[i][j];
      }
    }
    sum_col.push_back(sum);
  }

  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 back(const int &pos, const int &lvl, const int &total) {
  if (lvl == R) {
    final_sol = partial_sol;
    column_sort(total);
    return;
  }
  for (int i = pos + 1; (i < M) && (M - i >= R - lvl); ++i) {
    if (total - sum_row[i] > MAX) {
      partial_sol[i] = 1;
      back(i, lvl + 1, total - sum_row[i]);
      partial_sol[i] = 0;
    }
  }
}

void solve() {
  back(-1, 0, total);
}

void output() {
  ofstream out("elimin.out");
  out<<MAX;
  out.close();
}

int main() {
  input();
  solve();
  output();
  return 0;
}