oj.uz

Submission #127786

# Submission time Handle Problem Language Result Execution time Memory
127786 2019-07-10T06:09:11 Z 윤교준(#3110) Dragon 2 (JOI17_dragon2) C++14
15 / 100
 4000 ms 6520 KB 
#pragma GCC optimize("Ofast")
#include <bits/stdc++.h>
#define rf(x) (x)=0;while(*p<48)im='-'==*p++;while(47<*p)(x)=((x)<<3)+((x)<<1)+(*p++&15);if(im)(x)=-(x);
#define eb emplace_back
#define allv(V) ((V).begin()),((V).end())
#define sorv(V) sort(allv(V))
#define univ(V) (V).erase(unique(allv(V)),(V).end())
using namespace std;
typedef unsigned int ui;
typedef long long ll;
static unsigned char str[5500055], *p=str;
static bool im=false;

const int MAXN = 30005;
const int MAXK = 30005;
const int MAXQ = 100005;

vector<ui> QV[MAXK], QRV[MAXK];
vector<ui> AV[MAXK];
ll DA[MAXN], DB[MAXN];
ui QI[MAXK], QRI[MAXK];
bitset<MAXK> haveToCal;
bitset<MAXK> chk;
bitset<MAXN> isInv;

int X[MAXN], Y[MAXN];
ui A[MAXN];
ui B[MAXQ], C[MAXQ];

static ui Ans[MAXQ];
int AX, AY, BX, BY;
int N, K, Q;

void process() {
	ll xa, ya, taa, tab, t, ta, tb;
	for(ui ga = K; ga; ga--) if(haveToCal[ga]) {
		chk.reset();
		memset(QI, 0, (K+1)<<2);
		memset(QRI, 0, (K+1)<<2);
		for(ui qi : QV[ga]) {
			ui gb = qi & 32767;
			QI[gb] = qi >> 15;
			chk[gb] = true;
		}
		for(ui qi : QRV[ga]) {
			ui gb = qi & 32767;
			QRI[gb] = qi >> 15;
			chk[gb] = true;
		}
		for(ui a : AV[ga]) {
			xa = X[a]; ya = Y[a]; taa = -DA[a]; tab = -DB[a];
			if(isInv[a]) {
				for(ui b = a-1, *gb = A+b; b; b--, gb--) if(chk[*gb]) {
					t = xa*Y[b] - ya*X[b];
					ta = taa + DA[b] - t;
					tb = tab + DB[b] - t;
					if((ta < 0 && tb < 0) || (ta > 0 && tb > 0)) continue;
					if(ta > 0) {
						Ans[QI[*gb]]++;
						if(!isInv[b]) Ans[QRI[*gb]]++;
					} else Ans[QRI[*gb]]++;
				}
			} else {
				for(ui b = a-1, *gb = A+b; b; b--, gb--) if(chk[*gb]) {
					t = xa*Y[b] - ya*X[b];
					ta = taa + DA[b] - t;
					tb = tab + DB[b] - t;
					if((ta < 0 && tb < 0) || (ta > 0 && tb > 0)) continue;
					if(tb > 0) {
						Ans[QI[*gb]]++;
						if(isInv[b]) Ans[QRI[*gb]]++;
					} else Ans[QRI[*gb]]++;
				}
			}
		}
	}
}

void precal() {
	for(ui i = N; i; i--) AV[A[i]].eb(i);
	for(int i = 1; i <= N; i++) {
		ll x = X[i], y = Y[i];
		DA[i] = y*AX - x*AY;
		DB[i] = y*BX - x*BY;

		ll t = (ll(AX)*BY - ll(AY)*BX) + (y*BX - x*BY) + (x*AY - y*AX);
		if(t > 0) isInv[i] = true;
	}
	for(ui i = Q; i; i--) {
		QV[B[i]].eb(i<<15 | C[i]);
		QRV[C[i]].eb(i<<15 | B[i]);
		haveToCal[B[i]] = haveToCal[C[i]] = true;
	}
}

int main() {
	fread(str, 1, 5500055, stdin);

	rf(N) rf(K)
	for(int i = 1; i <= N; i++) { rf(X[i]) rf(Y[i]) rf(A[i]) }
	rf(AX) rf(AY) rf(BX) rf(BY) rf(Q)
	for(int i = 1; i <= Q; i++) { rf(B[i]) rf(C[i]) }
	
	precal();
	process();

	for(int i = 1; i <= Q; i++)
		printf("%u\n", Ans[i]);
	return 0;
}

Compilation message

dragon2.cpp: In function 'int main()':
dragon2.cpp:97:7: warning: ignoring return value of 'size_t fread(void*, size_t, size_t, FILE*)', declared with attribute warn_unused_result [-Wunused-result]
  fread(str, 1, 5500055, stdin);
  ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~

Subtask #1
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 59 ms 2680 KB Output is correct
2 Correct 57 ms 2680 KB Output is correct
3 Correct 66 ms 2880 KB Output is correct
4 Correct 51 ms 5956 KB Output is correct
5 Correct 42 ms 6520 KB Output is correct
6 Correct 15 ms 2812 KB Output is correct
7 Correct 15 ms 2936 KB Output is correct
8 Correct 50 ms 2684 KB Output is correct
9 Correct 22 ms 2680 KB Output is correct
10 Correct 28 ms 2712 KB Output is correct

Subtask #2
0 / 45.0

# Verdict Execution time Memory Grader output
1 Execution timed out 4027 ms 4088 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #3
0 / 40.0

# Verdict Execution time Memory Grader output
1 Correct 59 ms 2680 KB Output is correct
2 Correct 57 ms 2680 KB Output is correct
3 Correct 66 ms 2880 KB Output is correct
4 Correct 51 ms 5956 KB Output is correct
5 Correct 42 ms 6520 KB Output is correct
6 Correct 15 ms 2812 KB Output is correct
7 Correct 15 ms 2936 KB Output is correct
8 Correct 50 ms 2684 KB Output is correct
9 Correct 22 ms 2680 KB Output is correct
10 Correct 28 ms 2712 KB Output is correct
11 Execution timed out 4027 ms 4088 KB Time limit exceeded
12 Halted 0 ms 0 KB -