HDU-5934-scc

http://acm.hdu.edu.cn/showproblem.php?pid=5934

题目大意：

有n个炸弹，给出他们的坐标(x,y)以及伤害半径r和引爆代价c。

在炸弹爆炸半径以内以及边缘上的炸弹会被连锁引爆。求最少需要多少的代价引爆所有的炸弹。  

思路：

求scc缩点，然后产生若干个DAG，找到每个DAG得起点，对每个起点内的所有点找一个最小的代价点贡献给ans，累加求解就是答案，

代码：

#include<bits/stdc++.h>
using namespace std;
#define pii pair<int,int>
#define mp make_pair
const int maxn=1005;
int T,N;
long long x[maxn],y[maxn],r[maxn],c[maxn];
vector<int>g[maxn];
int sum,scc_cnt;
int low[maxn],dfn[maxn],scc[maxn],in[maxn];
stack<int>S;
bool ok(int i,int j){
    return (y[i]-y[j])*(y[i]-y[j])+(x[i]-x[j])*(x[i]-x[j])<=r[i]*r[i];
}
void dfs(int u){
    low[u]=dfn[u]=++sum;
    S.push(u);
    for(int v:g[u]){
        if(!dfn[v]){
            dfs(v);
            low[u]=min(low[u],low[v]);
        }else if(!scc[v]) low[u]=min(low[u],dfn[v]);
    }
    if(low[u]==dfn[u]){
        ++scc_cnt;
        int x=-1;
        while(x!=u){
            x=S.top();
            S.pop();
            scc[x]=scc_cnt;
        }
    }
}
int main(){ 
    scanf("%d",&T);
    for(int cas=1;cas<=T;++cas){
        cin>>N;
        sum=scc_cnt=0;
        while(!S.empty())S.pop();
        for(int i=1;i<=N;++i)
        low[i]=dfn[i]=scc[i]=0,
        g[i].clear(),scanf("%lld%lld%lld%lld",&x[i],&y[i],&r[i],&c[i]);
        for(int i=1;i<=N;++i){
            for(int j=1;j<=N;++j){
                if(i==j)continue;
                if(ok(i,j))g[i].push_back(j);
            }
        }
        for(int i=1;i<=N;++i){
            if(!dfn[i])dfs(i);
        }//puts("---");
        for(int i=1;i<=scc_cnt;++i)in[i]=0;
        for(int i=1;i<=N;++i){
            for(int j:g[i]){
                if(scc[i]==scc[j])continue;
                in[scc[j]]++;
            }
        }
        long long ans=0;
        for(int i=1;i<=scc_cnt;++i){
            if(in[i]==0){
                long long minc=99999999;
                for(int j=1;j<=N;++j){
                    if(scc[j]==i){
                        minc=min(minc,c[j]);
                    }
                }
                ans+=minc;
            }
        }
        cout<<"Case #"<<cas<<": "<<ans<<endl;
    }
    return 0;
}

-