/** * author: otera * created: 03.08.2024 21:23:28 **/ #include #include #include #include #include #include namespace atcoder { namespace internal { #ifndef _MSC_VER template using is_signed_int128 = typename std::conditional::value || std::is_same::value, std::true_type, std::false_type>::type; template using is_unsigned_int128 = typename std::conditional::value || std::is_same::value, std::true_type, std::false_type>::type; template using make_unsigned_int128 = typename std::conditional::value, __uint128_t, unsigned __int128>; template using is_integral = typename std::conditional::value || is_signed_int128::value || is_unsigned_int128::value, std::true_type, std::false_type>::type; template using is_signed_int = typename std::conditional<(is_integral::value && std::is_signed::value) || is_signed_int128::value, std::true_type, std::false_type>::type; template using is_unsigned_int = typename std::conditional<(is_integral::value && std::is_unsigned::value) || is_unsigned_int128::value, std::true_type, std::false_type>::type; template using to_unsigned = typename std::conditional< is_signed_int128::value, make_unsigned_int128, typename std::conditional::value, std::make_unsigned, std::common_type>::type>::type; #else template using is_integral = typename std::is_integral; template using is_signed_int = typename std::conditional::value && std::is_signed::value, std::true_type, std::false_type>::type; template using is_unsigned_int = typename std::conditional::value && std::is_unsigned::value, std::true_type, std::false_type>::type; template using to_unsigned = typename std::conditional::value, std::make_unsigned, std::common_type>::type; #endif template using is_signed_int_t = std::enable_if_t::value>; template using is_unsigned_int_t = std::enable_if_t::value>; template using to_unsigned_t = typename to_unsigned::type; } // namespace internal } // namespace atcoder namespace atcoder { template struct fenwick_tree { using U = internal::to_unsigned_t; public: fenwick_tree() : _n(0) {} explicit fenwick_tree(int n) : _n(n), data(n) {} void add(int p, T x) { assert(0 <= p && p < _n); p++; while (p <= _n) { data[p - 1] += U(x); p += p & -p; } } T sum(int l, int r) { assert(0 <= l && l <= r && r <= _n); return sum(r) - sum(l); } private: int _n; std::vector data; U sum(int r) { U s = 0; while (r > 0) { s += data[r - 1]; r -= r & -r; } return s; } }; } // namespace atcoder namespace otera { template struct compressor { public: static constexpr int absent = -1; compressor() : _xs(std::vector{}) {} compressor(const std::vector &vs) { add(vs); build(); } void add(const std::vector &vs) { std::copy(vs.begin(), vs.end(), std::back_inserter(_xs)); } void add(const T &e) { _xs.emplace_back(e); } void build() { std::sort(_xs.begin(), _xs.end()); _xs.erase(std::unique(_xs.begin(), _xs.end()), _xs.end()); } int size() const { return _xs.size(); } bool has_key(const T &e) const { return std::binary_search(_xs.begin(), _xs.end(), e); } std::vector get(const std::vector &vs) const { std::vector res; std::transform(vs.begin(), vs.end(), std::back_inserter(res), [&](const T &e) { return get(e); }); return res; } int get(const T &e, int default_value = absent) const { const int id = min_geq_index(e); return id != size() and _xs[id] == e ? id : default_value; } const T &operator[](int compressed_index) const { return _xs[compressed_index]; } T min_gt(const T &e, const T &default_value) const { auto it = std::upper_bound(_xs.begin(), _xs.end(), e); return it == _xs.end() ? default_value : *it; } T min_geq(const T &e, const T &default_value) const { auto it = std::lower_bound(_xs.begin(), _xs.end(), e); return it == _xs.end() ? default_value : *it; } T max_lt(const T &e, const T &default_value) const { auto it = std::lower_bound(_xs.begin(), _xs.end(), e); return it == _xs.begin() ? default_value : *std::prev(it); } T max_leq(const T &e, const T &default_value) const { auto it = std::upper_bound(_xs.begin(), _xs.end(), e); return it == _xs.begin() ? default_value : *std::prev(it); } int min_gt_index(const T &e) const { return std::upper_bound(_xs.begin(), _xs.end(), e) - _xs.begin(); } int min_geq_index(const T &e) const { return std::lower_bound(_xs.begin(), _xs.end(), e) - _xs.begin(); } int max_lt_index(const T &e) const { return int(std::lower_bound(_xs.begin(), _xs.end(), e)- _xs.begin()) - 1; } int max_leq_index(const T &e) const { return int(std::upper_bound(_xs.begin(), _xs.end(), e)- _xs.begin()) - 1; } private: std::vector _xs; }; } // namespace otera namespace otera { template long long inversion_number(std::vector a) { int _n = (int)a.size(); compressor ca(a); int _sz = (int)ca.size(); atcoder::fenwick_tree bit(_sz); long long res = 0; for(int i = 0; i < _n; ++ i) { int _pos = ca.get(a[i]); res += (long long)bit.sum(_pos + 1, _sz); bit.add(_pos, 1); } return res; } } // namespace otera namespace otera {} using namespace std; using namespace otera; struct io_setup { io_setup(int precision = 20) { ios::sync_with_stdio(false); cin.tie(0); cout << fixed << setprecision(precision); } } io_setup_ {}; #define int long long using ll = long long; using ld = long double; using ull = unsigned long long; using i128 = __int128_t; using u128 = __uint128_t; #define repa(i, n) for(int i = 0; i < (n); ++ i) #define repb(i, a, b) for(int i = (a); i < (b); ++ i) #define repc(i, a, b, c) for(int i = (a); i < (b); i += (c)) #define overload4(a, b, c, d, e, ...) e #define overload3(a, b, c, d, ...) d #define rep(...) overload4(__VA_ARGS__, repc, repb, repa)(__VA_ARGS__) #define rep1a(i, n) for(int i = 1; i <= (n); ++ i) #define rep1b(i, a, b) for(int i = (a); i <= (b); ++ i) #define rep1c(i, a, b, c) for(int i = (a); i <= (b); i += (c)) #define rep1(...) overload4(__VA_ARGS__, rep1c, rep1b, rep1a)(__VA_ARGS__) #define rev_repa(i, n) for(int i=(n)-1;i>=0;i--) #define rev_repb(i, a, b) for(int i=(b)-1;i>=(a);i--) #define rev_rep(...) overload3(__VA_ARGS__, rev_repb, rev_repa)(__VA_ARGS__) #define rev_rep1a(i, n) for(int i=(n);i>=1;i--) #define rev_rep1b(i, a, b) for(int i=(b);i>=(a);i--) #define rev_rep1(...) overload3(__VA_ARGS__, rev_rep1b, rev_rep1a)(__VA_ARGS__) #define for_subset(t, s) for(ll t = (s); t >= 0; t = (t == 0 ? -1 : (t - 1) & (s))) typedef pair P; typedef pair LP; #define pb push_back #define pf push_front #define ppb pop_back #define ppf pop_front #define eb emplace_back #define mkp make_pair #define mkt make_tuple #define fr first #define sc second #define all(c) std::begin(c), std::end(c) #define rall(c) std::rbegin(c), std::rend(c) #define lb(c, x) distance(std::begin(c), lower_bound(all(c), (x))) #define ub(c, x) distance(std::begin(c), upper_bound(all(c), (x))) #define Sort(a) sort(all(a)) #define Rev(a) reverse(all(a)) #define Uniq(a) sort(all(a));(a).erase(unique(all(a)),std::end(a)) #define si(c) (int)(c).size() inline ll popcnt(ull a){ return __builtin_popcountll(a); } #define kth_bit(x, k) ((x>>k)&1) #define unless(A) if(!(A)) #define elif else if ll modulo(ll n, ll d){ return (n % d + d) % d; }; ll intpow(ll a, ll b){ ll ans = 1; while(b){ if(b & 1) ans *= a; a *= a; b /= 2; } return ans; } ll intpow(ll a, ll b, ll m) {ll ans = 1; while(b){ if(b & 1) (ans *= a) %= m; (a *= a) %= m; b /= 2; } return ans; } template inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } template constexpr inline int pow_m1(T n) { return -(n & 1) | 1; } template constexpr inline T fld(const T x, const T y) { return (x ^ y) >= 0 ? x / y : (x - (y + pow_m1(y >= 0))) / y; } template constexpr inline T cld(const T x, const T y) { return (x ^ y) <= 0 ? x / y : (x + (y + pow_m1(y >= 0))) / y; } #define INT(...) int __VA_ARGS__;in(__VA_ARGS__) #define LL(...) ll __VA_ARGS__;in(__VA_ARGS__) #define ULL(...) ull __VA_ARGS__;in(__VA_ARGS__) #define STR(...) string __VA_ARGS__;in(__VA_ARGS__) #define CHR(...) char __VA_ARGS__;in(__VA_ARGS__) #define DBL(...) double __VA_ARGS__;in(__VA_ARGS__) #define LD(...) ld __VA_ARGS__;in(__VA_ARGS__) #define vec(type,name,...) vectorname(__VA_ARGS__) #define VEC(type,name,size) vectorname(size);in(name) #define vv(type,name,h,...) vector>name(h,vector(__VA_ARGS__)) #define VV(type,name,h,w) vector>name(h,vector(w));in(name) #define vvv(type,name,h,w,...) vector>>name(h,vector>(w,vector(__VA_ARGS__))) template using vc = vector; template using vvc = vector>; template using vvvc = vector>; template using vvvvc = vector>; template using pq = priority_queue; template using pqg = priority_queue, greater>; template using umap = unordered_map; template void scan(T& a){ cin >> a; } template void scan(vector& a){ for(auto&& i : a) scan(i); } void in(){} template void in(Head& head, Tail&... tail){ scan(head); in(tail...); } void print(){ cout << ' '; } template void print(const T& a){ cout << a; } template void print(const vector& a){ if(a.empty()) return; print(a[0]); for(auto i = a.begin(); ++i != a.end(); ){ cout << ' '; print(*i); } } int out(){ cout << '\n'; return 0; } template int out(const T& t){ print(t); cout << '\n'; return 0; } template int out(const Head& head, const Tail&... tail){ print(head); cout << ' '; out(tail...); return 0; } #define CHOOSE(a) CHOOSE2 a #define CHOOSE2(a0,a1,a2,a3,a4,x,...) x #define debug_1(x1) cout<<#x1<<": "< bit(n); deque ans; rep(i, n) { int pre = bit.sum(0, p[i]); int post = bit.sum(p[i] + 1, n); if(pre == post) { if(ans.front() < p[i] + 1) { ans.push_back(p[i] + 1); } else { ans.push_front(p[i] + 1); } } else if(pre < post) { ans.push_front(p[i] + 1); } else { ans.push_back(p[i] + 1); } bit.add(p[i], 1); } vc res; for(int x: ans) { res.eb(x); } int sum = inversion_number(res); out(sum); out(res); } signed main() { int testcase = 1; in(testcase); while(testcase--) solve(); return 0; }