#include <bits/stdc++.h>
using namespace std;
#pragma GCC target ("avx2")
#pragma GCC optimization ("O3")
#pragma GCC optimization ("unroll-loops")
#define fast ios::sync_with_stdio(0);cin.tie(0);
#define s second
#define f first
typedef long long ll;
const ll MOD = 998244353;
const ll LOGN = 20;
const ll MAXN = 5e6 + 101;
struct Edge {
int to;
short cost;
char cap;
Edge(int _to, short _cost, char _cap) : to(_to), cost(_cost), cap(_cap) { }
};
struct Flow {
vector<Edge> edges;
vector<vector<int>> graph;
vector<int> from;
vector<short> dist;
int n_, Q, source, terminal;
Flow(int n) {
Q = -1;
n_ = n;
graph = vector<vector<int>>(n + 10);
dist = vector<short>(n + 10);
from = vector<int>(n + 10);
}
void addEdge(int from, int to, int cost, int cap) {
Q++;
edges.push_back(Edge(to, cost, cap));
graph[from].push_back(Q);
Q++;
edges.push_back(Edge(from, -cost, 0));
graph[to].push_back(Q);
}
bool SPFA() {
queue<int> q;
vector<bool> que(n_ + 10, 0);
for (int i = 0; i < dist.size(); i++)
dist[i] = 10000;
dist[source] = 0;
que[source] = 1;
q.push(source);
while (!q.empty()) {
int node = q.front();
q.pop();
que[node] = false;
for (auto u : graph[node]) {
if (edges[u].cap && dist[edges[u].to] > edges[u].cost + dist[node]) {
dist[edges[u].to] = edges[u].cost + dist[node];
from[edges[u].to] = u;
if (!que[edges[u].to]) {
q.push(edges[u].to);
que[edges[u].to] = true;
}
}
}
}
return (dist[terminal] != 10000);
}
ll minCost() {
ll cost = 0;
while (SPFA()) {
cost += dist[terminal];
int now = terminal;
while (now != source) {
int edge_id = from[now];
int from = edges[edge_id ^ 1].to;
now = from;
edges[edge_id].cap--;
edges[edge_id ^ 1].cap++;
}
}
return cost;
}
};
short grid[2005][2005];
int main() {
fast
int N, K;
cin >> N >> K;
ll ans = 0;
for (int i = 0; i < N; i++) {
for (int j = 1; j <= N; j++) {
cin >> grid[i][j];
ans += grid[i][j];
}
}
Flow F = Flow(N * N + 10);
F.addEdge(N * N + 2, 0, 0, K);
for (int i = 0; i < N; i++) {
for (int j = 1; j <= N; j++) {
if ((i + 1 + j) % 2)
F.addEdge(0, i * N + j, -grid[i][j], 1);
}
}
for (int i = 0; i < N; i++) {
for (int j = 1; j <= N; j++) {
if ((i + 1 + j) % 2) {
if (j != N)
F.addEdge(i * N + j, i * N + j + 1, 0, 1);
if (i != N - 1)
F.addEdge(i * N + j, (i+1) * N + j, 0, 1);
if (j != 1)
F.addEdge(i * N + j, i * N + j - 1, 0, 1);
if (i != 0)
F.addEdge(i * N + j, (i-1) * N + j, 0, 1);
}
}
}
for (int i = 0; i < N; i++) {
for (int j = 1; j <= N; j++) {
if ((i + 1 + j) % 2 == 0)
F.addEdge(i * N + j, N * N + 1, -grid[i][j], 1);
}
}
F.source = N * N + 2, F.terminal = N * N + 1;
ans += F.minCost();
cout << ans << "\n";
}
Compilation message
domino.cpp:5: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
5 | #pragma GCC optimization ("O3")
|
domino.cpp:6: warning: ignoring '#pragma GCC optimization' [-Wunknown-pragmas]
6 | #pragma GCC optimization ("unroll-loops")
|
domino.cpp: In member function 'bool Flow::SPFA()':
domino.cpp:47:21: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<short int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
47 | for (int i = 0; i < dist.size(); i++)
| ~~^~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
62 ms |
37084 KB |
Output is correct |
2 |
Correct |
65 ms |
36820 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
2908 KB |
Output is correct |
2 |
Correct |
1 ms |
2908 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
642 ms |
524288 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
641 ms |
301500 KB |
Output is correct |
2 |
Correct |
602 ms |
300928 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
309 ms |
137124 KB |
Output is correct |
2 |
Correct |
286 ms |
136928 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
660 ms |
524288 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
2840 KB |
Output is correct |
2 |
Correct |
1 ms |
2908 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
633 ms |
524288 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
4 ms |
4048 KB |
Output is correct |
2 |
Correct |
4 ms |
3828 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
1 ms |
348 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
633 ms |
524288 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
3288 KB |
Output is correct |
2 |
Correct |
2 ms |
3288 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
415 ms |
137412 KB |
Output is correct |
2 |
Correct |
386 ms |
136892 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
635 ms |
524288 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |