#761842 (C++17) No.1473 おでぶなおばけさん

提出ソース

結果

問題	No.1473 おでぶなおばけさん
ユーザー	h-izu
提出日時	2022-05-18 00:08:49
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	265 ms / 2,000 ms
コード長	6,736 bytes
コンパイル時間	5,256 ms
コンパイル使用メモリ	284,380 KB
実行使用メモリ	16,640 KB
最終ジャッジ日時	2024-09-15 23:01:30
合計ジャッジ時間	11,320 ms
ジャッジサーバーID （参考情報）	judge4 / judge6

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	AC	2 ms 5,376 KB
testcase_02	AC	142 ms 12,672 KB
testcase_03	AC	92 ms 11,008 KB
testcase_04	AC	97 ms 8,064 KB
testcase_05	AC	39 ms 5,376 KB
testcase_06	AC	170 ms 12,288 KB
testcase_07	AC	148 ms 11,904 KB
testcase_08	AC	265 ms 11,904 KB
testcase_09	AC	132 ms 11,776 KB
testcase_10	AC	48 ms 7,424 KB
testcase_11	AC	61 ms 8,960 KB
testcase_12	AC	46 ms 8,320 KB
testcase_13	AC	27 ms 6,528 KB
testcase_14	AC	22 ms 5,504 KB
testcase_15	AC	47 ms 7,680 KB
testcase_16	AC	38 ms 6,940 KB
testcase_17	AC	4 ms 5,376 KB
testcase_18	AC	7 ms 5,376 KB
testcase_19	AC	53 ms 8,064 KB
testcase_20	AC	74 ms 10,368 KB
testcase_21	AC	91 ms 10,880 KB
testcase_22	AC	63 ms 13,568 KB
testcase_23	AC	50 ms 10,368 KB
testcase_24	AC	51 ms 11,776 KB
testcase_25	AC	258 ms 13,184 KB
testcase_26	AC	254 ms 13,824 KB
testcase_27	AC	53 ms 6,968 KB
testcase_28	AC	162 ms 12,032 KB
testcase_29	AC	96 ms 11,648 KB
testcase_30	AC	132 ms 13,056 KB
testcase_31	AC	100 ms 13,056 KB
testcase_32	AC	103 ms 12,288 KB
testcase_33	AC	72 ms 11,136 KB
testcase_34	AC	46 ms 7,552 KB
testcase_35	AC	39 ms 7,424 KB
testcase_36	AC	78 ms 10,112 KB
testcase_37	AC	116 ms 10,368 KB
testcase_38	AC	17 ms 5,376 KB
testcase_39	AC	39 ms 6,940 KB
testcase_40	AC	37 ms 6,912 KB
testcase_41	AC	37 ms 6,640 KB
testcase_42	AC	36 ms 6,640 KB
testcase_43	AC	46 ms 9,088 KB
testcase_44	AC	46 ms 9,088 KB
testcase_45	AC	46 ms 9,092 KB
testcase_46	AC	95 ms 14,336 KB
testcase_47	AC	117 ms 16,640 KB
testcase_48	AC	114 ms 16,384 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#ifdef __LOCAL
#define _GLIBCXX_DEBUG
#endif
#include <bits/stdc++.h>

#include <atcoder/all>
#include <vector>
using namespace std;
using namespace atcoder;

typedef unsigned long long ull;
typedef long long ll;

const double PI = 3.14159265358979323846;

#define rep(i, n) for (ll i = 0; i < (ll)(n); i++)

// a^p
// 2^3 = 2 * 2^2
// 2^2 = 2 * (2^1)
// 2^1 = 2
ll modpow(ll a, ll p, ll mod) {
  if (p == 0) return 1;

  if (p % 2 == 0) {
    ll half = modpow(a, p / 2, mod) % mod;
    return half * half % mod;
  } else {
    return a * modpow(a, p - 1, mod) % mod;
  }
}

// a^p
ll powpow(ll a, ll p) {
  if (p == 0) return 1;

  if (p % 2 == 0) {
    ll half = pow(a, p / 2);
    return half * half;
  } else {
    return a * pow(a, p - 1);
  }
}

// a/b
// https://qiita.com/drken/items/3b4fdf0a78e7a138cd9a
ll moddiv(ll a, ll b, ll mod) { return a * modpow(b, mod - 2, mod); }

// nCa を求める
ll modCombination(ll n, ll a, ll mod) {
  if (n < 0 || a < 0 || n < a) return 0;
  if (n - a < a) {
    return modCombination(n, n - a, mod);
  }

  ll denominator = 1;  // 分母
  ll numerator = 1;    // 分子

  for (ll i = 0; i < a; i++) {
    denominator *= a - i;
    numerator *= n - i;
    denominator %= mod;
    numerator %= mod;
  }

  return numerator * modpow(denominator, mod - 2, mod) % mod;
}

vector<vector<ll>> combination(ll n) {
  vector<vector<ll>> C(n + 1, vector<ll>(n + 1));
  C[0][0] = 1;
  rep(i, n) rep(j, i + 1) {
    C[i + 1][j + 1] += C[i][j];
    C[i + 1][j] += C[i][j];
  }

  return C;
}

// ref. https://drken1215.hatenablog.com/entry/2018/06/08/210000
class ModCombinationTale {
 private:
  ll n;
  ll mod;
  vector<ll> fac, finv, inv;

 public:
  ModCombinationTale(ll n, ll mod) : n(n), mod(mod) {
    fac.resize(n + 1);
    finv.resize(n + 1);
    inv.resize(n + 1);

    fac[0] = fac[1] = 1;
    finv[0] = finv[1] = 1;
    inv[1] = 1;

    for (ll i = 2; i <= n; i++) {
      fac[i] = fac[i - 1] * i % mod;
      inv[i] = mod - inv[mod % i] * (mod / i) % mod;
      finv[i] = finv[i - 1] * inv[i] % mod;
    }
  }

  ll operator()(ll n, ll k) {
    if (n < k) return 0;
    if (n < 0 || k < 0) return 0;
    return fac[n] * (finv[k] * finv[n - k] % mod) % mod;
  }
};

class UnionFind {
 private:
  vector<ll> parents;

 public:
  UnionFind(ll n) : parents(n, -1) {}

  bool issame(ll x, ll y) { return root(x) == root(y); }

  bool merge(ll x, ll y) {
    if (issame(x, y)) return false;

    ll rx = root(x);
    ll ry = root(y);
    if (parents[rx] > parents[ry]) swap(rx, ry);
    // サイズ情報を更新
    parents[rx] += parents[ry];
    // yの親を更新
    parents[ry] = rx;

    return true;
  }

  ll size(ll x) { return -parents[root(x)]; }

  ll root(ll x) {
    if (parents[x] < 0) return x;
    // 根の親の値に木の(-)サイズの情報を入れる
    return parents[x] = root(parents[x]);
  }
};

// cf. https://qiita.com/drken/items/a14e9af0ca2d857dad23
vector<ll> enum_divisors(ll n) {
  vector<ll> res;
  // sqrt(n)まで試し割り
  for (ll i = 1; i * i <= n; i++) {
    if (n % i == 0) {
      res.push_back(i);
      // 重複しないならばiの相方であるn/iも約数
      // e.g. n=25のときのi=5は重複
      if (n / i != i) res.push_back(n / i);
    }
  }

  sort(res.begin(), res.end());
  return res;
}

// cf. https://qiita.com/drken/items/a14e9af0ca2d857dad23
map<ll, ll> prime_factors(ll n) {
  map<ll, ll> res;
  // sqrt(n)まで試し割り
  for (ll a = 2; a * a <= n; a++) {
    if (n % a != 0) continue;

    // nで割れる限り割る
    while (n % a == 0) {
      res[a]++;
      n /= a;
    }
  }
  if (n != 1) res[n]++;

  return res;
}

ll gcd(ll a, ll b) {
  if (b == 0)
    return a;
  else
    return gcd(b, a % b);
}

ll lcm(ll a, ll b) { return (a / gcd(a, b)) * b; }

// cf. https://qiita.com/drken/items/56a6b68edef8fc605821
class AccumSum2D {
 private:
  vector<vector<ll>> sum;
  ll H;
  ll W;

 public:
  AccumSum2D(vector<vector<ll>> &A) {
    H = (ll)A.size();
    W = (ll)A[0].size();
    sum.resize(H + 1, vector<ll>(W + 1));

    for (ll i = 0; i < H; i++) {
      for (ll j = 0; j < W; j++) {
        sum[i + 1][j + 1] = sum[i][j + 1] + sum[i + 1][j] - sum[i][j] + A[i][j];
      }
    }
  }

  // クエリ [x1, x2) × [y1, y2) の長方形区域の和
  ll Sum(ll x1, ll x2, ll y1, ll y2) {
    return sum[x2][y2] - sum[x1][y2] - sum[x2][y1] + sum[x1][y1];
  }
};

// p/q
struct fraction {
  ll p, q;
  fraction(ll _p = 0, ll _q = 1) : p(_p), q(_q) {
    if (q == 0) {
      p = 1;
      return;
    }

    ll g = gcd(p, q);
    p /= g;
    q /= g;
  }

  bool operator<(const fraction &other) const {
    return p * other.q < q * other.p;
  }

  bool operator<=(const fraction &other) const {
    return p * other.q <= q * other.p;
  }

  bool operator==(const fraction &other) const {
    return p == other.p && q == other.q;
  }
};

// res[i][c] := i 文字目以降で最初に文字 c が登場する index (存在しないときは n)
vector<vector<ll>> calcNext(const string &S) {
  ll n = (ll)S.size();
  vector<vector<ll>> res(n + 1, vector<ll>(26, n));
  for (ll i = n - 1; i >= 0; --i) {
    for (ll j = 0; j < 26; ++j) res[i][j] = res[i + 1][j];
    res[i][S[i] - 'a'] = i;
  }
  return res;
}

struct Edge {
  ll to, d;
  Edge(ll to, ll d) : to(to), d(d) {}
};

int main() {
  cin.tie(nullptr);
  ios_base::sync_with_stdio(false);
  cout << fixed << setprecision(15);

  ll N, M;
  cin >> N >> M;

  vector G(N, vector<Edge>());
  rep(i, M) {
    ll s, t, d;
    cin >> s >> t >> d;
    s--;
    t--;
    G[s].emplace_back(t, d);
    G[t].emplace_back(s, d);
  }

  auto is_ok = [&](ll x) -> bool {
    vector<bool> visited(N);
    bool ret = false;
    auto dfs = [&](auto f, ll node) -> void {
      if (node == N - 1) {
        ret = true;
        // visited[node] = false;
        // cout << node << endl;
        return;
      }
      // cout << node << " -> ";
      visited[node] = true;
      for (auto e : G[node]) {
        if (visited[e.to]) continue;
        if (e.d < x) continue;
        f(f, e.to);
      }
      // visited[node] = false;
    };

    dfs(dfs, 0);

    return ret;
  };

  // cout << is_ok(4) << endl;

  ll ng = ll(1e9) + 1;
  ll ok = 0;
  while (abs(ok - ng) > 1) {
    ll mid = (ok + ng) / 2;
    if (is_ok(mid)) {
      ok = mid;
    } else {
      ng = mid;
    }
  }

  ll INF = 1e18;
  vector<ll> dist(N, INF);
  dist[0] = 0;
  queue<ll> q;
  q.push(0);
  while (q.size()) {
    ll node = q.front();
    q.pop();
    for (auto e : G[node]) {
      if (e.d < ok) continue;
      if (dist[e.to] > dist[node] + 1) {
        dist[e.to] = dist[node] + 1;
        q.push(e.to);
      }
    }
  }
  cout << ok << " " << dist[N - 1] << endl;

  return 0;
}

yukicoder

結果

テストケース

ソースコード