【UOJ388】配对树(dsu on tree+线段树)
题意:
给出一颗含有\(n\)个结点的无根树,之后给出一个长度为\(m\)的序列,每个元素在\([1,n]\)之间。
现在序列中每个长度为偶数的区间的完成时间定义为树上最小配对方法中每对匹配点间距离的总和。
现在要求所有长度为偶数的区间的完成时间的和。
思路:
- 首先不妨将这颗树转化为有根树,最终不会影响答案。
- 注意到性质:偶数个点的两两匹配方式是唯一的,都是最深的两个点相互匹配,这样才能保证没有重复计算的边。
- 在子树内部直接计算不好算,要考虑很多东西(一开始就想偏了QAQ)。因为匹配方式唯一,所以子树中若全部匹配完成,那么至多只会剩下一个点。显然,若剩下一个点,此时子树及其父亲这条边肯定会算上。
- 我们将问题转化为单独考虑一条边的贡献,一条边的贡献次数即为序列元素在该子树中出现奇数次的偶数区间个数。
- 考虑暴力计算:对于当前的子树,暴力给序列上打上标记,然后做一个前缀和,若存在\(i,j,i<j\),满足\(i\equiv j(mod\ 2)\)且\(s_j\equiv s_i+1(mod\ 2)\)即可。
- 如果用线段树维护,那么只需要维护奇数位置、偶数位置上面前缀模\(2\)意义下\(1\)的个数即可快速求得当前子树内部的答案。
- 时间复杂度\(O(n^2logn)\)。
- 因为这涉及到子树问题,且不带修改,所以可以直接施展\(dsu\ on\ tree\),最终复杂度为\(O(nlog^2n)\)。
感觉这种将子树内部的问题转化为子树+一条边的问题似乎是一种套路?(上次atcoder有个题也是这样)。
以后这种计算贡献的题可以尝试转换一下思路,单独考虑每个元素的贡献。
这个题还可以直接线段树合并来写,复杂度是\(O(nlogn)\),不过(不会写)懒得写了。。
代码如下:
/*
* Author: heyuhhh
* Created Time: 2019/11/16 8:54:36
*/
#include <bits/stdc++.h>
#define MP make_pair
#define fi first
#define se second
#define sz(x) (int)(x).size()
#define all(x) (x).begin(), (x).end()
#define INF 0x3f3f3f3f
//#define Local
#ifdef Local
#define dbg(args...) do { cout << #args << " -> "; err(args); } while (0)
void err() { std::cout << '\n'; }
template<typename T, typename...Args>
void err(T a, Args...args) { std::cout << a << ' '; err(args...); }
#else
#define dbg(...)
#endif
void pt() {std::cout << '\n'; }
template<typename T, typename...Args>
void pt(T a, Args...args) {std::cout << a << ' '; pt(args...); }
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
//head
const int N = 1e5 + 5, MOD = 998244353;
template <class T>
inline void read(T& x) {
static char c;
x = 0;
bool sign = 0;
while (!isdigit(c = getchar()))
if (c == '-')
sign = 1;
for (; isdigit(c); x = x * 10 + c - '0', c = getchar())
;
if (sign)
x = -x;
}
int n, m;
int a[N];
struct Edge{
int v, w, next;
}e[N << 1];
int head[N], tot;
void adde(int u, int v, int w) {
e[tot].v = v; e[tot].w = w; e[tot].next = head[u]; head[u] = tot++;
}
int pre[N], last[N];
int sz[N], bson[N], son;
void dfs(int u, int fa) {
int mx = 0; sz[u] = 1;
for(int i = head[u]; i != -1; i = e[i].next) {
int v = e[i].v;
if(v != fa) {
a[v] = e[i].w;
dfs(v, u);
sz[u] += sz[v];
if(sz[v] > mx) mx = sz[v], bson[u] = v;
}
}
}
int odd[N << 2], even[N << 2], rev[N << 2];
int ans;
void Reverse(int o, int l, int r) {
if(l == r) {
if(l & 1) odd[o] ^= 1;
else even[o] ^= 1;
} else {
odd[o] = (r + 1) / 2 - l / 2 - odd[o];
even[o] = r / 2 - (l - 1) / 2 - even[o];
}
rev[o] ^= 1;
}
void push_up(int o) {
odd[o] = odd[o << 1] + odd[o << 1|1];
even[o] = even[o << 1] + even[o << 1|1];
}
void push_down(int o, int l, int r) {
if(rev[o]) {
int mid = (l + r) >> 1;
Reverse(o << 1, l, mid);
Reverse(o << 1|1, mid + 1, r);
rev[o] = 0;
}
}
void upd(int o, int l, int r, int L, int R) {
if(L <= l && r <= R) {
Reverse(o, l, r);
return;
}
push_down(o, l, r);
int mid = (l + r) >> 1;
if(L <= mid) upd(o << 1, l, mid, L, R);
if(R > mid) upd(o << 1|1, mid + 1, r, L, R);
push_up(o);
}
int Get() {
int res = (1ll * (m / 2 + 1 - even[1]) * even[1] % MOD + 1ll * ((m + 1) / 2 - odd[1]) * odd[1] % MOD) % MOD;
return res;
}
void calc(int u, int fa) {
for(int i = last[u]; i; i = pre[i]) upd(1, 1, m, i, m);
for(int i = head[u]; i != -1; i = e[i].next) {
int v = e[i].v;
if(v == fa || v == son) continue;
calc(v, u);
}
}
void dfs2(int u, int fa, int op) {
for(int i = head[u]; i != -1; i = e[i].next) {
int v = e[i].v;
if(v != fa && v != bson[u]) dfs2(v, u, 0);
}
if(bson[u]) dfs2(bson[u], u, 1);
son = bson[u];
calc(u, fa);
son = 0;
dbg(u, get());
ans = (ans + 1ll * a[u] * Get() % MOD) % MOD;
if(!op) calc(u, fa);
}
void run(){
memset(head, -1, sizeof(head));
read(n), read(m);
for(int i = 1; i < n; i++) {
int u, v, w;
read(u), read(v), read(w);
adde(u, v, w); adde(v, u, w);
}
for(int i = 1; i <= m; i++) {
int x; read(x);
pre[i] = last[x];
last[x] = i;
}
dfs(1, 0);
dfs2(1, 0, 1);
cout << ans << '\n';
}
int main() {
run();
return 0;
}
【UOJ388】配对树(dsu on tree+线段树)的更多相关文章
- BZOJ.3307.雨天的尾巴(dsu on tree/线段树合并)
BZOJ 洛谷 \(dsu\ on\ tree\).(线段树合并的做法也挺显然不写了) 如果没写过\(dsu\)可以看这里. 对修改操作做一下差分放到对应点上,就成了求每个点子树内出现次数最多的颜色, ...
- 【cf375】D. Tree and Queries(dsu on tree+线段树)
传送门 题意: 给出一颗以\(1\)为根的有根树,每个结点有个颜色\(c_i\). 之后要回答\(m\)组询问,每组询问包含\(v_i,k_i\),要回答以\(v_i\)为根的子树中,颜色出现次数不小 ...
- poj 3237 Tree(树链剖分,线段树)
Tree Time Limit: 5000MS Memory Limit: 131072K Total Submissions: 7268 Accepted: 1969 Description ...
- HDU 3333 Turing Tree 线段树+离线处理
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=3333 Turing Tree Time Limit: 6000/3000 MS (Java/Othe ...
- 树(一)——线段树
问题 现在有1~30这30个数,数N被抽上的概率正比于1/sqrt(N+1),求满足这个概率分布的随机数发生器. 思路 第一,如何解决这个"概率正比"问题. 第二,如何产生满足条件 ...
- bzoj 4034 [HAOI2015] T2(树链剖分,线段树)
4034: [HAOI2015]T2 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 1536 Solved: 508[Submit][Status] ...
- bzoj 1036 [ZJOI2008]树的统计Count(树链剖分,线段树)
1036: [ZJOI2008]树的统计Count Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 10677 Solved: 4313[Submit ...
- bzoj 3626 [LNOI2014]LCA(离线处理+树链剖分,线段树)
3626: [LNOI2014]LCA Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1272 Solved: 451[Submit][Status ...
- HDU 5877 2016大连网络赛 Weak Pair(树状数组,线段树,动态开点,启发式合并,可持久化线段树)
Weak Pair Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others) Tota ...
随机推荐
- selenium时间等待方法
在UI自动化测试中,必然会遇到环境不稳定.网络慢等情况.当你觉得定位没有问题,但程序却直接报了元素不可见时,那你就需要思考是否因为程序运行太快或者页面加载太慢而造成了元素不可见,必须要再等待直至元素可 ...
- python 错误信息是:sudo :apt-get:command not found
1.问题描述 错误信息是:sudo :apt-get:command not found 2.问题原因及解决 在centos下用yum install xxx yum和apt-get的区别一般来说著名 ...
- 运行springboot项目报错:Field userMapper in XX required a bean of type 'xx' that could not be found.
运行springboot项目报错: *************************** APPLICATION FAILED TO START ************************** ...
- Fiddler安装
1.在网上搜一下fiddler的安装包,下载.下载完成能看到一个exe文件. 2.点击文件安装,同意,选择一个目录进行安装即可. 3.安装完成,打开应用是一个这样的界面. 需要安装包的下面评论.... ...
- LG1344 「USACO4.4」Pollutant Control 最小割
问题描述 LG1344 题解 我太菜了,我一开始竟然没有看出这是个最小割裸题... 两个询问. 第一个询问,直接跑最小割就好了. 第二个询问,建图的时候边权建 \(1\) ,代表割掉这条边需要 \(1 ...
- Nginx与keepalived实现高可用
主keepalived设置 #安装keepalived [root@localhost ~]# yum -y install keepalived #安装nginx [root@localhost ~ ...
- IDEA创建maven项目慢的不行
方法二 下载archetype-catalog.xml文件,在maven的VM Options加上-DarchetypeCatalog=local 默认情况下,创建maven项目是从网络下载catal ...
- Java类加载机制以及双亲委派模型
一.Java类加载机制 1.概述 Class文件由类装载器装载后,在JVM中将形成一份描述Class结构的元信息对象,通过该元信息对象可以获知Class的结构信息:如构造函数,属性和方法等,Java允 ...
- NopCommerce 4.2的安装与运行
一.关于NopCommerce NopCommerce是国外ASP.Net领域一个高质量的B2C开源电商项目,最新版本4.2基于ASP.NET Core MVC 2.2和EF Core 2.2开发,其 ...
- Jenkins的安装后配置
配置Jenkins server 1.安装插件:选择自定义插件安装 Locale plugin Localization: Chinese (Simplified) Git Subversion HT ...