Algoogle

Algorithm for Programming Contest

AOJ 1023 Amazing Graze

Category: AOJ Tag: implementation

Amazing Graze

問題概要


2次元空間上にAN個の戦闘機とBN個の弾がある.
戦闘機の戦闘力は半径4R以内にある弾の数になる. 戦闘力の総和を求めよ.

解法


戦闘機ごとに全ての弾と判定していてはだめ.
弾をx座標ごとに分けて, 戦闘機の前後4Rのx座標にある弾のうち, 上下4Rのy座標にある弾と判定する.

コード


(1023.cpp) download
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
#include <bits/stdc++.h>
using namespace std;

int n, m, r, x[100010], y[100010];
vector<int> ys[10010];

inline int dist(int a, int b, int c, int d) { return (a-c)*(a-c)+(d-b)*(d-b);}

int solve()
{
        for (int i = 0; i < 10010; i++) sort(begin(ys[i]),end(ys[i]));

        int ans = 0;
        for (int i = 0; i < n; i++) {
                for(int xx = max(0,x[i]-4*r); xx <= min(10000, x[i]+4*r); xx++) {
                        auto lb = lower_bound(begin(ys[xx]),end(ys[xx]), y[i]-4*r),
                             ub = upper_bound(begin(ys[xx]),end(ys[xx]), y[i]+4*r);
                        for(auto it = lb; it != ub; it++)
                                if(dist(x[i], y[i], xx, *it) <= 16*r*r) ans++;
                }
        }
        return ans;
}

bool input()
{
        scanf("%d%d%d", &n, &m, &r);
        for (int i = 0; i < n; i++) scanf("%d%d", x+i, y+i);
        for (int i = 0; i < 10010; i++) ys[i].clear();
        for (int i = 0; i < m; i++) {
                int xx, yy;
                scanf("%d%d", &xx, &yy);
                ys[xx].push_back(yy);
        }
        return n or m;
}

int main()
{
        while(input()) printf("%d\n", solve());
        return 0;
}

Comments