#include <bits/stdc++.h> using namespace std; #define sim template < class c #define ris return * this #define dor > debug & operator << #define eni(x) sim > typename \ enable_if<sizeof dud<c>(0) x 1, debug&>::type operator<<(c i) { sim > struct rge {c b, e; }; sim > rge<c> range(c i, c j) { return rge<c>{i, j}; } sim > auto dud(c* x) -> decltype(cerr << *x, 0); sim > char dud(...); struct debug { #ifdef LOCAL ~debug() { cerr << endl; } eni(!=) cerr << boolalpha << i; ris; } eni(==) ris << range(begin(i), end(i)); } sim, class b dor(pair < b, c > d) { ris << "(" << d.first << ", " << d.second << ")"; } sim dor(rge<c> d) { *this << "["; for (auto it = d.b; it != d.e; ++it) *this << ", " + 2 * (it == d.b) << *it; ris << "]"; } #else sim dor(const c&) { ris; } #endif }; #define imie(...) " [" << #__VA_ARGS__ ": " << (__VA_ARGS__) << "] " typedef long long ll; const int nax = 1e6 + 5; int g[nax]; vector<int> inv[nax]; ll sum[nax], leftmost[nax]; ll would[nax]; vector<pair<int,int>> events[nax]; void addEvent(int a, int b, ll when) { assert(when >= 0); if(when < nax) events[when].push_back({a, b}); } ll sufit(ll a, ll b) { return (a + b - 1) / b; } int main() { int n, q; scanf("%d%d", &n, &q); vector<ll> in(n+1, 0LL); ++n; set<int> starts; for(int i = 0; i < n; ++i) { if(i) scanf("%lld", &in[i]); leftmost[i] = in[i]; g[i] = i; inv[i] = {i}; if(i) addEvent(i - 1, i, in[i] - in[i-1]); starts.insert(i); } ll A = 0, B = 0; auto pairs = [&] (ll x) { return x * (x - 1) / 2; }; auto add = [&] (int x, int mul) { assert(g[x] == x); A += mul * pairs(inv[x].size()); B += mul * (leftmost[x] * inv[x].size() - sum[x]); }; for(int single = 0; single < nax; ++single) { for(int i = 0; i < (int) events[single].size(); ++i) { pair<int,int> p = events[single][i]; int a = g[p.first], b = g[p.second]; if(a == b) continue; debug() << "start"; //~ debug() << imie(a) imie(b) imie(starts); add(a, -1); add(b, -1); if(inv[a].size() > inv[b].size()) swap(a, b); for(int x : inv[a]) { inv[b].push_back(x); g[x] = b; } inv[a].clear(); leftmost[b] = min(leftmost[b], leftmost[a]); sum[b] += sum[a]; add(b, 1); starts.erase(a); debug() << "a"; auto it = starts.find(b); assert(it != starts.end()); if(it != starts.begin()) { auto it2 = prev(it); addEvent(*it2, *it, sufit(leftmost[*it] - leftmost[*it2], inv[*it2].size())); } auto it2 = next(it); if(it2 != starts.end()) addEvent(*it, *it2, sufit(leftmost[*it2] - leftmost[*it], inv[*it].size())); debug() << "b"; } //~ debug() << imie(A) imie(B); would[single] = A * single + B; } while(q--) { int single; scanf("%d", &single); printf("%lld\n", would[single]); } }
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 | #include <bits/stdc++.h> using namespace std; #define sim template < class c #define ris return * this #define dor > debug & operator << #define eni(x) sim > typename \ enable_if<sizeof dud<c>(0) x 1, debug&>::type operator<<(c i) { sim > struct rge {c b, e; }; sim > rge<c> range(c i, c j) { return rge<c>{i, j}; } sim > auto dud(c* x) -> decltype(cerr << *x, 0); sim > char dud(...); struct debug { #ifdef LOCAL ~debug() { cerr << endl; } eni(!=) cerr << boolalpha << i; ris; } eni(==) ris << range(begin(i), end(i)); } sim, class b dor(pair < b, c > d) { ris << "(" << d.first << ", " << d.second << ")"; } sim dor(rge<c> d) { *this << "["; for (auto it = d.b; it != d.e; ++it) *this << ", " + 2 * (it == d.b) << *it; ris << "]"; } #else sim dor(const c&) { ris; } #endif }; #define imie(...) " [" << #__VA_ARGS__ ": " << (__VA_ARGS__) << "] " typedef long long ll; const int nax = 1e6 + 5; int g[nax]; vector<int> inv[nax]; ll sum[nax], leftmost[nax]; ll would[nax]; vector<pair<int,int>> events[nax]; void addEvent(int a, int b, ll when) { assert(when >= 0); if(when < nax) events[when].push_back({a, b}); } ll sufit(ll a, ll b) { return (a + b - 1) / b; } int main() { int n, q; scanf("%d%d", &n, &q); vector<ll> in(n+1, 0LL); ++n; set<int> starts; for(int i = 0; i < n; ++i) { if(i) scanf("%lld", &in[i]); leftmost[i] = in[i]; g[i] = i; inv[i] = {i}; if(i) addEvent(i - 1, i, in[i] - in[i-1]); starts.insert(i); } ll A = 0, B = 0; auto pairs = [&] (ll x) { return x * (x - 1) / 2; }; auto add = [&] (int x, int mul) { assert(g[x] == x); A += mul * pairs(inv[x].size()); B += mul * (leftmost[x] * inv[x].size() - sum[x]); }; for(int single = 0; single < nax; ++single) { for(int i = 0; i < (int) events[single].size(); ++i) { pair<int,int> p = events[single][i]; int a = g[p.first], b = g[p.second]; if(a == b) continue; debug() << "start"; //~ debug() << imie(a) imie(b) imie(starts); add(a, -1); add(b, -1); if(inv[a].size() > inv[b].size()) swap(a, b); for(int x : inv[a]) { inv[b].push_back(x); g[x] = b; } inv[a].clear(); leftmost[b] = min(leftmost[b], leftmost[a]); sum[b] += sum[a]; add(b, 1); starts.erase(a); debug() << "a"; auto it = starts.find(b); assert(it != starts.end()); if(it != starts.begin()) { auto it2 = prev(it); addEvent(*it2, *it, sufit(leftmost[*it] - leftmost[*it2], inv[*it2].size())); } auto it2 = next(it); if(it2 != starts.end()) addEvent(*it, *it2, sufit(leftmost[*it2] - leftmost[*it], inv[*it].size())); debug() << "b"; } //~ debug() << imie(A) imie(B); would[single] = A * single + B; } while(q--) { int single; scanf("%d", &single); printf("%lld\n", would[single]); } } |