#include using namespace std; using int128 = __int128_t; using int64 = long long; using int32 = int; using uint128 = __uint128_t; using uint64 = unsigned long long; using uint32 = unsigned int; #define ALL(obj) (obj).begin(),(obj).end() template using priority_queue_reverse = priority_queue,greater>; constexpr int64 MOD = 1'000'000'000LL + 7; //' constexpr int64 MOD2 = 998244353; constexpr int64 HIGHINF = 1'000'000'000'000'000'000LL; constexpr int64 LOWINF = 1'000'000'000'000'000LL; //' constexpr long double PI = 3.1415926535897932384626433L; template vector multivector(size_t N,T init){return vector(N,init);} template auto multivector(size_t N,T... t){return vector(N,multivector(t...));} template void corner(bool flg, T hoge) {if (flg) {cout << hoge << endl; exit(0);}} template ostream &operator<<(ostream &o, const map&obj) {o << "{"; for (auto &x : obj) o << " {" << x.first << " : " << x.second << "}" << ","; o << " }"; return o;} template ostream &operator<<(ostream &o, const set&obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr) o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;} template ostream &operator<<(ostream &o, const multiset&obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr) o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;} template ostream &operator<<(ostream &o, const vector&obj) {o << "{"; for (int i = 0; i < (int)obj.size(); ++i)o << (i > 0 ? ", " : "") << obj[i]; o << "}"; return o;} template ostream &operator<<(ostream &o, const deque&obj) {o << "{"; for (int i = 0; i < (int)obj.size(); ++i)o << (i > 0 ? ", " : "") << obj[i]; o << "}"; return o;} template ostream &operator<<(ostream &o, const pair&obj) {o << "{" << obj.first << ", " << obj.second << "}"; return o;} void print(void) {cout << endl;} template void print(Head&& head) {cout << head;print();} template void print(Head&& head, Tail&&... tail) {cout << head << " ";print(forward(tail)...);} template void chmax(T& a, const T b){a=max(a,b);} template void chmin(T& a, const T b){a=min(a,b);} vector split(const string &str, const char delemiter) {vector res;stringstream ss(str);string buffer; while( getline(ss, buffer, delemiter) ) res.push_back(buffer); return res;} inline constexpr int msb(int x) {return x?31-__builtin_clz(x):-1;} inline constexpr int64 ceil_div(const int64 a,const int64 b) {return (a+(b-1))/b;}// return ceil(a/b) void YN(bool flg) {cout << (flg ? "YES" : "NO") << endl;} void Yn(bool flg) {cout << (flg ? "Yes" : "No") << endl;} void yn(bool flg) {cout << (flg ? "yes" : "no") << endl;} /* * @title BinaryIndexedTree * @docs md/segment/BinaryIndexedTree.md */ template class BinaryIndexedTree { using TypeNode = typename Operator::TypeNode; size_t depth; size_t length; size_t num; vector node; public: //[0,N) constructed, inplace [0,1) + [1,N+1) //you can ignore inplace offset BinaryIndexedTree(const size_t num) : num(num) { for (depth = 1,length = 1; length < num; depth++,length *= 2); node.resize(length+1, Operator::unit_node); } //[idx,idx+1) update void update(size_t idx, TypeNode var) { assert(0 <= idx && idx < length); for (++idx; idx <= length; idx += idx & -idx) node[idx] = Operator::func_node(node[idx],var); } //[0,idx) get TypeNode get(size_t idx) { TypeNode ret = Operator::unit_node; for (idx = min(length,idx); idx > 0; idx -= idx & -idx) ret = Operator::func_node(ret,node[idx]); return ret; } //return [0,length] int binary_search(TypeNode var) { if(!Operator::func_check(node.back(),var)) return num; TypeNode ret = Operator::unit_node; size_t off = 0; for(size_t idx = length; idx; idx >>= 1){ if(off + idx <= length && !Operator::func_check(Operator::func_node(ret,node[off+idx]),var)) { off += idx; ret = Operator::func_node(ret,node[off]); } } return min(off,num); } void print() { cout << "{ " << get(1); for(int i = 1; i < length; ++i) cout << ", " << get(i+1); cout << " }" << endl; } }; template struct NodePrefixSumPointAdd { using TypeNode = T; inline static constexpr TypeNode unit_node = 0; inline static constexpr TypeNode func_node(const TypeNode& l,const TypeNode& r){return l+r;} // Binary Search for first index where func_check is true inline static constexpr bool func_check(const TypeNode nodeVal,const TypeNode& val){return val <= nodeVal;} }; /* * @title RangeInversionQuery - 区間転倒数 * @docs md/segment/RangeModeQuery.md */ template class RangeInversionQuery { vector compressed; vector prefix_inv; vector suffix_inv; vector> sqrt_bucket_freq; vector> sqrt_bucket_inv; vector> sqrt_bucket_sort_index; vector sqrt_bucket_size; size_t N,B; public: RangeInversionQuery(const vector& ar, T pre=-1) : compressed(ar.size()),prefix_inv(ar.size()),suffix_inv(ar.size()) { N = ar.size(); B = sqrt(N) + 1; //zarts { vector> ord(N); for(size_t i=0;i freq(N+1,0); for(size_t i=0;i> bit(N+1); for(size_t i=0;i(B,0)); for(size_t i=0;i br + 1) { inv += prefix_inv[r-1]; if(l) inv -= prefix_inv[l-1]; bl = l / B; long long sum = 0; for(size_t i: sqrt_bucket_sort_index[bl]) { if(r <= i) continue; if(l <= i) sum++; else inv -= sum; } } else { if(bl<=br) inv += sqrt_bucket_inv[bl][br]; inv += suffix_inv[l]; inv += prefix_inv[r-1]; size_t ml = bl*B; for(size_t i=l;i ql,qr; for(size_t i: sqrt_bucket_sort_index[bl-1]) ql.push(i); for(size_t i: sqrt_bucket_sort_index[br+1]) qr.push(i); long long sum = 0; while(ql.size()) { auto xl = ql.front(); if(xl < l) { ql.pop(); continue; } if(qr.empty()) { inv += sum; ql.pop(); continue; } auto xr = qr.front(); if(xr >= r) { qr.pop(); continue; } if(compressed[xl] > compressed[xr]) { sum++; qr.pop(); continue; } else { inv += sum; ql.pop(); continue; } } } return inv; } }; /** * @url * @est */ int main() { int N; cin >> N; vector a(N),b(N),c(N+1); for(int i = 0; i < N; ++i) { cin >> a[i]; c[a[i]]=i; } for(int i = 0; i < N; ++i) { cin >> b[i]; b[i]=c[b[i]]; } vector idx(N); iota(ALL(idx),0); sort(ALL(idx),[&](int l,int r){return b[l] tmp(M+N,0); for(int i=0;i riq(tmp); cout << riq.get(M,M+N) << endl; return 0; }