HDU Traffic Real Time Query System
题目大意是：对于（n， m）的图，给定边a， b查询从a到b要经过的割点的最少数目。
先tarjan算法求双连通然后缩点，即对于每个割点将周围的每个双连通看成一个点与之相连。然后求解LCA即可，距离dis[u]表示从根出发到u的遍历过程中经过的割顶的数目，利用
tarjan离线算法， 最后答案是：dis[u] + dis[v] - 2*dis[findset(v)] + （findset(v) > bcc_cnt)。注意findset（v） > bcc_cnt表示当LCA（u，v） 为割顶时的判断，此时必须加1。
代码：
#include <iostream>
#include <sstream>
#include <cstdio>
#include <climits>
#include <cstring>
#include <cstdlib>
#include <string>
#include <stack>
#include <map>
#include <cmath>
#include <vector>
#include <queue>
#include <algorithm>
#define esp 1e-6
#define pi acos(-1.0)
#define pb push_back
#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#define mp(a, b) make_pair((a), (b))
#define in freopen("in.txt", "r", stdin);
#define out freopen("out.txt", "w", stdout);
#define print(a) printf("%d\n",(a));
#define bug puts("********))))))");
#define stop system("pause");
#define Rep(i, c) for(__typeof(c.end()) i = c.begin(); i != c.end(); i++)
#define inf 0x0f0f0f0f using namespace std;
typedef long long LL;
typedef vector<int> VI;
typedef pair<int, int> pii;
typedef vector<pii> VII;
typedef vector<pii, int> VIII;
typedef VI:: iterator IT;
#define eid first
#define vtx second const int maxn = + ;
const int maxm = + ;
int pre[maxn], low[maxn], bccno[maxn], ebccno[maxm], iscut[maxn], vis[maxm], ans[maxn], pa[maxn], bcc_cnt, dfs_clock;
VI g[maxn], bcc[maxn];
VII adj[maxn], query[maxn];
struct edge{
int u, v;
} ee[maxm];
stack<int> S; int dfs(int u, int fa)
{
int lowu = pre[u] = ++dfs_clock;
int child = ;
for(int i = ; i < adj[u].size(); i++)
{
int v = adj[u][i].vtx;
if(!pre[v])
{
S.push(adj[u][i].eid);
child++;
vis[adj[u][i].eid] = ;
int lowv = dfs(v, u);
lowu = min(lowu, lowv);
if(lowv >= pre[u])
{
iscut[u] = ;
bcc_cnt++;
bcc[bcc_cnt].clear();
for(;;)
{
int e = S.top(); S.pop();
bcc[bcc_cnt].pb(e);
ebccno[e] = bcc_cnt;
bccno[ee[e].u] = bccno[ee[e].v] = bcc_cnt;
if(e == adj[u][i].eid) break;
}
}
}
else if(pre[v] < pre[u])
{
if(v != fa)
{
lowu = min(lowu, pre[v]);
S.push(adj[u][i].eid);
vis[adj[u][i].eid] = ;
}
else {
if( !vis[adj[u][i].eid])
{
lowu = min(lowu, pre[v]);
S.push(adj[u][i].eid);
vis[adj[u][i].eid] = ;
}
}
}
}
if(fa < && child == )
{
iscut[u] = ;
}
return low[u] = lowu;
}
void find_bcc(int n)
{
memset(pre, , sizeof(pre));
memset(bccno, , sizeof(bccno));
memset(iscut, , sizeof(iscut));
memset(ebccno, , sizeof(ebccno));
memset(vis, , sizeof(vis)); while(!S.empty()) S.pop();
dfs_clock = bcc_cnt = ;
for(int i = ; i < n ;i++)
if(!pre[i])
dfs(i, -);
}
int flag[maxn];
int dis[maxn];
int vv[maxn];
int findset(int x)
{
return pa[x] == x ? x : pa[x] = findset(pa[x]);
}
void tarjan(int u, int d)
{
vv[u] = ;
dis[u] = d;
pa[u] = u;
if(u > bcc_cnt)
dis[u]++;
for(int i = ; i < query[u].size(); i++)
{
int v = query[u][i].first;
int id = query[u][i].second;
if(vis[v])
ans[id] = dis[u] + dis[v] - *dis[findset(v)] + (findset(v) > bcc_cnt);
}
for(int i = ; i < g[u].size(); i++)
{
int v = g[u][i];
if(!vv[v])
{
tarjan(v, dis[u]);
pa[v] = u;
}
}
}
int main(void)
{
int n, m;
while(scanf("%d%d", &n, &m), n||m)
{
for(int i = ; i < maxn; i++)
adj[i].clear(), g[i].clear(), query[i].clear();
memset(flag, -, sizeof(flag));
for(int i = ; i <= m; i++)
{
int u, v;
scanf("%d%d", &u, &v);
u--, v--;
adj[u].pb(mp(i, v));
adj[v].pb(mp(i, u));
ee[i].u = u, ee[i].v = v;
}
find_bcc(n);
int cut_cnt = ;
for(int i = ; i < n; i++)
if(iscut[i])
{
cut_cnt++;
int u = cut_cnt+bcc_cnt;
for(int k = ; k < adj[i].size(); k++)
{
int v = ebccno[adj[i][k].first];
if(flag[v] != u )
{
flag[v] = u;
g[v].pb(u);
g[u].pb(v);
}
}
}
int q;
for(int i = scanf("%d", &q); i <= q; i++)
{
int s, t;
scanf("%d%d", &s, &t);
s = ebccno[s], t = ebccno[t];
query[s].pb(mp(t, i));
query[t].pb(mp(s, i));
}
memset(vv, , sizeof(vv));
for(int i = ; i <= bcc_cnt+cut_cnt; i++)
if(!vv[i])
tarjan(i, );
for(int i = ; i <= q; i++)
printf("%d\n", ans[i]);
}
return ;
}
这个也可用Sparse Table（ST算法）结合RMQ求解LCA
下面是代码：
#include <iostream>
#include <sstream>
#include <cstdio>
#include <climits>
#include <cstring>
#include <cstdlib>
#include <string>
#include <stack>
#include <map>
#include <cmath>
#include <vector>
#include <queue>
#include <algorithm>
#define esp 1e-6
#define pi acos(-1.0)
#define pb push_back
#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#define mp(a, b) make_pair((a), (b))
#define in freopen("in.txt", "r", stdin);
#define out freopen("out.txt", "w", stdout);
#define print(a) printf("%d\n",(a));
#define bug puts("********))))))");
#define stop system("pause");
#define Rep(i, c) for(__typeof(c.end()) i = c.begin(); i != c.end(); i++)
#define inf 0x0f0f0f0f using namespace std;
typedef long long LL;
typedef vector<int> VI;
typedef pair<int, int> pii;
typedef vector<pii> VII;
typedef vector<pii, int> VIII;
typedef VI:: iterator IT;
#define eid first
#define vtx second const int maxn = + ;
const int maxm = + ;
int pre[maxn], low[maxn], bccno[maxn], ebccno[maxm], iscut[maxn], vis[maxm], ans[maxn], pa[maxn], bcc_cnt, dfs_clock;
int dep[maxn], E[maxn], R[maxn];
int dp[maxn*][];
VI g[maxn], bcc[maxn];
VII adj[maxn], query[maxn];
struct edge
{
int u, v;
} ee[maxm];
stack<int> S;
int cnt;
int dfs(int u, int fa)
{
int lowu = pre[u] = ++dfs_clock;
int child = ;
for(int i = ; i < adj[u].size(); i++)
{
int v = adj[u][i].vtx;
if(!pre[v])
{
S.push(adj[u][i].eid);
child++;
vis[adj[u][i].eid] = ;
int lowv = dfs(v, u);
lowu = min(lowu, lowv);
if(lowv >= pre[u])
{
iscut[u] = ;
bcc_cnt++;
bcc[bcc_cnt].clear();
for(;;)
{
int e = S.top();
S.pop();
bcc[bcc_cnt].pb(e);
ebccno[e] = bcc_cnt;
bccno[ee[e].u] = bccno[ee[e].v] = bcc_cnt;
if(e == adj[u][i].eid) break;
}
}
}
else if(pre[v] < pre[u])
{
if(v != fa)
{
lowu = min(lowu, pre[v]);
S.push(adj[u][i].eid);
vis[adj[u][i].eid] = ;
}
else
{
if( !vis[adj[u][i].eid])
{
lowu = min(lowu, pre[v]);
S.push(adj[u][i].eid);
vis[adj[u][i].eid] = ;
}
}
}
}
if(fa < && child == )
{
iscut[u] = ;
}
return low[u] = lowu;
}
void find_bcc(int n)
{
memset(pre, , sizeof(pre));
memset(bccno, , sizeof(bccno));
memset(iscut, , sizeof(iscut));
memset(ebccno, , sizeof(ebccno));
memset(vis, , sizeof(vis)); while(!S.empty()) S.pop();
dfs_clock = bcc_cnt = ;
for(int i = ; i < n ; i++)
if(!pre[i])
dfs(i, -);
}
int flag[maxn];
int dis[maxn];
int vv[maxn];
int findset(int x)
{
return pa[x] == x ? x : pa[x] = findset(pa[x]);
}
void ST(int u, int d, int sum)
{
vv[u] = ;
R[u] = ++cnt;
dep[cnt] = d;
E[cnt] = u;
dis[u] = sum;
if(u > bcc_cnt)
dis[u]++;
for(int i = ; i < g[u].size(); i++)
{
int v = g[u][i];
if(!vv[v])
{
ST(v, d+, dis[u]);
E[++cnt] = u;
dep[cnt] = d;
}
}
}
void Init(void)
{
for(int i = ; i <= cnt; i++)
dp[i][] = i;
for(int j = ; j <= (int)(log(cnt)/log(2.0)); j++)
for(int i = ; i + (<<j) - <= cnt; i++)
{
if(dep[dp[i][j-]] < dep[dp[i+(<<(j-))][j-]])
dp[i][j] = dp[i][j-];
else dp[i][j] = dp[i + (<<(j-))][j-];
}
}
int RMQ(int u, int v)
{
int k = (int)(log(v - u + )/log(2.0));
if(dep[dp[u][k]] < dep[dp[v - (<<k) + ][k]])
return dp[u][k];
return dp[v - (<<k) + ][k];
}
int main(void)
{ int n, m;
while(scanf("%d%d", &n, &m), n||m)
{
for(int i = ; i < maxn; i++)
adj[i].clear(), g[i].clear(), query[i].clear();
memset(flag, -, sizeof(flag));
for(int i = ; i <= m; i++)
{
int u, v;
scanf("%d%d", &u, &v);
u--, v--;
adj[u].pb(mp(i, v));
adj[v].pb(mp(i, u));
ee[i].u = u, ee[i].v = v;
}
find_bcc(n);
int cut_cnt = ;
for(int i = ; i < n; i++)
if(iscut[i])
{
cut_cnt++;
int u = cut_cnt+bcc_cnt;
for(int k = ; k < adj[i].size(); k++)
{
int v = ebccno[adj[i][k].first];
if(flag[v] != u )
{
flag[v] = u;
g[v].pb(u);
g[u].pb(v);
}
}
}
memset(vv, , sizeof(vv));
cnt = ;
for(int i = ; i <= bcc_cnt+cut_cnt; i++)
if(!vv[i])
ST(i, , );
int q;
Init();//记得初始化啊 。。。。。。。。
for(int i = scanf("%d", &q); i <= q; i++)
{
int s, t;
scanf("%d%d", &s, &t);
s = ebccno[s], t = ebccno[t];
int ss , tt;
ss = min(R[s], R[t]);
tt = max(R[s], R[t]);
int lca = RMQ(ss, tt);
lca = E[lca];
int ans = dis[s] + dis[t] - * dis[lca] + (lca > bcc_cnt);
printf("%d\n", ans);
}
}
return ;
}
两者时间上差不多，而且ST算法需要初始化，更容易出错
HDU Traffic Real Time Query System的更多相关文章
- HDU 3686 Traffic Real Time Query System (图论)
HDU 3686 Traffic Real Time Query System 题目大意 给一个N个点M条边的无向图,然后有Q个询问X,Y,问第X边到第Y边必需要经过的点有多少个. solution ...
- CH#24C 逃不掉的路 和 HDU3686 Traffic Real Time Query System
逃不掉的路 CH Round #24 - 三体杯 Round #1 题目描述 现代社会,路是必不可少的.任意两个城镇都有路相连,而且往往不止一条.但有些路连年被各种XXOO,走着很不爽.按理说条条大路 ...
- HDU 3686 Traffic Real Time Query System（双连通分量缩点+LCA）（2010 Asia Hangzhou Regional Contest）
Problem Description City C is really a nightmare of all drivers for its traffic jams. To solve the t ...
- HDU3686 Traffic Real Time Query System 题解
题目 City C is really a nightmare of all drivers for its traffic jams. To solve the traffic problem, t ...
- Traffic Real Time Query System 圆方树+LCA
题目描述 City C is really a nightmare of all drivers for its traffic jams. To solve the traffic problem, ...
- hdu 3686 Traffic Real Time Query System 点双两通分量 + LCA。这题有重边！！！
http://acm.hdu.edu.cn/showproblem.php?pid=3686 我要把这题记录下来. 一直wa. 自己生成数据都是AC的.现在还是wa.留坑. 我感觉我现在倒下去床上就能 ...
- Traffic Real Time Query System HDU - 3686
https://vjudge.net/problem/HDU-3686 点双啊,就是在求割顶的时候,另外用一个栈来存一些边 在遍历u点出发的边时,遇到树边或反向边(u,v)就把此边加入栈(可能要记一下 ...
- HDU 3686 Traffic Real Time Query System（点双连通）
题意 给定一张 \(n\) 个点 \(m\) 条边的无向图,\(q\) 次询问,每次询问两边之间的必经之点个数. 思路 求两点之间必经之边的个数用的是边双缩点,再求树上距离.而对比边双和点双之 ...
- UVALive-4839 HDU-3686 Traffic Real Time Query System 题解
题目大意: 有一张无向连通图,问从一条边走到另一条边必定要经过的点有几个. 思路: 先用tarjan将双连通分量都并起来,剩下的再将割点独立出来,建成一棵树,之后记录每个点到根有几个割点,再用RMQ求 ...
随机推荐
- 【算法】A*改进算法
目的:我这里希望实现一个java A* 游戏里的战斗寻径 定义部分: 这个定义引用自 http://www.cnblogs.com/kanego/archive/2011/08/30/2159070. ...
- vs2010开发android的准备工作
安装 Mono for Android for Visual Studio 2010 需要下面4个步骤: 安装 JDK 安装 Android SDK 配置模拟器 安装 Mono for Android ...
- Dr.com校园网客户端故障解决方法
一,登录客户端的时候提示“登录超时失败” 解决办法:单击桌面右下角的联网图标然后右键,打开网络和共享中心,在左侧导航栏中找到更改网络适配器,找到本地连接(或者无线连接),点击右键找到详细信息,看看ip ...
- Android系统简介(上):历史渊源
上个月,看到微信的一系列文章,讲到Linux的鼻祖-李纳斯的传记<Just for Fun>, 其人神乎其能, 其人生过程非常有趣,值得每个程序员细细品味. 而实际上,对我而已,虽然做软件 ...
- iOS App Transport Security
网络请求提示:Application Transport Security has blocked a cleartext HTTP (http://) resource load since it ...
- iOS夯实：ARC时代的内存管理
iOS夯实:ARC时代的内存管理 文章转自 ARC时代的内存管理 什么是ARC Automatic Reference Counting (ARC) is a compiler feature tha ...
- iOS真机调试之我见
入职20多天,以前一直以为iOS真机调试是多么复杂的事情,但在公司大牛的帮助下:终于理清头绪,由于公司证书已申请,文章中免不了旁征博引. 1.首先,得有苹果开发者账号,如果在公司,公司会提供:不在公 ...
- asp.net 备份和恢复数据库
观看了Insus的视频写下来的,代码可能有点冗长,如有好的想法的,可以多交流. 前台: <form id="form1" runat="server"&g ...
- JavaScript学习笔记(6)——JavaScript语法之对象
JavaScript 中的所有事物都是对象:字符串.数字.数组.日期,等等. 在 JavaScript 中,对象是拥有属性和方法的数据. 属性是与对象相关的值. 方法是能够在对象上执行的动作. 提示: ...
- [Guava官方文档翻译] 2.使用和避免使用null (Using And Avoiding Null Explained)
本文地址:http://www.cnblogs.com/hamhog/p/3536647.html "null很恶心." -Doug Lea "这是一个令我追悔莫及的错误 ...