#pragma GCC optimize("Ofast,unroll-loops") #include // #include // #include // #include // #include // #include // #include // #include // #include // #include #include // #include // #include // #include // #include // #include // #include // #include // #include // #include #include // #include // #include // #include // #include // #include // #include // #include // #include #include // #include // #include // #include // #include // #include // #include // #include using namespace std; #define BIT(i,j) (((i)>>(j))&1LL) #define MASK(i) (1LL<<(i)) #define ALL(x) (x).begin(),(x).end() #define SZ(x) (int)(x).size() #define fi first #define se second #define ull ungsigned long long #define ll long long #define ld long double #define vi vector #define vvi vector #define vvvi vector #define pii pair #define vpii vector #define vvpii vector #define REPDIS(i,be,en,j) for(int i = (be); i<=(en); i+=j) #define REPD(i,be,en) for(int i = (be); i>=(en); i--) #define REP(i,be,en) for(int i = (be); i<=(en); i++) #define endl "\n" #define MP make_pair // #define int ll //-----------------------------------------------------------------------------------------------// inline void scan(){} template inline void scan(F &f,R&... r){cin>>f;scan(r...);} inline void print(){} template inline void print(F f,R... r){cout< void __print(const pair &x) {cerr << '{'; __print(x.first); cerr << ','; __print(x.second); cerr << '}';} template void __print(const T &x) {int f = 0; cerr << '{'; for (auto &i: x) cerr << (f++ ? "," : ""), __print(i); cerr << "}";} void _print() {cerr << "]\n";} template void _print(T t, V... v) {__print(t); if (sizeof...(v)) cerr << ", "; _print(v...);} #ifndef ONLINE_JUDGE #define debug(x...) cerr << "[" << #x << "] = ["; _print(x) #else #define debug(x...) #endif //------------------------------------------------------------------------------------------------// const ll LINF = 1e18; const int INF = 1e9; const int LOG = 20; const int MAXN = 2e5+7; const int N = 1e2+3; const int MOD = 1e9+7; const int BASE = 1e5; const ld EPS = 1e-9; // const ld PI = acos(-1); const int OFFSET = 1e3; //------------------------------------------------------------------------------------------------// template bool mini(T1 &a, T2 b){if(a>b){a=b;return true;}return false;} template bool maxi(T1 &a, T2 b){if(a T gcd(T a, T b) { while(b) { a %= b; swap(a,b); } return a; } template T lcm(T a, T b) { return a/gcd(a,b)*b; } //------------------------------------------------------------------------------------------------// /* ---------------------------------------------------------------- END OF TEMPLATE ---------------------------------------------------------------- */ int BUF=100000; int buf_readed; int buf_index; char buffer[100000]; void my_init(){ buf_index=buf_readed=0; } /*****************************************/ char my_getchar(){ char x; if(buf_index==buf_readed){ buf_index=0; buf_readed=fread(buffer,1,BUF,stdin); } x=buffer[buf_index]; buf_index++; return x; } /*************************************/ int readInt(){ int r=0; char c; while(1){ c=my_getchar(); if((c>='0' && c<='9')){ r=(c-'0'); break; } } while(1){ c=my_getchar(); if(!(c>='0' && c<='9')) break; r=r*10+(c-'0'); } return r; } /***************************************/ vi adj[MAXN]; int a[MAXN],res[MAXN],del[MAXN],child[MAXN],n; multiset tam[MAXN]; int st[MAXN*4]; vpii b; void seg_update(int id, int l, int r, int pos, int val){ if(l==r){ st[id] = val; return ; } int mid = (l+r)>>1; if(pos<=mid) seg_update(id<<1,l,mid,pos,val); if(pos>mid) seg_update(id<<1|1,mid+1,r,pos,val); st[id] = max(st[id<<1],st[id<<1|1]); } int seg_get(int id, int l, int r, int lo, int hi){ if(l>hi || r=lo && r<=hi) return st[id]; int mid = (l+r)>>1; return max(seg_get(id<<1,l,mid,lo,hi),seg_get(id<<1|1,mid+1,r,lo,hi)); } void seg_update(int pos, int val){ tam[pos].insert(val); seg_update(1,1,n,pos,*tam[pos].rbegin()); } void seg_rem(int pos, int val){ tam[pos].erase(tam[pos].find(val)); int x = -INF; if(SZ(tam[pos])>0) x = *tam[pos].rbegin(); seg_update(1,1,n,pos,x); } int seg_get(int l, int r){ return seg_get(1,1,n,l,r); } int dfs1(int u, int par){ child[u] = 1; for(int &v:adj[u]){ if(del[v] || v==par) continue; child[u] += dfs1(v,u); } return child[u]; } int dfs2(int u, int par, int siz){ for(int &v:adj[u]){ if(v==par || del[v]) continue; if(child[v]*2>siz) return dfs2(v,u,siz); } return u; } void update(int u, int par, int dis){ seg_update(a[u],dis); for(int &v:adj[u]){ if(del[v] || v==par) continue; update(v,u,dis+1); } } void rem(int u, int par, int dis){ seg_rem(a[u],dis); for(int &v:adj[u]){ if(del[v] || v==par) continue; rem(v,u,dis+1); } } void rem1(int u, int par){ if(SZ(tam[a[u]])>0){ tam[a[u]].clear(); seg_update(1,1,n,a[u],-INF); } for(int &v:adj[u]){ if(del[v] || v==par) continue; rem1(v,u); } } void solve(int u, int par, int dis){ maxi(res[u],dis+seg_get(a[u],n)); for(int &v:adj[u]){ if(del[v] || v==par) continue; solve(v,u,dis+1); } } void solve1(int u, int par, int dis){ b.emplace_back(u,dis); seg_rem(a[u],dis); for(int &v:adj[u]){ if(del[v] || v==par) continue; solve1(v,u,dis+1); } } void calc(int u){ int siz = dfs1(u,0); u = dfs2(u,0,siz); del[u] = 1; update(u,0,0); maxi(res[u],seg_get(a[u],n)); for(int &v:adj[u]){ if(del[v]) continue; b.clear(); solve1(v,u,1); for(auto &&[x,dis]:b) maxi(res[x],dis+seg_get(a[x],n)); for(auto &&[x,dis]:b) seg_update(a[x],dis); } rem1(u,0); // rem(u,0,0); for(int &v:adj[u]){ if(del[v]) continue; calc(v); } } void solve(){ n = readInt(); vi compress; for(int i = 1; i<=n; i++){ a[i] = readInt(); compress.emplace_back(a[i]); } sort(ALL(compress)); compress.resize(unique(ALL(compress))-compress.begin()); for(int i = 1; i<=n; i++){ a[i] = upper_bound(ALL(compress),a[i])-compress.begin(); } for(int i = 1; iTL \>AC **/