#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <fstream>
#include <sstream>
#include <set>
#include <map>
#include <vector>
#include <list>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <string>
#include <queue>
#include <bitset> //UWAGA - w czasie kompilacji musi byc znany rozmiar wektora - nie mozna go zmienic
#include <cassert>
#include <iomanip> //do setprecision
#include <ctime>
#include <complex>
using namespace std;
#define FOR(i,b,e) for(int i=(b);i<(e);++i)
#define FORQ(i,b,e) for(int i=(b);i<=(e);++i)
#define FORD(i,b,e) for(int i=(b)-1;i>=(e);--i)
#define REP(x, n) for(int x = 0; x < (n); ++x)
#define ST first
#define ND second
#define PB push_back
#define MP make_pair
#define LL long long
#define ULL unsigned LL
#define LD long double
const double pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342;
#define MR 500010
LL pre[MR];
int a[MR];
struct magic
{
LL base, prevD;
int ind;
};
vector <magic> V;
int main()
{
magic wartownik;
wartownik.base = wartownik.prevD = wartownik.ind = 0;
V.push_back(wartownik);
int n, m;
scanf("%d%d", &n, &m);
FORQ(i,1,n) scanf("%d", &a[i]);
sort(a+1, a+n+1);
FORQ(i,1,n) pre[i] = a[i] + pre[i-1];
wartownik.ind = n;
V.PB(wartownik);
REP(c,m)
{
LL b, d;
scanf("%lld%lld", &d, &b);
LL res = 0;
while(!V.empty())
{
if(V.back().base + (d-V.back().prevD)*a[V.back().ind] > b) // cos rosnie
{
assert(V.size() > 1);
// spr czy rosnie w calej grupie
int pos;
if(b <= V.back().base)
{
pos = V[V.size()-2].ind+1;
}
else
{
int pom = (b-V.back().base)/(d-V.back().prevD);
assert(V.back().base + (d-V.back().prevD)*pom <= b);
assert(V.back().base + (d-V.back().prevD)*(pom+1) > b);
pos = upper_bound(a+V[V.size()-2].ind+1, a+V.back().ind+1, pom) - a;
}
if(pos == V[V.size()-2].ind+1) // w calej grupie rosnie
{
assert(pos > 0);
pos--;
int ile = V.back().ind - pos;
res += ile*V.back().base + (pre[V.back().ind]-pre[pos])*(d-V.back().prevD) - ile*b;
V.pop_back();
}
else // rosnie tylko w czesci
{
pos--;
int ile = V.back().ind - pos;
res += ile*V.back().base + (pre[V.back().ind]-pre[pos])*(d-V.back().prevD) - ile*b;
V.back().ind = pos;
wartownik.base = b;
wartownik.prevD = d;
assert(wartownik.ind == n);
V.PB(wartownik);
break;
}
}
else //nic nie rosnie
{
if(V.back().ind < n) // dodaj nowa grupe
{
wartownik.base = b;
wartownik.prevD = d;
assert(wartownik.ind == n);
V.PB(wartownik);
}
break;
}
}
printf("%lld\n", res);
}
return 0;
}
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 123 124 125 126 127 128 | #include <cstdio> #include <cstdlib> #include <iostream> #include <fstream> #include <sstream> #include <set> #include <map> #include <vector> #include <list> #include <algorithm> #include <cstring> #include <cmath> #include <string> #include <queue> #include <bitset> //UWAGA - w czasie kompilacji musi byc znany rozmiar wektora - nie mozna go zmienic #include <cassert> #include <iomanip> //do setprecision #include <ctime> #include <complex> using namespace std; #define FOR(i,b,e) for(int i=(b);i<(e);++i) #define FORQ(i,b,e) for(int i=(b);i<=(e);++i) #define FORD(i,b,e) for(int i=(b)-1;i>=(e);--i) #define REP(x, n) for(int x = 0; x < (n); ++x) #define ST first #define ND second #define PB push_back #define MP make_pair #define LL long long #define ULL unsigned LL #define LD long double const double pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342; #define MR 500010 LL pre[MR]; int a[MR]; struct magic { LL base, prevD; int ind; }; vector <magic> V; int main() { magic wartownik; wartownik.base = wartownik.prevD = wartownik.ind = 0; V.push_back(wartownik); int n, m; scanf("%d%d", &n, &m); FORQ(i,1,n) scanf("%d", &a[i]); sort(a+1, a+n+1); FORQ(i,1,n) pre[i] = a[i] + pre[i-1]; wartownik.ind = n; V.PB(wartownik); REP(c,m) { LL b, d; scanf("%lld%lld", &d, &b); LL res = 0; while(!V.empty()) { if(V.back().base + (d-V.back().prevD)*a[V.back().ind] > b) // cos rosnie { assert(V.size() > 1); // spr czy rosnie w calej grupie int pos; if(b <= V.back().base) { pos = V[V.size()-2].ind+1; } else { int pom = (b-V.back().base)/(d-V.back().prevD); assert(V.back().base + (d-V.back().prevD)*pom <= b); assert(V.back().base + (d-V.back().prevD)*(pom+1) > b); pos = upper_bound(a+V[V.size()-2].ind+1, a+V.back().ind+1, pom) - a; } if(pos == V[V.size()-2].ind+1) // w calej grupie rosnie { assert(pos > 0); pos--; int ile = V.back().ind - pos; res += ile*V.back().base + (pre[V.back().ind]-pre[pos])*(d-V.back().prevD) - ile*b; V.pop_back(); } else // rosnie tylko w czesci { pos--; int ile = V.back().ind - pos; res += ile*V.back().base + (pre[V.back().ind]-pre[pos])*(d-V.back().prevD) - ile*b; V.back().ind = pos; wartownik.base = b; wartownik.prevD = d; assert(wartownik.ind == n); V.PB(wartownik); break; } } else //nic nie rosnie { if(V.back().ind < n) // dodaj nowa grupe { wartownik.base = b; wartownik.prevD = d; assert(wartownik.ind == n); V.PB(wartownik); } break; } } printf("%lld\n", res); } return 0; } |
polski