aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAnthony Wang2021-03-12 20:59:30 -0600
committerAnthony Wang2021-03-12 20:59:30 -0600
commit334ce6d602b385dabafd8d12374efcf827b7e318 (patch)
tree079fa6ffe9236aa303794444fd1fc7323a27e92e
parent1dcc5b4273c85a9453af0ecfbab365e1f1fb9a16 (diff)
Add minimizing.cpp
-rw-r--r--21/feb/plat/minimizing.cpp78
1 files changed, 78 insertions, 0 deletions
diff --git a/21/feb/plat/minimizing.cpp b/21/feb/plat/minimizing.cpp
new file mode 100644
index 0000000..9d5bf95
--- /dev/null
+++ b/21/feb/plat/minimizing.cpp
@@ -0,0 +1,78 @@
+#include <bits/stdc++.h>
+#define f first
+#define s second
+using namespace std;
+using ll = long long;
+using ii = pair<int, int>;
+constexpr int MX = 1e5+5;
+
+int ans, D[2*MX];
+vector<int> G[MX], H[2*MX];
+
+int main() {
+ if (fopen("in", "r")) freopen("in", "r", stdin), freopen("out", "w", stdout);
+ cin.tie(0)->sync_with_stdio(0);
+
+ int T; cin >> T;
+ while (T--) {
+ int N, M; cin >> N >> M;
+ for (int u = 1; u <= N; ++u) G[u].clear();
+ for (int u = 1; u <= 2*N; ++u) H[u].clear();
+ for (int i = 0; i < M; ++i) {
+ int x, y; cin >> x >> y;
+ G[x].push_back(y), G[y].push_back(x);
+ H[x].push_back(y+N), H[y+N].push_back(x);
+ H[y].push_back(x+N), H[x+N].push_back(y);
+ }
+
+ for (int u = 1; u <= 2*N; ++u) D[u] = 1e9;
+ queue<int> q;
+ D[1] = 0, q.push(1);
+ while (q.size()) {
+ int u = q.front(); q.pop();
+ for (int v : H[u]) if (D[v] == 1e9)
+ D[v] = D[u]+1, q.push(v);
+ }
+
+ if (*max_element(D+1, D+2*N+1) == 1e9) {
+ cout << N-1 << '\n';
+ continue;
+ }
+
+ /*vector<ii> C;
+ C.emplace_back(0, 0);
+ for (int u = 1; u <= N; ++u) {
+ C.emplace_back(D[u], D[u+N]);
+ }*/
+
+ ans = 0;
+ map<ii, int> f;
+ map<ii, vector<int>> b;
+ for (int u = 1; u <= N; ++u) {
+ ii p(min(D[u], D[u+N]), max(D[u], D[u+N]));
+ ++f[p], b[p].push_back(u);
+ }
+ map<ii, int> ea;
+ for (auto [p, c] : f) {
+ int pr = ea[ii(p.f-1, p.s+1)];
+ if (p.f+1 == p.s) {
+ if (p.f == 0) ans += (c+1)/2;
+ else if (f[ii(p.f-1, p.s-1)]) ans += max((c-pr)+(pr+1)/2, (c+1)/2);
+ else {
+ if (pr < c) ans += c-pr;
+ ans += (c+1)/2;
+ }
+ }
+ else {
+ ans += c;
+ if (p.f == 0) ea[p] = c;
+ else if (f[ii(p.f-1, p.s-1)]) ea[p] += min(c, pr);
+ else {
+ if (pr < c) ans += c-pr;
+ ea[p] = c;
+ }
+ }
+ }
+ cout << ans << '\n';
+ }
+}