/** * author: otera **/ #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 = 0; 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) assert((a) > (b));for(int i=(a);i>(b);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) assert((a) >= (b));for(int i=(a);i>=(b);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; } #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<<": "< a; vc cnt(n, 0); int res = m; rep(i, n) { cnt[i] = l[i]; res -= l[i]; } auto check = [&](int x) -> bool { int sum = 0; rep(i, n) { sum += min(r[i] - l[i], max(x - l[i], 0LL)); } return sum <= res; }; int ok = 0, ng = m; while(ng - ok > 1) { int mid = (ok + ng) / 2; if(check(mid)) ok = mid; else ng = mid; } debug(ok); dump(cnt); rep(i, n) { res -= min(r[i] - l[i], max(ok - l[i], 0LL)); cnt[i] += min(r[i] - l[i], max(ok - l[i], 0LL)); } dump(cnt); rep(i, n) { if(res > 0 and cnt[i] < ok + 1 and ok + 1 <= r[i]) { cnt[i] ++; -- res; } } dump(cnt); ll ans = 0; ll sum = 0; rev_rep(i, n) { ans += sum * cnt[i]; sum += cnt[i]; } debug(sum); out(ans); } signed main() { int testcase = 1; in(testcase); while(testcase--) solve(); return 0; }