This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include <bits/stdc++.h>
//#include <ext/pb_ds/assoc_container.hpp>
//#include <ext/pb_ds/tree_policy.hpp>
//#pragma GCC optimize("unroll-loops")
//#pragma GCC optimize("-O3")
//#pragma GCC optimize("Ofast")
//#pragma GCC optimize("fast-math")
//#pragma GCC optimize("no-stack-protector")
#define F first
#define S second
#define sz(x) int(x.size())
#define pb push_back
#define N 101001
#define M ll(1e9 + 7)
#define inf 1e9 + 1e9
using namespace std;
//using namespace __gnu_pbds;
typedef long double ld;
typedef long long ll;
typedef short int si;
typedef array <int, 3> a3;
//typedef tree <int, null_type, less_equal<int>, rb_tree_tag, tree_order_statistics_node_update> ordered_set;
const ld E = 1e-9;
vector <ll> g[2005];
ll A[2005], B[2005], T[2005], C[2005];
ll opr(ll i, ll b) {return (A[i] == b ? B[i] : A[i]);}
set <pair <ld, ll> > se;
ll ans[2005];
pair <ld, array <ll, 2> > dst[2005];
int main()
{
//freopen("input.txt", "r", stdin); //freopen("output4.txt", "w", stdout);
ios_base::sync_with_stdio(0); istream::sync_with_stdio(0); cin.tie(0); cout.tie(0);
int n, m;
cin >> n >> m;
for (int i = 0; i < m; i++)
{
cin >> A[i] >> B[i] >> T[i] >> C[i];
A[i]--; B[i]--;
g[A[i]].pb(i);
g[B[i]].pb(i);
}
ld ca = 1, cb = 0;
while (1)
{
for (int i = 1; i < n; i++) dst[i] = {1e18, {ll(1e9), ll(1e9)}};
se.insert({0.0, 0});
while (sz(se))
{
int v = (*se.begin()).S;
se.erase(se.begin());
for (auto i : g[v])
{
int nxt = opr(i, v);
ld val = dst[v].F + ca * ld(T[i]) + cb * ld(C[i]);
if (val + E < dst[nxt].F)
{
se.erase({dst[nxt].F, nxt});
dst[nxt] = {val, {dst[v].S[0] + T[i], dst[v].S[1] + C[i]}};
se.insert({dst[nxt].F, nxt});
}
}
}
ld nw = 0;
for (int v = 0; v < n; v++)
{
if (dst[v].F == 1e18) continue;
if (ans[v] == 0 || ans[v] > dst[v].S[0] * dst[v].S[1]) ans[v] = dst[v].S[0] * dst[v].S[1];
for (auto i : g[v])
{
int nxt = opr(i, v);
ld difa = (dst[v].S[0] + T[i]) - dst[nxt].S[0], difb = (dst[v].S[1] + C[i]) - dst[nxt].S[1];
if (difa == difb) continue;
ld Z = difb / (difb - difa);
if (Z + E < ca && Z - E> nw) nw = Z;
}
}
if (nw + E >= ca) break;
ca = nw;
cb = 1.0 - ca;
}
for (int i = 1; i < n; i++) if (ans[i] == 0) cout << -1 << '\n'; else cout << ans[i] << '\n';
}
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |