#include "plants.h"
#include <bits/stdc++.h>
#define sz(x) (int)(x.size())
#define debug printf
#define lp(i,a,b) for(int i = a ; i < b; i++)
#define pb push_back
#define ff first
#define ss second
#define mk make_pair
#define pii pair<int,int>
#define ll long long
#define all(x) x.begin(),x.end()
const int MAX = 2e5+10 ;
const int LOG = 20 ;
const int inf = 1e8+10 ;
using namespace std ;
struct Seg
{
ll lz[MAX*4] , mn[MAX*4];
int idx[MAX*4] ;
void clean()
{
memset(lz, 0, sizeof lz) ;
memset(mn, 0, sizeof mn) ;
}
int m(int l, int r) { return (l+r)>>1 ; }
void refresh(int pos, int l, int r)
{
if(!lz[pos]) return ;
mn[pos] += lz[pos] ;
if(l == r) return (void)(lz[pos] = 0 ) ;
lz[pos<<1] += lz[pos] ;
lz[pos<<1|1] += lz[pos] ;
lz[pos] = 0 ;
}
void upd(int pos, int l, int r, int beg, int en , int val )
{
refresh(pos,l,r) ;
if(l > en || r < beg ) return ;
if(l >= beg && r <= en)
{
lz[pos] += (ll)val ;
refresh(pos,l,r) ;
if(l == r) idx[pos] = l ;
return ;
}
upd(pos<<1 , l, m(l,r) , beg, en, val ) ;
upd(pos<<1|1 , m(l,r)+1,r, beg, en, val ) ;
if(mn[pos<<1] <= mn[pos<<1|1])
{
idx[pos] = idx[pos<<1] ;
mn[pos] = mn[pos<<1] ;
}
else
{
idx[pos] = idx[pos<<1|1] ;
mn[pos] = mn[pos<<1|1] ;
}
}
pii qry(int pos, int l, int r, int beg, int en )
{
refresh(pos,l,r) ;
if( l > en || r < beg ) return mk(inf, -1) ;
if(l >= beg && r <= en) return mk( mn[pos] , idx[pos] ) ;
pii al = qry(pos<<1 , l, m(l,r) , beg, en ) ;
pii ar = qry(pos<<1|1 , m(l,r)+1,r, beg, en ) ;
if( al.ff <= ar.ff ) return al ;
return ar ;
}
} seg ;
int n , last , k , mx ;
int dp[2][LOG][MAX] ;
vector<int> h , lef , rig ;
vector<pii> intervalLeft[MAX] , intervalRight[MAX] ;
void extract(int x)
{
while(true)
{
pii p = mk(inf, -1) ;
for(auto e : intervalLeft[x] )
p = min(p, seg.qry(1,0,n-1, e.ff, e.ss) ) ;
if(p.ff == 0)
{
extract(p.ss) ;
continue ;
}
h[x] = last-- ;
seg.upd(1,0,n-1, x , x , inf ) ;
for(auto e : intervalLeft[x] ) seg.upd(1,0,n-1, e.ff, e.ss, -1) ;
return ;
}
}
void init(int K, vector<int> r)
{
k = K ;
n = sz(r) ;
h.resize(n,0) ;
last = n ;
//Find intervals
intervalLeft[0].push_back( mk(n-k+1,n-1) ) ;
lp(i,1,n)
{
intervalLeft[i].push_back( mk(max(0,i-k+1) , i-1) ) ;
if( i >= k-1 ) continue ;
int falta = k-1 - i ;
intervalLeft[i].pb( mk( n-falta , n-1 ) ) ;
}
intervalRight[n-1].pb( mk( 0 , k-2 ) ) ;
lp(i,0,n-1)
{
intervalRight[i].pb( mk( i+1, min(n-1, i+k-1) ) ) ;
if( i+k-1 <= n-1 ) continue ;
int lef = k-1 - (n-i-1) ;
intervalRight[i].pb( mk(0,lef-1)) ;
}
for(int i = 0 ; i < n ; i++ )
seg.upd(1,0,n-1, i , i , r[i] );
while(true)
{
pii p = seg.qry(1,0,n-1,0,n-1) ;
if(p.ff != 0) break ;
extract(p.ss) ;
}
if(2*k > n ) return ;
seg.clean() ;
vector<int> idx(n) ;
iota(all(idx) , 0 ) ;
sort(all(idx), [&](int i, int j){ return h[i] < h[j]; } ) ;
memset(dp, -1, sizeof dp ) ;
for(auto e : idx )
{
pii p = mk(0 , -1 ) ;
for(auto e : intervalLeft[e] )
p = min(p , seg.qry(1,0,n-1, e.ff, e.ss) ) ;
dp[0][0][e] = p.ss ;
p = mk(0, -1 ) ;
for(auto e : intervalRight[e] )
p = min(p, seg.qry(1,0,n-1, e.ff, e.ss) ) ;
dp[1][0][e] = p.ss ;
seg.upd(1,0,n-1,e,e, -h[e] ) ;
}
lp(i,0,2)
lp(j,1,LOG)
lp(g,0,n)
{
if(dp[i][j-1][g] == -1) continue ;
dp[i][j][g] = dp[i][j-1][ dp[i][j-1][g] ] ;
}
//for(auto e : h ) cout << e << " " ;
// cout << endl ;
}
bool isOn(int m, int i) { return ((1<<i)&m) != 0 ; }
int getDist( bool clockwise , int a, int b )
{
if(!clockwise)
{
if(b < a) return a - b ;
return a + n - b ;
}
if( b > a ) return b-a ;
return n-a+b ;
}
int compare_plants(int x, int y)
{
if( 2*k > n || getDist(false, x, y) < k || getDist(true, x, y) < k)
return h[x] > h[y] ? 1 : -1 ;
for(int cur : {x,y} )
for(int type = 0 ; type < 2 ; type++ )
{
int filho = cur , _filho ;
int other = (cur == x ) ? y : x ;
if(h[cur] < h[other]) continue ;
for(int j = LOG - 1 ; j >= 0 ; j-- )
{
filho = cur ;
for(int i = j ; i >= 0 ; i-- )
{
_filho = dp[type][i][filho] ;
if(_filho == -1) continue ;
int distToCur = getDist(type, _filho, cur ) ;
int distToOther = getDist(type, _filho, other ) ;
if(distToCur < distToOther) continue ;
if(distToOther < k )
{
if(h[other] <= h[_filho] )
return (cur == x) ? 1 : -1 ;
}
else if(h[_filho] > h[other] )filho = _filho ;
}
if(h[other] <= h[filho] && getDist(type, filho, other) < k )
return (cur == x) ? 1 : -1 ;
}
}
return 0 ;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
7 ms |
9728 KB |
Output is correct |
2 |
Correct |
7 ms |
9728 KB |
Output is correct |
3 |
Correct |
30 ms |
53496 KB |
Output is correct |
4 |
Correct |
30 ms |
53504 KB |
Output is correct |
5 |
Correct |
32 ms |
53504 KB |
Output is correct |
6 |
Correct |
195 ms |
56440 KB |
Output is correct |
7 |
Correct |
369 ms |
58232 KB |
Output is correct |
8 |
Correct |
828 ms |
73592 KB |
Output is correct |
9 |
Correct |
920 ms |
73592 KB |
Output is correct |
10 |
Correct |
1097 ms |
73720 KB |
Output is correct |
11 |
Correct |
1367 ms |
73720 KB |
Output is correct |
12 |
Correct |
1622 ms |
73592 KB |
Output is correct |
13 |
Correct |
956 ms |
73720 KB |
Output is correct |
14 |
Correct |
774 ms |
73720 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
8 ms |
9728 KB |
Output is correct |
2 |
Correct |
6 ms |
9680 KB |
Output is correct |
3 |
Correct |
6 ms |
9728 KB |
Output is correct |
4 |
Correct |
7 ms |
9728 KB |
Output is correct |
5 |
Correct |
7 ms |
9728 KB |
Output is correct |
6 |
Correct |
11 ms |
9856 KB |
Output is correct |
7 |
Correct |
90 ms |
13432 KB |
Output is correct |
8 |
Correct |
8 ms |
9856 KB |
Output is correct |
9 |
Correct |
10 ms |
9856 KB |
Output is correct |
10 |
Correct |
88 ms |
13304 KB |
Output is correct |
11 |
Correct |
90 ms |
13176 KB |
Output is correct |
12 |
Correct |
91 ms |
13432 KB |
Output is correct |
13 |
Correct |
87 ms |
13276 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
8 ms |
9728 KB |
Output is correct |
2 |
Correct |
6 ms |
9680 KB |
Output is correct |
3 |
Correct |
6 ms |
9728 KB |
Output is correct |
4 |
Correct |
7 ms |
9728 KB |
Output is correct |
5 |
Correct |
7 ms |
9728 KB |
Output is correct |
6 |
Correct |
11 ms |
9856 KB |
Output is correct |
7 |
Correct |
90 ms |
13432 KB |
Output is correct |
8 |
Correct |
8 ms |
9856 KB |
Output is correct |
9 |
Correct |
10 ms |
9856 KB |
Output is correct |
10 |
Correct |
88 ms |
13304 KB |
Output is correct |
11 |
Correct |
90 ms |
13176 KB |
Output is correct |
12 |
Correct |
91 ms |
13432 KB |
Output is correct |
13 |
Correct |
87 ms |
13276 KB |
Output is correct |
14 |
Correct |
135 ms |
15364 KB |
Output is correct |
15 |
Correct |
813 ms |
37368 KB |
Output is correct |
16 |
Correct |
136 ms |
15352 KB |
Output is correct |
17 |
Correct |
785 ms |
37288 KB |
Output is correct |
18 |
Correct |
450 ms |
37368 KB |
Output is correct |
19 |
Correct |
456 ms |
37240 KB |
Output is correct |
20 |
Correct |
734 ms |
37288 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
30 ms |
53504 KB |
Output is correct |
2 |
Correct |
30 ms |
53496 KB |
Output is correct |
3 |
Correct |
268 ms |
56568 KB |
Output is correct |
4 |
Correct |
1251 ms |
77432 KB |
Output is correct |
5 |
Correct |
1506 ms |
73976 KB |
Output is correct |
6 |
Incorrect |
1618 ms |
73720 KB |
Output isn't correct |
7 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
7 ms |
9728 KB |
Output is correct |
2 |
Correct |
6 ms |
9728 KB |
Output is correct |
3 |
Correct |
31 ms |
53496 KB |
Output is correct |
4 |
Correct |
31 ms |
53504 KB |
Output is correct |
5 |
Correct |
30 ms |
53496 KB |
Output is correct |
6 |
Correct |
38 ms |
53632 KB |
Output is correct |
7 |
Correct |
106 ms |
54136 KB |
Output is correct |
8 |
Correct |
50 ms |
54264 KB |
Output is correct |
9 |
Incorrect |
98 ms |
54168 KB |
Output isn't correct |
10 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
6 ms |
9728 KB |
Output is correct |
2 |
Correct |
6 ms |
9728 KB |
Output is correct |
3 |
Correct |
30 ms |
53504 KB |
Output is correct |
4 |
Correct |
30 ms |
53496 KB |
Output is correct |
5 |
Correct |
33 ms |
53632 KB |
Output is correct |
6 |
Correct |
1080 ms |
73808 KB |
Output is correct |
7 |
Correct |
1396 ms |
73720 KB |
Output is correct |
8 |
Correct |
1329 ms |
73848 KB |
Output is correct |
9 |
Correct |
1375 ms |
73736 KB |
Output is correct |
10 |
Correct |
851 ms |
74104 KB |
Output is correct |
11 |
Correct |
988 ms |
73720 KB |
Output is correct |
12 |
Correct |
1094 ms |
78460 KB |
Output is correct |
13 |
Correct |
1283 ms |
74104 KB |
Output is correct |
14 |
Correct |
1431 ms |
73720 KB |
Output is correct |
15 |
Correct |
1402 ms |
73760 KB |
Output is correct |
16 |
Correct |
827 ms |
73848 KB |
Output is correct |
17 |
Correct |
1060 ms |
73700 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
7 ms |
9728 KB |
Output is correct |
2 |
Correct |
7 ms |
9728 KB |
Output is correct |
3 |
Correct |
30 ms |
53496 KB |
Output is correct |
4 |
Correct |
30 ms |
53504 KB |
Output is correct |
5 |
Correct |
32 ms |
53504 KB |
Output is correct |
6 |
Correct |
195 ms |
56440 KB |
Output is correct |
7 |
Correct |
369 ms |
58232 KB |
Output is correct |
8 |
Correct |
828 ms |
73592 KB |
Output is correct |
9 |
Correct |
920 ms |
73592 KB |
Output is correct |
10 |
Correct |
1097 ms |
73720 KB |
Output is correct |
11 |
Correct |
1367 ms |
73720 KB |
Output is correct |
12 |
Correct |
1622 ms |
73592 KB |
Output is correct |
13 |
Correct |
956 ms |
73720 KB |
Output is correct |
14 |
Correct |
774 ms |
73720 KB |
Output is correct |
15 |
Correct |
8 ms |
9728 KB |
Output is correct |
16 |
Correct |
6 ms |
9680 KB |
Output is correct |
17 |
Correct |
6 ms |
9728 KB |
Output is correct |
18 |
Correct |
7 ms |
9728 KB |
Output is correct |
19 |
Correct |
7 ms |
9728 KB |
Output is correct |
20 |
Correct |
11 ms |
9856 KB |
Output is correct |
21 |
Correct |
90 ms |
13432 KB |
Output is correct |
22 |
Correct |
8 ms |
9856 KB |
Output is correct |
23 |
Correct |
10 ms |
9856 KB |
Output is correct |
24 |
Correct |
88 ms |
13304 KB |
Output is correct |
25 |
Correct |
90 ms |
13176 KB |
Output is correct |
26 |
Correct |
91 ms |
13432 KB |
Output is correct |
27 |
Correct |
87 ms |
13276 KB |
Output is correct |
28 |
Correct |
135 ms |
15364 KB |
Output is correct |
29 |
Correct |
813 ms |
37368 KB |
Output is correct |
30 |
Correct |
136 ms |
15352 KB |
Output is correct |
31 |
Correct |
785 ms |
37288 KB |
Output is correct |
32 |
Correct |
450 ms |
37368 KB |
Output is correct |
33 |
Correct |
456 ms |
37240 KB |
Output is correct |
34 |
Correct |
734 ms |
37288 KB |
Output is correct |
35 |
Correct |
30 ms |
53504 KB |
Output is correct |
36 |
Correct |
30 ms |
53496 KB |
Output is correct |
37 |
Correct |
268 ms |
56568 KB |
Output is correct |
38 |
Correct |
1251 ms |
77432 KB |
Output is correct |
39 |
Correct |
1506 ms |
73976 KB |
Output is correct |
40 |
Incorrect |
1618 ms |
73720 KB |
Output isn't correct |
41 |
Halted |
0 ms |
0 KB |
- |