HDU 5934 强连通分量

点击打开链接

题意：给了n个炸弹，然后让你引爆其中的炸弹，然后使得所有的炸弹都爆炸，一个炸弹爆炸会使得它半径内的炸弹爆炸，以此类推

思路：强连通缩点，然后找入度为0的所有的联通块，然后在联通块内找一个花费最小的即可

#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=1010;
int V;
vector<int>G[maxn],rG[maxn],vs,ans[maxn];
bool used[maxn],vis[maxn][maxn],fvis[maxn];
int cmp[maxn],in[maxn];
void add_edge(int from,int to){
    G[from].push_back(to);
    rG[to].push_back(from);
}
void dfs(int v){
    used[v]=1;
    for(unsigned int i=0;i<G[v].size();i++){
        if(!used[G[v][i]]) dfs(G[v][i]);
    }
    vs.push_back(v);
}
void rdfs(int v,int k){
    used[v]=1;
    cmp[v]=k;
    for(unsigned int i=0;i<rG[v].size();i++){
        if(!used[rG[v][i]]) rdfs(rG[v][i],k);
    }
}
int scc(){
    memset(used,0,sizeof(used));
    vs.clear();
    for(int v=0;v<V;v++){
        if(!used[v]) dfs(v);
    }
    memset(used,0,sizeof(used));
    int sum=0;
    for(int i=vs.size()-1;i>=0;i--){
        if(!used[vs[i]]) rdfs(vs[i],sum++);
    }
    return sum;
}
ll x[maxn],y[maxn],r[maxn],c[maxn];
void slove(){
    for(int i=0;i<V;i++){
        for(int j=0;j<V;j++){
            if(i==j) continue;
            if(r[i]*r[i]>=(y[i]-y[j])*(y[i]-y[j])+(x[i]-x[j])*(x[i]-x[j])) add_edge(i,j),vis[i][j]=1;
        }
    }
}
int main(){
    int T,n,cas=1;
    scanf("%d",&T);
    while(T--){
        scanf("%d",&V);
        for(int i=0;i<maxn;i++) G[i].clear(),rG[i].clear(),ans[i].clear();
        for(int i=0;i<V;i++) for(int j=0;j<V;j++) vis[i][j]=0;
        memset(in,0,sizeof(in));
        memset(fvis,0,sizeof(fvis));
        for(int i=0;i<V;i++) scanf("%I64d%I64d%I64d%I64d",&x[i],&y[i],&r[i],&c[i]);
        slove();
        scc();
        int cnt=0;
        for(int i=0;i<V;i++){
            for(int j=0;j<V;j++){
                if(i==j||vis[i][j]==0) continue;
                if(cmp[i]==cmp[j]) continue;
                in[cmp[j]]++;
            }
        }
        for(int i=0;i<V;i++) ans[cmp[i]].push_back(i);
        for(int i=0;i<V;i++){
            if(fvis[cmp[i]]) continue;
            if(in[cmp[i]]==0){
                ll tmp1=inf;
                fvis[cmp[i]]=1;
                for(unsigned int j=0;j<ans[cmp[i]].size();j++){
                    int ttt=ans[cmp[i]][j];
                    tmp1=min(tmp1,c[ttt]);
                }
                if(tmp1!=inf) cnt+=tmp1;
            }
        }
        printf("Case #%d: %d\n",cas++,cnt);
    }
    return 0;
}