#include using namespace std; struct SCC{ int N; vector used; vector vs; vector cmp; void dfs(int v, const vector G[]){ used[v] = true; for(auto u : G[v]){ if(!used[u]) dfs(u, G); } vs.push_back(v); } void rdfs(int v, int k, const vector>& rG){ used[v] = true; cmp[v] = k; for(auto u : rG[v]){ if(!used[u]) rdfs(u, k, rG); } } SCC(int N, const vector G[]):N(N){ vector> rG(N); for(int i=0; i=0; i--){ if(!used[vs[i]]) rdfs(vs[i], k++, rG); } } int dump_graph(const vector G[], vector edges2[]){ for(int i=0; i> N; vector S(N), T(N), U(N); for(int i=0; i> S[i]; for(int i=0; i> T[i]; for(int i=0; i> U[i]; static vector edges[500*500*2]; auto encode = [&](int i, int j, int b){ return b*N*N+i*N+j; }; for(int i=0; i