HDU 5934-Bomb(强连通缩点)

来源:互联网 发布:淘宝怎么选择发货地 编辑:程序博客网 时间:2024/05/29 10:52


题意:给你n个炸弹的坐标和爆炸半径,已经引爆每个炸弹的花费,若其他炸弹在某炸弹的爆炸范围内,

则其他炸弹也会被引爆,并之后牵连其他炸弹,问你最少花费多少代价可以引爆所有炸弹。

题解:我们可以考虑将所有能牵连的炸弹之间建边,跑一发强连通,然后将每个强连通分量缩成

一个点(因为互相都能炸到啊),中间保存这个强连通分量的最小花费。对于缩成的点我们考虑

建立新图,当两个炸弹能够互相影响到并且不属于一个强两桶分量的话,对其建边,然后找入度为0

的点开始累加答案,记住,一定要找当前花费最少的入度为0的点。。。。

#include<stack>#include<vector>#include<stdio.h>#include<string.h>#include<algorithm>using namespace std;typedef long long ll;struct node{ll x,y,c,r;}a[1005];stack<int>s;vector<int>q[1005],q1[1005];   ll dis(node a,node b){     return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);}                  int in[1005],dfn[1005],low[1005],t,cnt;int vis[1005],num,belong[1005],cost[1005];void tarjan(int u){int i,now;dfn[u]=low[u]=++t;s.push(u);in[u]=1;vis[u]=1;for(i=0;i<q[u].size();i++){int v=q[u][i];if(vis[v]==0){tarjan(v);low[u]=min(low[u],low[v]);continue;}if(belong[v]==0)low[u]=min(low[u],dfn[v]);}if(dfn[u]==low[u]){num++;while(s.empty()==0){            now=s.top();              s.pop();              in[u]=0;              belong[now]=num;  cost[num]=min(cost[num],(int)a[now].c);            if(now==u)                  break;  }}}void dfs(int x){vis[x]=1;for(int i=0;i<q1[x].size();i++){int v=q1[x][i];if(vis[v]==0)dfs(v);}}int main(void){int T,n,i,j,cases=0;scanf("%d",&T);while(T--){num=t=0;memset(in,0,sizeof(in));memset(vis,0,sizeof(vis));memset(dfn,0,sizeof(dfn));memset(low,0,sizeof(low));memset(cost,127,sizeof(cost));memset(belong,0,sizeof(belong));scanf("%d",&n);for(i=1;i<=n;i++)q[i].clear(),q1[i].clear();for(i=1;i<=n;i++)scanf("%lld%lld%lld%lld",&a[i].x,&a[i].y,&a[i].r,&a[i].c);for(i=1;i<=n;i++)for(j=1;j<=n;j++){if(i==j)continue;ll d=dis(a[i],a[j]);if(d<=a[i].r*a[i].r)q[i].push_back(j);}for(i=1;i<=n;i++)if(vis[i]==0)  tarjan(i);  for(i=1;i<=n;i++)  for(j=0;j<q[i].size();j++){int v=q[i][j];if(belong[i]!=belong[v])q1[belong[i]].push_back(belong[v]);}ll ans=0;memset(vis,0,sizeof(vis));for(i=num;i>0;i--)if(vis[i]==0){ans+=cost[i];dfs(i);}printf("Case #%d: %lld\n",++cases,ans);}return 0;}