| 16087 | 2025-03-30 20:29:55 | szil | Autópálya infláció | cpp17 | Accepted 100/100 | 1.914s | 5748 KiB |
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const ll MOD = 1e9+7;
const int MAXN = 3001;
const int MAXM = 6001;
const int BLOCK_CNT = 50;
const int BLOCK_SIZE = 60;
const int TH = 30;
const ll MULTI = (1ll<<(1+BLOCK_SIZE))%MOD;
struct bigint {
vector<uint64_t> data;
bigint() {
data.assign(BLOCK_CNT, ~0ull);
}
void clear() {
data.assign(BLOCK_CNT, 0ull);
}
bool operator<(const bigint &rhs) {
for (int i = 0; i < BLOCK_CNT; i++) {
if (data[i] != rhs.data[i]) {
return data[i] < rhs.data[i];
}
}
return 0;
}
bigint operator+(const bigint &rhs) {
bigint res;
res.clear();
bool carry = 0;
for (int i = BLOCK_CNT-1; i >= 0; i--) {
res.data[i] = data[i] + rhs.data[i] + carry;
carry = (res.data[i] >> (BLOCK_SIZE + 1)) & 1;
res.data[i] &= ~(1ull << (BLOCK_SIZE + 1));
}
return res;
}
void inflate() {
bool carry = 0;
for (int i = BLOCK_CNT-1; i >= 0; i--) {
bool next_carry = (data[i] >> BLOCK_SIZE) & 1;
data[i] &= ~(1ull << BLOCK_SIZE);
data[i] <<= 1;
data[i] |= carry;
carry = next_carry;
}
}
ll calc_ans() {
ll ans = 0;
ll mul = 1;
for (int i = BLOCK_CNT-1; i >= 0; i--) {
data[i] %= MOD;
ans += data[i] * mul;
ans %= MOD;
mul *= MULTI;
mul %= MOD;
}
return ans;
}
};
int u[MAXM], v[MAXM];
void solve() {
int n, m; cin >> n >> m;
vector<bigint> w(m+1);
for (int i = 1; i <= m; i++) {
w[i].clear();
cin >> u[i] >> v[i] >> w[i].data.back();
}
vector<bigint> dist_last(n+1);
vector<int> first_seen(n+1, n);
vector<bool> vis_last(n+1);
dist_last[1].clear();
vis_last[1] = 1;
int is_ok = -1;
first_seen[1] = 0;
for (int i = 1; i <= n; i++) {
vector<bigint> dist_curr = dist_last;
vector<bool> vis_curr = vis_last;
for (int j = 1; j <= m; j++) {
if (first_seen[u[j]] + TH >= i && first_seen[v[j]] + TH >= i && vis_last[u[j]] && dist_last[u[j]] + w[j] < dist_curr[v[j]]) {
dist_curr[v[j]] = dist_last[u[j]] + w[j];
vis_curr[v[j]] = 1;
first_seen[v[j]] = min(first_seen[v[j]], i);
}
}
for (int j = 1; j <= m; j++) {
if (first_seen[u[j]] + TH >= i && first_seen[v[j]] + TH >= i && vis_last[v[j]] && dist_last[v[j]] + w[j] < dist_curr[u[j]]) {
dist_curr[u[j]] = dist_last[v[j]] + w[j];
vis_curr[u[j]] = 1;
first_seen[u[j]] = min(first_seen[u[j]], i);
}
}
dist_last = dist_curr;
vis_last = vis_curr;
for (int j = 1; j <= m; j++) {
if (first_seen[u[j]] + TH >= i && first_seen[v[j]] + TH >= i) w[j].inflate();
}
}
for (int i = 2; i <= n; i++) {
cout << dist_last[i].calc_ans() << "\n";
}
cout << "\n";
}
int main() {
ios::sync_with_stdio(0); cin.tie(0);
int t = 1;
// cin >> t;
while (t--) solve();
return 0;
}| Subtask | Sum | Test | Verdict | Time | Memory | ||
|---|---|---|---|---|---|---|---|
| subtask1 | 0/0 | ||||||
| 1 | Accepted | 1ms | 316 KiB | ||||
| 2 | Accepted | 1.437s | 5660 KiB | ||||
| subtask2 | 8/8 | ||||||
| 3 | Accepted | 1.269s | 4440 KiB | ||||
| 4 | Accepted | 1.299s | 4336 KiB | ||||
| 5 | Accepted | 1.307s | 4148 KiB | ||||
| 6 | Accepted | 1.299s | 4340 KiB | ||||
| 7 | Accepted | 1.633s | 4336 KiB | ||||
| subtask3 | 15/15 | ||||||
| 8 | Accepted | 1ms | 316 KiB | ||||
| 9 | Accepted | 1ms | 316 KiB | ||||
| 10 | Accepted | 1ms | 316 KiB | ||||
| 11 | Accepted | 1ms | 316 KiB | ||||
| 12 | Accepted | 1ms | 508 KiB | ||||
| 13 | Accepted | 1ms | 316 KiB | ||||
| 14 | Accepted | 1ms | 316 KiB | ||||
| 15 | Accepted | 1ms | 316 KiB | ||||
| 16 | Accepted | 1ms | 316 KiB | ||||
| 17 | Accepted | 1ms | 316 KiB | ||||
| 18 | Accepted | 1ms | 316 KiB | ||||
| subtask4 | 34/34 | ||||||
| 19 | Accepted | 1.411s | 5660 KiB | ||||
| 20 | Accepted | 1.472s | 5660 KiB | ||||
| 21 | Accepted | 1.478s | 5672 KiB | ||||
| 22 | Accepted | 1.427s | 5660 KiB | ||||
| 23 | Accepted | 1.48s | 5528 KiB | ||||
| 24 | Accepted | 1.544s | 5528 KiB | ||||
| 25 | Accepted | 1.402s | 5432 KiB | ||||
| 26 | Accepted | 1.424s | 5660 KiB | ||||
| 27 | Accepted | 1.774s | 5652 KiB | ||||
| 28 | Accepted | 1.347s | 4204 KiB | ||||
| 29 | Accepted | 1.35s | 4332 KiB | ||||
| 30 | Accepted | 1.348s | 3892 KiB | ||||
| 31 | Accepted | 1.322s | 4148 KiB | ||||
| 32 | Accepted | 1.623s | 5644 KiB | ||||
| 33 | Accepted | 1.646s | 5428 KiB | ||||
| subtask5 | 21/21 | ||||||
| 34 | Accepted | 61ms | 1076 KiB | ||||
| 35 | Accepted | 71ms | 1588 KiB | ||||
| 36 | Accepted | 70ms | 1772 KiB | ||||
| 37 | Accepted | 70ms | 1588 KiB | ||||
| 38 | Accepted | 68ms | 1740 KiB | ||||
| 39 | Accepted | 68ms | 1736 KiB | ||||
| 40 | Accepted | 78ms | 1588 KiB | ||||
| 41 | Accepted | 79ms | 1732 KiB | ||||
| 42 | Accepted | 37ms | 1076 KiB | ||||
| 43 | Accepted | 43ms | 1076 KiB | ||||
| 44 | Accepted | 43ms | 1076 KiB | ||||
| 45 | Accepted | 43ms | 1080 KiB | ||||
| 46 | Accepted | 39ms | 1076 KiB | ||||
| 47 | Accepted | 37ms | 820 KiB | ||||
| 48 | Accepted | 50ms | 1076 KiB | ||||
| 49 | Accepted | 50ms | 1276 KiB | ||||
| subtask6 | 22/22 | ||||||
| 50 | Accepted | 1.61s | 5532 KiB | ||||
| 51 | Accepted | 1.493s | 5748 KiB | ||||
| 52 | Accepted | 1.475s | 5688 KiB | ||||
| 53 | Accepted | 1.416s | 5664 KiB | ||||
| 54 | Accepted | 1.411s | 5684 KiB | ||||
| 55 | Accepted | 1.442s | 5684 KiB | ||||
| 56 | Accepted | 1.82s | 5684 KiB | ||||
| 57 | Accepted | 1.776s | 5656 KiB | ||||
| 58 | Accepted | 158ms | 1588 KiB | ||||
| 59 | Accepted | 1.365s | 4332 KiB | ||||
| 60 | Accepted | 1.376s | 4148 KiB | ||||
| 61 | Accepted | 1.373s | 4336 KiB | ||||
| 62 | Accepted | 1.378s | 4148 KiB | ||||
| 63 | Accepted | 1.33s | 4148 KiB | ||||
| 64 | Accepted | 1.605s | 5428 KiB | ||||
| 65 | Accepted | 1.618s | 5428 KiB | ||||
| 66 | Accepted | 1.914s | 5660 KiB | ||||