#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #pragma warning(disable:4996) typedef long long ll; #define INF 2140000000 using namespace std; // // G[from]={toの集合} というグラフと始点 s を入力。 // 始点sから各点への距離を求める。pprは最短距離のときのprevを記録 // vector bfs(int s, const vector>& G, vector* ppr) { queue que; vector d(G.size(), INF); if(ppr) ppr->resize(G.size(), -1); d[s] = 0; que.push(s); while (!que.empty()) { int curr = que.front(); que.pop(); int i; for (i = 0; i < (int)G[curr].size(); i++) { int next = G[curr][i]; if (d[curr] + 1 < d[next]) { d[next] = d[curr] + 1; if(ppr) (*ppr)[next] = curr; que.push(next); } } } return d; } void solve() { int n, m; cin >> n >> m; vector> g(n*2); // g[from]={to} の形のグラフ。ただし頂点を倍加(始点からの移動距離が偶数/奇数を区別) vector cntin(n); // 入次数をカウント map, int> ed; // (from,to)->edgeid int i, j; for (i = 0; i < m; i++) { int a, b; cin >> a >> b; a--; b--; ed[make_pair(a, b)] = i; cntin[b]++; g[a*2].push_back(b*2+1); g[a*2+1].push_back(b*2); } vector ans; // 長さnの0,1からなる配列。0,1それぞれR,Gを表す for (i = 0; i < n; i++) { vector pr; vector dist = bfs(2 * i, g, &pr); for (j = 0; j < n; j++) { if (dist[j * 2] == INF && dist[j * 2 + 1] == INF) { printf("-1"); return; } } if (ans.empty() && cntin[i] > 1 && dist[i * 2 + 1] < INF) { // 奇数閉路で、かつ i に他からも有向辺が来ている vector flag(n); // loopに使われたnodeに印をつける vector loop; int curr = i * 2 + 1; int bad = 0; while (curr != i * 2 ) { if (flag[curr/2]) { bad = 1; break; // 同じnodeを2回通るループは除外 } else { flag[curr / 2] = 1; } int prev = pr[curr]; int edge = ed[make_pair(prev / 2, curr / 2)]; loop.push_back(edge); curr = prev; } if (bad == 0) { vector ans0(m); // 長さmの0,1からなる配列。0,1それぞれR,Gを表す for(j=0; j<=(int)loop.size()/2; j++){ // ループの長さを 2k+1 とするとき、後半の(k+1)個だけ1をセット。それ以外は全て0 ans0[loop[j]]=1; } ans.swap(ans0); } } } if (ans.empty()) { cout << "-1" << endl; return; } for(i=0;i