結果

問題 No.3508 OR Mapping
コンテスト
ユーザー Nzt3
提出日時 2026-03-17 02:00:48
言語 C++23
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
WA  
実行時間 -
コード長 3,344 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,825 ms
コンパイル使用メモリ 355,948 KB
実行使用メモリ 226,312 KB
最終ジャッジ日時 2026-04-17 19:39:14
合計ジャッジ時間 28,329 ms
ジャッジサーバーID
(参考情報)
judge2_1 / judge1_1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 60 WA * 5
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code

#include <bits/stdc++.h>
using namespace std;

using ll = long long;
using pii = pair<int, int>;
using pll = pair<ll, ll>;
using vi = vector<int>;
using vl = vector<ll>;
#define rep3(i, a, b, c) for (ll i = (a); i < (b); i += (c))
#define rep2(i, a, b) rep3(i, a, b, 1)
#define rep1(i, n) rep2(i, 0, n)
#define rep0(n) rep1(aaaaa, n)
#define ov4(a, b, c, d, name, ...) name
#define rep(...) ov4(__VA_ARGS__, rep3, rep2, rep1, rep0)(__VA_ARGS__)
#define per(i, a, b) for (ll i = (a) - 1; i >= (b); i--)
#define fore(e, v) for (auto &&e : v)
#define all(a) begin(a), end(a)
#define sz(a) (int)(size(a))
#define lb(v, x) (lower_bound(all(v), x) - begin(v))
#define eb emplace_back

template <typename T, typename S> bool chmin(T &a, const S &b) {
  return a > b ? a = b, 1 : 0;
}
template <typename T, typename S> bool chmax(T &a, const S &b) {
  return a < b ? a = b, 1 : 0;
}
const int INF = 1e9 + 100;
const ll INFL = 3e18 + 100;
#define i128 __int128_t
struct _ {
  _() { cin.tie(0)->sync_with_stdio(0), cout.tie(0); }
} __;

template <typename G> struct SCC {
  G g; // vector<vi> or vector<vl>
  vector<vi> rg;
  vi comp, ord, used;
  int num; // 連結成分の数
  SCC(G g) : g(g), rg(sz(g)), comp(sz(g), -1), ord(sz(g)), used(sz(g)) {
    rep(i, sz(g)) fore(e, g[i]) rg[e].eb(i);
    build();
  }; // 4587f363
  int operator[](int k) { return comp[k]; }
  void dfs(int x) {
    if (used[x])
      return;
    used[x] = true;
    fore(e, g[x]) if (!used[e]) dfs(e);
    ord.eb(x);
  } // 14b1a925
  void rdfs(int x, int cnt) {
    if (comp[x] != -1)
      return;
    comp[x] = cnt;
    fore(e, rg[x]) if (comp[e] == -1) rdfs(e, cnt);
  } // de35b0e0
  void build() {
    rep(i, g.size()) dfs(i);
    reverse(all(ord));
    num = 0;
    fore(i, ord) if (comp[i] == -1) { rdfs(i, num), num++; }
  } // d3ed3d5a
};

int main() {
  int N, M, K;
  cin >> N >> M >> K;
  vector G(N, vi());
  rep(M) {
    int U, V;
    cin >> U >> V;
    --U, --V;
    G[U].push_back(V);
  }
  SCC G_SCC(G);
  auto comp = G_SCC.comp;
  int L = G_SCC.num;
  vector<set<int>> G_DAG(L), revG_DAG(L);

  auto ng = [&]() {
    cout << "No\n";
    exit(0);
  };
  vi dist(N, INF), oddcycle(L);

  // SCC Graphの構成
  rep(i, N) {
    if (dist[i] < INF)
      continue;
    int c = comp[i];
    queue<int> bfs;
    bfs.push(i);
    dist[i] = 0;
    while (bfs.size()) {
      int v = bfs.front();
      bfs.pop();
      for (auto u : G[v]) {
        if (comp[u] == c) {
          if (dist[u] > dist[v] + 1) {
            bfs.push(u);
            dist[u] = dist[v] + 1;
          } else if (dist[u] % 2 == dist[v] % 2) {
            oddcycle[c] = 1;
          }
        } else {
          G_DAG[c].insert(comp[u]);
          revG_DAG[comp[u]].insert(c);
        }
      }
    }
  }
  vi cnt(L);
  rep(i, N) cnt[comp[i]]++;
  // for(auto i:comp)cerr<<i<<' ';cerr<<endl;
  // for(auto i:cnt)cerr<<i<<' ';cerr<<endl;
  if (cnt[comp[0]] == 1)
    ng();
  rep(i, L - 1) if (cnt[i] == cnt[i + 1] && cnt[i] == 1) ng();
  rep(i, L) if (cnt[i] != 1 && !oddcycle[i]) ng();

  queue<int> bfs;
  vi G_DAG_outdeg(L);
  rep(i, L) G_DAG_outdeg[i] = sz(G_DAG[i]);
  rep(i, L - 1) if (G_DAG_outdeg[i] == 0) ng();
  per(i, L, 0) {
    for (auto u : revG_DAG[i]) {
      G_DAG_outdeg[u] -= 1;
      if (G_DAG_outdeg[u] == 0 && u != i - 1)
        ng();
    }
  }
  cout << "Yes\n";
}
0