#include <bits/stdc++.h> using namespace std; #define REP(i,a,b) for (int i = (a); i <= (b); ++i) #define REPD(i,a,b) for (int i = (a); i >= (b); --i) #define FORI(i,n) REP(i,1,n) #define FOR(i,n) REP(i,0,int(n)-1) #define mp make_pair #define pb push_back #define pii pair<int,int> #define vi vector<int> #define ll long long #define SZ(x) int((x).size()) #define DBG(v) cerr << #v << " = " << (v) << endl; #define FOREACH(i,t) for (typeof(t.begin()) i=t.begin(); i!=t.end(); i++) #define fi first #define se second const int maxn = 1000100; int a[maxn]; long long pref[maxn]; long long d[maxn], b[maxn]; int last[maxn], pos[maxn]; long long getsum(int low, int high) { return pref[high]-pref[low]; } long long add(int prv, int cur, int left, int right) { long long res = (b[prv]-b[cur])*(right-left) + (d[cur]-d[prv])*getsum(left,right); //printf("dodaje %lld do wyniku\n", res); return res; } int main () { int n,m,it=0; scanf("%d%d", &n, &m); FOR(i,n) scanf("%d", &a[i]); sort(a,a+n); pref[0] = 0; FORI(i,n) pref[i] = pref[i-1] + a[i-1]; last[it] = 0; pos[it] = 0; d[0] = 0; b[0] = 0; //m=10; //printf("najwyzszy ma %d\n", a[n-1]); FORI(i,m) { scanf("%lld%lld", &d[i], &b[i]); //printf("dzien %lld, wysokosc %lld\n", d[i], b[i]); long long res = 0; int right = n; while (it >= 0) { //printf("wyciagam przedzial [%d %d] z dnia %lld (wysokosc %lld)\n", pos[it], right, d[last[it]], b[last[it]]); long long left = b[last[it]] + (d[i]-d[last[it]])*a[pos[it]]; //printf("lewy rosnie o %d\n", a[pos[it]]); if (left >= b[i]) { //printf("caly bedzie sciety (%lld vs %lld)\n", left, b[i]); res += add(last[it], i, pos[it], right); right = pos[it]; it--; } else { //printf("nie bedzie sciety (%lld vs %lld)\n", left, b[i]); break; } } if (it >= 0) { //printf("dziele przedzial [%d %d] z dnia %lld (wysokosc %lld, i=%d)\n", pos[it], right, d[last[it]], b[last[it]], last[it]); long long thr = b[i] - b[last[it]]; long long days = d[i] - d[last[it]]; //printf("szukam granicy dla %lld w %lld dni (%lld na dzien)\n", thr, days, (thr+days-1)/days); int q = lower_bound(a+pos[it], a+right, (thr+days)/days) - a; //printf("punkt podzialu %d\n", q); if (q < n) { res += add(last[it], i, q, right); it++; last[it] = i; pos[it] = q; } } else { //printf("skonczyly sie przedzialy\n"); it++; last[it] = i; pos[it] = 0; } printf("%lld\n", res); } }
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 | #include <bits/stdc++.h> using namespace std; #define REP(i,a,b) for (int i = (a); i <= (b); ++i) #define REPD(i,a,b) for (int i = (a); i >= (b); --i) #define FORI(i,n) REP(i,1,n) #define FOR(i,n) REP(i,0,int(n)-1) #define mp make_pair #define pb push_back #define pii pair<int,int> #define vi vector<int> #define ll long long #define SZ(x) int((x).size()) #define DBG(v) cerr << #v << " = " << (v) << endl; #define FOREACH(i,t) for (typeof(t.begin()) i=t.begin(); i!=t.end(); i++) #define fi first #define se second const int maxn = 1000100; int a[maxn]; long long pref[maxn]; long long d[maxn], b[maxn]; int last[maxn], pos[maxn]; long long getsum(int low, int high) { return pref[high]-pref[low]; } long long add(int prv, int cur, int left, int right) { long long res = (b[prv]-b[cur])*(right-left) + (d[cur]-d[prv])*getsum(left,right); //printf("dodaje %lld do wyniku\n", res); return res; } int main () { int n,m,it=0; scanf("%d%d", &n, &m); FOR(i,n) scanf("%d", &a[i]); sort(a,a+n); pref[0] = 0; FORI(i,n) pref[i] = pref[i-1] + a[i-1]; last[it] = 0; pos[it] = 0; d[0] = 0; b[0] = 0; //m=10; //printf("najwyzszy ma %d\n", a[n-1]); FORI(i,m) { scanf("%lld%lld", &d[i], &b[i]); //printf("dzien %lld, wysokosc %lld\n", d[i], b[i]); long long res = 0; int right = n; while (it >= 0) { //printf("wyciagam przedzial [%d %d] z dnia %lld (wysokosc %lld)\n", pos[it], right, d[last[it]], b[last[it]]); long long left = b[last[it]] + (d[i]-d[last[it]])*a[pos[it]]; //printf("lewy rosnie o %d\n", a[pos[it]]); if (left >= b[i]) { //printf("caly bedzie sciety (%lld vs %lld)\n", left, b[i]); res += add(last[it], i, pos[it], right); right = pos[it]; it--; } else { //printf("nie bedzie sciety (%lld vs %lld)\n", left, b[i]); break; } } if (it >= 0) { //printf("dziele przedzial [%d %d] z dnia %lld (wysokosc %lld, i=%d)\n", pos[it], right, d[last[it]], b[last[it]], last[it]); long long thr = b[i] - b[last[it]]; long long days = d[i] - d[last[it]]; //printf("szukam granicy dla %lld w %lld dni (%lld na dzien)\n", thr, days, (thr+days-1)/days); int q = lower_bound(a+pos[it], a+right, (thr+days)/days) - a; //printf("punkt podzialu %d\n", q); if (q < n) { res += add(last[it], i, q, right); it++; last[it] = i; pos[it] = q; } } else { //printf("skonczyly sie przedzialy\n"); it++; last[it] = i; pos[it] = 0; } printf("%lld\n", res); } } |