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求 ...
随机推荐
- JavaScript设计模式(一)
使用JavaScript框架和库过程中, 我遇到过很多感觉上'奇形怪状'的代码. 大多数情况下, 按照惯例编写代码也能够写出很多出色的功能. 但是如果不从根本上理解它们实现的方法, 就没办法完全充分发 ...
- 通过GitHub部署网站到Azure WebSite
一.将代码上传到GitHub存储库中 第一步:创建Github新账户 第二步:新建仓库 第三步:安装Github shell程序,地址:http://windows.github.com/ 第四步:打 ...
- Android 使用BaseAdapter 插入不同类型数据
在使用过程有时要在listview中插入不同类型的数据,比如说position＝0的位置插入,广告,其它列表显示数据的情况. 一定要重写两种方法 @Override public int getIte ...
- JAVA操作LDAP总结
一.LDAP概念 LDAP的全称为Lightweight Directory Access Protocol(轻量级目录访问协议), 基于X.500标准, 支持 TCP/IP. LDAP目录为数据库, ...
- C# ACM poj1004
水题.. public static void acm1004(float[] a) { ; foreach (var item in a) { sum += item; } Console.Writ ...
- 通过正则获取url参数
1.通过正则来获取url地址栏的参数: ---------------------------我是分割线-------------------------------- var reg1=/([^?& ...
- Oracle PL/SQL 异常处理
Oracle数据库中的异常:没有异常的转移,因为没有受检异常和非受检异常得区分. 1.异常的产生: 2.异常的处理: declare --变量定义,初始化赋值. begin --变量的赋值,函数调用, ...
- Codevs 1684 垃圾陷阱
1684 垃圾陷阱 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 黄金 Gold 题目描述 Description 卡门--农夫约翰极其珍视的一条Holsteins奶牛--已经落了 ...
- strcpy函数导致release版程序崩溃
最近在写一个读取模型文件的小程序.很随意的使用了strcpy函数进行char字符数组的拷贝,这个数组是需要传递给PostMessage作为WPARAM的参数.代码部分如下: char pStrCurr ...
- 九度OJ 1024 畅通工程 -- 并查集、贪心算法(最小生成树)
题目地址:http://ac.jobdu.com/problem.php?pid=1024 题目描述: 省政府"畅通工程"的目标是使全省任何两个村庄间都可以实现公路交通(但 ...