「6月雅礼集训 2017 Day2」A
【题目大意】
给出一棵树,求有多少对点(u,v)满足其路径上不存在两个点a,b满足(a,b)=1
n<=10^5
【题解】
考虑找出所有不符合的点对,共有n*ln(n)对,他们要么是祖先->儿子边,要么是不是。
考虑祖先->儿子边,那么一个点在祖先以上,一个点在儿子以下的点对全部无法访问。
考虑另外一种边,就是LCA不是两个端点的,这就比较好统计了,两个点在这两棵子树的点对无法访问。
考虑用DFS序,这样子树就是连续的一段(祖先以上是连续两段)
然后就是一个二维覆盖问题,用扫描线+线段树即可解决。
复杂度O(nln(n)logn)
注意。。扫描线数组要开到 4 * n * ln(n) 不然。。会奇怪的WA/RE。。。
# include <stdio.h>
# include <string.h>
# include <iostream>
# include <algorithm> using namespace std; typedef long long ll;
typedef unsigned long long ull;
typedef long double ld; # define RG register
# define ST static const int M = 2e5 + , N = 1e5 + , Max = * M;
const int mod = ; int n, head[N], nxt[M], to[M], tot = ;
inline void add(int u, int v) {
++tot; nxt[tot] = head[u]; head[u] = tot; to[tot] = v;
}
inline void adde(int u, int v) {
add(u, v), add(v, u);
} int in[N], out[N], DFN = ;
int dep[N], fa[N][];
inline void dfs(int x, int fat = ) {
in[x] = ++DFN; dep[x] = dep[fat] + ;
fa[x][] = fat;
for (int i=; i<=; ++i) fa[x][i] = fa[fa[x][i-]][i-];
for (int i=head[x]; i; i=nxt[i]) {
if(to[i] == fat) continue;
dfs(to[i], x);
}
out[x] = DFN;
} inline int lca(int u, int v) {
if(dep[u] < dep[v]) swap(u, v);
for (int i=; ~i; --i)
if((dep[u] - dep[v]) & (<<i)) u = fa[u][i];
if(u == v) return u;
for (int i=; ~i; --i)
if(fa[u][i] != fa[v][i]) u = fa[u][i], v = fa[v][i];
return fa[u][];
} inline int jump(int u, int anc) {
for (int i=; ~i; --i)
if(dep[fa[u][i]] > dep[anc]) u = fa[u][i];
return u;
} struct pa {
int x, yl, yr, d;
pa() {}
pa(int x, int yl, int yr, int d) : x(x), yl(yl), yr(yr), d(d) {}
friend bool operator < (pa a, pa b) {
return a.x < b.x;
}
}p[Max * ]; int pn = ; inline void ADD(int xl, int xr, int yl, int yr) {
p[++pn] = pa(xl, yl, yr, );
p[++pn] = pa(xr+, yl, yr, -);
} inline void doit(int x, int y) {
int par = lca(x, y);
// if(par == -1) cout << x << ' ' << y << endl;
if(dep[x] > dep[y]) swap(x, y);
if(x == par) {
int pars = jump(y, par);
ADD(, in[pars] - , in[y], out[y]);
ADD(in[y], out[y], out[pars] + , n);
} else {
if(in[x] > in[y]) swap(x, y);
ADD(in[x], out[x], in[y], out[y]);
}
} struct SMT {
int w[Max], tag[Max];
# define ls (x<<)
# define rs (x<<|)
inline void set() {
memset(w, , sizeof w);
memset(tag, , sizeof tag);
}
inline int gs(int x, int l, int r) {
if(tag[x]) return r-l+;
else return w[x];
}
inline void edt(int x, int l, int r, int L, int R, int d) {
if(L > R) return ;
if(L <= l && r <= R) {tag[x] += d; return ;}
int mid = l+r>>;
if(L <= mid) edt(ls, l, mid, L, R, d);
if(R > mid) edt(rs, mid+, r, L, R, d);
w[x] = gs(ls, l, mid) + gs(rs, mid+, r);
}
inline int sum(int x, int l, int r, int L, int R) {
if(L > R) return ;
if(tag[x]) return min(R, r) - max(L, l) + ;
if(L <= l && r <= R) return gs(x, l, r);
int mid = l+r>>, ret = ;
if(L <= mid) ret += sum(ls, l, mid, L, R);
if(R > mid) ret += sum(rs, mid+, r, L, R);
return ret;
}
# undef ls
# undef rs
}T; int main() {
// freopen("A.in", "r", stdin);
// freopen("A.out", "w", stdout);
cin >> n;
for (int i=, u, v; i<n; ++i) {
scanf("%d%d", &u, &v);
adde(u, v);
}
dfs(, );
for (int i=; i<=n; ++i)
for (int j=i+i; j<=n; j+=i) doit(i, j); sort(p+, p+pn+); T.set();
ll ans = (ll)n * (n-) / ;
for (int i=, j=; i<=n; ++i) {
while(j<=pn && p[j].x == i) T.edt(, , n, p[j].yl, p[j].yr, p[j].d), ++j;
ans -= T.sum(, , n, i+, n);
}
cout << ans;
return ;
}
「6月雅礼集训 2017 Day2」A的更多相关文章
- 「6月雅礼集训 2017 Day2」C
[题目大意] 有一棵n个点的完全二叉树,边权均为1,每个点有小鸟容量c[i] 依次来了m只小鸟,第i只小鸟初始位置在pos[i]上,问来了x只小鸟的时候,怎样安排小鸟的路线可以使得小鸟移动的边权和最小 ...
- 「6月雅礼集训 2017 Day2」B
[题目大意] 求n*n的棋盘,每行每列都有2个黑格子的方案数. n<=10^7 [题解] zzq的做法好神奇啊 行列建点,二分图 左边有i个点,右边有j个点的方案数 f[i,j] 左边有i个点, ...
- 「6月雅礼集训 2017 Day10」quote
[题目大意] 一个合法的引号序列是空串:如果引号序列合法,那么在两边加上同一个引号也合法:或是把两个合法的引号序列拼起来也是合法的. 求长度为$n$,字符集大小为$k$的合法引号序列的个数.多组数据. ...
- 「6月雅礼集训 2017 Day4」qyh(bzoj2687 交与并)
原题传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=2687 [题目大意] 给出若干区间,求一个区间的大于等于2的子集,使得 |区间并| 和 | ...
- 「6月雅礼集训 2017 Day11」delight
[题目大意] 有$n$天,每天能吃饭.睡觉.什么事也不干 每天吃饭的愉悦值为$e_i$,睡觉的愉悦值为$s_i$,什么都不干愉悦值为0. 要求每连续$k$天都要有至少$E$天吃饭,$S$天睡觉. 求最 ...
- 「6月雅礼集训 2017 Day11」jump
[题目大意] 有$n$个位置,每个位置有一个数$x_i$,代表从$i$经过1步可以到达的点在$[\max(1, i-x_i), \min(i+x_i, n)]$中. 定义$(i,j)$的距离表示从$i ...
- 「6月雅礼集训 2017 Day11」tree
[题目大意] 给出一棵带权树,有两类点,一类黑点,一类白点. 求切断黑点和白点间路径的最小代价. $n \leq 10^5$ [题解] 直接最小割能过..但是树形dp明显更好写 设$f_{x,0/1/ ...
- 「6月雅礼集训 2017 Day10」perm(CodeForces 698F)
[题目大意] 给出一个$n$个数的序列$\{a_n\}$,其中有些地方的数为0,要求你把这个序列填成一个1到$n$的排列,使得: $(a_i, a_j) = 1$,当且仅当$(i, j) = 1$.多 ...
- 「6月雅礼集训 2017 Day8」route
[题目大意] 给出平面上$n$个点,求一条连接$n$个点的不相交的路径,使得转换的方向符合所给长度为$n-2$的字符串. $n \leq 5000$ [题解] 考虑取凸包上一点,然后如果下一个是‘R' ...
随机推荐
- BeyondCompare3 提示许可证密钥已被撤销解决方法
今天对比文件提示 许可证密钥已被撤销:3281-0350! 找过了几个注册码还是不行. 正确简单的解决方法: 1.找到 BCState.xml 文件 2.编辑器打开,删除<TCheckFor ...
- localStorage简析
声明:引用自http://www.cnblogs.com/st-leslie/p/5617130.html 一.什么是localStorage.sessionStorage 在HTML5中,新加入了一 ...
- C++学习002-C++代码中插入汇编语句
在C++中我们有时会遇到使用汇编语言的情况,这时可以在前面加上关键字“_asm”宏. 如下示例 编写环境 :vs2015 int main() { __asm mov al, 0x20; __asm ...
- Python request 简单使用
Requests 是用Python语言编写,基于 urllib,采用 Apache2 Licensed 开源协议的 HTTP 库.它比 urllib 更加方便,可以节约我们大量的工作,完全满足 HTT ...
- 在Kotlin编写RecyclerView适配器(KAD 16)
作者:Antonio Leiva 时间:Mar 14, 2017 原文链接:https://antonioleiva.com/recyclerview-adapter-kotlin/ 通过创建Recy ...
- debian8.2安装kubernetes
master上通过kubeadm安装Kubernetes 添加国内阿里源后安装kubeadm: deb https://mirrors.aliyun.com/kubernetes/apt/ kuber ...
- mongodb数据库高级操作
1.创建索引 2.索引名称 3.其他索引 4.explain 5.操作索引 6.高级特性 7.固定集合 8.导入导出 9.上锁 10.添加用户 11.主从复制
- HDU 4468 Spy(KMP+贪心)(2012 Asia Chengdu Regional Contest)
Description “Be subtle! Be subtle! And use your spies for every kind of business. ”― Sun Tzu“A spy w ...
- Python中运算符"=="和"is"的差别分析
前言 在讲is和==这两种运算符区别之前,首先要知道Python中对象包含的三个基本要素,分别是:id(身份标识).python type()(数据类型)和value(值).is和==都是对对象进行比 ...
- Mininet实验 动态改变转发规则
介绍 拓扑如下: 在该环境下,假设H1 ping H4,初始的路由规则是S1-S2-S5,一秒后,路由转发规则变为S1-S3-S5,再过一秒,规则变为S1-S4-S5,然后再回到最初的转发规则S1-S ...