#include <bits/stdc++.h>
#include "garden.h"
#include "gardenlib.h"
using namespace std;
#define type(x) __typeof((x).begin())
#define foreach(it, x) for(type(x) it = (x).begin(); it != (x).end(); it++)
typedef long long ll;
typedef pair < int, int > ii;
const int inf = 1e9 + 333;
const ll linf = 1e18 + inf;
const int N = 3e5 + 5;
int n, m, p;
bool was[N], h[N];
int go[N];
vector < int > come[N];
vector < ii > v[N];
vector < int > cycle;
int findCycle(int x) {
h[x] = 1;
cycle.push_back(x);
if(!h[go[x]])
return findCycle(go[x]);
return go[x];
}
int curOut, curSize, cycleCnt, tme;
int out[N], dep[N], cycleSize[N], whichCycle[N], st[N], nd[N], onCycle[N];
void dfs(int x, int d) {
h[x] = 1;
st[x] = ++tme;
whichCycle[x] = cycleCnt;
out[x] = curOut;
dep[x] = d;
cycleSize[x] = curSize;
foreach(it, come[x]) {
int u = *it;
if(!h[u])
dfs(u, d + 1);
}
nd[x] = tme;
}
void init() {
for(int i = 1; i <= n; i++) {
for(int it = 0; it < v[i].size(); it++) {
int u = v[i][it].first;
int e = v[i][it].second;
was[i + (it * n)] = 1;
if(e == v[u][0].second and v[u].size() > 1) {
go[i + (it * n)] = u + n;
come[u + n].push_back(i + (it * n));
}
else {
go[i + (it * n)] = u;
come[u].push_back(i + (it * n));
}
}
}
for(int i = 1; i <= 2 * n; i++) {
//printf("go[%d] = %d\n", i, go[i]);
}
for(int i = 1; i <= 2 * n; i++)
onCycle[i] = -1;
for(int i = 1; i <= 2 * n; i++) {
if(was[i] and !h[i]) {
cycleCnt++;
cycle.clear();
int cyc = findCycle(i);
foreach(it, cycle) {
int x = *it;
h[x] = 0;
}
reverse(cycle.begin(), cycle.end());
while(cycle.back() != cyc)
cycle.pop_back();
reverse(cycle.begin(), cycle.end());
for(int it = 0; it < cycle.size(); it++) {
onCycle[cycle[it]] = it;
}
foreach(it, cycle) {
int x = *it;
h[x] = 1;
}
curSize = cycle.size();
foreach(it, cycle) {
int x = *it;
curOut = x;
dfs(x, 0);
}
}
}
}
bool check(int x, int y, int k) {
if(whichCycle[x] != whichCycle[y])
return 0;
if(onCycle[y] == -1) {
if(out[x] == out[y] and st[y] <= st[x] and st[x] <= nd[y])
return dep[x] - dep[y] == k;
return 0;
}
int dist = dep[x], size = cycleSize[x];
x = out[x];
x = onCycle[x];
y = onCycle[y];
if(x <= y)
dist += y - x;
else
dist += size - x + y;
if(dist > k)
return 0;
return (k - dist) % size == 0;
}
int solve(int x) {
int ans = 0;
for(int i = 1; i <= n; i++) {
bool flag = 0;
flag |= check(i, p, x);
if(was[p + n])
flag |= check(i, p + n, x);
ans += flag;
}
return ans;
}
void count_routes(int N, int M, int P, int R[][2], int Q, int G[]) {
n = N;
m = M;
p = P + 1;
for(int i = 0; i < M; i++) {
int x = R[i][0] + 1;
int y = R[i][1] + 1;
if(v[x].size() < 2)
v[x].push_back(ii(y, i));
if(v[y].size() < 2)
v[y].push_back(ii(x, i));
}
init();
for(int i = 0; i < Q; i++) {
answer(solve(G[i]));
}
}
Compilation message
garden.cpp: In function 'void init()':
garden.cpp:54:22: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for(int it = 0; it < v[i].size(); it++) {
~~~^~~~~~~~~~~~~
garden.cpp:86:23: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for(int it = 0; it < cycle.size(); it++) {
~~~^~~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
16 ms |
14704 KB |
Output is correct |
2 |
Correct |
15 ms |
14688 KB |
Output is correct |
3 |
Correct |
15 ms |
14712 KB |
Output is correct |
4 |
Correct |
19 ms |
14452 KB |
Output is correct |
5 |
Correct |
15 ms |
14512 KB |
Output is correct |
6 |
Correct |
17 ms |
14808 KB |
Output is correct |
7 |
Correct |
15 ms |
14468 KB |
Output is correct |
8 |
Correct |
12 ms |
14660 KB |
Output is correct |
9 |
Correct |
18 ms |
14656 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
16 ms |
14704 KB |
Output is correct |
2 |
Correct |
15 ms |
14688 KB |
Output is correct |
3 |
Correct |
15 ms |
14712 KB |
Output is correct |
4 |
Correct |
19 ms |
14452 KB |
Output is correct |
5 |
Correct |
15 ms |
14512 KB |
Output is correct |
6 |
Correct |
17 ms |
14808 KB |
Output is correct |
7 |
Correct |
15 ms |
14468 KB |
Output is correct |
8 |
Correct |
12 ms |
14660 KB |
Output is correct |
9 |
Correct |
18 ms |
14656 KB |
Output is correct |
10 |
Correct |
16 ms |
14480 KB |
Output is correct |
11 |
Correct |
31 ms |
18012 KB |
Output is correct |
12 |
Correct |
61 ms |
20604 KB |
Output is correct |
13 |
Correct |
69 ms |
30616 KB |
Output is correct |
14 |
Correct |
213 ms |
34792 KB |
Output is correct |
15 |
Correct |
226 ms |
35000 KB |
Output is correct |
16 |
Correct |
168 ms |
29304 KB |
Output is correct |
17 |
Correct |
145 ms |
27264 KB |
Output is correct |
18 |
Correct |
56 ms |
20540 KB |
Output is correct |
19 |
Correct |
214 ms |
34792 KB |
Output is correct |
20 |
Correct |
214 ms |
35008 KB |
Output is correct |
21 |
Correct |
167 ms |
29340 KB |
Output is correct |
22 |
Correct |
151 ms |
27520 KB |
Output is correct |
23 |
Correct |
220 ms |
36692 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
16 ms |
14704 KB |
Output is correct |
2 |
Correct |
15 ms |
14688 KB |
Output is correct |
3 |
Correct |
15 ms |
14712 KB |
Output is correct |
4 |
Correct |
19 ms |
14452 KB |
Output is correct |
5 |
Correct |
15 ms |
14512 KB |
Output is correct |
6 |
Correct |
17 ms |
14808 KB |
Output is correct |
7 |
Correct |
15 ms |
14468 KB |
Output is correct |
8 |
Correct |
12 ms |
14660 KB |
Output is correct |
9 |
Correct |
18 ms |
14656 KB |
Output is correct |
10 |
Correct |
16 ms |
14480 KB |
Output is correct |
11 |
Correct |
31 ms |
18012 KB |
Output is correct |
12 |
Correct |
61 ms |
20604 KB |
Output is correct |
13 |
Correct |
69 ms |
30616 KB |
Output is correct |
14 |
Correct |
213 ms |
34792 KB |
Output is correct |
15 |
Correct |
226 ms |
35000 KB |
Output is correct |
16 |
Correct |
168 ms |
29304 KB |
Output is correct |
17 |
Correct |
145 ms |
27264 KB |
Output is correct |
18 |
Correct |
56 ms |
20540 KB |
Output is correct |
19 |
Correct |
214 ms |
34792 KB |
Output is correct |
20 |
Correct |
214 ms |
35008 KB |
Output is correct |
21 |
Correct |
167 ms |
29340 KB |
Output is correct |
22 |
Correct |
151 ms |
27520 KB |
Output is correct |
23 |
Correct |
220 ms |
36692 KB |
Output is correct |
24 |
Correct |
17 ms |
14564 KB |
Output is correct |
25 |
Correct |
292 ms |
18144 KB |
Output is correct |
26 |
Correct |
463 ms |
20616 KB |
Output is correct |
27 |
Correct |
2861 ms |
30688 KB |
Output is correct |
28 |
Correct |
2229 ms |
34956 KB |
Output is correct |
29 |
Correct |
3443 ms |
35320 KB |
Output is correct |
30 |
Correct |
1905 ms |
29464 KB |
Output is correct |
31 |
Correct |
3206 ms |
27516 KB |
Output is correct |
32 |
Correct |
1071 ms |
20576 KB |
Output is correct |
33 |
Correct |
2242 ms |
34936 KB |
Output is correct |
34 |
Correct |
3540 ms |
35304 KB |
Output is correct |
35 |
Correct |
3352 ms |
29476 KB |
Output is correct |
36 |
Correct |
3053 ms |
27664 KB |
Output is correct |
37 |
Correct |
2185 ms |
36976 KB |
Output is correct |
38 |
Correct |
4623 ms |
40484 KB |
Output is correct |