#include<bits/stdc++.h>
#define fast ios_base::sync_with_stdio(0);cin.tie(NULL);cout.tie(NULL)
#define pb push_back
#define F first
#define S second
#define int long long
using namespace std;
const int N = 1e6 + 5, MOD = 1e9 + 7;
vector<pair<int, int> > g[N];
int n, par[N], dis[N], dp[N], s[N], d[N];
bool vis[N];
struct Line{
int m, c;
int val(int x){
return m*x + c;
}
};
deque<Line> hulls;
bool check(Line a, Line b, Line cc){
int x = (b.c - a.c);
x = x * (a.m - cc.m);
int y = (cc.c - a.c);
y = y * (a.m - b.m);
if(x > y)
return true;
return false;
}
void insert(int mm, int cc){
Line cur;
cur.m = mm;
cur.c = cc;
if(!hulls.empty() and hulls[0].m == mm){
if(hulls[0].c > cc)
hulls.pop_front();
else
return;
}
while(hulls.size() > 1){
if(check(hulls[1], hulls[0], cur))
hulls.pop_front();
else
break;
}
hulls.push_front(cur);
}
int get(int x){
if(hulls.empty())
return 1e15;
int low = 0, high = (int) hulls.size() - 2, ans = hulls.size() - 1;
while(low <= high){
int mid = (low + high) / 2;
if(hulls[mid].val(x) <= hulls[mid + 1].val(x)){
ans = mid;
high = mid - 1;
}
else
low = mid + 1;
}
return hulls[ans].val(x);
}
void dfs(int node, int p, int sum){
dis[node] = sum;
par[node] = p;
for(auto k : g[node]){
if(k.F == p) continue;
dfs(k.F, node, sum + k.S);
}
}
void bfs(int node){
queue<int> q;
q.push(node);
vis[node] = true;
while(!q.empty()){
int node = q.front();
q.pop();
if(node != 1 and n <= 2500){
// calculate dp[node]
int cur = par[node], mn = 1e15;
while(cur != 0){
mn = min(mn, dp[cur] - dis[cur] * d[node]);
cur = par[cur];
}
dp[node] = mn + s[node] + dis[node] * d[node];
}
for(auto k : g[node]){
if(!vis[k.F]){
vis[k.F] = true;
q.push({k.F});
}
}
}
}
void dfs1(int node){
int mn = s[node] + dis[node] * d[node];
mn = min(mn, get(-d[node]) + s[node] + dis[node] * d[node]);
dp[node] = mn;
if(node != 1)
insert(dis[node], dp[node]);
for(auto k : g[node]){
if(k.F == par[node]) continue;
dfs1(k.F);
}
if(hulls.front().m == dis[node] and hulls.front().c == dp[node])
hulls.pop_front();
}
void solve(){
cin >> n;
for(int i = 1; i < n; ++i){
int u, v, c;
cin >> u >> v >> c;
g[u].pb({v, c});
g[v].pb({u, c});
}
for(int i = 1; i < n; ++i)
cin >> s[i + 1] >> d[i + 1];
dfs(1, 0, 0);
memset(vis, 0, sizeof vis);
if(n <= 2500)
bfs(1);
else
dfs1(1);
for(int i = 2; i <= n; ++i)
cout << dp[i] << " ";
}
int32_t main(){
fast;
int t = 1;
// cin >> t;
while(t--)
solve();
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
13 ms |
24796 KB |
Output is correct |
2 |
Correct |
27 ms |
25192 KB |
Output is correct |
3 |
Incorrect |
48 ms |
31832 KB |
Output isn't correct |
4 |
Runtime error |
67 ms |
35068 KB |
Memory limit exceeded |
5 |
Runtime error |
116 ms |
38356 KB |
Memory limit exceeded |
6 |
Runtime error |
126 ms |
41384 KB |
Memory limit exceeded |
7 |
Runtime error |
61 ms |
33556 KB |
Memory limit exceeded |
8 |
Runtime error |
105 ms |
37996 KB |
Memory limit exceeded |
9 |
Runtime error |
107 ms |
38916 KB |
Memory limit exceeded |
10 |
Runtime error |
91 ms |
38004 KB |
Memory limit exceeded |