#pragma region Macros #include #define rep(i, n) for(int(i) = 0; (i) < (n); (i)++) #define rrep(i, n) for(int(i) = (n)-1; (i) >= 0; (i)--) #define FOR(i, m, n) for(int(i) = (m); (i) < (n); (i)++) #define ROF(i, m, n) for(int(i) = (n)-1; (i) >= (m); (i)--) #define ALL(v) (v).begin(), (v).end() #define LLA(v) (v).rbegin(), (v).rend() #define SZ(v) (int)(v).size() #define INT(...) \ int __VA_ARGS__; \ read(__VA_ARGS__) #define LL(...) \ ll __VA_ARGS__; \ read(__VA_ARGS__) #define DOUBLE(...) \ double __VA_ARGS__; \ read(__VA_ARGS__) #define CHAR(...) \ char __VA_ARGS__; \ read(__VA_ARGS__) #define STRING(...) \ string __VA_ARGS__; \ read(__VA_ARGS__) #define VEC(type, name, size) \ vector name(size); \ read(name) #define VEC2(type, name, height, width) \ vector> name(height, vector(width)); \ read(name) #define DVEC(type, name1, name2, size) \ vector name1(size), name2(size); \ read(name1, name2) #define TVEC(type, name1, name2, name3, size) \ vector name1(size), name2(size), name3(size); \ read(name1, name2, name3) using namespace std; using ll = long long; using pii = pair; using pll = pair; const int INF = 1 << 30; const ll LINF = 1LL << 60; const int MOD = 1e9 + 7; const char newl = '\n'; const int dx[] = {1, 0, -1, 0}; const int dy[] = {0, 1, 0, -1}; template inline bool between(T x, T l, T r) { return l <= x && x < r; } template inline vector make_vec(size_t a, T val) { return vector(a, val); } template inline auto make_vec(size_t a, Ts... ts) { return vector(a, make_vec(ts...)); } void read() {} template inline void read(T &a) { cin >> a; } template inline void read(pair &p) { read(p.first), read(p.second); } template inline void read(vector &v) { for(auto &&a : v) read(a); } template inline void read(vector &a, vector &b) { for(int i = 0; i < a.size(); i++) { read(a[i]); read(b[i]); } } template inline void read(vector &a, vector &b, vector &c) { for(int i = 0; i < a.size(); i++) { read(a[i]); read(b[i]); read(c[i]); } } template inline void read(Head &head, Tail &...tail) { read(head), read(tail...); } template void _write(const T &a) { cout << a; } template void _write(const std::pair &a) { cout << a.first << ' ' << a.second; } template void write(const T &a) { _write(a); cout << newl; } template void write(const vector &a) { for(int i = 0; i < a.size(); i++) { _write(a[i]); cout << (i + 1 == a.size() ? newl : ' '); } } template void write(const Head &head, const Tail &...tail) { _write(head); cout << ' '; write(tail...); } template void writel(const T &a) { cout << a << '\n'; } template void writel(const vector &a) { for(int i = 0; i < a.size(); i++) { _write(a[i]); cout << newl; } } template void writel(const Head &head, const Tail &...tail) { _write(head); cout << newl; write(tail...); } template void _debug(const T &a) { cerr << a; } template void _debug(const std::pair &a) { cerr << a.first << ' ' << a.second; } template void debug(const T &a) { _debug(a); cerr << newl; } template void debug(const vector &a) { for(int i = 0; i < a.size(); i++) { _debug(a[i]); cerr << (i + 1 == a.size() ? newl : ' '); } } template void debug(const Head &head, const Tail &...tail) { _debug(head); cerr << ' '; debug(tail...); } template void debugl(const T &a) { cerr << a << '\n'; } template void debugl(const vector &a) { for(int i = 0; i < a.size(); i++) { _debug(a[i]); cerr << newl; } } template void debugl(const Head &head, const Tail &...tail) { _debug(head); cerr << newl; debug(tail...); } template auto sum(const vector &a) { return accumulate(ALL(a), T(0)); } template auto min(const vector &a) { return *min_element(ALL(a)); } template auto max(const vector &a) { return *max_element(ALL(a)); } template void msort(vector &a, vector &b) { assert(a.size() == b.size()); vector> ab(a.size()); for(int i = 0; i < a.size(); i++) ab[i] = {a[i], b[i]}; sort(ALL(ab)); for(int i = 0; i < a.size(); i++) { a[i] = ab[i].first; b[i] = ab[i].second; } } template void msort(vector &a, vector &b, vector &c) { assert(a.size() == b.size() && b.size() == c.size()); vector> abc(a.size()); for(int i = 0; i < a.size(); i++) abc[i] = {a[i], b[i], c[i]}; sort(ALL(abc)); for(int i = 0; i < a.size(); i++) { a[i] = get<0>(abc[i]); b[i] = get<1>(abc[i]); c[i] = get<2>(abc[i]); } } template inline bool chmax(T &a, U b) { if(a < b) { a = b; return 1; } return 0; } template inline bool chmin(T &a, U b) { if(a > b) { a = b; return 1; } return 0; } int digit(ll a) { ll ret = 0; while(a && ++ret) a /= 10; return ret; } int digit_sum(ll a) { ll ret = 0; while(a) { ret += a % 10; a /= 10; } return ret; } ll llpow(ll a, ll n) { ll ret = 1; rep(i, n) ret *= a; return ret; } inline int bsf(int v) { return __builtin_ctz(v); } // 最下位の1が下から何番目か inline int bsf(ll v) { return __builtin_ctzll(v); } inline int bsr(int v) { return 31 - __builtin_clz(v); } // 最上位の1が下から何番目か inline int bsr(ll v) { return 63 - __builtin_clzll(v); } inline int lsb(int v) { return v & -v; } // 最下位の1だけ残す inline ll lsb(ll v) { return v & -v; } inline int msb(int v) { return 1 << bsr(v); } // 最上位の1だけ残す inline ll msb(ll v) { return 1LL << bsr(v); } struct IO { IO() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(10); } } io; #pragma endregion /** * @brief 0-indexed (Internally 1-indexed) Fenwick Tree * * @tparam T Universal set. */ template struct FenwickTree { int n; std::vector dat; FenwickTree(int n) : n(++n), dat(n, 0) {} void add(int i, T x) { for(++i; i < n; i += i & -i) dat[i] = dat[i] + x; } /** * @brief Sum values in [0,i] * * @param i Right end of interval. Note that closed-interval [0,i]. * @return T Summation values in [0,i]. */ T sum(int i) { T res = 0; for(++i; i > 0; i -= i & -i) { res = res + dat[i]; } return res; } /** * @brief Sum values in [l,r] * * @param l Left end of interval. * @param r Right end of interval. * @return T Summation values in [l,r] */ T sum(int l, int r) { return sum(r) - sum(l - 1); } /** * @brief Sum all values. * * @return T Summation values. */ T sumAll() { return sum(n - 1); } /** * @brief If regard as histogram of elements, get k-th element. * * @param k Get k-th element. * @return int Index to which k-th element belongs. */ int get(T k) { ++k; int res = 0; int N = 1; while(N < n) N *= 2; for(int i = N / 2; i > 0; i /= 2) { if(res + i < n && dat[res + i] < k) { k = k - dat[res + i]; res = res + i; } } return res + 1; } }; void solve() { INT(n); VEC(ll, a, n); VEC(ll, b, n); VEC(ll, c, n); sort(ALL(a)); sort(ALL(b)); sort(ALL(c)); vector rn(3001), gn(3001), bn(3001); rep(i, n) { rn[a[i]]++; gn[b[i]]++; bn[c[i]]++; } vector sum(3002); rep(i, 3001) { sum[i + 1] = sum[i] + rn[i]; } FenwickTree bit(3001); rep(i, n) bit.add(a[i], 1); ll res = 0; FOR(i, 1, 3001) FOR(j, 1, 3001) { // res += gn[i] * bn[j] * (sum[min(i + j, 3001)] - sum[max(i, j)]); res += gn[i] * bn[j] * bit.sum(max(i, j), min(3000, i + j - 1)); } write(res); } int main() { int t = 1; // cin >> t; while(t--) { solve(); } }