HDU - 5934

tarjan 視頻講解

/**
 *      題目連接：https://vjudge.net/problem/HDU-5934
 * 題意：給你n個炸彈，引爆每一個炸彈會有必定的花費。每一個炸彈給出座標x，y，半徑r，引爆花費；
 * 引爆一個炸彈會把範圍內的炸彈引爆，連鎖反應。 如今想把全部炸彈引爆的最小花費。
 * 
 * 解題思路：強連通縮點。根據a可以引爆b，能夠在建一條a到b的單向邊。若是是一個強連通（這一部分的圖，
 * 任意兩點均可以相互到達）那麼就把這個強連通份量變成一個點，值最份量的最小值。這樣圖就變成有向無環圖了。
 * 考慮到每一個點的花費都是大於0的，因此引爆開始點最划算，即爲入度爲0的點。
 * 
 * 前置技能 tarjan 縮點。
*/

#include <bits/stdc++.h>
using namespace std;

const int maxn=1000+10;
const int INF=2e9+1e8;
vector<int>E[maxn];
struct Point 
{
    int x,y,r,cost;
}boom[maxn];
bool judge(Point a,Point b)
{
    if( 1ll*(a.x-b.x)*(a.x-b.x)+1ll*(a.y-b.y)*(a.y-b.y)<=1ll*a.r*a.r ) return true;
    return false;
}
int dfn[maxn],low[maxn],id,vis[maxn],ans,deg[maxn];
int num[maxn],cnt,cost[maxn];//對點進行從新編號，（數組num），按照聯通份量進行編號
stack<int>S;
void init()
{
    id=cnt=0;
    memset(deg,0,sizeof(deg));
    memset(num,0,sizeof(num));
    memset(vis,0,sizeof(vis));
    memset(dfn,0,sizeof(dfn));
}
void tarjan(int x)
{
    low[x]=dfn[x]=++id;
    S.push(x);
    vis[x]=1;
    for(int i=0;i<(int)E[x].size();i++)
    {
        int to=E[x][i];
        if(!dfn[to]) 
        {
            tarjan(to);
            low[x]=min(low[x],low[to]);
        }
        else if(vis[to]) low[x]=min(low[x],dfn[to]);
    }
    if(low[x]==dfn[x])
    {
        int mincost=INF,in=0;
        cnt++;
        while(1)
        {
            int now=S.top();
            S.pop();
            vis[now]=0;
            num[now]=cnt;
            mincost=min(mincost,boom[now].cost);
            if(now==x) break;
        }
        cost[cnt]=mincost;
    }
}
int main()
{
    int T,cas=1;
    scanf("%d",&T);
    while(T--)
    {
        int n;
        scanf("%d",&n);
        init();
        for(int i=1;i<=n;i++)
        {
            E[i].clear();
            scanf("%d%d%d%d",&boom[i].x,&boom[i].y,&boom[i].r,&boom[i].cost);
        }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                if(i==j) continue;
                if(judge(boom[i],boom[j])) E[i].push_back(j);
            }
        }
        for(int i=1;i<=n;i++)
            if(!dfn[i]) tarjan(i);
        for(int i=1;i<=n;i++)
        {
            for(int j=0;j<(int)E[i].size();j++)
            {
                int to=E[i][j];
                if(num[i]!=num[to]) deg[num[to]]++;
            }
        }
        ans=0;
        for(int i=1;i<=cnt;i++) if(deg[i]==0) 
            ans+=cost[i];
        printf("Case #%d: %d\n",cas++,ans);
    }
    return 0;
}