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求 ...
随机推荐
- 抽象（abstract）升级版变接口（interface） 继承（extends）升级版叫实现（implements） 升级版啊升级版 接口可以多继承
Client --------------------------------------------------- public class Client{ public static void m ...
- RabbitMQ 原文译04--路由
在前一篇文章中我们构建了一个简单的日志系统,我们可以向多个接受者广播消息. 在这篇文章我,我们将要添加一些功能使得针对部分消息的接受成为可能,例如我们只对错误的消息进行磁盘记录,同时又可以把所有的消息 ...
- JAVA的StringBuffer类(转载整理)____非常重要的一个类，线程安全，不用每次创建一个对象，以及和String的区别
核心部分转载自:http://www.cnblogs.com/springcsc/archive/2009/12/03/1616330.html StringBuffer类和String一样,也用来代 ...
- WebClient.UploadValues Post中文乱码的解决方法
//using (System.Net.WebClient wc = new System.Net.WebClient()) //{ // wc.Encoding = Encoding.GetEnco ...
- CPrintDialog
CPrintDialog 封装windows为打印提供服务的通用窗体. BOOL GetDefaults(); //获取默认设备,不显示对话框 // Helpers for parsing infor ...
- Java实战之03Spring-02Spring的核心之IoC
二.Spring的核心之IoC(Inversion of Control 控制反转) 1.IoC是什么? 回顾:连接池原理中的使用. 注意:我们在连接池原理中是使用构造函数注入的,当然还可以通过提供s ...
- 11_Jaxws常用注解
[不使用注解] 默认namespace是服务类包名的倒序 默认portType是服务类的类名 ............... 注解的所起的作用: Jaxws提供的注解可以对WebService的接口规 ...
- 回溯算法之n皇后问题
今天在看深度优先算法的时候,联想到DFS本质不就是一个递归回溯算法问题,只不过它是应用在图论上的.OK,写下这篇博文也是为了回顾一下回溯算法设计吧. 学习回溯算法问题,最为经典的问题我想应该就是八皇后 ...
- 在windows下装2个mysql数据库的办法
如下 记录下以免找不到 http://blog.chinaunix.net/uid-77311-id-3450734.html 然后接下来是 1045的解决方案 可以在随便中找
- VS2010 EntityFramework Database First
本文演练介绍如何使用实体框架进行 Database First 开发.通过 Database First,可以从现有数据库对模型进行反向工程处理.模型存储在一个 EDMX 文件(扩展名为 .edmx) ...