/**
* author: wxhtzdy
* created: 24.06.2022 12:30:23
**/
#include <bits/stdc++.h>
using namespace std;
string to_string(string s) {
return '"' + s + '"';
}
string to_string(const char* s) {
return to_string((string) s);
}
string to_string(bool b) {
return (b ? "true" : "false");
}
template <typename A, typename B>
string to_string(pair<A, B> p) {
return "(" + to_string(p.first) + ", " + to_string(p.second) + ")";
}
template <typename A>
string to_string(A v) {
bool first = true;
string res = "{";
for (const auto &x : v) {
if (!first) {
res += ", ";
}
first = false;
res += to_string(x);
}
res += "}";
return res;
}
void debug_out() { cerr << endl; }
template <typename Head, typename... Tail>
void debug_out(Head H, Tail... T) {
cerr << " " << to_string(H);
debug_out(T...);
}
#ifdef LOCAL
#define debug(...) cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__)
#else
#define debug(...) 42
#endif
int f[1 << 24][9];
int main() {
ios::sync_with_stdio(false);
cin.tie(0);
int n, k;
cin >> n >> k;
vector<vector<int>> a(n, vector<int>(n));
for (int i = 0; i < n; i++) {
for (int j = 0; j < n; j++) {
cin >> a[i][j];
}
}
vector<array<int, 3>> b;
auto Id = [&](int i, int j) {
return i * n + j;
};
for (int i = 0; i < n; i++) {
for (int j = 0; j < n; j++) {
if (i + 1 < n) {
b.push_back({a[i][j] + a[i + 1][j], Id(i, j), Id(i + 1, j)});
}
if (j + 1 < n) {
b.push_back({a[i][j] + a[i][j + 1], Id(i, j), Id(i, j + 1)});
}
}
}
sort(b.rbegin(), b.rend());
while ((int) b.size() > 51) {
b.pop_back();
}
int sz = (int) b.size();
int R = min(23, sz / 2);
int L = sz - R;
for (int s = 0; s < (1 << R); s++) {
for (int j = 0; j <= k; j++) {
f[s][j] = 0;
}
}
for (int s = 0; s < (1 << R); s++) {
for (int i = 0; i < R; i++) {
for (int j = 0; j <= k; j++) {
if (s >> i & 1) {
f[s][j] = max(f[s][j], f[s ^ (1 << i)][j]);
}
}
}
for (int p = 1; p <= k; p++) {
f[s][p] = max(f[s][p], f[s][p - 1]);
}
if (__builtin_popcount(s) >= k) {
continue;
}
for (int i = 0; i < R; i++) {
for (int p = 0; p < k; p++) {
if (s >> i & 1) {
continue;
}
bool ok = true;
for (int j = 0; j < R; j++) {
if (!(s >> j & 1)) {
continue;
}
if (b[L + i][1] == b[L + j][1]) {
ok = false;
}
if (b[L + i][1] == b[L + j][2]) {
ok = false;
}
if (b[L + i][2] == b[L + j][1]) {
ok = false;
}
if (b[L + i][2] == b[L + j][2]) {
ok = false;
}
}
if (ok) {
f[s | (1 << i)][p + 1] = max(f[s | (1 << i)][p], f[s][p] + b[L + i][0]);
}
}
}
for (int p = 1; p <= k; p++) {
f[s][p] = max(f[s][p], f[s][p - 1]);
}
}
int ans = 0;
for (int s = 0; s < (1 << R); s++) {
for (int j = 0; j <= k; j++) {
ans = max(ans, f[s][j]);
}
}
for (int s = 0; s < (1 << L); s++) {
if (__builtin_popcount(s) > k) {
continue;
}
int t = (1 << R) - 1;
bool ok = true;
for (int i = 0; i < L; i++) {
if (s >> i & 1) {
for (int j = i + 1; j < L; j++) {
if (s >> j & 1) {
if (b[i][1] == b[j][1]) {
ok = false;
}
if (b[i][1] == b[j][2]) {
ok = false;
}
if (b[i][2] == b[j][1]) {
ok = false;
}
if (b[i][2] == b[j][2]) {
ok = false;
}
}
}
}
}
if (!ok) {
continue;
}
for (int i = 0; i < L; i++) {
if (s >> i & 1) {
for (int j = 0; j < R; j++) {
if (t >> j & 1) {
bool found = false;
if (b[i][1] == b[L + j][1]) {
found = true;
}
if (b[i][1] == b[L + j][2]) {
found = true;
}
if (b[i][2] == b[L + j][1]) {
found = true;
}
if (b[i][2] == b[L + j][2]) {
found = true;
}
if (found) {
t ^= (1 << j);
}
}
}
}
}
int ft = 0;
for (int i = 0; i < L; i++) {
if (s >> i & 1) {
ft += b[i][0];
}
}
ans = max(ans, ft + f[t][k - __builtin_popcount(s)]);
}
int sum = 0;
for (int i = 0; i < n; i++) {
for (int j = 0; j < n; j++) {
sum += a[i][j];
}
}
cout << sum - ans << '\n';
return 0;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1882 ms |
302732 KB |
Output is correct |
2 |
Correct |
1881 ms |
302692 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2099 ms |
295964 KB |
Output is correct |
2 |
Correct |
2028 ms |
295964 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3850 ms |
405620 KB |
Output is correct |
2 |
Correct |
3420 ms |
405504 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2442 ms |
295904 KB |
Output is correct |
2 |
Correct |
2302 ms |
295804 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3344 ms |
357592 KB |
Output is correct |
2 |
Correct |
3188 ms |
357496 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3058 ms |
323292 KB |
Output is correct |
2 |
Correct |
2995 ms |
323292 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
4022 ms |
405504 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3773 ms |
295964 KB |
Output is correct |
2 |
Correct |
3593 ms |
295968 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
4118 ms |
405464 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
4078 ms |
296136 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
5 ms |
468 KB |
Output is correct |
2 |
Correct |
6 ms |
340 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
4107 ms |
405512 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
4081 ms |
295976 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
4066 ms |
323328 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
897 ms |
37236 KB |
Output is correct |
2 |
Correct |
902 ms |
37240 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
4122 ms |
405460 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |