【bzoj3772】精神污染 STL+LCA+主席树
题目描述
输入
输出
样例输入
5 3
1 2
2 3
3 4
2 5
3 5
2 5
1 4
样例输出
1/3
题解
Orz PoPoQQQ
由题可知所求为能够覆盖的路径对数/总路径对数。
总对数=m*(m-1)/2,所以只要求能够覆盖的对数即可。
首先我们易知如果某条路径B的两个端点都在另一条路径A上,那么A必然会覆盖B,反之亦然。
那么问题就转化为求A中所有点在多少条B上。
我们设Bi为(xi,yi),则问题就是求对于所有xi,有多少个yi在A上。
可以先把所有xi对应的所有yi存起来,这里使用vector。
用主席树维护每个xi对应的yi,使用Dfs序,入栈位置+1,出栈位置-1。
那么直接查询[xA,yA]范围内数的个数即可。
当然,每条A不能直接使用Dfs序计算,需要拆成[xA,lca]和[yA,lca](用倍增lca解决,树剖也可以),再减去重复计算的[lca,lca]。
由于主席树是建立在每个节点的父亲上的,所以每棵树代表着[1,x],所求路径即为[1,xA]+[1,yA]-[1,lca]-[1,fa[lca]]。
当然自己包含自己是不行的,需要减去。
这样就可以求出能够覆盖的路径对数。
需要说明的是题目中已经声明没有重复路径,所以不需要去重。
最后求gcd再输出即可。
数组大小略坑。
#include <cstdio>
#include <vector>
#include <algorithm>
#define N 100010
using namespace std;
struct data
{
int x , y , lca;
}a[N];
vector<int> v[N];
int fa[N][20] , deep[N] , log[N] , head[N] , to[N * 2] , next[N * 2] , cnt;
int lp[N] , rp[N] , pl , q[N] , top;
int root[N * 2] , ls[N * 39] , rs[N * 39] , si[N * 39] , tot;
bool cmp(data a , data b)
{
return a.x == b.x ? a.y < b.y : a.x < b.x;
}
void add(int x , int y)
{
to[++cnt] = y;
next[cnt] = head[x];
head[x] = cnt;
}
void dfs(int x)
{
int i;
lp[x] = ++pl;
q[++top] = x;
for(i = head[x] ; i ; i = next[i])
{
if(to[i] != fa[x][0])
{
fa[to[i]][0] = x;
deep[to[i]] = deep[x] + 1;
dfs(to[i]);
}
}
rp[x] = ++pl;
}
int getlca(int x , int y)
{
int i;
if(deep[x] < deep[y]) swap(x , y);
for(i = log[deep[x] - deep[y]] ; i >= 0 ; i -- )
if(deep[x] - (1 << i) >= deep[y])
x = fa[x][i];
for(i = log[deep[x]] ; i >= 0 ; i -- )
if(fa[x][i] != fa[y][i])
x = fa[x][i] , y = fa[y][i];
return x == y ? x : fa[x][0];
}
void update(int x , int &y , int l , int r , int p , int a)
{
y = ++tot;
if(l == r)
{
si[y] = si[x] + a;
return;
}
int mid = (l + r) >> 1;
if(p <= mid) rs[y] = rs[x] , update(ls[x] , ls[y] , l , mid , p , a);
else ls[y] = ls[x] , update(rs[x] , rs[y] , mid + 1 , r , p , a);
si[y] = si[ls[y]] + si[rs[y]];
}
int query(int w , int x , int y , int z , int l , int r , int b , int e)
{
if(b <= l && r <= e)
return si[w] + si[x] - si[y] - si[z];
int mid = (l + r) >> 1 , sum = 0;
if(b <= mid) sum += query(ls[w] , ls[x] , ls[y] , ls[z] , l , mid , b , e);
if(e > mid) sum += query(rs[w] , rs[x] , rs[y] , rs[z] , mid + 1 , r , b , e);
return sum;
}
long long gcd(long long a , long long b)
{
return b ? gcd(b , a % b) : a;
}
int main()
{
int n , m , i , j , x , y;
long long A = 0 , B = 0 , T;
scanf("%d%d" , &n , &m);
for(i = 1 ; i < n ; i ++ )
scanf("%d%d" , &x , &y) , add(x , y) , add(y , x);
dfs(1);
for(i = 2 ; i <= n ; i ++ ) log[i] = log[i >> 1] + 1;
for(i = 1 ; i <= log[n] ; i ++ )
for(j = 1 ; j <= n ; j ++ )
if(deep[j] >= (1 << i))
fa[j][i] = fa[fa[j][i - 1]][i - 1];
for(i = 1 ; i <= m ; i ++ )
{
scanf("%d%d" , &a[i].x , &a[i].y);
a[i].lca = getlca(a[i].x , a[i].y);
v[a[i].x].push_back(a[i].y);
}
sort(a + 1 , a + m + 1 , cmp);
for(i = 1 ; i <= top ; i ++ )
{
x = q[i];
root[x] = root[fa[x][0]];
for(j = 0 ; j < (int)v[x].size() ; j ++ )
{
update(root[x] , root[x] , 1 , pl , lp[v[x][j]] , 1);
update(root[x] , root[x] , 1 , pl , rp[v[x][j]] , -1);
}
}
for(i = 1 ; i <= m ; i ++ )
{
A += query(root[a[i].x] , root[a[i].y] , root[a[i].lca] , root[fa[a[i].lca][0]] , 1 , pl , lp[a[i].lca] , lp[a[i].x]);
A += query(root[a[i].x] , root[a[i].y] , root[a[i].lca] , root[fa[a[i].lca][0]] , 1 , pl , lp[a[i].lca] , lp[a[i].y]);
A -= query(root[a[i].x] , root[a[i].y] , root[a[i].lca] , root[fa[a[i].lca][0]] , 1 , pl , lp[a[i].lca] , lp[a[i].lca]);
A -- ;
}
B = (long long)m * (m - 1) / 2;
T = gcd(A , B);
printf("%lld/%lld\n" , A / T , B / T);
return 0;
}
【bzoj3772】精神污染 STL+LCA+主席树的更多相关文章
- bzoj3772 精神污染 dfs 序+主席树
题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=3772 题解 很简单的一道题目. 上午研究一个题目的时候发现了这个题目是一个弱化版,所以来写了一 ...
- 【BZOJ3772】精神污染 DFS序+主席树
[BZOJ3772]精神污染 Description 兵库县位于日本列岛的中央位置,北临日本海,南面濑户内海直通太平洋,中央部位是森林和山地,与拥有关西机场的大阪府比邻而居,是关西地区面积最大的县,是 ...
- BZOJ3772精神污染——可持久化线段树+出栈入栈序
题目描述 兵库县位于日本列岛的中央位置,北临日本海,南面濑户内海直通太平洋,中央部位是森林和山地,与拥有关西机场的大阪府比邻而居,是关西地区面积最大的县,是集经济和文化于一体的一大地区,是日本西部门户 ...
- BZOJ 3772: 精神污染(dfs序+主席树)
传送门 解题思路 比较神仙的一道题.首先计算答案时可以每条路径所包含的路径数,对于\(x,y\)这条路径,可以在\(x\)这处开个\(vector\)存\(y\),然后计算时只需要算这个路径上每个点的 ...
- LCA+主席树 (求树上路径点权第k大)
SPOJ 10628. Count on a tree (树上第k大,LCA+主席树) 10628. Count on a tree Problem code: COT You are given ...
- 【bzoj3123】[Sdoi2013]森林 倍增LCA+主席树+启发式合并
题目描述 输入 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数.第三行包含N个非负 ...
- [BZOJ3772]精神污染 主席树上树+欧拉序
3772: 精神污染 Time Limit: 10 Sec Memory Limit: 64 MB Description 兵库县位于日本列岛的中央位置,北临日本海,南面濑户内海直通太平洋,中央部位 ...
- BZOJ3772 精神污染 【主席树 + dfs序】
题目 兵库县位于日本列岛的中央位置,北临日本海,南面濑户内海直通太平洋,中央部位是森林和山地,与拥有关西机场的大阪府比邻而居,是关西地区面积最大的县,是集经济和文化于一体的一大地区,是日本西部门户,海 ...
- BZOJ3772 精神污染 主席树 dfs序
欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ3772 题意概括 给出一个树,共n个节点. 有m条互不相同的树上路径. 现在让你随机选择2条路径,问 ...
随机推荐
- 20155335 俞昆 2016-2017-2 《Java程序设计》第九周学习总结
学号 2016-2017-2 <Java程序设计>第九周学习总结 ##JDBC入门 在正式介绍JDBC前,已知JDBC是用来执行SQL的解决方案,开发人员使用JDBC的标准接口,开发人员不 ...
- 人脸检测库libfacedetection介绍
libfacedetection是于仕琪老师放到GitHub上的二进制库,没有源码,它的License是MIT,可以商用.目前只提供了windows 32和64位的release动态库,主页为http ...
- php常用几个数组的区别
本文主要介绍的php数组函数主要有:sort.rsort.asort.arsort.ksort.krsort 测试数据定义一个关联数组如下: $data=[ 'f'=>123, 'b'=> ...
- 敏捷开发学习笔记-Agile development(AM)
以人为核心,迭代,循序渐进 项目被切分为多个子项目,每个子项目都经过测试,具备集成和可运行的特征 5个价值观:沟通.简单.反馈.勇气.谦逊 敏捷模型与瀑布模型的区别 相对于瀑布模型,提高开发效率和 ...
- 关于java中的“error: bad operand types for binary operator ”
今天做这个题目的时候(142. O(1) Check Power of 2),遇到一个错误“ bad operand types for binary operator '&' ”. 先贴一下 ...
- chrome编辑器与截图
在地址栏中输入 data:text/html,<html contenteditable>即可使用编辑功能,打开控制台,ctrl + shift + p 调用命令面板,输入 capture ...
- 【20180807模拟测试】t1 function
low逼的我也只能写这样的水题... 题面 对于一个整数,定义 f(x)为他的每个数位的阶乘的乘积.例如 f(135)=1! * 3! * 5! =720.给出一个数 a(可以包含前缀零),a 满足他 ...
- python3 SQLAlchemy模块使用
更详细的操作介绍:https://www.imooc.com/article/22343 定义: SQLAlchemy是Python编程语言下的一款ORM框架,该框架建立在数据库API之上,使用关系对 ...
- leetcode7_C++整数反转
给出一个 32 位的有符号整数,你需要将这个整数中每位上的数字进行反转. 示例 1: 输入: 输出: 示例 2: 输入: - 输出: - 示例 3: 输入: 输出: 注意: 假设我们的环境只能存 ...
- def语句和参数
如果调用print()或len()函数,你会传入一些值,放在括号内,在这里成为“参数”.也可以自己定义接受参数的函数.在文件编辑器中输入这个例子: def hello(name): print('He ...