結果
問題 | No.1660 Matrix Exponentiation |
ユーザー | tokusakurai |
提出日時 | 2021-08-27 22:08:33 |
言語 | C++17 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 112 ms / 2,000 ms |
コード長 | 3,759 bytes |
コンパイル時間 | 2,578 ms |
コンパイル使用メモリ | 218,644 KB |
実行使用メモリ | 26,728 KB |
最終ジャッジ日時 | 2023-08-13 09:07:59 |
合計ジャッジ時間 | 4,569 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge15 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
4,380 KB |
testcase_01 | AC | 2 ms
4,384 KB |
testcase_02 | AC | 2 ms
4,384 KB |
testcase_03 | AC | 1 ms
4,384 KB |
testcase_04 | AC | 1 ms
4,384 KB |
testcase_05 | AC | 1 ms
4,384 KB |
testcase_06 | AC | 2 ms
4,380 KB |
testcase_07 | AC | 1 ms
4,380 KB |
testcase_08 | AC | 1 ms
4,384 KB |
testcase_09 | AC | 11 ms
12,268 KB |
testcase_10 | AC | 11 ms
12,432 KB |
testcase_11 | AC | 9 ms
11,412 KB |
testcase_12 | AC | 2 ms
4,388 KB |
testcase_13 | AC | 2 ms
4,384 KB |
testcase_14 | AC | 1 ms
4,380 KB |
testcase_15 | AC | 2 ms
4,380 KB |
testcase_16 | AC | 2 ms
4,384 KB |
testcase_17 | AC | 1 ms
4,380 KB |
testcase_18 | AC | 2 ms
4,380 KB |
testcase_19 | AC | 56 ms
14,312 KB |
testcase_20 | AC | 40 ms
12,992 KB |
testcase_21 | AC | 46 ms
9,236 KB |
testcase_22 | AC | 8 ms
8,980 KB |
testcase_23 | AC | 25 ms
10,172 KB |
testcase_24 | AC | 112 ms
26,728 KB |
testcase_25 | AC | 47 ms
16,996 KB |
testcase_26 | AC | 108 ms
26,624 KB |
testcase_27 | AC | 40 ms
17,760 KB |
testcase_28 | AC | 80 ms
17,976 KB |
testcase_29 | AC | 95 ms
21,424 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep(i, n) for(int i = 0; i < n; i++) #define rep2(i, x, n) for(int i = x; i <= n; i++) #define rep3(i, x, n) for(int i = x; i >= n; i--) #define each(e, v) for(auto &e: v) #define pb push_back #define eb emplace_back #define all(x) x.begin(), x.end() #define rall(x) x.rbegin(), x.rend() #define sz(x) (int)x.size() using ll = long long; using pii = pair<int, int>; using pil = pair<int, ll>; using pli = pair<ll, int>; using pll = pair<ll, ll>; const int MOD = 1000000007; //const int MOD = 998244353; const int inf = (1<<30)-1; const ll INF = (1LL<<60)-1; template<typename T> bool chmax(T &x, const T &y) {return (x < y)? (x = y, true) : false;}; template<typename T> bool chmin(T &x, const T &y) {return (x > y)? (x = y, true) : false;}; struct io_setup{ io_setup(){ ios_base::sync_with_stdio(false); cin.tie(NULL); cout << fixed << setprecision(15); } } io_setup; template<bool directed = true> struct Graph{ struct edge{ int to, id; edge(int to, int id) : to(to), id(id) {} }; vector<vector<edge>> es; const int n; int m; vector<int> deg; Graph(int n) : es(n), n(n), m(0), deg(n, 0) {} void add_edge(int from, int to){ es[from].emplace_back(to, m); if(!directed) es[to].emplace_back(from, m); deg[to]++; m++; } void solve(){ vector<int> dp(n, 0); queue<int> que; rep(i, n){ if(deg[i] == 0) que.emplace(i); } while(!empty(que)){ int i = que.front(); que.pop(); each(e, es[i]){ chmax(dp[e.to], dp[i]+1); if(--deg[e.to] == 0){ que.emplace(e.to); } } } cout << *max_element(all(dp))+1 << '\n'; } }; template<bool directed = true> struct Strongly_Connected_Components{ struct edge{ int to, id; edge(int to, int id) : to(to), id(id) {} }; vector<vector<edge>> es, rs; vector<int> vs, comp; vector<bool> used; const int n; int m; Strongly_Connected_Components(int n) : es(n), rs(n), vs(n), comp(n), used(n), n(n), m(0) {} void add_edge(int from, int to){ es[from].emplace_back(to, m), rs[to].emplace_back(from, m); if(!directed) es[to].emplace_back(from, m), rs[from].emplace_back(to, m); m++; } void _dfs(int now){ used[now] = true; for(auto &e: es[now]){ if(!used[e.to]) _dfs(e.to); } vs.push_back(now); } void _rdfs(int now, int cnt){ used[now] = true, comp[now] = cnt; for(auto &e: rs[now]){ if(!used[e.to]) _rdfs(e.to, cnt); } } Graph<true> decompose(){ fill(begin(used), end(used), false); for(int i = 0; i < n; i++){ if(!used[i]) _dfs(i); } fill(begin(used), end(used), false), reverse(begin(vs), end(vs)); int cnt = 0; for(auto &e: vs){ if(!used[e]) _rdfs(e, cnt++); } Graph<true> G(cnt); for(int i = 0; i < n; i++){ for(auto &e: es[i]){ int u = comp[i], v = comp[e.to]; if(u != v) G.add_edge(u, v); } } return G; } int operator [] (int k) const {return comp[k];} }; int main(){ int N, K; cin >> N >> K; Strongly_Connected_Components G(N); bool flag = true; rep(i, K){ int u, v; cin >> u >> v; u--, v--; if(u == v) flag = false; G.add_edge(u, v); } Graph H = G.decompose(); if(!flag || H.n != N){ cout << "-1\n"; return 0; } H.solve(); }