/* Arkadiusz Wróbel * * Konkurs: Potyczki Algorytmiczne 2015, runda 1A * Zadanie: Siano */ #include <cstdio> #include <iostream> #include <algorithm> #include <cmath> #include <vector> #include <set> using namespace std; typedef long long LL; typedef pair<int, int> PII; typedef vector<int> VI; #define REP(I, N) for(int I=0; I<(N); ++I) #define FOR(I, M, N) for(int I=(M); I<=(N); ++I) #define FORD(I, M, N) for(int I=(M); I>=(N); --I) #define FOREACH(IT, CON) for(__typeof(CON.begin()) IT=CON.begin(); IT!=CON.end(); ++IT) #define ST first #define ND second #define MP make_pair #define PB push_back #define SIZE(CON) ((int)CON.size()) const int INF=1000000000; const LL INFLL=1000000000000000000LL; struct Koszenie { LL d, b; int odkad; Koszenie(LL d_, LL b_, int odkad_) : d(d_), b(b_), odkad(odkad_) {} Koszenie() : d(0), b(0), odkad(0) {} }; int n, m; LL t[1000000]; // numerowane od 1 (!!!) LL s[1000000]; // sumy prefiksowe LL d, b; vector<Koszenie> stos; LL wyn; inline LL sumuj(const int aa, const int bb) { // standardowy przedział return s[bb - 1] - s[aa - 1]; } inline int binsearch(int pocz, int kon, Koszenie kos) { while (pocz + 1 < kon) { int sred = (pocz + kon) / 2; if (kos.b + t[sred] * (d - kos.d) <= b) { pocz = sred; } else { kon = sred; } } return kon; } void brut() { LL lastD = 0; REP(nrKoszenia, m) { scanf("%lld%lld", &d, &b); LL tt = d - lastD; wyn = 0; FOR(i, 1, n) { s[i] += tt * t[i]; if (s[i] > b) { wyn += s[i] - b; s[i] = b; } } printf("%lld\n", wyn); lastD = d; } } int main() { scanf("%d%d", &n, &m); FOR(i, 1, n) { scanf("%lld", &t[i]); } if (n*m <= 2300000) { brut(); return 0; } sort(t+1, t+n+1); FOR(i, 1, n) { s[i] = t[i] + s[i - 1]; } stos.PB(Koszenie(0, 0, 0)); REP(nrKoszenia, m) { scanf("%lld%lld", &d, &b); wyn = 0; Koszenie kos; // bieżące koszenie int pocz; int kon = n + 1; while (true) { kos = stos.back(); pocz = kos.odkad; if (kos.b + t[pocz] * (d - kos.d) <= b) { break; } // printf(" %d %d\n", pocz, kon); /* Zliczanie */ wyn += (d - kos.d)*sumuj(pocz, kon) + (kos.b - b)*(kon - pocz); kon = pocz; stos.pop_back(); } // printf("binsearch %d %d\n", pocz, kon); pocz = binsearch(pocz, kon, kos); // Żeby kosić tylko to, co przekracza b // printf("x: %d\n", pocz); /* Zliczanie */ wyn += (d - kos.d)*sumuj(pocz, kon) + (kos.b - b)*(kon - pocz); if (pocz <= n) { stos.PB(Koszenie(d, b, pocz)); } printf("%lld\n", wyn); // REP(i, SIZE(stos)) {printf("%d: %lld %lld %d\n", i, stos[i].d, stos[i].b, stos[i].odkad);} } 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 129 130 131 132 133 134 135 136 137 | /* Arkadiusz Wróbel * * Konkurs: Potyczki Algorytmiczne 2015, runda 1A * Zadanie: Siano */ #include <cstdio> #include <iostream> #include <algorithm> #include <cmath> #include <vector> #include <set> using namespace std; typedef long long LL; typedef pair<int, int> PII; typedef vector<int> VI; #define REP(I, N) for(int I=0; I<(N); ++I) #define FOR(I, M, N) for(int I=(M); I<=(N); ++I) #define FORD(I, M, N) for(int I=(M); I>=(N); --I) #define FOREACH(IT, CON) for(__typeof(CON.begin()) IT=CON.begin(); IT!=CON.end(); ++IT) #define ST first #define ND second #define MP make_pair #define PB push_back #define SIZE(CON) ((int)CON.size()) const int INF=1000000000; const LL INFLL=1000000000000000000LL; struct Koszenie { LL d, b; int odkad; Koszenie(LL d_, LL b_, int odkad_) : d(d_), b(b_), odkad(odkad_) {} Koszenie() : d(0), b(0), odkad(0) {} }; int n, m; LL t[1000000]; // numerowane od 1 (!!!) LL s[1000000]; // sumy prefiksowe LL d, b; vector<Koszenie> stos; LL wyn; inline LL sumuj(const int aa, const int bb) { // standardowy przedział return s[bb - 1] - s[aa - 1]; } inline int binsearch(int pocz, int kon, Koszenie kos) { while (pocz + 1 < kon) { int sred = (pocz + kon) / 2; if (kos.b + t[sred] * (d - kos.d) <= b) { pocz = sred; } else { kon = sred; } } return kon; } void brut() { LL lastD = 0; REP(nrKoszenia, m) { scanf("%lld%lld", &d, &b); LL tt = d - lastD; wyn = 0; FOR(i, 1, n) { s[i] += tt * t[i]; if (s[i] > b) { wyn += s[i] - b; s[i] = b; } } printf("%lld\n", wyn); lastD = d; } } int main() { scanf("%d%d", &n, &m); FOR(i, 1, n) { scanf("%lld", &t[i]); } if (n*m <= 2300000) { brut(); return 0; } sort(t+1, t+n+1); FOR(i, 1, n) { s[i] = t[i] + s[i - 1]; } stos.PB(Koszenie(0, 0, 0)); REP(nrKoszenia, m) { scanf("%lld%lld", &d, &b); wyn = 0; Koszenie kos; // bieżące koszenie int pocz; int kon = n + 1; while (true) { kos = stos.back(); pocz = kos.odkad; if (kos.b + t[pocz] * (d - kos.d) <= b) { break; } // printf(" %d %d\n", pocz, kon); /* Zliczanie */ wyn += (d - kos.d)*sumuj(pocz, kon) + (kos.b - b)*(kon - pocz); kon = pocz; stos.pop_back(); } // printf("binsearch %d %d\n", pocz, kon); pocz = binsearch(pocz, kon, kos); // Żeby kosić tylko to, co przekracza b // printf("x: %d\n", pocz); /* Zliczanie */ wyn += (d - kos.d)*sumuj(pocz, kon) + (kos.b - b)*(kon - pocz); if (pocz <= n) { stos.PB(Koszenie(d, b, pocz)); } printf("%lld\n", wyn); // REP(i, SIZE(stos)) {printf("%d: %lld %lld %d\n", i, stos[i].d, stos[i].b, stos[i].odkad);} } return 0; } |