#include<bits/stdc++.h>
using namespace std;
#define int long long
#define ii pair <int, int>
#define app push_back
#define all(a) a.begin(), a.end()
#define bp __builtin_popcountll
#define ll long long
#define mp make_pair
#define f first
#define s second
#define Time (double)clock()/CLOCKS_PER_SEC
const int N = 30007;
struct Point {
int x, y, c;
Point operator - (Point p) {
return {x - p.x, y - p.y};
}
Point operator + (Point p) {
return {x + p.x, y + p.y};
}
int operator *(Point p) {
return x * p.y - y * p.x;
}
} a[N];
int sign(int x) {
if (x < 0)
return -1;
else if (x == 0)
return 0;
else
return 1;
}
vector <Point> cA[N][2], cB[N][2];
vector <Point> al[N];
signed main() {
#ifdef HOME
freopen("input.txt", "r", stdin);
#else
#define endl '\n'
ios_base::sync_with_stdio(0); cin.tie(0);
#endif
int n, m;
cin >> n >> m;
for (int i = 1; i <= n; ++i) {
cin >> a[i].x >> a[i].y >> a[i].c;
}
Point A, B;
cin >> A.x >> A.y;
cin >> B.x >> B.y;
auto half = [&](Point p) {
int ans = (p - A) * (B - A) > 0;
return ans;
};
for (int i = 1; i <= n; ++i) {
cA[a[i].c][half(a[i])].app(a[i]);
cB[a[i].c][half(a[i])].app(a[i]);
al[a[i].c].app(a[i]);
}
auto compA = [&](Point a, Point b) {
return (a - A) * (b - A) > 0;
};
auto compB = [&](Point a, Point b) {
return (a - B) * (b - B) > 0;
};
for (int i = 1; i <= m; ++i) {
for (int j = 0; j < 2; ++j) {
sort(all(cA[i][j]), compA);
sort(all(cB[i][j]), compB);
}
}
int q;
cin >> q;
while (q--) {
int c1, c2;
cin >> c1 >> c2;
swap(c1, c2);
int ans = 0;
for (auto a : al[c1]) {
int t = half(a);
{
int l = -1, r = cA[c2][t^1].size();
auto sim = A + (A - a);
while (l < r - 1) {
int m = (l + r) >> 1;
if (compA(cA[c2][t^1][m], sim))
l = m;
else
r = m;
}
if ((A - a) * (B - a) > 0) {
ans -= r;
}
else {
ans += r;
}
}
{
int l = -1, r = cB[c2][t^1].size();
auto sim = B + (B - a);
while (l < r - 1) {
int m = (l + r) >> 1;
if (compB(cB[c2][t^1][m], sim))
l = m;
else
r = m;
}
if ((A - a) * (B - a) > 0) {
ans += r;
}
else {
ans -= r;
}
}
if ((A - a) * (B - a) > 0) {
for (auto b : cA[c2][t]) {
if (compA(a, b) && compB(b, a)) {
++ans;
}
}
}
else {
for (auto b : cA[c2][t]) {
if (compA(b, a) && compB(a, b)) {
++ans;
}
}
}
}
cout << ans << endl;
}
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
11 ms |
4224 KB |
Output is correct |
2 |
Correct |
16 ms |
4352 KB |
Output is correct |
3 |
Correct |
25 ms |
4224 KB |
Output is correct |
4 |
Correct |
50 ms |
4616 KB |
Output is correct |
5 |
Correct |
38 ms |
4472 KB |
Output is correct |
6 |
Correct |
9 ms |
4224 KB |
Output is correct |
7 |
Correct |
9 ms |
4224 KB |
Output is correct |
8 |
Correct |
8 ms |
4224 KB |
Output is correct |
9 |
Correct |
11 ms |
4224 KB |
Output is correct |
10 |
Correct |
8 ms |
4224 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
278 ms |
7272 KB |
Output is correct |
2 |
Correct |
548 ms |
7928 KB |
Output is correct |
3 |
Correct |
39 ms |
7936 KB |
Output is correct |
4 |
Correct |
25 ms |
7808 KB |
Output is correct |
5 |
Correct |
32 ms |
7416 KB |
Output is correct |
6 |
Correct |
472 ms |
6868 KB |
Output is correct |
7 |
Correct |
443 ms |
6868 KB |
Output is correct |
8 |
Correct |
27 ms |
6876 KB |
Output is correct |
9 |
Correct |
295 ms |
7016 KB |
Output is correct |
10 |
Correct |
23 ms |
6888 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
11 ms |
4224 KB |
Output is correct |
2 |
Correct |
16 ms |
4352 KB |
Output is correct |
3 |
Correct |
25 ms |
4224 KB |
Output is correct |
4 |
Correct |
50 ms |
4616 KB |
Output is correct |
5 |
Correct |
38 ms |
4472 KB |
Output is correct |
6 |
Correct |
9 ms |
4224 KB |
Output is correct |
7 |
Correct |
9 ms |
4224 KB |
Output is correct |
8 |
Correct |
8 ms |
4224 KB |
Output is correct |
9 |
Correct |
11 ms |
4224 KB |
Output is correct |
10 |
Correct |
8 ms |
4224 KB |
Output is correct |
11 |
Correct |
278 ms |
7272 KB |
Output is correct |
12 |
Correct |
548 ms |
7928 KB |
Output is correct |
13 |
Correct |
39 ms |
7936 KB |
Output is correct |
14 |
Correct |
25 ms |
7808 KB |
Output is correct |
15 |
Correct |
32 ms |
7416 KB |
Output is correct |
16 |
Correct |
472 ms |
6868 KB |
Output is correct |
17 |
Correct |
443 ms |
6868 KB |
Output is correct |
18 |
Correct |
27 ms |
6876 KB |
Output is correct |
19 |
Correct |
295 ms |
7016 KB |
Output is correct |
20 |
Correct |
23 ms |
6888 KB |
Output is correct |
21 |
Correct |
273 ms |
7396 KB |
Output is correct |
22 |
Correct |
536 ms |
7800 KB |
Output is correct |
23 |
Correct |
711 ms |
7808 KB |
Output is correct |
24 |
Correct |
460 ms |
8184 KB |
Output is correct |
25 |
Correct |
77 ms |
7776 KB |
Output is correct |
26 |
Correct |
72 ms |
7672 KB |
Output is correct |
27 |
Correct |
30 ms |
7424 KB |
Output is correct |
28 |
Correct |
32 ms |
7552 KB |
Output is correct |
29 |
Execution timed out |
4080 ms |
7340 KB |
Time limit exceeded |
30 |
Halted |
0 ms |
0 KB |
- |