题目没有什么难理解的,就是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 的作用: