#include "walk.h"
#include<bits/stdc++.h>
using namespace std;
using ll = long long;
#define pll pair <ll,ll>
#define fi first
#define se second
#define MP make_pair
#define sz(a) (ll((a).size()))
#define BIT(mask,i) (((mask) >> (i))&1)
#define MASK(i) (1LL << (i))
struct walk{
ll l,r,y;
};
vector <walk> all;
const ll MAXN = 1e5 + 100;
const ll INF = 1e18;
vector <pll> g[MAXN*20];
ll dist[MAXN*20];
ll dijkstra(ll s,ll t){
memset(dist,0x3f,sizeof dist);
dist[s] = 0;
priority_queue <pll,vector <pll> ,greater <> > q;
q.push(MP(dist[s],s));
while (!q.empty()){
auto u = q.top().se;
ll val = q.top().fi;
q.pop();
if (dist[u] != val)continue;
for (auto tmp:g[u]){
ll v = tmp.fi,w = tmp.se;
if (dist[u] + w < dist[v]){
dist[v] = dist[u] + w;
q.push(MP(dist[v],v));
}
}
}
if (dist[t] >=INF)dist[t] = -1;
return dist[t];
}
ll SUSSYBAKA;
struct point{
ll x,y,id;
point(ll x1=-1,ll y1=-1):x(x1),y(y1),id(x1==-1?-1:SUSSYBAKA++){
}
bool operator < (const point &p)const {
return MP(x,y)<MP(p.x,p.y);
}
bool operator == (const point &p)const {
return MP(x,y)==MP(p.x,p.y);
}
};
vector <point> vertices;
long long min_distance(std::vector<int> X, std::vector<int> h, std::vector<int> l, std::vector<int> r, std::vector<int> y, int S, int G){
ll m = sz(l);
ll n = sz(X);
for (ll i = 0;i < m;i ++){
all.push_back({l[i],r[i],y[i]});
}
sort(all.begin(),all.end(),[](walk a1,walk a2){return a1.y<a2.y;});
vector <ll> order;
order.resize(n);
iota(order.begin(),order.end(),0);
sort(order.begin(),order.end(),[&](ll x,ll y){return h[x] > h[y];});
vertices.emplace_back(X[S],0);
vertices.emplace_back(X[G],0);
for (auto x:all){
vertices.emplace_back(X[x.l],x.y);
vertices.emplace_back(X[x.r],x.y);
}
auto add_vertex = [&](bool inv){
vector <ll> bd;
for (ll i = inv?n-1:0;inv?i>=0:i<n;inv?i--:i++){
// cout<<i<<endl;
if (sz(bd) && h[bd.back()] < h[i])bd.pop_back();
bd.push_back(i);
if (i==S||i==G){
ll ptr = sz(bd)-1;
for (auto x:all){
while (ptr>=0&&h[bd[ptr]] < x.y)ptr--;
if (x.l <= ptr && ptr <= x.r){
vertices.emplace_back(X[bd[ptr]],x.y);
}
}
}
}
};
add_vertex(0);
add_vertex(1);
// return -1;
vector <pair <pll,bool> > event;
for (auto x:all){
event.emplace_back(MP(X[x.l],x.y),1);
event.emplace_back(MP(X[x.r]+1,x.y),0);
}
sort(event.begin(),event.end());
sort(vertices.begin(),vertices.end());
vertices.resize(unique(vertices.begin(),vertices.end())-vertices.begin());
{
ll ptr = 0;
multiset <ll> ms;
vector <pll> Tm;
for (auto x:vertices){
while (ptr<sz(event) && event[ptr].fi.fi <= x.x){
if (event[ptr].se)ms.insert(event[ptr].fi.se);
else ms.erase(ms.find((event[ptr].fi.se)));
ptr++;
}
auto tmp = ms.upper_bound(x.y);
if (tmp != ms.end())Tm.emplace_back(x.x,*(tmp));
tmp = ms.lower_bound(x.y);
if (tmp != ms.begin())Tm.emplace_back(x.x,*prev(tmp));
}
for (auto x:Tm)vertices.emplace_back(x.fi,x.se);
}
sort(vertices.begin(),vertices.end());
vertices.resize(unique(vertices.begin(),vertices.end())-vertices.begin());
auto add_edge = [&](point a1,point a2){
ll w = abs(a1.x-a2.x) + abs(a1.y-a2.y);
g[a1.id].push_back(MP(a2.id,w));
g[a2.id].push_back(MP(a1.id,w));
};
for (ll j = 0;j + 1 < sz(vertices);j ++){
if (vertices[j].x == vertices[j+1].x)add_edge(vertices[j],vertices[j+1]);
}
auto cmp_y = [](point a1,point a2){return MP(a1.y,a1.x) < MP(a2.y,a2.x);};
sort(vertices.begin(),vertices.end(),cmp_y);
auto fi = [&](pll x){
return lower_bound(vertices.begin(),vertices.end(),x,
[](point a1,pll val){return MP(a1.y,a1.x) < MP(val.se,val.fi);})-vertices.begin();
};
vector <ll> cnt(sz(vertices));
for (auto x:all){
cnt[fi(MP(X[x.l],x.y))] ++;
cnt[fi(MP(X[x.r],x.y))] --;
}
for (ll i = 0;i < sz(vertices);i ++){
if (i)cnt[i]+=cnt[i-1];
if (cnt[i]){
add_edge(vertices[i],vertices[i+1]);
}
}
for (ll i = 0;i < sz(vertices);i ++){
if (MP(vertices[i].x,vertices[i].y) == MP(ll(X[S]),0LL))S = vertices[i].id;
if (MP(vertices[i].x,vertices[i].y) == MP(ll(X[G]),0LL))G = vertices[i].id;
}
return dijkstra(S,G);
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
10 ms |
63068 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
11 ms |
63068 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
78 ms |
76244 KB |
Output is correct |
2 |
Correct |
450 ms |
138384 KB |
Output is correct |
3 |
Correct |
525 ms |
141596 KB |
Output is correct |
4 |
Correct |
588 ms |
146016 KB |
Output is correct |
5 |
Correct |
660 ms |
149040 KB |
Output is correct |
6 |
Correct |
612 ms |
145640 KB |
Output is correct |
7 |
Correct |
266 ms |
106520 KB |
Output is correct |
8 |
Correct |
232 ms |
97544 KB |
Output is correct |
9 |
Correct |
561 ms |
145388 KB |
Output is correct |
10 |
Correct |
270 ms |
116644 KB |
Output is correct |
11 |
Correct |
18 ms |
64328 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
78 ms |
76244 KB |
Output is correct |
2 |
Correct |
450 ms |
138384 KB |
Output is correct |
3 |
Correct |
525 ms |
141596 KB |
Output is correct |
4 |
Correct |
588 ms |
146016 KB |
Output is correct |
5 |
Correct |
660 ms |
149040 KB |
Output is correct |
6 |
Correct |
612 ms |
145640 KB |
Output is correct |
7 |
Correct |
266 ms |
106520 KB |
Output is correct |
8 |
Correct |
232 ms |
97544 KB |
Output is correct |
9 |
Correct |
561 ms |
145388 KB |
Output is correct |
10 |
Correct |
270 ms |
116644 KB |
Output is correct |
11 |
Correct |
18 ms |
64328 KB |
Output is correct |
12 |
Correct |
531 ms |
143760 KB |
Output is correct |
13 |
Correct |
409 ms |
146396 KB |
Output is correct |
14 |
Correct |
693 ms |
174252 KB |
Output is correct |
15 |
Correct |
532 ms |
129120 KB |
Output is correct |
16 |
Correct |
455 ms |
136944 KB |
Output is correct |
17 |
Correct |
527 ms |
149176 KB |
Output is correct |
18 |
Correct |
455 ms |
130472 KB |
Output is correct |
19 |
Correct |
472 ms |
137644 KB |
Output is correct |
20 |
Correct |
285 ms |
107020 KB |
Output is correct |
21 |
Correct |
36 ms |
66200 KB |
Output is correct |
22 |
Correct |
425 ms |
134936 KB |
Output is correct |
23 |
Correct |
419 ms |
131496 KB |
Output is correct |
24 |
Correct |
335 ms |
121344 KB |
Output is correct |
25 |
Correct |
382 ms |
126596 KB |
Output is correct |
26 |
Correct |
309 ms |
123980 KB |
Output is correct |
27 |
Correct |
731 ms |
167300 KB |
Output is correct |
28 |
Correct |
426 ms |
148292 KB |
Output is correct |
29 |
Correct |
672 ms |
156884 KB |
Output is correct |
30 |
Correct |
291 ms |
110576 KB |
Output is correct |
31 |
Correct |
655 ms |
155320 KB |
Output is correct |
32 |
Correct |
237 ms |
111168 KB |
Output is correct |
33 |
Correct |
268 ms |
119240 KB |
Output is correct |
34 |
Correct |
288 ms |
120000 KB |
Output is correct |
35 |
Correct |
332 ms |
120320 KB |
Output is correct |
36 |
Correct |
271 ms |
110480 KB |
Output is correct |
37 |
Correct |
268 ms |
105676 KB |
Output is correct |
38 |
Correct |
273 ms |
117692 KB |
Output is correct |
39 |
Correct |
338 ms |
124468 KB |
Output is correct |
40 |
Correct |
287 ms |
115500 KB |
Output is correct |
41 |
Correct |
239 ms |
109728 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
10 ms |
63068 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |