oj.uz

Submission #125872

# Submission time Handle Problem Language Result Execution time Memory
125872 2019-07-06T13:20:55 Z eriksuenderhauf Nuclearia (CEOI15_nuclearia) C++11
50 / 100
 865 ms 534508 KB 
#pragma GCC optimize("O3")
#include <bits/stdc++.h>
using namespace std;
const int MAXN = 2500005;
typedef long long ll;
vector<vector<ll>> d1, d2;
vector<vector<ll>> offC, slC, offSlC, slSlC;
ll sm[MAXN], sl[MAXN], off[MAXN];
int w, h, n, q;

int get(int y, int x) {
	if (y == 0 || x == 0)
		return 0;
	return (y-1) * w + x-1;
}
 
int main() {
	scanf("%d %d %d", &w, &h, &n);
	d1.resize(h+3), d2.resize(h+3), offC.resize(h+3), slC.resize(h+3);
	offSlC.resize(h+3), slSlC.resize(h+3);
	for (int i = 0; i <= h+2; i++) {
		d1[i].resize(w+3), d2[i].resize(w+3);
		offC[i].resize(w+3), slC[i].resize(w+3);
		offSlC[i].resize(w+3), slSlC[i].resize(w+3);
	}
	for (int i = 0; i < n; i++) {
		ll x, y, a, b;
		scanf("%lld %lld %lld %lld", &x, &y, &a, &b);
		ll l = max(1ll, x - a / b), r = min((ll)w, x + a / b);
		ll c = max(1ll, y - a / b), d = min((ll)h, y + a / b);
		ll loc = x - abs(y-c), hic = x + abs(y-c);
		ll lod = x - abs(y-d), hid = x + abs(y-d);
		if (loc < 1ll) {
			d1[y-x+1][2] -= b;
			offSlC[c][2] -= b;
			offSlC[c][l] -= b;
			offSlC[c][l+1] += b;
			offSlC[y-x+1][2] += b;
			offSlC[y-x+1][l] += b;
			offSlC[y-x+1][l+1] -= b;
			slSlC[y-x+1][l] -= b;
			slSlC[y-x+1][l+1] += b;
		} else {
			d1[c][loc+1] -= b;
			slSlC[c][l] -= loc * b;
			slSlC[c][l+1] += loc * b;
			slSlC[c+1][l] += loc * b - b;
			slSlC[c+1][l+1] -= loc * b - b;
		}
		if (lod < 1ll) {
			d2[x+y-1][2] -= b;
			offSlC[x+y][2] -= b;
			offSlC[x+y][l] -= b;
			offSlC[x+y][l+1] += b;
			offSlC[d+1][2] += b;
			offSlC[d+1][l] += b;
			offSlC[d+1][l+1] -= b;
			slSlC[x+y+1][l] -= b;
			slSlC[x+y+1][l+1] += b;
		} else {
			d2[d][lod+1] -= b;
			slSlC[d+1][l] += lod * b - b;
			slSlC[d+1][l+1] -= lod * b - b;
			slSlC[d+2][l] -= lod * b;
			slSlC[d+2][l+1] += lod * b;
		}
		if (hid > (ll)w) {
			offSlC[w+y-x+1][r+1] -= b * w;
			offSlC[w+y-x+1][r+2] += b * w;
			offSlC[d+1][r+1] += b * w;
			offSlC[d+1][r+2] -= b * w;
			slSlC[w+y-x+1][r+1] += (ll)w * b + b;
			slSlC[w+y-x+1][r+2] -= (ll)w * b + b;
			slSlC[w+y-x+2][r+1] -= (ll)w * b;
			slSlC[w+y-x+2][r+2] += (ll)w * b;
		} else {
			d1[d+1][hid+2] += b;
			slSlC[d+1][r+1] += hid * b + b;
			slSlC[d+1][r+2] -= hid * b + b;
			slSlC[d+2][r+1] -= hid * b;
			slSlC[d+2][r+2] += hid * b;
		}
		if (hic > (ll)w) {
			offSlC[c][r+1] -= b * w;
			offSlC[c][r+2] += b * w;
			offSlC[x+y-w][r+1] += b * w;
			offSlC[x+y-w][r+2] -= b * w;
			slSlC[x+y-w][r+1] -= (ll)w * b;
			slSlC[x+y-w][r+2] += (ll)w * b;
			slSlC[x+y-w+1][r+1] += (ll)w*b + b;
			slSlC[x+y-w+1][r+2] -= (ll)w*b + b;
		} else {
			d2[c-1][hic+2] += b;
			slSlC[c][r+1] -= hic * b;
			slSlC[c][r+2] += hic * b;
			slSlC[c+1][r+1] += hic * b + b;
			slSlC[c+1][r+2] -= hic * b + b;
		}
		slSlC[y+1][l] += 2 * b;
		slSlC[y+1][l+1] -= 2 * b;
		slSlC[y+1][r+1] -= 2 * b;
		slSlC[y+1][r+2] += 2 * b;
		offC[c][l] += a;
		offC[c][r+1] -= a;
		offC[d+1][l] -= a;
		offC[d+1][r+1] += a;
		slC[c][l] -= (hic - x) * b;
		slC[c][r+1] += (hic - x) * b;
		slC[c+1][l] += (hic - x) * b + b;
		slC[c+1][r+1] -= (hic - x) * b + b;
		slC[y+1][l] -= 2 * b;
		slC[y+1][r+1] += 2 * b;
		slC[d+1][l] += (hid - x) * b + b;
		slC[d+1][r+1] -= (hid - x) * b + b;
		slC[d+2][l] -= (hid - x) * b;
		slC[d+2][r+1] += (hid - x) * b;
		offSlC[c][l] += b * l;
		offSlC[c][l+1] -= b * (l-1);
		offSlC[c][r+1] += b * (r+1);
		offSlC[c][r+2] -= b * r;
		offSlC[d+1][l] -= b * l;
		offSlC[d+1][l+1] += b * (l-1);
		offSlC[d+1][r+1] -= b * (r+1);
		offSlC[d+1][r+2] += b * r;
	}
	for (int i = 1; i <= h; i++) {
		for (int j = 1; j <= w; j++) {
			slC[i][j] += slC[i-1][j];
			slSlC[i][j] += slSlC[i-1][j];
			d1[i][j] += d1[i-1][j-1];
			offC[i][j] += offC[i-1][j];
			offSlC[i][j] += offSlC[i-1][j];
		}
	}
	for (int i = h-1; i > 0; i--)
		for (int j = 1; j <= w; j++)
			d2[i][j] += d2[i+1][j-1];
	for (int i = 1; i <= h; i++) {
		for (int j = 1; j <= w; j++) {
			slC[i][j] += slC[i-1][j];
			slSlC[i][j] += slSlC[i-1][j];
		}
		for (int j = 1; j <= w; j++) {
			sl[j] = d1[i][j] + d2[i][j] + offSlC[i][j] + slSlC[i][j];
			sl[j] += sl[j-1];
			off[j] = offC[i][j] + slC[i][j];
			off[j] += off[j-1];
		}
		for (int j = 1; j <= w; j++) {
			sl[j] += sl[j-1];
			sm[get(i,j)] = sm[get(i-1,j)] + sm[get(i,j-1)] - sm[get(i-1,j-1)] + sl[j] + off[j];
		}
		for (int j = 1; j <= w; j++)
			sl[j] = off[j] = 0;
	}
	scanf("%d", &q);
	while (q--) {
		ll x1, y1, x2, y2;
		scanf("%lld %lld %lld %lld", &x1, &y1, &x2, &y2);
		ll ans = sm[get(y2,x2)] - sm[get(y1-1,x2)] - sm[get(y2,x1-1)] + sm[get(y1-1,x1-1)];
		x1 = (x2-x1+1) * (y2-y1+1);
		if ((ans % x1) * 2 >= x1)
			printf("%lld\n", ans / x1 + 1);
		else
			printf("%lld\n", ans / x1);
	}
    return 0;
}

Compilation message

nuclearia.cpp: In function 'int main()':
nuclearia.cpp:18:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
  scanf("%d %d %d", &w, &h, &n);
  ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~
nuclearia.cpp:28:8: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   scanf("%lld %lld %lld %lld", &x, &y, &a, &b);
   ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
nuclearia.cpp:156:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
  scanf("%d", &q);
  ~~~~~^~~~~~~~~~
nuclearia.cpp:159:8: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   scanf("%lld %lld %lld %lld", &x1, &y1, &x2, &y2);
   ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Subtask #1
0 / 3.0

# Verdict Execution time Memory Grader output
1 Correct 538 ms 528712 KB Output is correct
2 Incorrect 100 ms 4668 KB Output isn't correct
3 Halted 0 ms 0 KB -

Subtask #2
0 / 2.0

# Verdict Execution time Memory Grader output
1 Correct 504 ms 528816 KB Output is correct
2 Incorrect 100 ms 4728 KB Output isn't correct

Subtask #3
0 / 3.0

# Verdict Execution time Memory Grader output
1 Correct 169 ms 138232 KB Output is correct
2 Incorrect 98 ms 2936 KB Output isn't correct

Subtask #4
0 / 2.0

# Verdict Execution time Memory Grader output
1 Correct 328 ms 178492 KB Output is correct
2 Incorrect 99 ms 2936 KB Output isn't correct

Subtask #5
0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 666 ms 534508 KB Output is correct
2 Incorrect 106 ms 5496 KB Output isn't correct

Subtask #6
0 / 4.0

# Verdict Execution time Memory Grader output
1 Incorrect 345 ms 214264 KB Output isn't correct
2 Halted 0 ms 0 KB -

Subtask #7
0 / 6.0

# Verdict Execution time Memory Grader output
1 Incorrect 315 ms 141388 KB Output isn't correct
2 Halted 0 ms 0 KB -

Subtask #8
0 / 4.0

# Verdict Execution time Memory Grader output
1 Correct 282 ms 135896 KB Output is correct
2 Incorrect 98 ms 2680 KB Output isn't correct

Subtask #9
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 865 ms 532780 KB Output is correct

Subtask #10
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 864 ms 532700 KB Output is correct

Subtask #11
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 564 ms 140260 KB Output is correct
2 Correct 573 ms 140296 KB Output is correct

Subtask #12
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 566 ms 140296 KB Output is correct
2 Correct 569 ms 141176 KB Output is correct

Subtask #13
0 / 12.0

# Verdict Execution time Memory Grader output
1 Correct 667 ms 147220 KB Output is correct
2 Incorrect 455 ms 141568 KB Output isn't correct

Subtask #14
0 / 8.0

# Verdict Execution time Memory Grader output
1 Incorrect 451 ms 140664 KB Output isn't correct
2 Halted 0 ms 0 KB -