ccf_无线网络

目前在一个很大的平面房间里有 n 个无线路由器,每个无线路由器都固定在某个点上。任何两个无线路由器只要距离不超过 r 就能互相建立网络连接。
　　除此以外,另有 m 个可以摆放无线路由器的位置。你可以在这些位置中选择至多 k 个增设新的路由器。
　　你的目标是使得第 1 个路由器和第 2 个路由器之间的网络连接经过尽量少的中转路由器。请问在最优方案下中转路由器的最少个数是多少?

输入格式

　　第一行包含四个正整数 n,m,k,r。(2 ≤ n ≤ 100,1 ≤ k ≤ m ≤ 100, 1 ≤ r ≤ 108)。
　　接下来 n 行,每行包含两个整数 xi 和 yi,表示一个已经放置好的无线 路由器在 (xi, yi) 点处。输入数据保证第 1 和第 2 个路由器在仅有这 n 个路由器的情况下已经可以互相连接(经过一系列的中转路由器)。
　　接下来 m 行,每行包含两个整数 xi 和 yi,表示 (xi, yi) 点处可以增设 一个路由器。
　　输入中所有的坐标的绝对值不超过 108,保证输入中的坐标各不相同。

输出格式

　　输出只有一个数,即在指定的位置中增设 k 个路由器后,从第 1 个路 由器到第 2 个路由器最少经过的中转路由器的个数。

样例输入

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

样例输出

2

解题思路：

     dijkstra算法

     dijkstra基本思路：它能求从某一个点出发到所有点的最短路径

     该题要求，从1号点出发到2号点尽可能花少的点，且不能超过k个增发路由器

三个数组:

    f1[maxn]; f1[i] = j 指1号路由器最少经过j个中继路由器连上i个路由器

    f2[maxn]; f2[i] = j 指1号路由器最少经过增设j个中继路由器连上i个路由器

    flag[maxn]; flag[i] = 1 指i号路由器已确定最少路线

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
using namespace std;
const int maxn = 205;
struct Node{
    int x,y;
}node[maxn];
int f1[maxn]; /// f1[i] = j 指1号路由器最少经过j个中继路由器连上i个路由器
int f2[maxn]; /// f2[i] = j 指1号路由器最少经过增设j个中继路由器连上i个路由器
int flag[maxn]; /// flag[i] = 1 指i号路由器已确定最少路线     
bool check(int p_one,int p_two,int r)
{
    ///判断两个路由器是否能连通
    return  sqrt((long long)(node[p_one].x - node[p_two].x)*(long long)(node[p_one].x - node[p_two].x)
    + (long long)(node[p_one].y - node[p_two].y)*(long long)(node[p_one].y - node[p_two].y)) <= r;
}
void init(int n,int m)
{
    for(int i=1;i<=n+m;++i)
    {
        f1[i] = maxn;
        f2[i] = maxn;
    }
    f1[1] = 0;
    f2[1] = 0;
}
int dij(int n,int m,int r,int Nchose)
{
    // 开始n-1次循环
    memset(flag,0,sizeof(flag));
    flag[1] = 1;
    for(int i=1;i<n+m;++i)
    {
        int k=i,mi_num = maxn,ami_num = maxn;
        /***找到最小值***/
        for(int j=2;j<=n+m;j++)
        {
            if(flag[j]) continue;
            if(f1[j] < mi_num ||(f1[j] == mi_num && f2[j] < ami_num))
            {
                mi_num = f1[j];
                ami_num = f2[j];
                k = j;
            }
        }
        //printf("-%d %d %d\n",k,f1[k],f2[k]);
        /*** 在最小值向周围更新权值 ***/
        flag[k] = 1;
        for(int j=2;j<=n+m;j++)
        {
            int t = j>n?1:0;
            //printf("%d %d %d %d %d %d %d %d\n",j,f1[j],f2[j],k,f1[k],f2[k],check(k,j,r),Nchose);
            if(check(k,j,r) && j!=k && f2[k]+t <= Nchose)
            {
                if(f1[k]+1<f1[j]){
                    f1[j] = f1[k]+1;
                    f2[j] = f2[k]+t;
                }
                else if(f1[k]+1 == f1[j] && f2[k]+t < f2[j])
                {
                    f1[j] = f1[k]+1;
                    f2[j] = f2[k]+t;
                }
            }
        }
    }
    return f1[2] -1;
}
int main()
{
    int n,m,k,r;
    scanf("%d%d%d%d",&n,&m,&k,&r);
    init(n,m);
    for(int i=1;i<=n+m;i++){
        scanf("%d%d",&node[i].x,&node[i].y);
        if(i>1 && check(1,i,r))
        {
            f1[i] = 1;
            f2[i] = 0;
            if(i>n) f2[i] = 1;
            if(i == 2) f1[i] = 0;
        }
    }
    int res = dij(n,m,r,k);
    printf("%d\n",res);
    return 0;
}