HDU3832 Earth Hour 2011 Multi-University Training Contest 1 - Host by HNU

题目地址：http://acm.hdu.edu.cn/showproblem.php?pid=3832

题意：有n块区域，求保证前三个区域联通的情况下最多可以删掉几个区域；

可以想象最优连通方案不是一条链就是一个三叉。因为是无权的，最少点数实际上就是最短路。当分叉点在特定点上的时候，三叉就退化为一条链，所以我们只需要枚举分叉点即可。

先3次bfs求出各个点到0，1，2三个区域的最短路，然后枚举点求联通最小量

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
#include<stack>
using namespace std;
#define eps 1e-6
#define INF 1000000
struct Point
{
    double x,y,r;
}p[230];
int n,d[3][230],map[230][230];
int vis[230];
double dis(Point a,Point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
void init()
{
    scanf("%d",&n);
    memset(map,0,sizeof(map));
    for(int i=0;i<3;i++)
    for(int j=0;j<n;j++)
    d[i][j]=INF;
    d[0][0]=0;
    d[1][1]=0;
    d[2][2]=0;
    for(int i=0;i<n;i++)
    {
        scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].r);
    }
    for(int i=0;i<n;i++)
    for(int j=i+1;j<n;j++)
    {
        if(dis(p[i],p[j])-p[i].r-p[j].r<=eps)
        {
            map[i][j]=1;
            map[j][i]=1;                                   // 在有公共区域的圆之间建边
        }
    }

}
void bfs(int k)
{
    memset(vis,0,sizeof(vis));
    queue<int> q;
    q.push(k);
    vis[k]=1;
    d[k][k]=0;
    while(!q.empty())
    {
         int x=q.front();q.pop();

         for(int i=0;i<n;i++)
         if(!vis[i]&&map[x][i])
         {
             vis[i]=1;
             d[k][i]=d[k][x]+1;
             q.push(i);
         }
    }
}

int main()
{
    int Case;
    scanf("%d",&Case);
    while(Case--)
    {
        init();
        int ans=INF;
        bfs(0);
        bfs(1);
        bfs(2);
        for(int i=0;i<n;i++)                           //枚举分叉点（也包括链的情况）
        {
            if(d[1][i]+d[2][i]+d[0][i]+1<ans)
            ans=d[1][i]+d[2][i]+d[0][i]+1;
        }
        if(ans==INF)printf("-1\n");
        else printf("%d\n",n-ans);
    }
return 0;
}