Submission #1248972

#TimeUsernameProblemLanguageResultExecution timeMemory
1248972eysbutnoTropical Garden (IOI11_garden)C++20
100 / 100
1303 ms21124 KiB
#include <bits/stdc++.h>
#include "garden.h"
#include "gardenlib.h"

using ll = long long;

constexpr int INF = 1e9;

void count_routes(int n, int m, int p, int r[][2], int q, int g[]) {
	std::vector<std::vector<std::array<int, 2>>> adj(n), rev(n);
	for (int i = 0; i < m; i++) {
		adj[r[i][0]].push_back({r[i][1], i});
		adj[r[i][1]].push_back({r[i][0], i});
	}

	for (int i = 0; i < n; i++) {
		if (adj[i].size() > 2) adj[i].resize(2);
		for (const auto &[j, w] : adj[i]) {
			rev[j].push_back({i, w});
		}
	}

	std::vector<std::array<int, 2>> dist(n, {INF, INF});
	std::vector<std::array<int, 2>> to_p(n, {INF, INF});
	for (int tt = 0; tt < 2; tt++) {
		// t = 0 means we take best edge from p
		// t = 1 means we don't take best edge 
		if (tt == 1 && adj[p].size() == 1) break;

		std::queue<std::array<int, 3>> bfs;
		bfs.push({0, p, tt});
		while (!bfs.empty()) {
			const auto [t, u, f] = bfs.front();
			bfs.pop();
			if (dist[u][f] != INF) continue;
			dist[u][f] = t;

			int wt = adj[u][0][1];
			for (const auto &[j, w] : rev[u]) {
				if (f == 0 && w == wt && adj[u].size() > 1) continue;
				if (f == 1 && w != wt && adj[u].size() > 1) continue;
				int flag = w != adj[j][0][1];
				bfs.push({t + 1, j, flag});
			}
		}

		for (int i = 0; i < n; i++) {
			to_p[i][tt] = dist[i][0];
		}

		dist.assign(n, {INF, INF});
	}

    const auto get_cycle = [&](int f) -> std::array<int, 2> {
        int cycle_len = 0, saw = 0;
        std::vector<std::array<bool, 2>> vis(n);

        int node = p, flag = f, trav = 0;
        while (!vis[node][flag]) {
            vis[node][flag] = true;
            if (node == p && flag == !f) saw = trav;

            const auto [nxt, wt] = adj[node][flag];
            if (adj[nxt].size() == 1) {
                node = nxt, flag = 0;
            } else {
                node = nxt, flag = (adj[nxt][0][1] == wt);
            }

            trav++;
        }

        if (node == p) cycle_len = trav;
	    else cycle_len = INT_MAX, saw = 0;
        return {cycle_len, saw};
    };

	std::array<int, 2> s1 = get_cycle(0);
    std::array<int, 2> s2 = {INT_MAX, 0};
    if (adj[p].size() > 1) s2 = get_cycle(1);

	for (int t = 0; t < q; t++) {
		int k = g[t], res = 0;
		for (int i = 0; i < n; i++) {
            if (to_p[i][0] <= k && ((k - to_p[i][0]) % s1[0] == 0 || (k - to_p[i][0]) % s1[0] == s1[1])) {
                res++;
            } else if (to_p[i][1] <= k && ((k - to_p[i][1]) % s2[0] == 0 || (k - to_p[i][1]) % s2[0] == s2[1])) {
                res++;
            }
		}

		// std::cout << res << '\n';
		answer(res);
	}
}
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...