oj.uz

Submission #127766

# Submission time Handle Problem Language Result Execution time Memory
127766 2019-07-10T05:45:09 Z 윤교준(#3110) Dragon 2 (JOI17_dragon2) C++14
15 / 100
 4000 ms 6520 KB 
#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> chkQ, chkQR;
bitset<MAXN> isInv;

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

static unsigned int 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]) {
		chkQ.reset(); chkQR.reset();
		for(ui qi : QV[ga]) {
			ui gb = qi & 32767;
			QI[gb] = qi >> 15;
			chkQ[gb] = true;
		}
		for(ui qi : QRV[ga]) {
			ui gb = qi & 32767;
			QRI[gb] = qi >> 15;
			chkQR[gb] = true;
		}
		for(ui a : AV[ga]) {
			xa = X[a]; ya = Y[a]; taa = -DA[a]; tab = -DB[a];
			bool flag = isInv[a];
			for(ui b = a-1, gb; b; b--) {
				gb = A[b];
				if(!chkQ[gb] && !chkQR[gb]) continue;
				t = xa*Y[b] - ya*X[b];
				ta = taa + DA[b] - t;
				tb = tab + DB[b] - t;
				if(chkQ[gb]) {
					if((!flag && ta < 0 && tb > 0) || (flag && ta > 0 && tb < 0)) {
						Ans[QI[gb]]++;
					}
				}
				if(chkQR[gb]) {
					if(isInv[b]) {
						if(ta < 0 && tb > 0) {
							Ans[QRI[gb]]++;
						}
					} else {
						if(ta > 0 && tb < 0) {
							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:95: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 52 ms 2808 KB Output is correct
2 Correct 73 ms 2704 KB Output is correct
3 Correct 90 ms 2856 KB Output is correct
4 Correct 53 ms 6008 KB Output is correct
5 Correct 43 ms 6520 KB Output is correct
6 Correct 16 ms 2936 KB Output is correct
7 Correct 15 ms 2936 KB Output is correct
8 Correct 26 ms 2680 KB Output is correct
9 Correct 21 ms 2680 KB Output is correct
10 Correct 22 ms 2680 KB Output is correct

Subtask #2
0 / 45.0

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

Subtask #3
0 / 40.0

# Verdict Execution time Memory Grader output
1 Correct 52 ms 2808 KB Output is correct
2 Correct 73 ms 2704 KB Output is correct
3 Correct 90 ms 2856 KB Output is correct
4 Correct 53 ms 6008 KB Output is correct
5 Correct 43 ms 6520 KB Output is correct
6 Correct 16 ms 2936 KB Output is correct
7 Correct 15 ms 2936 KB Output is correct
8 Correct 26 ms 2680 KB Output is correct
9 Correct 21 ms 2680 KB Output is correct
10 Correct 22 ms 2680 KB Output is correct
11 Execution timed out 4097 ms 4216 KB Time limit exceeded
12 Halted 0 ms 0 KB -