Problem p02703 (C++)

Buggy submission (status: Runtime Error):
#include <algorithm>
#include <bitset>
#include <cmath>
#include <functional>
#include <iomanip>
#include <iostream>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <unordered_set>
#include <vector>

using namespace std;
using ll = long long;
using ull = unsigned long long;
// const ll mod =  998244353 ;
const ll mod = 1000000007;
const ll inf = 1e18;
const int dx[] = {1, -1, 0, 0};
const int dy[] = {0, 0, 1, -1};

#define REP(i, n) for (ll i = 0; i < (n); ++i)
#define REP_FROM(i, j, n) for (ll i = (j); i < (n); ++i)
#define REP_REV(i, n) for (ll i = n - 1; i >= 0; --i)
#define all(x) (x).begin(), (x).end()
#define sz(x) int(x.size())

template <typename T> void print_vec(vector<T> a) {
  int n = a.size();
  REP(i, n) {
    cout >> a[i];
    if (i != n - 1)
      cout << " ";
    else
      cout << endl;
  }
}
ll power(ll base, ll exponent) {
  if (exponent % 2) {
    return power(base, exponent - 1) * base % mod;
  } else if (exponent) {
    ll root_ans = power(base, exponent / 2);
    return root_ans * root_ans % mod;
  } else {
    return 1;
  }
}
ll inverse(ll x) { return power(x, mod - 2); }
ll gcd(ll a, ll b) {
  if (a < b)
    gcd(b, a);
  ll r;
  while (r = a % b) {
    a = b;
    b = r;
  }
  return b;
}
template <typename T> ll sum(T begin, T end) {
  return accumulate(begin, end, 0ll);
}
struct combination {
  vector<ll> fact, inv;
  combination(int sz) : fact(sz + 1), inv(sz + 1) {
    fact[0] = 1;
    for (int i = 1; i <= sz; i++) {
      fact[i] = fact[i - 1] * i % mod;
    }
    inv[sz] = power(fact[sz], mod - 2);
    for (int i = sz - 1; i >= 0; i--) {
      inv[i] = inv[i + 1] * (i + 1) % mod;
    }
  }
  ll C(int p, int q) const {
    if (q < 0 || p < q)
      return 0;
    return (fact[p] * inv[q] % mod * inv[p - q] % mod);
  }
};

using Pair = pair<ll, ll>;

template <ll Modulus> struct ModInt {
  ll a;
  constexpr ModInt(const ll x = 0) noexcept : a((x % mod + mod) % mod) {}
  constexpr ll &value() noexcept { return a; }
  constexpr const ll &value() const noexcept { return a; }
  constexpr ModInt operator+(const ModInt x) const noexcept {
    return ModInt(*this) += x;
  }
  constexpr ModInt operator-(const ModInt x) const noexcept {
    return ModInt(*this) -= x;
  }
  constexpr ModInt operator*(const ModInt x) const noexcept {
    return ModInt(*this) *= x;
  }
  constexpr ModInt operator/(const ModInt x) const noexcept {
    return ModInt(*this) /= x;
  }
  constexpr ModInt &operator+=(const ModInt x) noexcept {
    a += x.a;
    if (a >= Modulus) {
      a -= Modulus;
    }
    return *this;
  }
  constexpr ModInt &operator-=(const ModInt x) noexcept {
    if (a < x.a) {
      a += Modulus;
    }
    a -= x.a;
    return *this;
  }
  constexpr ModInt &operator*=(const ModInt x) noexcept {
    a = a * x.a % Modulus;
    return *this;
  }
  constexpr ModInt &operator/=(ModInt x) noexcept {
    ll exp = Modulus - 2;
    while (exp) {
      if (exp % 2) {
        *this *= x;
      }
      x *= x;
      exp /= 2;
    }
    return *this;
  }
  constexpr ModInt operator-() noexcept { return ModInt(-a); }
  friend ostream &operator<<(ostream &os, const ModInt &m) {
    os << m.a;
    return os;
  }
};

using mint = ModInt<mod>;

signed main() {
  ios::sync_with_stdio(false);
  cout << fixed << setprecision(15);
  ll n, m, s;
  cin >> n >> m >> s;
  vector<vector<int>> a(n), b(n), g(n);
  REP(i, m) {
    int u, v, x, y;
    cin >> u >> v >> x >> y;
    u--;
    v--;
    g[u].push_back(v);
    a[u].push_back(x);
    b[u].push_back(y);

    g[v].push_back(u);
    a[v].push_back(x);
    b[v].push_back(y);
  }
  vector<int> c(n), d(n);
  REP(i, n) cin >> c[i] >> d[i];
  vector<vector<ll>> dp(n, vector<ll>(2501, inf));
  using Triple = tuple<ll, int, ll>;
  priority_queue<Triple, vector<Triple>, greater<>> q;
  q.emplace(0, 0, s);
  while (q.size()) {
    auto x = q.top();
    q.pop();
    ll time = get<0>(x);
    int node = get<1>(x);
    ll gold = get<2>(x);
    if (dp[node][gold] <= time)
      continue;
    dp[node][gold] = time;
    q.emplace(time + d[node], node, min(gold + c[node], 2500ll));
    REP(i, sz(g[node])) {
      if (gold < a[node][i])
        continue;
      q.emplace(time + b[node][i], g[node][i], gold - a[node][i]);
    }
  }
  REP_FROM(i, 1, n) { cout << *min_element(all(dp[i])) << endl; }
  return 0;
}


Fixed submission:
#include <algorithm>
#include <bitset>
#include <cmath>
#include <functional>
#include <iomanip>
#include <iostream>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <unordered_set>
#include <vector>

using namespace std;
using ll = long long;
using ull = unsigned long long;
// const ll mod =  998244353 ;
const ll mod = 1000000007;
const ll inf = 1e18;
const int dx[] = {1, -1, 0, 0};
const int dy[] = {0, 0, 1, -1};

#define REP(i, n) for (ll i = 0; i < (n); ++i)
#define REP_FROM(i, j, n) for (ll i = (j); i < (n); ++i)
#define REP_REV(i, n) for (ll i = n - 1; i >= 0; --i)
#define all(x) (x).begin(), (x).end()
#define sz(x) int(x.size())

template <typename T> void print_vec(vector<T> a) {
  int n = a.size();
  REP(i, n) {
    cout >> a[i];
    if (i != n - 1)
      cout << " ";
    else
      cout << endl;
  }
}
ll power(ll base, ll exponent) {
  if (exponent % 2) {
    return power(base, exponent - 1) * base % mod;
  } else if (exponent) {
    ll root_ans = power(base, exponent / 2);
    return root_ans * root_ans % mod;
  } else {
    return 1;
  }
}
ll inverse(ll x) { return power(x, mod - 2); }
ll gcd(ll a, ll b) {
  if (a < b)
    gcd(b, a);
  ll r;
  while (r = a % b) {
    a = b;
    b = r;
  }
  return b;
}
template <typename T> ll sum(T begin, T end) {
  return accumulate(begin, end, 0ll);
}
struct combination {
  vector<ll> fact, inv;
  combination(int sz) : fact(sz + 1), inv(sz + 1) {
    fact[0] = 1;
    for (int i = 1; i <= sz; i++) {
      fact[i] = fact[i - 1] * i % mod;
    }
    inv[sz] = power(fact[sz], mod - 2);
    for (int i = sz - 1; i >= 0; i--) {
      inv[i] = inv[i + 1] * (i + 1) % mod;
    }
  }
  ll C(int p, int q) const {
    if (q < 0 || p < q)
      return 0;
    return (fact[p] * inv[q] % mod * inv[p - q] % mod);
  }
};

using Pair = pair<ll, ll>;

template <ll Modulus> struct ModInt {
  ll a;
  constexpr ModInt(const ll x = 0) noexcept : a((x % mod + mod) % mod) {}
  constexpr ll &value() noexcept { return a; }
  constexpr const ll &value() const noexcept { return a; }
  constexpr ModInt operator+(const ModInt x) const noexcept {
    return ModInt(*this) += x;
  }
  constexpr ModInt operator-(const ModInt x) const noexcept {
    return ModInt(*this) -= x;
  }
  constexpr ModInt operator*(const ModInt x) const noexcept {
    return ModInt(*this) *= x;
  }
  constexpr ModInt operator/(const ModInt x) const noexcept {
    return ModInt(*this) /= x;
  }
  constexpr ModInt &operator+=(const ModInt x) noexcept {
    a += x.a;
    if (a >= Modulus) {
      a -= Modulus;
    }
    return *this;
  }
  constexpr ModInt &operator-=(const ModInt x) noexcept {
    if (a < x.a) {
      a += Modulus;
    }
    a -= x.a;
    return *this;
  }
  constexpr ModInt &operator*=(const ModInt x) noexcept {
    a = a * x.a % Modulus;
    return *this;
  }
  constexpr ModInt &operator/=(ModInt x) noexcept {
    ll exp = Modulus - 2;
    while (exp) {
      if (exp % 2) {
        *this *= x;
      }
      x *= x;
      exp /= 2;
    }
    return *this;
  }
  constexpr ModInt operator-() noexcept { return ModInt(-a); }
  friend ostream &operator<<(ostream &os, const ModInt &m) {
    os << m.a;
    return os;
  }
};

using mint = ModInt<mod>;

signed main() {
  ios::sync_with_stdio(false);
  cout << fixed << setprecision(15);
  ll n, m, s;
  cin >> n >> m >> s;
  vector<vector<int>> a(n), b(n), g(n);
  REP(i, m) {
    int u, v, x, y;
    cin >> u >> v >> x >> y;
    u--;
    v--;
    g[u].push_back(v);
    a[u].push_back(x);
    b[u].push_back(y);

    g[v].push_back(u);
    a[v].push_back(x);
    b[v].push_back(y);
  }
  vector<int> c(n), d(n);
  REP(i, n) cin >> c[i] >> d[i];
  vector<vector<ll>> dp(n, vector<ll>(2501, inf));
  using Triple = tuple<ll, int, ll>;
  priority_queue<Triple, vector<Triple>, greater<>> q;
  q.emplace(0, 0, min(s, 2500ll));
  while (q.size()) {
    auto x = q.top();
    q.pop();
    ll time = get<0>(x);
    int node = get<1>(x);
    ll gold = get<2>(x);
    if (dp[node][gold] <= time)
      continue;
    dp[node][gold] = time;
    q.emplace(time + d[node], node, min(gold + c[node], 2500ll));
    REP(i, sz(g[node])) {
      if (gold < a[node][i])
        continue;
      q.emplace(time + b[node][i], g[node][i], gold - a[node][i]);
    }
  }
  REP_FROM(i, 1, n) { cout << *min_element(all(dp[i])) << endl; }
  return 0;
}