HDU3832 Earth Hour 斯坦纳树（三点连通）

题目链接

http://acm.hdu.edu.cn/showproblem.php?pid=3832

题意

背景：为了响应停电一小时活动，某校决定关掉除了图书馆，教室，学生宿舍(编号为1，2，3)三处其他所有的灯，且要求这三处是相连的。每处的灯有一个圆形的照亮范围，如果两处的灯的照亮范围存在相交区域，那么称之为相连的。学校想在符合上述限制下，尽可能的关闭多的灯，求出这个数量，如果这3处就不相连，那么输出-1

思路

据说是斯坦纳树，就直接查结论了：以1，2，3点为起点跑三次最短路，边的权值设为1,依次获得他们到其他点的最短路d1[],d2[],d3[]。要使得保持连通，那么至少得留下
ans=min(d1[]+d2[]+d3[],ans) 再加个1这么多盏灯。也即找到一点，这一点到三点的路径中，包含的点数是最少的，因为边权为1，跑过多少边相当于跑过多少点，然后加上这个点本身，就是至少得留下的点数。

#include<cstdio>
#include<queue>
#include<iostream>
#include<vector>
#include<map>
#include<cstring>
#include<string>
#include<set>
#include<stack>
#include<algorithm>
#define cle(a) memset(a,0,sizeof(a))
#define inf(a) memset(a,0x3f,sizeof(a))
#define ll long long
#define Rep(i,a,n) for(int i=a;i<=n;i++)
using namespace std;
#define INF2 9223372036854775807ll
const int INF = ( 2e9 ) + 2;
const ll maxn = 210;
const int maxm = maxn*maxn/2;
vector<int> g[maxn];
struct P
{
    int x,y,r;
}p[maxn];
int d[4][maxn],vis[maxn];
bool check(int i,int j)
{
    int d=(p[i].x-p[j].x)*(p[i].x-p[j].x)+(p[i].y-p[j].y)*(p[i].y-p[j].y);
    return d<=(p[i].r+p[j].r)*(p[i].r+p[j].r);
}
void BFS(int s,int n,int *d)
{
    memset(vis,0,sizeof(vis));
    for(int i=1;i<=n;i++)d[i]=INF;
    d[s]=0;
    queue<int> q;
    q.push(s);
    vis[s]=1;
    while(!q.empty())
    {
        int u=q.front();q.pop();
        for(int i=0,L=g[u].size();i<L;i++)
        {
            int v=g[u][i];
            if(!vis[v])
            {
                d[v]=d[u]+1;
                vis[v]=1;
                q.push(v);
            }
        }
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n;
        scanf("%d",&n);

        for(int i=1;i<=n;i++)g[i].clear();
        for(int i=1;i<=n;i++)
        {
            scanf("%d%d%d",&p[i].x,&p[i].y,&p[i].r);
        }
        for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++)
        if(check(i,j))
        {
            g[i].push_back(j);
            g[j].push_back(i);
        }
        for(int i=1;i<=3;i++)
        BFS(i,n,d[i]);
        int ans=INF;
        for(int i=1;i<=n;i++)
        {
            if(d[1][i]==INF||d[2][i]==INF||d[3][i]==INF)continue;
            ans=min(d[1][i]+d[2][i]+d[3][i],ans);
        }
        if(ans<INF)
        printf("%d\n",n-ans-1);
        else
        printf("-1\n");
    }
}