题目描述

在 H 国的小 w 决定到从城市 \(u\) 到城市 \(v\) 旅行,但是此时小 c 由于各种原因不在城市 \(u\),但是小 c 决定到在中途与小 w 相遇

由于 H 国道路的原因,小 w 从城市 \(u\) 到城市 \(v\) 的路线不是固定的,为了合理分配时间,小 c 想知道从城市 \(u\) 到城市 \(v\) 有多少个城市小 w 一定会经过,特别地,\(u,v\) 也必须被算进去,也就是说无论如何答案不会小于 2

由于各种特殊的原因,小 c 并不知道小 w 的起点和终点,但是小 c 知道小 w 的起点和终点只有 \(M\) 种可能,所以对于这 \(M\) 种可能,小 c 都想知道小 w 一定会经过的城市数

H 国所有的边都是无向边,两个城市之间最多只有一条道路直接相连,没有一条道路连接相同的一个城市

任何时候,H 国不存在城市 \(u\) 和城市 \(v\) 满足从 \(u\) 无法到达 \(v\)

输入输出格式

输入格式:

输入第 1 行两个正整数 \(N,E\),表示 H 国的城市数,以及道路数

输入第 2 行至第 \(E+1\) 行,每行两个不同的正整数 \(u,v\),表示城市 \(u\) 到城市 \(v\) 之间有一条边

输入第 \(E+2\) 行一个正整数 \(M\)

输入第 \(E+3\) 行到第 \(E+M+2\) 行每行两个正整数 \(u,v\) 表示小 w 旅行的一种可能的路线

输出格式:

输出共 \(M\) 行,每行一个正整数

输入输出样例

输入样例#1:

5 6
1 2
1 3
2 3
3 4
4 5
3 5
1
1 5

输出样例#1:

3

说明

从城市 1 到城市 5 总共有 4 种可能 :

1 -> 2 -> 3 -> 4 -> 5

1 -> 2 -> 3 -> 5

1 -> 3 -> 4 -> 5

1 -> 3 -> 5

可以发现小 w 总会经过城市 1, 3, 5,所以答案为 3

你可以认为小 w 不会经过相同的城市两次,当然,如果你认为可以经过相同的城市两次也不会影响答案

subtask1 : 15分,\(N=5,M=50\)

subtask2 : 15分,\(N=100,M=5000\)

subtask3 : 20分,\(N=3000,M=500000\)

subtask4 : 20分,\(N=499999,M=500000,E=N−1\)

subtask5 : 30分,\(N=500000,M=500000\)

对于所有数据 : \(1\leq N\leq 500000, 1\leq M\leq 500000, 1\leq E\leq \min(\frac{N(N-1)}{2}, 1000000)\)

题解

建出圆方树,题目所求的就是 \(u\) 到 \(v\) 的路径上有多少圆点

树上差分一下就好了

#include<bits/stdc++.h>
#define ui unsigned int
#define ll long long
#define db double
#define ld long double
#define ull unsigned long long
const int MAXN=1000000+10;
int n,k,m,e,cnt,val[MAXN],beg[MAXN],nex[MAXN<<1],to[MAXN<<1],out[MAXN<<1],DFN[MAXN],LOW[MAXN],Visit_Num,vis[MAXN],d[MAXN],dep[MAXN],Jie[20][MAXN],s[MAXN<<1],top;
std::vector<int> point[MAXN];
namespace IO
{
const ui Buffsize=1<<15,Output=1<<23;
static char Ch[Buffsize],*S=Ch,*T=Ch;
inline char getc()
{
return((S==T)&&(T=(S=Ch)+fread(Ch,1,Buffsize,stdin),S==T)?0:*S++);
}
static char Out[Output],*nowps=Out;
inline void flush(){fwrite(Out,1,nowps-Out,stdout);nowps=Out;}
template<typename T>inline void read(T&x)
{
x=0;static char ch;T f=1;
for(ch=getc();!isdigit(ch);ch=getc())if(ch=='-')f=-1;
for(;isdigit(ch);ch=getc())x=(x<<3)+(x<<1)+(ch^'0');
x*=f;
}
template<typename T>inline void write(T x,char ch='\n')
{
if(!x)*nowps++='0';
if(x<0)*nowps++='-',x=-x;
static ui sta[111],tp;
for(tp=0;x;x/=10)sta[++tp]=x%10;
for(;tp;*nowps++=sta[tp--]^48);
*nowps++=ch;
}
}
using namespace IO;
template<typename T> inline void chkmin(T &x,T y){x=(y<x?y:x);}
template<typename T> inline void chkmax(T &x,T y){x=(y>x?y:x);}
template<typename T> inline T min(T x,T y){return x<y?x:y;}
template<typename T> inline T max(T x,T y){return x>y?x:y;}
inline void insert(int x,int y)
{
to[++e]=y;
nex[e]=beg[x];
out[e]=x;
beg[x]=e;
}
inline void Tarjan(int x,int f)
{
DFN[x]=LOW[x]=++Visit_Num;
for(register int i=beg[x];i;i=nex[i])
if(to[i]==f)continue;
else if(!DFN[to[i]])
{
s[++top]=i;
Tarjan(to[i],x);
chkmin(LOW[x],LOW[to[i]]);
if(LOW[to[i]]>=DFN[x])
{
static int temp;++cnt;
do{
temp=s[top--];
if(vis[out[temp]]!=cnt)
{
vis[out[temp]]=cnt;
point[cnt].push_back(out[temp]);
}
if(vis[to[temp]]!=cnt)
{
vis[to[temp]]=cnt;
point[cnt].push_back(to[temp]);
}
}while(out[temp]!=x||to[temp]!=to[i]);
}
}
else if(DFN[to[i]]<DFN[x])s[++top]=i,chkmin(LOW[x],DFN[to[i]]);
}
inline void dfs(int x,int f)
{
d[x]=d[f]+val[x];Jie[0][x]=f;dep[x]=dep[f]+1;
for(register int i=beg[x];i;i=nex[i])
if(to[i]!=f)dfs(to[i],x);
}
inline int LCA(int u,int v)
{
if(dep[u]<dep[v])std::swap(u,v);
int tmp=dep[u]-dep[v];
for(register int i=19;i>=0;--i)
if(tmp>>i&1)u=Jie[i][u];
if(u==v)return u;
for(register int i=19;i>=0;--i)
if(Jie[i][u]^Jie[i][v])u=Jie[i][u],v=Jie[i][v];
return Jie[0][u];
}
int main()
{
read(n);read(k);
for(register int i=1;i<=k;++i)
{
int u,v;read(u);read(v);
insert(u,v);insert(v,u);
}
Tarjan(1,0);e=0;
for(register int i=1;i<=n;++i)val[i]=1,beg[i]=0;
for(register int i=1;i<=cnt;++i)
for(register int j=0,lt=point[i].size();j<lt;++j)insert(i+n,point[i][j]),insert(point[i][j],i+n);
dfs(1,0);
for(register int j=1;j<20;++j)
for(register int i=1;i<=n+cnt;++i)Jie[j][i]=Jie[j-1][Jie[j-1][i]];
read(m);
static int u,v,lca;
while(m--)
{
read(u);read(v);lca=LCA(u,v);
write(d[u]+d[v]-d[lca]-d[Jie[0][lca]],'\n');
}
flush();
return 0;
}

【刷题】洛谷 P4320 道路相遇的更多相关文章

  1. [洛谷P4320]道路相遇

    题目大意:基本同上一题[bzoj5329][Sdoi2018]战略游戏,只是每个点集内只有两个点,且只有一组询问而已.(双倍经验?我反正就直接改了一下代码就交了) 题解:同上一题(链接见“题目大意”) ...

  2. 【题解】洛谷P1070 道路游戏(线性DP)

    次元传送门:洛谷P1070 思路 一开始以为要用什么玄学优化 没想到O3就可以过了 我们只需要设f[i]为到时间i时的最多金币 需要倒着推回去 即当前值可以从某个点来 那么状态转移方程为: f[i]= ...

  3. 洛谷 P3905 道路重建

    题目描述 从前,在一个王国中,在n个城市间有m条道路连接,而且任意两个城市之间至多有一条道路直接相连.在经过一次严重的战争之后,有d条道路被破坏了.国王想要修复国家的道路系统,现在有两个重要城市A和B ...

  4. 2018.10.30 一题 洛谷4660/bzoj1168 [BalticOI 2008]手套——思路!问题转化与抽象!+单调栈

    题目:https://www.luogu.org/problemnew/show/P4660 https://www.lydsy.com/JudgeOnline/problem.php?id=1168 ...

  5. 【洛谷 P4320】 道路相遇 (圆方树,LCA)

    题目链接 题意:给一张无向图和\(M\)个询问,问\(u,v\)之间的路径的必经之点的个数. 对图建出圆方树,然后必经之点就是两点路径经过的原点个数,用\((dep[u]+dep[v]-dep[LCA ...

  6. Solution -「洛谷 P4320」道路相遇

    \(\mathcal{Description}\)   Link.   给定一个 \(n\) 个点 \(m\) 条边的连通无向图,并给出 \(q\) 个点对 \((u,v)\),询问 \(u\) 到 ...

  7. 洛谷 P1070 道路游戏(noip 2009 普及组 第四题)

    题目描述 小新正在玩一个简单的电脑游戏. 游戏中有一条环形马路,马路上有 nn个机器人工厂,两个相邻机器人工厂之间由一小段马路连接.小新以某个机器人工厂为起点,按顺时针顺序依次将这 nn个机器人工厂编 ...

  8. 洛谷P1070 道路游戏(dp+优先队列优化)

    题目链接:传送门 题目大意: 有N条相连的环形道路.在1-M的时间内每条路上都会出现不同数量的金币(j时刻i工厂出现的金币数量为val[i][j]).每条路的起点处都有一个工厂,总共N个. 可以从任意 ...

  9. 洛谷 P1070 道路游戏 解题报告

    P1070 道路游戏 题目描述 小新正在玩一个简单的电脑游戏. 游戏中有一条环形马路,马路上有\(n\)个机器人工厂,两个相邻机器人工厂之间由一小段马路连接.小新以某个机器人工厂为起点,按顺时针顺序依 ...

随机推荐

  1. 在main函数前后执行的函数之 C语言

    在gcc中,可以使用attribute关键字,声明constructor和destructor,来指定了函数在main之前或之后运行,代码如下: #include <stdio.h> __ ...

  2. Android 网络编程之---HttpClient 与 HttpURLConnection 共用cookie

    HttpClient 与 HttpURLConnection 共用 SessionId HttpClient 与 HttpUrlConnection 是Android 中HTTP操作最常见的訪问方式. ...

  3. JSON传输数组的基本操作

    目标JSON结果如下: 生成JSON的的过程如下: Document document; Document::AllocatorType& allocator = document.GetAl ...

  4. centos 6.5 搭建开源堡垒机 Teleport 遇到的问题解决

    几款开源的堡垒机 下面进行 teleport 的安装: https://docs.tp4a.com/install/#11 异常1:libc.so.6: version `GLIBC_2.14' no ...

  5. datatables.js 简单使用--弹出编辑框或添加数据框

    内容:选中某一条记录,弹出编辑框 环境:asp.net mvc ,  bootstrap 显示效果: 代码: 至于怎么弄多选框,在上一篇博客里已经有说明. 主要用到了bootstrap的模态窗,下面代 ...

  6. 微服务监控zipkin、skywalking以及日志ELK监控系列

    0.整体架构 整体架构目录:ASP.NET Core分布式项目实战-目录 一.目录 1.zipkin监控 2.skywalking监控 3.ELK日志监控 asp.net Core 交流群:78746 ...

  7. JQ_五星级评分特效

    代码如下:<!DOCTYPE HTML><html><head><meta http-equiv="Content-Type" conte ...

  8. Jq_网站顶部定时折叠广告

    <!DOCTYPE html><html xmlns="http://www.w3.org/1999/xhtml"><head><titl ...

  9. MatCap冰冻效果Shader

    MatCap方案 使用说明 制作合适的MatCap贴图 这张图决定冰像不像,网上找.Vray渲个球.ASE或者ShaderForge连,甚至直接手绘,总之只要一张长得像下面的图 注意MatCap图只有 ...

  10. Unity特殊路径

    Resources: Resources文件可以在根目录下,也可以在子目录下,只要叫Resources就好.Resources目录下所有资源将被打包进游戏存放资源的archive中,Resources ...