#include <bits/stdc++.h>
using namespace std;
#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 from, to, cost;
char cap;
Edge() { }
Edge(int _from, int _to, int _cost, char _cap) : from(_from), to(_to), cost(_cost), cap(_cap) { }
};
struct Flow {
vector<Edge> edges;
vector<vector<int>> graph;
vector<short> dist, from;
int Q, source, terminal;
Flow(int n) {
Q = -1;
graph = vector<vector<int>>(n + 10);
dist = vector<short>(n + 10);
from = vector<short>(n + 10);
}
void addEdge(int from, int to, int cost, int cap) {
Q++;
edges.push_back(Edge(from, to, cost, cap));
graph[from].push_back(Q);
Q++;
edges.push_back(Edge(to, from, -cost, 0));
graph[to].push_back(Q);
}
bool shortestPath() {
for (int i = 0; i < dist.size(); i++)
dist[i] = 10000;
dist[source] = 0;
priority_queue<pair<int,int>, vector<pair<int,int>>, greater<pair<int,int>>> pq;
pq.push({0, source});
while (!pq.empty()) {
int node = pq.top().s;
int d = pq.top().f;
pq.pop();
for (auto u : graph[node]) {
if (edges[u].cap && dist[edges[u].to] > edges[u].cost + d) {
dist[edges[u].to] = edges[u].cost + d;
from[edges[u].to] = u;
pq.push({dist[edges[u].to], edges[u].to});
}
}
}
return (dist[terminal] != 10000);
}
ll minCost() {
ll cost = 0;
while (shortestPath()) {
cost += dist[terminal];
int now = terminal;
while (now != source) {
int edge_id = from[now];
int from = edges[edge_id].from;
now = from;
edges[edge_id].cap--;
edges[edge_id ^ 1].cap++;
}
}
return cost;
}
};
int 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: In member function 'bool Flow::shortestPath()':
domino.cpp:39:21: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<short int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
39 | for (int i = 0; i < dist.size(); i++)
| ~~^~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
190 ms |
54924 KB |
Output is correct |
2 |
Correct |
97 ms |
53424 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
2 ms |
3000 KB |
Output is correct |
2 |
Correct |
2 ms |
3032 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
763 ms |
524288 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
344 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
4167 ms |
524288 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
2597 ms |
204168 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
720 ms |
524288 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3 ms |
3032 KB |
Output is correct |
2 |
Correct |
2 ms |
3028 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
650 ms |
524288 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
9 ms |
8656 KB |
Execution killed with signal 11 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
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 |
693 ms |
524288 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
7 ms |
3536 KB |
Output is correct |
2 |
Correct |
5 ms |
3540 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
4065 ms |
205760 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
348 KB |
Output is correct |
2 |
Correct |
1 ms |
356 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Runtime error |
713 ms |
524288 KB |
Execution killed with signal 9 |
2 |
Halted |
0 ms |
0 KB |
- |