結果
| 問題 | 
                            No.1660 Matrix Exponentiation
                             | 
                    
| コンテスト | |
| ユーザー | 
                             Nachia
                         | 
                    
| 提出日時 | 2021-08-27 21:47:49 | 
| 言語 | C++17(gcc12)  (gcc 12.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 88 ms / 2,000 ms | 
| コード長 | 1,837 bytes | 
| コンパイル時間 | 3,805 ms | 
| コンパイル使用メモリ | 129,424 KB | 
| 最終ジャッジ日時 | 2025-01-24 02:49:43 | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge5 / judge1 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 27 | 
ソースコード
#include <iostream>
#include <vector>
#include <algorithm>
#include <cmath>
using namespace std;
using i32 = int32_t;
using u32 = uint32_t;
using i64 = int64_t;
using u64 = uint64_t;
#define rep(i,n) for(int i=0; i<(n); i++)
#include <vector>
vector<vector<int>> SCC(const vector<vector<int>>& E){
  int n = E.size();
  vector<int> F(n, 0);
  vector<int> St;
  vector<int> I;
  for(int i=0; i<n; i++) if(F[i] == 0){
    St.push_back(i);
    while(St.size()){
      int p = St.back();
      St.pop_back();
      if(p < 0){
        I.push_back(-1-p);
        continue;
      }
      if(F[p] != 0) continue;
      F[p] = 1;
      St.push_back(-1-p);
      for(int e : E[p]) St.push_back(e);
    }
  }
  vector<vector<int>> R(n);
  vector<vector<int>> G;
  for(int i=0; i<n; i++) for(int e : E[i]){
    R[e].push_back(i);
  }
  for(int i=n-1; i>=0; i--){
    int s = I[i];
    if(F[s] == 0) continue;
    vector<int> g;
    St.push_back(s);
    while(St.size()){
      int p = St.back();
      St.pop_back();
      if(F[p] == 0) continue;
      F[p] = 0;
      g.push_back(p);
      for(int e : R[p]) St.push_back(e);
    }
    G.emplace_back(move(g));
  }
  for(int i=0; i<G.size(); i++){
    for(int p : G[i]) F[p] = i;
  }
  return move(G);
}
int main(){
  int N; cin >> N;
  int K; cin >> K;
  vector<vector<int>> E(N);
  rep(i,K){
    int u,v; cin >> u >> v; u--; v--;
    if(u == v){ cout << "-1\n"; return 0; }
    E[u].push_back(v);
  }
  auto G = SCC(E);
  if(G.size() != N){ cout << "-1\n"; return 0; }
  vector<int> D(N,1);
  int ans = 1;
  for(auto& A : G) for(auto p : A) for(int e : E[p]) D[e] = max(D[e], D[p]+1);
  rep(i,N) ans = max(ans,D[i]);
  cout << ans << endl;
  return 0;
}
struct ios_do_not_sync{
  ios_do_not_sync(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
  }
} ios_do_not_sync_instance;
            
            
            
        
            
Nachia