// https://oj.uz/problem/view/IOI11_garden > p326
#include "garden.h"
#include "gardenlib.h"
#include <bits/stdc++.h>
using namespace std;
const int MAXN = 150000;
const int MAXLOG = 30;
const int MAXM = 150000;
vector<vector<int>> adj;
int pai[2*MAXN];
vector<int> pai_inv[2*MAXN];
bool in_cycle[2*MAXN];
int cycle_pai[2*MAXN];
int nivel[2*MAXN];
int cycle_sz[2*MAXN];
int in_grau[2*MAXN];
int state[2*MAXN];
int cycle_idx[2*MAXN];
int idx_in_cycle[2*MAXN];
bool mark[2*MAXN];
int qty_cycles = 0;
int n, m;
void map_cycle(int src) {
int cycle_id = qty_cycles++;
int v = src;
int v_idx_in_cycle = 0;
while (state[v] == 0) {
state[v] = 2;
in_cycle[v] = true;
cycle_pai[v] = v;
nivel[v] = 0;
cycle_idx[v] = cycle_id;
idx_in_cycle[v] = v_idx_in_cycle;
cycle_sz[cycle_id]++;
v = pai[v];
v_idx_in_cycle++;
}
}
void top_sort() {
queue<int> rm; vector<int> top_order;
for (int i=0; i<2*n; i++) {
if (in_grau[i] == 0) rm.push(i);
}
while (!rm.empty()) {
int v = rm.front(); rm.pop();
top_order.push_back(v);
in_grau[pai[v]]--;
if (in_grau[pai[v]] == 0) rm.push(pai[v]);
state[v] = 3; in_cycle[v] = false;
}
for (int i=0; i<2*n; i++) {
if (state[i] == 0) {
map_cycle(i);
}
}
for (int i=(int)top_order.size()-1; i>=0; i--) {
int v = top_order[i];
cycle_pai[v] = cycle_pai[pai[v]];
nivel[v] = nivel[pai[v]]+1;
cycle_idx[v] = cycle_idx[pai[v]];
idx_in_cycle[v] = -1;
}
}
pair<int,int> calc_pai(pair<int,int> u) {
int e = 0;
if (u.second == 1 and (int)adj[u.first].size() > 1) { // já veio de grande
e = 1;
}
int v = adj[u.first][e];
bool is_v_biggest = false;
if (adj[v][0] == u.first) is_v_biggest = true;
return {v, is_v_biggest};
};
int cycle_dist(int u, int v) { // u até v
int sz = cycle_sz[cycle_idx[u]];
int idx_u = idx_in_cycle[u];
int idx_v = idx_in_cycle[v];
if (idx_u <= idx_v) return idx_v - idx_u;
else return (sz - idx_u + idx_v);
}
void BFS(int src, int nivel_target) {
bool vis[2*n];
for (int i=0; i<2*n; i++) vis[i] = false;
queue<int> visit;
vis[src] = true;
visit.push(src);
while (!visit.empty()) {
int v = visit.front(); visit.pop();
if (nivel[v] == nivel_target and (v&1)==0) {
mark[v] = true;
continue;
}
for (int viz : pai_inv[v]) {
if (vis[viz]) continue;
vis[viz] = true;
visit.push(viz);
}
}
}
void mark_ans(int dest, int steps) {
if (in_cycle[dest]) {
for (int i=0; i<2*n; i++) {
if (i & 1) continue;
if (cycle_idx[i] != cycle_idx[dest]) continue;
int first = nivel[i] + cycle_dist(cycle_pai[i], dest);
if (steps < first) continue;
if (((steps - first) % cycle_sz[cycle_idx[dest]]) == 0) {
mark[i] = true;
}
}
} else {
BFS(dest, steps+nivel[dest]);
}
}
void count_routes(int N, int M, int P, int R[][2], int Q, int G[])
{
/*
• P – destino que conta se chegar
• R[M][2] – arestas ordenadas de forma decrescente pela beleza (são distintas)
• G[Q] qnts arestas cada grupo passará.
*/
n = N; m = M;
for (int i=0; i<2*n; i++) {
pai_inv[i].clear();
in_grau[i] = 0;
state[i] = 0;
cycle_sz[i] = 0;
}
adj = vector<vector<int>>(n);
qty_cycles = 0;
for (int e=0; e<m; e++) {
// já ficam ordenadas
int u = R[e][0], v = R[e][1];
adj[u].push_back(v);
adj[v].push_back(u);
}
for (int i=0; i<n; i++) {
for (int stat=0; stat<=1; stat++) {
int idx = 2*i + stat;
pair<int,int> pai_aux = calc_pai({i, stat});
pai[idx] = 2*pai_aux.first + pai_aux.second;
in_grau[pai[idx]]++;
pai_inv[pai[idx]].push_back(idx);
}
}
top_sort();
for(int g=0; g<Q; g++) {
// quantos lugares saem e chegam em P
int ans = 0;
mark_ans(2*P, G[g]);
mark_ans(2*P+1, G[g]);
for (int i=0; i<2*n; i++) {
ans += mark[i];
mark[i] = false;
}
answer(ans);
}
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
7516 KB |
Output is correct |
2 |
Correct |
3 ms |
7516 KB |
Output is correct |
3 |
Correct |
3 ms |
7516 KB |
Output is correct |
4 |
Correct |
3 ms |
7516 KB |
Output is correct |
5 |
Correct |
3 ms |
7516 KB |
Output is correct |
6 |
Correct |
3 ms |
7580 KB |
Output is correct |
7 |
Correct |
3 ms |
7516 KB |
Output is correct |
8 |
Correct |
3 ms |
7516 KB |
Output is correct |
9 |
Correct |
4 ms |
7776 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
7516 KB |
Output is correct |
2 |
Correct |
3 ms |
7516 KB |
Output is correct |
3 |
Correct |
3 ms |
7516 KB |
Output is correct |
4 |
Correct |
3 ms |
7516 KB |
Output is correct |
5 |
Correct |
3 ms |
7516 KB |
Output is correct |
6 |
Correct |
3 ms |
7580 KB |
Output is correct |
7 |
Correct |
3 ms |
7516 KB |
Output is correct |
8 |
Correct |
3 ms |
7516 KB |
Output is correct |
9 |
Correct |
4 ms |
7776 KB |
Output is correct |
10 |
Correct |
3 ms |
7516 KB |
Output is correct |
11 |
Correct |
10 ms |
11612 KB |
Output is correct |
12 |
Correct |
21 ms |
14212 KB |
Output is correct |
13 |
Correct |
27 ms |
23856 KB |
Output is correct |
14 |
Correct |
75 ms |
31168 KB |
Output is correct |
15 |
Correct |
76 ms |
31496 KB |
Output is correct |
16 |
Correct |
59 ms |
24272 KB |
Output is correct |
17 |
Correct |
49 ms |
21776 KB |
Output is correct |
18 |
Correct |
20 ms |
14288 KB |
Output is correct |
19 |
Correct |
95 ms |
31220 KB |
Output is correct |
20 |
Correct |
80 ms |
31428 KB |
Output is correct |
21 |
Correct |
63 ms |
24296 KB |
Output is correct |
22 |
Correct |
54 ms |
21764 KB |
Output is correct |
23 |
Correct |
87 ms |
33456 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
7516 KB |
Output is correct |
2 |
Correct |
3 ms |
7516 KB |
Output is correct |
3 |
Correct |
3 ms |
7516 KB |
Output is correct |
4 |
Correct |
3 ms |
7516 KB |
Output is correct |
5 |
Correct |
3 ms |
7516 KB |
Output is correct |
6 |
Correct |
3 ms |
7580 KB |
Output is correct |
7 |
Correct |
3 ms |
7516 KB |
Output is correct |
8 |
Correct |
3 ms |
7516 KB |
Output is correct |
9 |
Correct |
4 ms |
7776 KB |
Output is correct |
10 |
Correct |
3 ms |
7516 KB |
Output is correct |
11 |
Correct |
10 ms |
11612 KB |
Output is correct |
12 |
Correct |
21 ms |
14212 KB |
Output is correct |
13 |
Correct |
27 ms |
23856 KB |
Output is correct |
14 |
Correct |
75 ms |
31168 KB |
Output is correct |
15 |
Correct |
76 ms |
31496 KB |
Output is correct |
16 |
Correct |
59 ms |
24272 KB |
Output is correct |
17 |
Correct |
49 ms |
21776 KB |
Output is correct |
18 |
Correct |
20 ms |
14288 KB |
Output is correct |
19 |
Correct |
95 ms |
31220 KB |
Output is correct |
20 |
Correct |
80 ms |
31428 KB |
Output is correct |
21 |
Correct |
63 ms |
24296 KB |
Output is correct |
22 |
Correct |
54 ms |
21764 KB |
Output is correct |
23 |
Correct |
87 ms |
33456 KB |
Output is correct |
24 |
Correct |
4 ms |
7512 KB |
Output is correct |
25 |
Correct |
193 ms |
11612 KB |
Output is correct |
26 |
Correct |
340 ms |
15836 KB |
Output is correct |
27 |
Correct |
1071 ms |
25232 KB |
Output is correct |
28 |
Correct |
1427 ms |
31616 KB |
Output is correct |
29 |
Correct |
1717 ms |
31940 KB |
Output is correct |
30 |
Correct |
928 ms |
25292 KB |
Output is correct |
31 |
Correct |
1002 ms |
22728 KB |
Output is correct |
32 |
Correct |
184 ms |
15820 KB |
Output is correct |
33 |
Correct |
1364 ms |
31688 KB |
Output is correct |
34 |
Correct |
1681 ms |
32052 KB |
Output is correct |
35 |
Correct |
684 ms |
25544 KB |
Output is correct |
36 |
Correct |
3673 ms |
22728 KB |
Output is correct |
37 |
Correct |
1528 ms |
33716 KB |
Output is correct |
38 |
Correct |
3252 ms |
36680 KB |
Output is correct |