#pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") #pragma GCC optimize("inline") #include using namespace std; template struct cLtraits_identity{ using type = T; } ; template using cLtraits_try_make_signed = typename conditional< is_integral::value, make_signed, cLtraits_identity >::type; template struct cLtraits_common_type{ using tS = typename cLtraits_try_make_signed::type; using tT = typename cLtraits_try_make_signed::type; using type = typename common_type::type; } ; void*wmem; char memarr[96000000]; template inline auto min_L(S a, T b) -> typename cLtraits_common_type::type{ return (typename cLtraits_common_type::type) a <= (typename cLtraits_common_type::type) b ? a : b; } template inline auto max_L(S a, T b) -> typename cLtraits_common_type::type{ return (typename cLtraits_common_type::type) a >= (typename cLtraits_common_type::type) b ? a : b; } template inline void walloc1d(T **arr, int x, void **mem = &wmem){ static int skip[16] = {0, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1}; (*mem) = (void*)( ((char*)(*mem)) + skip[((unsigned long long)(*mem)) & 15] ); (*arr)=(T*)(*mem); (*mem)=((*arr)+x); } template inline void walloc1d(T **arr, int x1, int x2, void **mem = &wmem){ walloc1d(arr, x2-x1, mem); (*arr) -= x1; } template void sortA_L(int N, T1 a[], void *mem = wmem){ sort(a, a+N); } struct Rand{ unsigned x; unsigned y; unsigned z; unsigned w; Rand(void){ x=123456789; y=362436069; z=521288629; w=(unsigned)time(NULL); } Rand(unsigned seed){ x=123456789; y=362436069; z=521288629; w=seed; } inline unsigned get(void){ unsigned t; t = (x^(x<<11)); x=y; y=z; z=w; w = (w^(w>>19))^(t^(t>>8)); return w; } inline double getUni(void){ return get()/4294967296.0; } inline int get(int a){ return (int)(a*getUni()); } inline int get(int a, int b){ return a+(int)((b-a+1)*getUni()); } inline long long get(long long a){ return(long long)(a*getUni()); } inline long long get(long long a, long long b){ return a+(long long)((b-a+1)*getUni()); } inline double get(double a, double b){ return a+(b-a)*getUni(); } inline int getExp(int a){ return(int)(exp(getUni()*log(a+1.0))-1.0); } inline int getExp(int a, int b){ return a+(int)(exp(getUni()*log((b-a+1)+1.0))-1.0); } } ; inline int my_getchar_unlocked(){ static char buf[1048576]; static int s = 1048576; static int e = 1048576; if(s == e && e == 1048576){ e = fread_unlocked(buf, 1, 1048576, stdin); s = 0; } if(s == e){ return EOF; } return buf[s++]; } inline void my_putchar_unlocked(const int k){ putchar_unlocked(k); } inline void wt_L(char a){ my_putchar_unlocked(a); } inline void wt_L(int x){ int s=0; int m=0; char f[10]; if(x<0){ m=1; x=-x; } while(x){ f[s++]=x%10; x/=10; } if(!s){ f[s++]=0; } if(m){ my_putchar_unlocked('-'); } while(s--){ my_putchar_unlocked(f[s]+'0'); } } inline void wt_L(const char c[]){ int i=0; for(i=0;c[i]!='\0';i++){ my_putchar_unlocked(c[i]); } } template inline S chmin(S &a, T b){ if(a>b){ a=b; } return a; } template inline S chmax(S &a, T b){ if(a=0){ t=d[t]; } while(d[a]>=0){ k=d[a]; d[a]=t; a=k; } return a; } inline int connect(int a, int b){ if(d[a]>=0){ a=get(a); } if(d[b]>=0){ b=get(b); } if(a==b){ return 0; } if(d[a] < d[b]){ d[a] += d[b]; d[b] = a; } else{ d[b] += d[a]; d[a] = b; } return 1; } inline int operator()(int a){ return get(a); } inline int operator()(int a, int b){ return connect(a,b); } inline int& operator[](const int a){ return d[a]; } inline int size(int a){ a = get(a); return -d[a]; } inline int sizeList(int res[]){ int i; int sz=0; for(i=(0);i<(N);i++){ if(d[i]<0){ res[sz++] = -d[i]; } } return sz; } inline int comp(int res[], void *mem = wmem){ int i; int sz=0; int*cnt; walloc1d(&cnt, N, &mem); for(i=(0);i<(N);i++){ cnt[i] = 0; } for(i=(0);i<(N);i++){ cnt[get(i)] = 1; } for(i=(0);i<(N);i++){ if(cnt[i]){ cnt[i] = sz++; } } for(i=(0);i<(N);i++){ res[i] = cnt[get(i)]; } return sz; } } ; long long cReader_ll(long long mn, long long mx, char nx){ int i; int fg = 0; int m = 1; int f = -1; long long res = 0; double tmp = 0; for(;;){ i = my_getchar_unlocked(); if(fg==0 && i=='-'){ fg++; m = -1; } else if('0' <= i && i <= '9'){ fg++; if(f == -1){ f = i - '0'; } res = 10 * res + i - '0'; tmp = 10 * tmp + i - '0'; assert(tmp < 1e20); } else{ break; } } assert(tmp / 2 <= res); assert((m==1 && fg >= 1) || (m==-1 && fg >= 2)); assert(mn <= m * res && m * res <= mx); assert(!(res == 0 && m == -1)); assert(!(res != 0 && f == 0)); assert(!(res == 0 && fg >= 2)); assert(i == nx); return m * res; } void cReader_eof(){ int i; i = my_getchar_unlocked(); assert(i == EOF); } int N; int A[100000]; int B[100000]; int Z[100000]; unionFind uf; unionFind uf1; unionFind uf2; vector oobaka(){ int i; int j; int k; set> a; set> b; vector res; int mx; int tmp; for(k=(0);k<(N);k++){ vector v; queue> q; a.clear(); b.clear(); mx = 0; v.clear(); for(i=(0);i<(N);i++){ v.push_back(A[i]); } while(q.size()){ q.pop(); } q.push(v); a.insert(v); while(q.size()){ v = q.front(); q.pop(); for(i=(1);i<(N);i++){ if(v[i-1] <= k && v[i] <= k){ swap(v[i-1], v[i]); if(a.count(v)==0){ q.push(v); a.insert(v); } swap(v[i-1], v[i]); } } } v.clear(); for(i=(0);i<(N);i++){ v.push_back(B[i]); } while(q.size()){ q.pop(); } q.push(v); b.insert(v); while(q.size()){ v = q.front(); q.pop(); for(i=(1);i<(N);i++){ if(v[i-1] <= k && v[i] <= k){ swap(v[i-1], v[i]); if(b.count(v)==0){ q.push(v); b.insert(v); } swap(v[i-1], v[i]); } } } for(auto aa : a){ for(auto bb : b){ int x; int tmp = 0; for(x=(0);x<(N);x++){ for(i=(0);i<(N);i++){ if(aa[i]==x){ break; } } for(j=(0);j<(N);j++){ if(bb[j]==Z[x]){ break; } } if(i==j){ tmp++; } } chmax(mx, tmp); } } res.push_back(mx); } return res; } vector baka(){ int i; int j; int k; vector res; for(k=(0);k<(N);k++){ int m; uf1.init(N); uf2.init(N); int use[N] = {}; int tmp = 0; for(i=(1);i<(N);i++){ if(A[i-1] <= k && A[i] <= k){ uf1(i-1, i); } } for(i=(1);i<(N);i++){ if(B[i-1] <= k && B[i] <= k){ uf2(i-1, i); } } for(m=(0);m<(N);m++){ int x; for(i=(0);i<(N);i++){ if(A[i] == m){ break; } } for(j=(0);j<(N);j++){ if(B[j] == Z[m]){ break; } } for(x=(0);x<(N);x++){ if(!use[x] && uf1(i) == uf1(x) && uf2(j) == uf2(x)){ tmp++; use[x] = 1; break; } } } res.push_back(tmp); } return res; } struct RBMM{ int mn; int mx; } ; struct RBV{ int mn; int mx; int v; } ; RBMM rollbackUnionFind_md_func(RBMM a, RBMM b){ return {min_L(a.mn, b.mn),max_L(a.mx, b.mx)}; } RBV rollbackUnionFind_md_func(RBV a, RBV b){ return {min_L(a.mn, b.mn),max_L(a.mx, b.mx), a.v + b.v}; } template struct rollbackUnionFind_md{ int*d; int N; int M; int*h1; int*h2; int sz; int snapsz; T*v; T*h3; inline void malloc(const int n){ d = (int*)std::malloc(n*sizeof(int)); h1 = (int*)std::malloc(2*n*sizeof(int)); h2 = (int*)std::malloc(2*n*sizeof(int)); v = (T*)std::malloc(n*sizeof(T)); h3 = (T*)std::malloc(2*n*sizeof(T)); M = n; } inline void malloc(const int n, const int fg){ d = (int*)std::malloc(n*sizeof(int)); h1 = (int*)std::malloc(2*n*sizeof(int)); h2 = (int*)std::malloc(2*n*sizeof(int)); v = (T*)std::malloc(n*sizeof(T)); h3 = (T*)std::malloc(2*n*sizeof(T)); M = n; if(fg){ init(n); } } inline void free(void){ std::free(d); std::free(h1); std::free(h2); std::free(v); std::free(h3); } inline void walloc(const int n, void **mem=&wmem){ walloc1d(&d, n, mem); walloc1d(&h1, 2*n, mem); walloc1d(&h2, 2*n, mem); walloc1d(&v, n, mem); walloc1d(&h3, 2*n, mem); M = n; } inline void walloc(const int n, const int fg, void **mem=&wmem){ walloc1d(&d, n, mem); walloc1d(&h1, 2*n, mem); walloc1d(&h2, 2*n, mem); walloc1d(&v, n, mem); walloc1d(&h3, 2*n, mem); M = n; if(fg){ init(n); } } inline void init(const int n){ int i; N = n; for(i=(0);i<(n);i++){ d[i] = -1; } sz = 0; snapsz = 0; } inline void init(void){ init(M); } inline int get(int a){ while(d[a]>=0){ a=d[a]; } return a; } inline int connect(int a, int b){ if(d[a]>=0){ a=get(a); } if(d[b]>=0){ b=get(b); } if(a==b){ return 0; } h1[sz] = a; h2[sz] = d[a]; h3[sz] = v[a]; sz++; h1[sz] = b; h2[sz] = d[b]; h3[sz] = v[b]; sz++; if(d[a] < d[b]){ d[a] += d[b]; d[b] = a; v[a] = rollbackUnionFind_md_func(v[a], v[b]); } else{ d[b] += d[a]; d[a] = b; v[b] = rollbackUnionFind_md_func(v[a], v[b]); } return 1; } inline int operator()(int a){ return get(a); } inline int operator()(int a, int b){ return connect(a,b); } inline int& operator[](const int a){ return d[a]; } inline int size(int a){ a = get(a); return -d[a]; } inline int sizeList(int res[]){ int i; int sz=0; for(i=(0);i<(N);i++){ if(d[i]<0){ res[sz++] = -d[i]; } } return sz; } inline int undo(void){ if(sz==0){ return 0; } sz--; d[h1[sz]] = h2[sz]; v[h1[sz]] = h3[sz]; sz--; d[h1[sz]] = h2[sz]; v[h1[sz]] = h3[sz]; if(snapsz > sz){ snapsz = sz; } return 1; } inline void snapshot(void){ snapsz = sz; } inline void rollback(void){ while(snapsz < sz){ sz--; d[h1[sz]] = h2[sz]; v[h1[sz]] = h3[sz]; sz--; d[h1[sz]] = h2[sz]; v[h1[sz]] = h3[sz]; } } inline T getVal(int a){ a = get(a); return v[a]; } inline void setVal(int a, T val){ a = get(a); v[a] = val; } } ; rollbackUnionFind_md ruf1; rollbackUnionFind_md ruf2; rollbackUnionFind_md ruf; int invA[100000]; int invB[100000]; int ftm[100000]; int plc[100000]; int ind[100000]; int liss[20][100000]; set unused; int solve_uf_doit_sub(int i, int j){ int res = 0; RBV v; if(ruf(i) != ruf(j)){ v = ruf.getVal(i); res -=min_L(v.v, v.mx - v.mn + 1); v = ruf.getVal(j); res -=min_L(v.v, v.mx - v.mn + 1); ruf(i,j); v = ruf.getVal(i); res +=min_L(v.v, v.mx - v.mn + 1); } return res; } int solve_uf_doit(int i){ int x; int res = 0; RBV v; x = invA[i]; if(x-1 >= 0 && A[x-1] <= i){ ruf1(x-1,x); if(ruf2(x-1) == ruf2(x)){ res += solve_uf_doit_sub(x-1, x); } } if(x+1 < N && A[x+1] <= i){ ruf1(x,x+1); if(ruf2(x) == ruf2(x+1)){ res += solve_uf_doit_sub(x, x+1); } } x = invB[i]; if(x-1 >= 0 && B[x-1] <= i){ ruf2(x-1,x); if(ruf1(x-1) == ruf1(x)){ res += solve_uf_doit_sub(x-1, x); } } if(x+1 < N && B[x+1] <= i){ ruf2(x,x+1); if(ruf1(x) == ruf1(x+1)){ res += solve_uf_doit_sub(x, x+1); } } return res; } void solve_uf_doit2(int i){ int x; x = invA[i]; if(x-1 >= 0 && A[x-1] <= i){ ruf1(x-1,x); } if(x+1 < N && A[x+1] <= i){ ruf1(x,x+1); } x = invB[i]; if(x-1 >= 0 && B[x-1] <= i){ ruf2(x-1,x); } if(x+1 < N && B[x+1] <= i){ ruf2(x,x+1); } } void solve_sub(int s, int e, int sz, int *lis, int dep){ int i; int j; int k; int m; int sz1; int sz2; int x; int y; set::iterator it; if(sz == 0){ for(i=(s);i<(e);i++){ solve_uf_doit2(i); } return; } if(e-s==1){ solve_uf_doit2(s); for(i=(0);i<(sz);i++){ x =max_L(ruf1.getVal(invA[lis[i]]).mn, ruf2.getVal(invB[Z[lis[i]]]).mn); y =min_L(ruf1.getVal(invA[lis[i]]).mx, ruf2.getVal(invB[Z[lis[i]]]).mx); it = unused.lower_bound(x); if(it == unused.end() || *it > y){ continue; } ftm[lis[i]] = s; unused.erase(it); } return; } m = (s + e) / 2; for(i=(s);i<(m);i++){ solve_uf_doit2(i); } sz1 = sz2 = 0; for(i=(0);i<(sz);i++){ x =max_L(ruf1.getVal(invA[lis[i]]).mn, ruf2.getVal(invB[Z[lis[i]]]).mn); y =min_L(ruf1.getVal(invA[lis[i]]).mx, ruf2.getVal(invB[Z[lis[i]]]).mx); it = unused.lower_bound(x); if(x > y){ continue; } if(it == unused.end() || *it > y){ continue; } liss[dep+1][sz1++] = lis[i]; } if(sz1){ ruf1.rollback(); ruf2.rollback(); ruf.rollback(); chmin(sz1, 10 * (m - s)); solve_sub(s,m,sz1,liss[dep+1],dep+1); } for(i=(0);i<(sz);i++){ if(ftm[lis[i]] == -1){ lis[sz2++] = lis[i]; } } ruf1.snapshot(); ruf2.snapshot(); ruf.snapshot(); solve_sub(m,e,sz2,lis,dep+1); for(i=(m);i<(e);i++){ solve_uf_doit2(i); } } vector solve(){ int i; int j; int k; vector res(N); RBV v; int resval; ruf1.init(N); ruf2.init(N); ruf.init(N); for(i=(0);i<(N);i++){ ruf1.setVal(i, {i,i}); } for(i=(0);i<(N);i++){ ruf2.setVal(i, {i,i}); } for(i=(0);i<(N);i++){ ruf.setVal(i, {i,i,0}); } for(i=(0);i<(N);i++){ invA[A[i]] = i; } for(i=(0);i<(N);i++){ invB[B[i]] = i; } for(i=(0);i<(N);i++){ liss[0][i] = i; } for(i=(0);i<(N);i++){ ftm[i] = -1; } for(i=(0);i<(N);i++){ unused.insert(i); } solve_sub(0, N, N, liss[0], 0); for(i=(0);i<(N);i++){ ind[i] = i; } sortA_L(N,ftm); k = 0; resval = 0; for(i=(0);i<(N);i++){ resval += solve_uf_doit(i); while(k < N && ftm[k] == i){ resval++; k++; } res[i] = resval; } return res; } int cnt[100000]; int main(){ int i; wmem = memarr; Rand rnd; uf.walloc(100000); uf1.walloc(100000); uf2.walloc(100000); ruf1.walloc(100000); ruf2.walloc(100000); ruf.walloc(100000); N = cReader_ll(2, 100000, '\n'); for(i=(0);i<(N);i++){ A[i] = cReader_ll(1, N, i==N-1?'\n':' ') - 1; } for(i=(0);i<(N);i++){ B[i] = cReader_ll(1, N, i==N-1?'\n':' ') - 1; } for(i=(0);i<(N);i++){ Z[i] = cReader_ll(1, N, i==N-1?'\n':' ') - 1; } for(i=(0);i<(N);i++){ cnt[i] = 0; } for(i=(0);i<(N);i++){ cnt[A[i]]++; } for(i=(0);i<(N);i++){ assert(cnt[i]==1); } for(i=(0);i<(N);i++){ cnt[i] = 0; } for(i=(0);i<(N);i++){ cnt[B[i]]++; } for(i=(0);i<(N);i++){ assert(cnt[i]==1); } for(i=(0);i<(N);i++){ cnt[i] = 0; } for(i=(0);i<(N);i++){ cnt[Z[i]]++; } for(i=(0);i<(N);i++){ assert(cnt[i]==1); } cReader_eof(); vector res; res = solve(); { int x7rJoZt9; for(x7rJoZt9=(0);x7rJoZt9<(N);x7rJoZt9++){ wt_L(res[x7rJoZt9]); wt_L('\n'); } } return 0; puts(""); for(;;){ vector r1; vector r2; vector r3; N = rnd.get(2,6); for(i=(0);i<(N);i++){ A[i] = B[i] = Z[i] = i; } random_shuffle(A,A+N); random_shuffle(B,B+N); random_shuffle(Z,Z+N); wt_L(N); wt_L(' '); wt_L(":"); wt_L(' '); { int MsHTBJEM; for(MsHTBJEM=(0);MsHTBJEM<(N);MsHTBJEM++){ wt_L(A[MsHTBJEM]); wt_L(' '); } } wt_L(":"); wt_L(' '); { int dd01qmIa; for(dd01qmIa=(0);dd01qmIa<(N);dd01qmIa++){ wt_L(B[dd01qmIa]); wt_L(' '); } } wt_L(":"); wt_L(' '); { int Zs_a1Ftk; if(N==0){ wt_L('\n'); } else{ for(Zs_a1Ftk=(0);Zs_a1Ftk<(N-1);Zs_a1Ftk++){ wt_L(Z[Zs_a1Ftk]); wt_L(' '); } wt_L(Z[Zs_a1Ftk]); wt_L('\n'); } } r1 = oobaka(); r2 = baka(); r3 = solve(); { int B8qxtQZy; if(N==0){ wt_L('\n'); } else{ for(B8qxtQZy=(0);B8qxtQZy<(N-1);B8qxtQZy++){ wt_L(r1[B8qxtQZy]); wt_L(' '); } wt_L(r1[B8qxtQZy]); wt_L('\n'); } } { int pRrsnthG; if(N==0){ wt_L('\n'); } else{ for(pRrsnthG=(0);pRrsnthG<(N-1);pRrsnthG++){ wt_L(r2[pRrsnthG]); wt_L(' '); } wt_L(r2[pRrsnthG]); wt_L('\n'); } } { int Vdr_4aft; if(N==0){ wt_L('\n'); } else{ for(Vdr_4aft=(0);Vdr_4aft<(N-1);Vdr_4aft++){ wt_L(r3[Vdr_4aft]); wt_L(' '); } wt_L(r3[Vdr_4aft]); wt_L('\n'); } } assert(r1==r2); assert(r2==r3); } return 0; } // cLay version 20210904-1 // --- original code --- // int N, A[1d5], B[], Z[]; // unionFind uf, uf1, uf2; // // VI oobaka(){ // int i, j, k; // set a, b; // VI res; // int mx, tmp; // rep(k,N){ // VI v; // queue q; // a.clear(); // b.clear(); // // mx = 0; // // v.clear(); // rep(i,N) v.push_back(A[i]); // while(q.size()) q.pop(); // q.push(v); // a.insert(v); // while(q.size()){ // v = q.front(); // q.pop(); // rep(i,1,N) if(v[i-1] <= k && v[i] <= k){ // swap(v[i-1], v[i]); // if(a.count(v)==0) q.push(v), a.insert(v); // swap(v[i-1], v[i]); // } // } // // v.clear(); // rep(i,N) v.push_back(B[i]); // while(q.size()) q.pop(); // q.push(v); // b.insert(v); // while(q.size()){ // v = q.front(); // q.pop(); // rep(i,1,N) if(v[i-1] <= k && v[i] <= k){ // swap(v[i-1], v[i]); // if(b.count(v)==0) q.push(v), b.insert(v); // swap(v[i-1], v[i]); // } // } // // for(auto aa : a) for(auto bb : b){ // int tmp = 0; // rep(x,N){ // rep(i,N) if(aa[i]==x) break; // rep(j,N) if(bb[j]==Z[x]) break; // if(i==j) tmp++; // } // mx >?= tmp; // } // // res.push_back(mx); // } // return res; // } // // VI baka(){ // int i, j, k; // VI res; // rep(k,N){ // uf1.init(N); // uf2.init(N); // int use[N] = {}, tmp = 0; // // rep(i,1,N) if(A[i-1] <= k && A[i] <= k) uf1(i-1, i); // rep(i,1,N) if(B[i-1] <= k && B[i] <= k) uf2(i-1, i); // // rep(m,N){ // rep(i,N) if(A[i] == m) break; // rep(j,N) if(B[j] == Z[m]) break; // rep(x,N) if(!use[x] && uf1(i) == uf1(x) && uf2(j) == uf2(x)){ // tmp++; // use[x] = 1; // break; // } // } // // res.push_back(tmp); // } // return res; // } // // // // // struct RBMM {int mn, mx;}; // struct RBV {int mn, mx, v;}; // // RBMM rollbackUnionFind_md_func(RBMM a, RBMM b){ // return {min(a.mn, b.mn), max(a.mx, b.mx)}; // } // // RBV rollbackUnionFind_md_func(RBV a, RBV b){ // return {min(a.mn, b.mn), max(a.mx, b.mx), a.v + b.v}; // } // // template // struct rollbackUnionFind_md{ // int *d, N, M; // int *h1, *h2, sz, snapsz; // T *v, *h3; // inline void malloc(const int n){ // d = (int*)std::malloc(n*sizeof(int)); // h1 = (int*)std::malloc(2*n*sizeof(int)); // h2 = (int*)std::malloc(2*n*sizeof(int)); // v = (T*)std::malloc(n*sizeof(T)); // h3 = (T*)std::malloc(2*n*sizeof(T)); // M = n; // } // inline void malloc(const int n, const int fg){ // d = (int*)std::malloc(n*sizeof(int)); // h1 = (int*)std::malloc(2*n*sizeof(int)); // h2 = (int*)std::malloc(2*n*sizeof(int)); // v = (T*)std::malloc(n*sizeof(T)); // h3 = (T*)std::malloc(2*n*sizeof(T)); // M = n; // if(fg) init(n); // } // inline void free(void){ // std::free(d); // std::free(h1); // std::free(h2); // std::free(v); // std::free(h3); // } // inline void walloc(const int n, void **mem=&wmem){ // walloc1d(&d, n, mem); // walloc1d(&h1, 2*n, mem); // walloc1d(&h2, 2*n, mem); // walloc1d(&v, n, mem); // walloc1d(&h3, 2*n, mem); // M = n; // } // inline void walloc(const int n, const int fg, void **mem=&wmem){ // walloc1d(&d, n, mem); // walloc1d(&h1, 2*n, mem); // walloc1d(&h2, 2*n, mem); // walloc1d(&v, n, mem); // walloc1d(&h3, 2*n, mem); // M = n; // if(fg) init(n); // } // inline void init(const int n){ // int i; // N = n; // rep(i,n) d[i] = -1; // sz = 0; // snapsz = 0; // } // inline void init(void){ // init(M); // } // inline int get(int a){ // while(d[a]>=0) a=d[a]; // return a; // } // inline int connect(int a, int b){ // if(d[a]>=0) a=get(a); // if(d[b]>=0) b=get(b); // if(a==b) return 0; // h1[sz] = a; h2[sz] = d[a]; h3[sz] = v[a]; sz++; // h1[sz] = b; h2[sz] = d[b]; h3[sz] = v[b]; sz++; // if(d[a] < d[b]) d[a] += d[b], d[b] = a, v[a] = rollbackUnionFind_md_func(v[a], v[b]); // else d[b] += d[a], d[a] = b, v[b] = rollbackUnionFind_md_func(v[a], v[b]); // return 1; // } // inline int operator()(int a){ // return get(a); // } // inline int operator()(int a, int b){ // return connect(a,b); // } // inline int& operator[](const int a){ // return d[a]; // } // inline int size(int a){ // a = get(a); // return -d[a]; // } // inline int sizeList(int res[]){ // int i, sz=0; // rep(i,N) if(d[i]<0) res[sz++] = -d[i]; // return sz; // } // inline int undo(void){ // if(sz==0) return 0; // sz--; d[h1[sz]] = h2[sz]; v[h1[sz]] = h3[sz]; // sz--; d[h1[sz]] = h2[sz]; v[h1[sz]] = h3[sz]; // if(snapsz > sz) snapsz = sz; // return 1; // } // inline void snapshot(void){ // snapsz = sz; // } // inline void rollback(void){ // while(snapsz < sz){ // sz--; d[h1[sz]] = h2[sz]; v[h1[sz]] = h3[sz]; // sz--; d[h1[sz]] = h2[sz]; v[h1[sz]] = h3[sz]; // } // } // inline T getVal(int a){ // a = get(a); // return v[a]; // } // inline void setVal(int a, T val){ // a = get(a); // v[a] = val; // } // }; // // // rollbackUnionFind_md ruf1, ruf2; // rollbackUnionFind_md ruf; // // int invA[1d5], invB[1d5]; // int ftm[1d5], plc[1d5], ind[1d5]; // int liss[20][1d5]; // set unused; // // int solve_uf_doit_sub(int i, int j){ // int res = 0; // RBV v; // if(ruf(i) != ruf(j)){ // v = ruf.getVal(i); // res -= min(v.v, v.mx - v.mn + 1); // v = ruf.getVal(j); // res -= min(v.v, v.mx - v.mn + 1); // ruf(i,j); // v = ruf.getVal(i); // res += min(v.v, v.mx - v.mn + 1); // } // return res; // } // // int solve_uf_doit(int i){ // int x, res = 0; // RBV v; // x = invA[i]; // if(x-1 >= 0 && A[x-1] <= i){ // ruf1(x-1,x); // if(ruf2(x-1) == ruf2(x)) res += solve_uf_doit_sub(x-1, x); // } // if(x+1 < N && A[x+1] <= i){ // ruf1(x,x+1); // if(ruf2(x) == ruf2(x+1)) res += solve_uf_doit_sub(x, x+1); // } // x = invB[i]; // if(x-1 >= 0 && B[x-1] <= i){ // ruf2(x-1,x); // if(ruf1(x-1) == ruf1(x)) res += solve_uf_doit_sub(x-1, x); // } // if(x+1 < N && B[x+1] <= i){ // ruf2(x,x+1); // if(ruf1(x) == ruf1(x+1)) res += solve_uf_doit_sub(x, x+1); // } // return res; // } // // void solve_uf_doit2(int i){ // int x; // x = invA[i]; // if(x-1 >= 0 && A[x-1] <= i){ // ruf1(x-1,x); // } // if(x+1 < N && A[x+1] <= i){ // ruf1(x,x+1); // } // x = invB[i]; // if(x-1 >= 0 && B[x-1] <= i){ // ruf2(x-1,x); // } // if(x+1 < N && B[x+1] <= i){ // ruf2(x,x+1); // } // } // // void solve_sub(int s, int e, int sz, int *lis, int dep){ // int i, j, k, m, sz1, sz2, x, y; // set::iterator it; // // if(sz == 0){ // rep(i,s,e) solve_uf_doit2(i); // return; // } // // if(e-s==1){ // solve_uf_doit2(s); // rep(i,sz){ // x = max(ruf1.getVal(invA[lis[i]]).mn, ruf2.getVal(invB[Z[lis[i]]]).mn); // y = min(ruf1.getVal(invA[lis[i]]).mx, ruf2.getVal(invB[Z[lis[i]]]).mx); // it = unused.lower_bound(x); // if(it == unused.end() || *it > y) continue; // ftm[lis[i]] = s; // unused.erase(it); // // wt(lis[i],s,x,y); // } // return; // } // // m = (s + e) / 2; // rep(i,s,m) solve_uf_doit2(i); // // sz1 = sz2 = 0; // rep(i,sz){ // x = max(ruf1.getVal(invA[lis[i]]).mn, ruf2.getVal(invB[Z[lis[i]]]).mn); // y = min(ruf1.getVal(invA[lis[i]]).mx, ruf2.getVal(invB[Z[lis[i]]]).mx); // it = unused.lower_bound(x); // if(x > y) continue; // if(it == unused.end() || *it > y) continue; // liss[dep+1][sz1++] = lis[i]; // } // // if(sz1){ // ruf1.rollback(); // ruf2.rollback(); // ruf.rollback(); // // sz1