tarjan 视频讲解

/**
* 题目链接:https://vjudge.net/problem/HDU-5934
* 题意:给你n个炸弹,引爆每个炸弹会有一定的花费。每个炸弹给出坐标x,y,半径r,引爆花费;
* 引爆一个炸弹会把范围内的炸弹引爆,连锁反应。 现在想把所有炸弹引爆的最小花费。
*
* 解题思路:强连通缩点。根据a能够引爆b,可以在建一条a到b的单向边。如果是一个强连通(这一部分的图,
* 任意两点都可以相互到达)那么就把这个强连通分量变成一个点,值最分量的最小值。这样图就变成有向无环图了。
* 考虑到每个点的花费都是大于0的,所以引爆开始点最划算,即为入度为0的点。
*
* 前置技能 tarjan 缩点。
*/ #include <bits/stdc++.h>
using namespace std; const int maxn=1000+10;
const int INF=2e9+1e8;
vector<int>E[maxn];
struct Point
{
int x,y,r,cost;
}boom[maxn];
bool judge(Point a,Point b)
{
if( 1ll*(a.x-b.x)*(a.x-b.x)+1ll*(a.y-b.y)*(a.y-b.y)<=1ll*a.r*a.r ) return true;
return false;
}
int dfn[maxn],low[maxn],id,vis[maxn],ans,deg[maxn];
int num[maxn],cnt,cost[maxn];//对点进行重新编号,(数组num),按照联通分量进行编号
stack<int>S;
void init()
{
id=cnt=0;
memset(deg,0,sizeof(deg));
memset(num,0,sizeof(num));
memset(vis,0,sizeof(vis));
memset(dfn,0,sizeof(dfn));
}
void tarjan(int x)
{
low[x]=dfn[x]=++id;
S.push(x);
vis[x]=1;
for(int i=0;i<(int)E[x].size();i++)
{
int to=E[x][i];
if(!dfn[to])
{
tarjan(to);
low[x]=min(low[x],low[to]);
}
else if(vis[to]) low[x]=min(low[x],dfn[to]);
}
if(low[x]==dfn[x])
{
int mincost=INF,in=0;
cnt++;
while(1)
{
int now=S.top();
S.pop();
vis[now]=0;
num[now]=cnt;
mincost=min(mincost,boom[now].cost);
if(now==x) break;
}
cost[cnt]=mincost;
}
}
int main()
{
int T,cas=1;
scanf("%d",&T);
while(T--)
{
int n;
scanf("%d",&n);
init();
for(int i=1;i<=n;i++)
{
E[i].clear();
scanf("%d%d%d%d",&boom[i].x,&boom[i].y,&boom[i].r,&boom[i].cost);
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
if(i==j) continue;
if(judge(boom[i],boom[j])) E[i].push_back(j);
}
}
for(int i=1;i<=n;i++)
if(!dfn[i]) tarjan(i);
for(int i=1;i<=n;i++)
{
for(int j=0;j<(int)E[i].size();j++)
{
int to=E[i][j];
if(num[i]!=num[to]) deg[num[to]]++;
}
}
ans=0;
for(int i=1;i<=cnt;i++) if(deg[i]==0)
ans+=cost[i];
printf("Case #%d: %d\n",cas++,ans);
}
return 0;
}

最新文章

  1. Singleton in C++11 style
  2. Android设置AlertDialog点击按钮对话框不关闭(转)
  3. IEEE754测试-软件
  4. I Count Two Three---hdu5878(打表+二分)
  5. PL/SQL中LOOP循环控制语句
  6. 一起啃PRML - 1.2.4 The Gaussian distribution 高斯分布 正态分布
  7. fragment的实现与互相通信
  8. HTML5 JavaScript API
  9. javaEE-string家族三大流氓
  10. 解析JSON的两种方法eval()和JSON.parse()
  11. python3+arcface2.0 离线人脸识别 demo
  12. MongoDB系列:四、spring整合mongodb,带用户验证
  13. 利用airbase-ng建立Soft AP
  14. Yii2的mongodb的聚合操作
  15. nodejs中async使用
  16. hive vs hbase
  17. 公告栏添加时钟——利用canvas画出一个时钟
  18. postgresql----几何类型和函数
  19. 小刘的机器学习---SVM
  20. MySQL ACID及四种隔离级别的解释

热门文章

  1. 在Ubuntu下编译FFmpeg
  2. (转) 对svn分支合并类型和深度的理解
  3. PHP基础知识学习总结
  4. Android Studio报Error:Execution failed for task &amp;#39;:Companion:preDexDebug&amp;#39;.
  5. Unity3d NGUI 360度旋转
  6. css3中font-face属性的用法详解
  7. Spring 和 filter
  8. 【BZOJ2790】[Poi2012]Distance 筛素数+调和级数
  9. spring boot ajax post 前后端
  10. ubuntu 16.04 jenkins pipline的实现 最终docker启动服务