HDU 5934(强连通分量)

题目链接：点击打开链接

题目大意：给你一大堆炸弹的位置还有爆炸范围以及各自的花费，问你如何花最少的钱让所有的炸弹炸掉

题目思路：选拔赛的时候还以为是并查集..太菜了...其实这道题就是先套tarjan的模板(需要记录每个分量中最小花费的炸弹)，然后就可以将同一个分量里的点看成一个点，然后建新图，其实就是看一下入度，如果入度是0，就说明不可以用别的炸弹把他炸了，如果入度不是0，就说明别的分量炸的时候可以顺便把他炸了，所以就只用把入度为0的强连通分量里的最小分量加起来就可以了。

以下是代码：

#include <iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<math.h>
using namespace std;
#define MAXN 1005
#define MAXM 1000005
#define inf 0x3f3f3f3f
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define per(i,a,b) for(int i=a;i>=b;i--)
int ea=1;
struct node{
    ll x,y,r,c;
}a[MAXN];
ll cost[MAXN];
struct Edge{
    int to,next;
}edge[MAXM];
int head[MAXN],tot;
int low[MAXN],dfn[MAXN],Stack[MAXN],belong[MAXN],du[MAXN];
int index,top;
int scc;
bool instack[MAXN];
int num[MAXN];
void add(int u,int v){
    edge[tot].to=v;edge[tot].next=head[u];head[u]=tot++;
}
void tarjan(int u){
    int v;
    low[u]=dfn[u]=++index;
    Stack[top++]=u;
    instack[u]=true;
    for(int i=head[u];i!=-1;i=edge[i].next){
        v=edge[i].to;
        if(!dfn[v]){
            tarjan(v);
            if(low[u]>low[v])low[u]=low[v];
        }
        else if(instack[v]&&low[u]>dfn[v])low[u]=dfn[v];
    }
    if(low[u]==dfn[u]){
        scc++;
        do{
            v=Stack[--top];
            cost[scc]=min(cost[scc],a[v].c);
            instack[v]=false;
            belong[v]=scc;
            num[scc]++;
        }while(v!=u);
    }
}
void solve(int N){
    memset(dfn,0,sizeof(dfn));
    memset(instack,0,sizeof(instack));
    memset(num,0,sizeof(num));
    index=scc=top=0;
    rep(i,1,N){
        if(!dfn[i]){
            tarjan(i);
        }
    }
}
void init(){
    tot=0;
    memset(head,-1,sizeof(head));
    memset(du,0,sizeof(du));
    rep(i,1,MAXN){
        cost[i]=inf;
    }
}

int main()
{
    int t,n;
    scanf("%d",&t);
    while(t--){
        init();
        scanf("%d",&n);
        rep(i,1,n){
            scanf("%lld%lld%lld%lld",&a[i].x,&a[i].y,&a[i].r,&a[i].c);
        }
        rep(i,1,n){
            rep(j,1,n){
                if(i==j)continue;
                if((a[i].x-a[j].x)*(a[i].x-a[j].x)+(a[i].y-a[j].y)*(a[i].y-a[j].y)<=a[i].r*a[i].r){
                    add(i,j);
                }
            }
        }
        solve(n);
        rep(i,1,n){
            for(int j=head[i];j!=-1;j=edge[j].next){
                int v=edge[j].to;
                if(belong[i]!=belong[v]){
                    du[belong[v]]++;
                }
            }
        }
        ll ans=0;
        rep(i,1,scc){
            //cout<<cost[i]<<endl;
            if(!du[i]){
                ans+=cost[i];
            }
        }
        printf("Case #%d: %lld\n",ea++,ans);
    }
    return 0;
}