HDU-5934 BOMB (强连通分量）

---

题目没有什么难理解的，就是T次情况，每次有N个炸弹，告诉你炸弹的坐标，爆炸范围，消耗，在另一个爆炸范围内的炸弹，不需要单独引爆，输出最少的消耗；

---

我们首先需要计算炸弹之间的距离，如果A可以引爆B，那么有一条A指向B的边。
再对这个图进行联通分量的求解，找到几个联通分量
最后对已知的联通分量进行缩减（如果两个独立的联通分量，可以有一条边来连接，那么只需要引爆前一个就可以了。）

---

关于联通分量的一些事：

对于一个图，我们dfs遍历时有下图的顺序：
我们用数组dfn来存放这些初始时间顺序，在他访问到点的时候，我们还需要更新，看看能否通过其他路径，回到更早的点去。我们用low来存放更新后的顺序。相同low的是同一组联通分量
如图，我们这个有三组联通分量。

---

对于此题，e不需要单独引爆，可以通过上一个联通分量来引爆，所以不需要单独消耗。

AC代码：

//[HDU-5934] BOMB
//强连通分量tarjan算法
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <algorithm>
#define ll long long
const int maxn=1005;
const int inf = 0x3f3f3f3f;
using namespace std;
inline int read() {
	int x = 0, f = 1; char ch = getchar();
	while (!isdigit(ch)) {
		if (ch == '-') f = -1;
		ch = getchar();
	}
	while (isdigit(ch)) {
		x = (x << 1) + (x << 3) + (ch & 15);
		ch = getchar();
	}
	return x * f;
}
struct Bomb
{
    ll x, y, r, c;
}bomb[maxn];//存放每个炸弹的信息
struct Edge
{
    int to, nxt;
}edge[maxn*maxn];
ll cost[maxn];
int low[maxn],dfn[maxn],stk[maxn],vis[maxn],head[maxn],in[maxn],belong[maxn];
//stk模拟栈，dfn存放正常的顺序，low存放修改后回溯到更小的顺序，vis存放是不是访问过
//tim为时间戳，top为栈的操作数
int tim,top,cnt,cut;
void init()
{
    memset(head, -1, sizeof(head));
    cnt = tim = top = cut = 0;
    memset(vis, 0, sizeof(vis));
    memset(cost, inf, sizeof(cost));
    memset(dfn, 0, sizeof(dfn));
    memset(low, 0, sizeof(low));
    memset(in, 0, sizeof(in));
    memset(belong, 0, sizeof(belong));
}

void add(int u,int v){
    edge[cnt].to=v;
    edge[cnt].nxt=head[u];
    head[u]=cnt++;
}//head于edge组成一个链状数组，head为头节点，于u相关的链以head[u]开头
//edge的 to存放u指向的节点，nxt 为下一个节点。


void Tarjan(int u){
    low[u]=dfn[u]=++tim;
    stk[top++]=u;
    vis[u]=1;
    for(int i=head[u];i!=-1;i=edge[i].nxt){
        int x=edge[i].to;
        if(!dfn[x]){
            Tarjan(x);
            low[u]=min(low[u],low[x]);
        }
        else if(vis[x])
            low[u]=min(low[u],dfn[x]);
    }
    int v;
    if(low[u]==dfn[u]){
        cut++;//联通分量++；
        do{
            v=stk[--top];//出栈
            belong[v]=cut;//v属于联通分量cut。
            cost[cut]=min(cost[cut],bomb[v].c);//找到一个联通分量里最小的那个花费
            vis[v]=0;
        }while(u!=v);
    }
}

int main(){
    int T, N;
    scanf("%d",&T);
    for(int ca=1;ca<=T;ca++){
        init();
        scanf("%d",&N);
        for(int i=1;i<=N;i++){
           scanf("%lld%lld%lld%lld",&bomb[i].x,&bomb[i].y, &bomb[i].r, &bomb[i].c);
        }
        for(int i=1;i<=N;i++)
            for(int j=1;j<=N;j++){
                if(i==j) continue;
                if((bomb[i].x-bomb[j].x)*(bomb[i].x-bomb[j].x)+(bomb[i].y-bomb[j].y)*(bomb[i].y-bomb[j].y) <= bomb[i].r*bomb[i].r)
                    add(i, j);
            }
        for(int i=1;i<=N;i++)
            if(!dfn[i]) Tarjan(i);
        for(int i=1;i<=N;i++)
            for(int j=head[i];j!=-1;j=edge[j].nxt){
                int x=edge[j].to;
                if(belong[i]!=belong[x])//i可以到的点与i不属于同一个联通分量时
                    in[belong[x]]++;//把v的联通分量标记一下。
            }
        ll ans=0;
        for(int i=1;i<=cut ;i++)
            if(!in[i])//跳过被标记的联通分量
                ans+=cost[i];
        printf("Case #%d: %lld\n",ca,ans);  
    }
    return 0;
}

这里画图描述一下 head 和 edge 的作用：