結果
問題 | No.1115 二つの数列 / Two Sequences |
ユーザー | T101010101 |
提出日時 | 2023-07-14 21:28:52 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 111 ms / 2,000 ms |
コード長 | 11,361 bytes |
コンパイル時間 | 5,987 ms |
コンパイル使用メモリ | 336,716 KB |
実行使用メモリ | 20,092 KB |
最終ジャッジ日時 | 2024-09-16 06:17:59 |
合計ジャッジ時間 | 9,965 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 9 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,248 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 92 ms
18,704 KB |
testcase_04 | AC | 106 ms
19,604 KB |
testcase_05 | AC | 91 ms
18,824 KB |
testcase_06 | AC | 82 ms
18,256 KB |
testcase_07 | AC | 101 ms
19,052 KB |
testcase_08 | AC | 5 ms
5,376 KB |
testcase_09 | AC | 59 ms
10,800 KB |
testcase_10 | AC | 111 ms
17,276 KB |
testcase_11 | AC | 2 ms
5,376 KB |
testcase_12 | AC | 106 ms
20,092 KB |
testcase_13 | AC | 103 ms
17,148 KB |
testcase_14 | AC | 103 ms
20,092 KB |
testcase_15 | AC | 2 ms
5,376 KB |
testcase_16 | AC | 2 ms
5,376 KB |
testcase_17 | AC | 2 ms
5,376 KB |
testcase_18 | AC | 2 ms
5,376 KB |
testcase_19 | AC | 2 ms
5,376 KB |
testcase_20 | AC | 2 ms
5,376 KB |
testcase_21 | AC | 2 ms
5,376 KB |
testcase_22 | AC | 2 ms
5,376 KB |
testcase_23 | AC | 9 ms
5,376 KB |
testcase_24 | AC | 35 ms
10,364 KB |
testcase_25 | AC | 80 ms
17,716 KB |
testcase_26 | AC | 17 ms
6,772 KB |
testcase_27 | AC | 42 ms
10,692 KB |
testcase_28 | AC | 56 ms
12,228 KB |
testcase_29 | AC | 81 ms
17,876 KB |
testcase_30 | AC | 99 ms
18,768 KB |
testcase_31 | AC | 23 ms
7,464 KB |
testcase_32 | AC | 14 ms
5,376 KB |
testcase_33 | AC | 94 ms
16,272 KB |
testcase_34 | AC | 2 ms
5,376 KB |
testcase_35 | AC | 2 ms
5,376 KB |
testcase_36 | AC | 2 ms
5,376 KB |
testcase_37 | AC | 2 ms
5,376 KB |
testcase_38 | AC | 2 ms
5,376 KB |
testcase_39 | AC | 2 ms
5,376 KB |
ソースコード
#pragma region Macros // #pragma GCC optimize("O3,unroll-loops") // #pragma GCC target("sse,sse2,sse3,ssse3,sse4,fma,abm,mmx,avx,avx2") #include <bits/extc++.h> #include <atcoder/all> using namespace atcoder; using namespace std; using namespace __gnu_pbds; // #include <boost/multiprecision/cpp_dec_float.hpp> // #include <boost/multiprecision/cpp_int.hpp> // namespace mp = boost::multiprecision; // using Bint = mp::cpp_int; // using Bdouble = mp::number<mp::cpp_dec_float<128>>; #define pb emplace_back #define int ll #define endl '\n' #define sqrt __builtin_sqrtl using ll = long long; using ld = long double; const ld PI = acosl(-1); const int INF = 1 << 30; const ll INFL = 1LL << 61; const int MOD = 998244353; // const int MOD = 1000000007; const ld EPS = 1e-10; const bool equals(ld a, ld b) { return fabs((a) - (b)) < EPS; } const vector<int> dx = {0, 1, 0, -1, 1, 1, -1, -1}; // → ↓ ← ↑ ↘ ↙ ↖ ↗ const vector<int> dy = {1, 0, -1, 0, 1, -1, -1, 1}; struct Edge { int from, to; ll cost; Edge(int to, ll cost) : from(-1), to(to), cost(cost) {} Edge(int from, int to, ll cost) : from(from), to(to), cost(cost) {} }; __attribute__((constructor)) void constructor() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(12); } struct custom_hash { static uint64_t splitmix64(uint64_t x) { x += 0x9e3779b97f4a7c15; x = (x ^ (x >> 30)) * 0xbf58476d1ce4e5b9; x = (x ^ (x >> 27)) * 0x94d049bb133111eb; return x ^ (x >> 31); } size_t operator()(uint64_t x) const { static const uint64_t FIXED_RANDOM = chrono::steady_clock::now().time_since_epoch().count(); return splitmix64(x + FIXED_RANDOM); } }; int POW(int x, int n) { __int128_t ret = 1; // if (ret >= INFL) return INFL; if (n < 0) { cout << "error" << endl; return 0; } else if (x == 1 or n == 0) ret = 1; else if (x == -1 && n % 2 == 0) ret = 1; else if (x == -1) ret = -1; else if (n % 2 == 0) ret = POW(x * x, n / 2); // TODO:オーバーフロー対策 else ret = x * POW(x, n - 1); if (ret > 8e18) ret = 0; return ret; } int floor(int x, int y) { return (x > 0 ? x / y : (x - y + 1) / y); } int ceil(int x, int y) { return (x > 0 ? (x + y - 1) / y : x / y); } ll per(int x, int y) { if (y == 0) { cout << "error" << endl; return INFL; } if (x >= 0 && y > 0) return x / y; if (x >= 0 && y < 0) return x / y - (x % y < 0); if (x < 0 && y < 0) return x / y + (x % y < 0); // if (x < 0 && y > 0) return x / y - (x % y < 0); } ll mod(int x, int y) { if (y == 0) { cout << "error" << endl; return INFL; } if (x >= 0 && y > 0) return x % y; if (x >= 0 && y < 0) return x % y; if (x < 0 && y < 0) { __int128_t ret = x % y; ret += (__int128_t)abs(y) * INFL; ret %= abs(y); return ret; } // if (x < 0 && y > 0) { __int128_t ret = x % y; ret += (__int128_t)abs(y) * INFL; ret %= abs(y); return ret; // } } int modf(ld x) { if (equals(x, 0)) return 0; else if (x < 0) return ceill(x); else return floorl(x); } template <class T> bool chmax(T &a, const T& b) { if (a < b) { a = b; return true; } return false; } template <class T> bool chmin(T &a, const T& b) { if (a > b) { a = b; return true; } return false; } int countl_zero(int N) { return __builtin_clzll(N); } int countl_one(int N) { int ret = 0; while (N % 2) { N /= 2; ret++; } return ret; } int countr_zero(int N) { return __builtin_ctzll(N); } int countr_one(int N) { int ret = 0, k = 63 - __builtin_clzll(N); while (k != -1 && (N & (1LL << k))) { k--; ret++; } return ret; } int popcount(int N) { return __builtin_popcountll(N); } int unpopcount(int N) { return 64 - __builtin_clzll(N) - __builtin_popcountll(N); } int top_bit(int N) { return 63 - __builtin_clzll(N);} // 2^kの位 int bot_bit(int N) { return __builtin_ctz(N);} // 2^kの位 int MSB(int N) { return 1 << (63 - __builtin_clzll(N)); } // mask // LSBどこいった? int bit_width(int N) { return 64 - __builtin_clzll(N); } // 桁数 int ceil_log2(int N) { return 63 - __builtin_clzll(N); } int bit_floor(int N) { return 1 << (63 - __builtin_clzll(N)); } int floor_log2(int N) { return 64 - __builtin_clzll(N-1); } int bit_ceil(int N) { return 1 << (64 - __builtin_clzll(N-1)) - (N==1); } class UnionFind { public: UnionFind() = default; UnionFind(int N) : par(N), sz(N, 1) { iota(par.begin(), par.end(), 0); } int root(int x) { if (par[x] == x) return x; return (par[x] = root(par[x])); } bool unite(int x, int y) { int rx = root(x); int ry = root(y); if (rx == ry) return false; if (sz[rx] < sz[ry]) swap(rx, ry); sz[rx] += sz[ry]; par[ry] = rx; return true; } bool issame(int x, int y) { return (root(x) == root(y)); } int size(int x) { return sz[root(x)]; } vector<vector<int>> groups(int N) { vector<vector<int>> G(N); for (int x = 0; x < N; x++) { G[root(x)].push_back(x); } G.erase( remove_if(G.begin(), G.end(), [&](const vector<int>& V) { return V.empty(); }), G.end()); return G; } private: vector<int> par; vector<int> sz; }; template<int mod> class Modint{ public: int val = 0; Modint(int x = 0) { while (x < 0) x += mod; val = x % mod; } Modint(const Modint &r) { val = r.val; } Modint operator -() { return Modint(-val); } // 単項 Modint operator +(const Modint &r) { return Modint(*this) += r; } Modint operator +(const int &q) { Modint r(q); return Modint(*this) += r; } Modint operator -(const Modint &r) { return Modint(*this) -= r; } Modint operator -(const int &q) { Modint r(q); return Modint(*this) -= r; } Modint operator *(const Modint &r) { return Modint(*this) *= r; } Modint operator *(const int &q) { Modint r(q); return Modint(*this) *= r; } Modint operator /(const Modint &r) { return Modint(*this) /= r; } Modint operator /(const int &q) { Modint r(q); return Modint(*this) /= r; } Modint& operator ++() { val++; if (val >= mod) val -= mod; return *this; } // 前置 Modint operator ++(signed) { ++*this; return *this; } // 後置 Modint& operator --() { val--; if (val < 0) val += mod; return *this; } Modint operator --(signed) { --*this; return *this; } Modint &operator +=(const Modint &r) { val += r.val; if (val >= mod) val -= mod; return *this; } Modint &operator +=(const int &q) { Modint r(q); val += r.val; if (val >= mod) val -= mod; return *this; } Modint &operator -=(const Modint &r) { if (val < r.val) val += mod; val -= r.val; return *this; } Modint &operator -=(const int &q) { Modint r(q); if (val < r.val) val += mod; val -= r.val; return *this; } Modint &operator *=(const Modint &r) { val = val * r.val % mod; return *this; } Modint &operator *=(const int &q) { Modint r(q); val = val * r.val % mod; return *this; } Modint &operator /=(const Modint &r) { int a = r.val, b = mod, u = 1, v = 0; while (b) {int t = a / b; a -= t * b; swap(a, b); u -= t * v; swap(u, v);} val = val * u % mod; if (val < 0) val += mod; return *this; } Modint &operator /=(const int &q) { Modint r(q); int a = r.val, b = mod, u = 1, v = 0; while (b) {int t = a / b; a -= t * b; swap(a, b); u -= t * v; swap(u, v);} val = val * u % mod; if (val < 0) val += mod; return *this; } bool operator ==(const Modint& r) { return this -> val == r.val; } bool operator <(const Modint& r) { return this -> val < r.val; } bool operator >(const Modint& r) { return this -> val > r.val; } bool operator !=(const Modint& r) { return this -> val != r.val; } }; using mint = Modint<MOD>; istream &operator >>(istream &is, mint& x) { int t; is >> t; x = t; return (is); } ostream &operator <<(ostream &os, const mint& x) { return os << x.val; } mint modpow(const mint &x, int n) { // TODO: n <= -1 if (n == 0) return 1; mint t = modpow(x, n / 2); t = t * t; if (n & 1) t = t * x; return t; } int modpow(__int128_t x, int n, int mod) { // TODO: n <= -1 __int128_t ret = 1; while (n > 0) { if (n % 2 == 1) ret = ret * x % mod; x = x * x % mod; n /= 2; } return ret; } int modinv(__int128_t x, int mod) { if (x == 1) return 1; return mod - modinv(mod % x, mod) * (mod / x) % mod; } istream &operator >>(istream &is, __int128_t& x) { string S; is >> S; __int128_t ret = 0; int f = 1; if (S[0] == '-') f = -1; for (int i = 0; i < S.length(); i++) if ('0' <= S[i] && S[i] <= '9') ret = ret * 10 + S[i] - '0'; x = ret * f; return (is); } ostream &operator <<(ostream &os, __int128_t x) { ostream::sentry s(os); if (s) { __uint128_t tmp = x < 0 ? -x : x; char buffer[128]; char *d = end(buffer); do { --d; *d = "0123456789"[tmp % 10]; tmp /= 10; } while (tmp != 0); if (x < 0) { --d; *d = '-'; } int len = end(buffer) - d; if (os.rdbuf()->sputn(d, len) != len) { os.setstate(ios_base::badbit); } } return os; } __int128_t stoll(string &S) { __int128_t ret = 0; int f = 1; if (S[0] == '-') f = -1; for (int i = 0; i < S.length(); i++) if ('0' <= S[i] && S[i] <= '9') ret = ret * 10 + S[i] - '0'; return ret * f; } string to_string(__int128_t x) { string ret = ""; if (x < 0) { ret += "-"; x *= -1; } while (x) { ret += (char)('0' + x % 10); x /= 10; } reverse(ret.begin(), ret.end()); return ret; } string to_string(char c) { string s = {c}; return s; } vector<mint> _fac, _finv, _inv; void COMinit(int N) { _fac.resize(N + 1); _finv.resize(N + 1); _inv.resize(N + 1); _fac[0] = _fac[1] = 1; _finv[0] = _finv[1] = 1; _inv[1] = 1; for (int i = 2; i <= N; i++) { _fac[i] = _fac[i-1] * mint(i); _inv[i] = -_inv[MOD % i] * mint(MOD / i); _finv[i] = _finv[i - 1] * _inv[i]; } } mint FAC(int N) { if (N < 0) return 0; return _fac[N]; } mint COM(int N, int K) { if (N < K) return 0; if (N < 0 or K < 0) return 0; return _fac[N] * _finv[K] * _finv[N - K]; } mint PERM(int N, int K) { if (N < K) return 0; if (N < 0 or K < 0) return 0; return _fac[N] * _finv[N - K]; } mint NHK(int N, int K) { if (N == 0 && K == 0) return 1; return COM(N + K - 1, K); } #pragma endregion signed main() { int N; cin >> N; vector<int> A(N); gp_hash_table<int, int, custom_hash> mp; for (int i = 0; i < N; i++) { cin >> A[i]; mp[A[i]] = i; } vector<int> B(N); for (int i = 0; i < N; i++) { cin >> B[i]; B[i] = mp[B[i]]; } tree<int, null_type, less<int>, rb_tree_tag, tree_order_statistics_node_update> tr; int t = 0; for (int i = N - 1; i >= 0; i--) { tr.insert(B[i]); t += tr.order_of_key(B[i]); } cout << t << endl; }