CCF CSPコンテストの質問 - のWi-Fi（201403から4）

問題の説明

　　現在、大きな平面室、各無線ルータの無線ルータは、いくつかの点で固定されたn個あります。任意の2つの無線ルータは、長い距離がRを超えないように相互にネットワーク接続を確立することができます。
　　また、m番目の無線ルーターの別の可能な配置。あなたはこれらの場所に新しいルータのほとんどのk個の作成時に選択することができます。
　　あなたの目標は、第1のルータ及び第2のルータ間のネットワーク接続が可能な限りわずか中継ルータを介して接続されているようにすることです。最適解での中継ルータの最小数はどのくらいあるのでしょうか？

入力形式

　　最初の行は、4つの正の整数N、M、K、Rを含んでいます。（2≤N≤100,1≤K ≤M≤100、1≤R≤108）。
　　次のn行は、二つの整数XIとYIを含む各ラインは、良好な無線ルータ（XI、YI）の点に配置されている表します。のみこの状況で第1の入力データと第二のルータは、nルータ（中継ルータ一連の）相互接続することができることを確認しました。
　　続いて行をM、二つの整数Xi及びYIを含む各行は、追加地点で表す（XI、李）がルータであってもよいです。
　　すべての入力座標の絶対値は、入力が変化座標ことを確実にするために、より多く108以下です。

出力フォーマット

　　出力後すなわちそれを通して少なくとも1つの第2のルータにルータから、指定された場所に中継ルータの数を追加のルータをK、ほんのわずかです。

サンプル入力

5 3 1 3
0 0
5
0 3
0 5
3 5
3 3
4 4
3 0

サンプル出力

2

コード

#include <iostream>
#include <utility>
#include <vector>
#include <algorithm>
#include <queue>
using namespace std; 

int N;
int M;
int K;
int R;
vector<pair<int, int> > sites;

bool isConnected(const pair<int, int> &a, const pair<int, int> &b) {
    int w = a.first-b.first;
    int h = a.second-b.second;
    return (long long)R * R >= (long long)w*w + (long long)h*h;
}

int main() {
    cin >> N >> M >> K >> R;
    int x, y;
    for (int i = 0; i < N+M; i++) {
        cin >> x >> y;
        sites.push_back(pair<int, int>(x, y));
    }
    vector<bool> used(N+M, false);
    vector<int> newSitesUsed(N+M, 0);
    queue<int> myque;
    used[0] = true;
    myque.push(0);
    int numOfSites = 0; 
    bool flag = true;
    while (flag) {
        int sz = myque.size();
        int cur;
        for (int i = 0; i < sz; i++) {
            cur = myque.front();
            myque.pop();
            if (cur == 1) {
                flag = false;
                break;
            }
            for (int j = 0; j < N+M; j++) {
                if (j != cur && isConnected(sites[j], sites[cur])) {
                    if (!used[j]) {
                        int t = (j >= N ? newSitesUsed[cur]+1 : newSitesUsed[cur]);
                        if (t <= K) {
                            used[j] = true;
                            newSitesUsed[j] = t;
                            myque.push(j);
                        }
                    } else {
                        int t = (j >= N ? newSitesUsed[cur]+1 : newSitesUsed[cur]);
                        if (t < newSitesUsed[j]) {
                            newSitesUsed[j] = t;
                            myque.push(j);
                        }
                    }
                }
            }
        }
        numOfSites++;
    }
    cout << numOfSites-2 << endl;
    return 0;
}