POJ3114 Countries in War (强连通分量 + 缩点 + 最短路径 + 好题)
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <vector>
using namespace std;
const int Max = ;
const int INF = 0x3f3f3f3f;
int n, m, dfs_clock, scc_cnt, scnt;
int g[Max][Max], pre[Max], low[Max], Stack[Max], sccno[Max];
int G[Max][Max];
int head[Max], num;
struct Edge
{
int v, Next;
};
Edge edge[Max * Max];
void addEdge(int u, int v)
{
edge[num].v = v;
edge[num].Next = head[u];
head[u] = num++;
}
void init()
{
memset(head, -, sizeof(head));
memset(pre, , sizeof(pre));
//memset(low, 0, sizeof(low));
memset(sccno, , sizeof(sccno));
scnt = dfs_clock = scc_cnt = num = ;
for (int i = ; i <= n; i++)
for (int j = i; j <= n; j++)
{
if (i == j)
G[i][j] = g[i][j] = ;
else
{
g[i][j] = g[j][i] = INF;
G[i][j] = G[j][i] = INF;
}
}
}
void dfs(int u)
{
pre[u] = low[u] = ++dfs_clock;
Stack[scnt++] = u;
for (int i = head[u]; i != -; i = edge[i].Next)
{
int v = edge[i].v;
if (!pre[v])
{
dfs(v);
low[u] = min(low[u], low[v]);
}
else if (!sccno[v])
low[u] = min(low[u], pre[v]);
}
if (low[u] == pre[u])
{
scc_cnt++;
for (; ;)
{
int x = Stack[--scnt];
sccno[x] = scc_cnt;
if ( x == u)
break;
}
}
}
void find_scc()
{
for (int i = ; i <= n; i++)
{
if (!pre[i])
dfs(i);
}
}
void build_new_graphic()
{
for (int i = ; i <= n; i++)
{
for (int j = ; j <= n; j++)
{
if (i != j && sccno[i] != sccno[j] && g[i][j] != INF) // 不同的连通分量号建立一条有向边。
{
G[ sccno[i] ][ sccno[j] ] = min(g[i][j], G[ sccno[i] ][ sccno[j] ]);
}
}
}
}
int dist[Max], vis[Max];
void dijkstra(int start, int goal)
{
//利用起点start,终点goal来搞,以前做惯了,直接用起点是1来做了
for (int i = ; i <= scc_cnt; i++)
dist[i] = G[start][i];
memset(vis, , sizeof(vis));
dist[start] = ;
vis[start] = ;
for (int i = ; i <= scc_cnt; i++)
{
int minn = INF, pos = ; // 这里初始化pos为1,否则当下面的循环不满足条件是,执行vis[pos]会出错
for (int j = ; j <= scc_cnt; j++)
{
if (!vis[j] && minn > dist[j])
{
minn = dist[j];
pos = j;
}
}
vis[pos] = ;
for (int j = ; j <= scc_cnt; j++)
{
if (!vis[j] && dist[j] > dist[pos] + G[pos][j])
dist[j] = dist[pos] + G[pos][j];
}
}
if (dist[goal] != INF)
printf("%d\n", dist[goal]);
else
printf("Nao e possivel entregar a carta\n");
}
int main()
{
while (scanf("%d%d", &n, &m) != EOF)
{
if (n == && m == )
break;
init();
int u, v, c;
for (int i = ; i <= m; i++)
{
scanf("%d%d%d", &u, &v, &c);
if (g[u][v] > c)
{
g[u][v] = c; // 判断重边
}
addEdge(u, v);
}
find_scc(); // 找强连通分量
//cout << scc_cnt << endl;
build_new_graphic(); // 重新构图 int k;
scanf("%d", &k);
while (k--)
{
scanf("%d%d", &u, &v);
if (sccno[u] == sccno[v]) // 同一连通分量直接输出
printf("0\n");
else
{
dijkstra(sccno[u], sccno[v]);
}
}
printf("\n");
} return ;
}
POJ3114 Countries in War (强连通分量 + 缩点 + 最短路径 + 好题)的更多相关文章
- Countries in War(强连通分量及其缩点)
http://poj.org/problem?id=3114 题意:有n个城市,m条边,由a城市到b城市的通信时间为w,若a城市与b城市连通,b城市与a城市也连通,则a,b城市之间的通信时间为0,求出 ...
- POJ1236Network of Schools[强连通分量|缩点]
Network of Schools Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 16571 Accepted: 65 ...
- POJ1236Network of Schools(强连通分量 + 缩点)
题目链接Network of Schools 参考斌神博客 强连通分量缩点求入度为0的个数和出度为0的分量个数 题目大意:N(2<N<100)各学校之间有单向的网络,每个学校得到一套软件后 ...
- HD2767Proving Equivalences(有向图强连通分量+缩点)
题目链接 题意:有n个节点的图,现在给出了m个边,问最小加多少边是的图是强连通的 分析:首先找到强连通分量,然后把每一个强连通分量缩成一个点,然后就得到了一个DAG.接下来,设有a个节点(每个节点对应 ...
- UVa11324 The Largest Clique(强连通分量+缩点+记忆化搜索)
题目给一张有向图G,要在其传递闭包T(G)上删除若干点,使得留下来的所有点具有单连通性,问最多能留下几个点. 其实这道题在T(G)上的连通性等同于在G上的连通性,所以考虑G就行了. 那么问题就简单了, ...
- ZOJ3795 Grouping(强连通分量+缩点+记忆化搜索)
题目给一张有向图,要把点分组,问最少要几个组使得同组内的任意两点不连通. 首先考虑找出强连通分量缩点后形成DAG,强连通分量内的点肯定各自一组,两个强连通分量的拓扑序能确定的也得各自一组. 能在同一组 ...
- POJ2553 The Bottom of a Graph(强连通分量+缩点)
题目是问,一个有向图有多少个点v满足∀w∈V:(v→w)⇒(w→v). 把图的强连通分量缩点,那么答案显然就是所有出度为0的点. 用Tarjan找强连通分量: #include<cstdio&g ...
- uva 11324 The Largest Clique(强连通分量缩点+DAG动态规划)
http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=25&page=sh ...
- poj 2762 Going from u to v or from v to u?(强连通分量+缩点重构图+拓扑排序)
http://poj.org/problem?id=2762 Going from u to v or from v to u? Time Limit: 2000MS Memory Limit: ...
随机推荐
- 再记一次w3wp占用CPU过高的解决过程(Dictionary和线程安全)
在此之前项目有发生过两次类似的状况,都得以解决,但最近又会发现偶尔CPU会跑满,虽然之前使用过WinDbg解决过两次问题但人的记忆是不可靠的,今天处理同样问题的时候还是遇到了一些障碍,这一次希望可以记 ...
- 【转】赶集网mysql开发36军规
原文链接:https://segmentfault.com/a/1190000004711147 写在前面的话:总是在灾难发生后,才想起容灾的重要性:总是在吃过亏后,才记得曾经有人提醒过. (一)核心 ...
- Beta版本冲刺———第四天
会议照片: 项目燃尽图: 1.项目进展: 今天解决的进度:新增加了一个撤销按钮,实现对上一步操作的撤销. 仍在进行对排行榜分数变更的实现. 2.每个人每天做的事情 郭怡锋:汇总工作进度,对此总结,进行 ...
- SVN常见图标的含义
项目视图 The Package Explorer view - 已忽略版本控制的文件.可以通过Window → Preferences → Team → Ignored Resources.来忽 ...
- 揭秘PHP匿名函数
揭秘PHP匿名函数 定义:匿名函数就是没有名字的函数. 有2种形式的匿名函数: 形式1:将一个匿名函数"赋值"给一个变量--此时该变量就代表该匿名函数了! 形式2: 是直接将一个匿 ...
- 在Word2013中多次应用格式刷
顾名思义,格式刷是为了方便需要跨区域操作时候,能快速的应用格式到相应文本.那么怎么使用word进行格式刷的多次使用呢.我们先来看单次的,这个比较容易,只要在先需要的格式单击一次格式刷,再到需要的文本执 ...
- jQuery基础--样式篇(3)
1.jQuiery对象与DOM对象 对于刚刚接触jQuery的初学者,我们要清楚认识一点:jQuery对象与DOM对象是不一样的.可能一时半会分不清楚哪些是jQuery对象,哪些是DOM对象,下面重点 ...
- [转]理解JavaScript中的事件处理
什么是事件? 事件(Event)是JavaScript应用跳动的心脏 ,也是把所有东西粘在一起的胶水.当我们与浏览器中 Web 页面进行某些类型的交互时,事件就发生了.事件可能是用户在某些内容上的点击 ...
- Android获取屏幕宽度、高度的4种方法
记录学习之用,有相同的问题可以参考 方法一: WindowManager wm = (WindowManager) this .getSystemService(Context.WINDOW_SERV ...
- ES6新特性:Proxy代理器
ES6新特性:Proxy: 要使用的话, 直接在浏览器中执行即可, node和babel目前还没有Proxy的polyfill;,要使用的话,直接在浏览器中运行就好了, 浏览器的兼容性为:chrome ...