FZOJ 2176 easy problem ( 树链剖分 )

2024-05-29 18:38

本文主要是介绍FZOJ 2176 easy problem ( 树链剖分 ),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题目链接~~>

做题感悟:感觉做多了树链剖分的题目,有许多是树链剖分 + 想法。。

解题思路:

                这题很明显的一点就是 k 很小,那就是告诉你可以从 k 入手,怎样入手呢 ? 观察可以发现无非最多是 k 类点 ,0 ~ k-1 ,分别表示与根的距离模 k .这样就可以把点分类加权值,但是每个线段树里存的还是所有元素,查询的时候对应查找。

代码:

#include<iostream>
#include<sstream>
#include<map>
#include<cmath>
#include<fstream>
#include<queue>
#include<vector>
#include<sstream>
#include<cstring>
#include<cstdio>
#include<stack>
#include<bitset>
#include<ctime>
#include<string>
#include<cctype>
#include<iomanip>
#include<algorithm>
using namespace std  ;
#define INT long long int
#define L(x)  (x * 2)
#define R(x)  (x * 2 + 1)
const int INF = 0x3f3f3f3f ;
const double esp = 0.00000000001 ;
const double PI = acos(-1.0) ;
const INT mod = 1234567891 ;
const int MY = (1<<5) + 5 ;
const int MX = 50000 + 5 ;
const int S = 20 ;
int n ,m ,k ,num ,idx ;
int head[MX] ,top[MX] ,siz[MX] ,son[MX] ,dep[MX] ,ti[MX] ,to[MX] ,father[MX] ;
struct Edge
{int v ,next ;
}E[MX*2] ;
void addedge(int u ,int v)
{E[num].v = v ; E[num].next = head[u] ;head[u] = num++ ;E[num].v = u ; E[num].next = head[v] ;head[v] = num++ ;
}
void dfs_find(int u ,int fa)
{dep[u] = dep[fa] + 1 ;siz[u] = 1 ;son[u] = 0 ;father[u] = fa ;for(int i = head[u] ;i != -1 ;i = E[i].next){int v = E[i].v ;if(v == fa)  continue ;dfs_find(v ,u) ;siz[u] += siz[v] ;if(siz[son[u]] < siz[v])  son[u] = v ;}
}
void dfs_time(int u ,int fa)
{top[u] = fa ;ti[u] = idx++ ;if(son[u])   dfs_time(son[u] ,top[u]) ;for(int i = head[u] ;i != -1 ;i = E[i].next){int v = E[i].v ;if(v == son[u] || v == father[u])  continue ;dfs_time(v ,v) ;}to[u] = idx-1 ;
}
struct node
{int le ,rt ,add ,c ;
}T[6][MX*4] ;
void push_down(int cnt ,int i)
{int& temp = T[cnt][i].add ;if(temp){T[cnt][L(i)].add += temp ;T[cnt][L(i)].c += temp ;T[cnt][R(i)].add += temp ;T[cnt][R(i)].c += temp ;temp = 0 ;}
}
void build(int cnt ,int i ,int le ,int rt)
{T[cnt][i].le = le ; T[cnt][i].rt = rt ;T[cnt][i].add = T[cnt][i].c = 0 ;if(le == rt)  return ;int Mid = (le + rt)>>1 ;build(cnt ,L(i) ,le ,Mid) ;build(cnt ,R(i) ,Mid+1 ,rt) ;
}
void section(int cnt ,int i ,int le ,int rt ,int val)
{if(T[cnt][i].le == le && T[cnt][i].rt == rt){T[cnt][i].add += val ;T[cnt][i].c += val ;return ;}push_down(cnt ,i) ;int Mid = (T[cnt][i].le + T[cnt][i].rt)>>1 ;if(le > Mid)   section(cnt ,R(i) ,le ,rt ,val) ;else if(rt <= Mid)   section(cnt ,L(i) ,le ,rt ,val) ;else{section(cnt ,L(i) ,le ,Mid ,val) ;section(cnt ,R(i) ,Mid+1 ,rt ,val) ;}
}
int Query(int cnt ,int i ,int pos)
{if(T[cnt][i].le == T[cnt][i].rt)return  T[cnt][i].c ;push_down(cnt ,i) ;int Mid = (T[cnt][i].le + T[cnt][i].rt)>>1 ;if(pos <= Mid)    return Query(cnt ,L(i) ,pos) ; // 在左边else    return Query(cnt ,R(i) ,pos) ;
}
int main()
{//freopen("input.txt" ,"r" ,stdin) ;int Tx ,u ,v ,cnt ,x ,dt ,pos ,val ,cse = 1 ;scanf("%d" ,&Tx) ;while(Tx--){scanf("%d%d%d" ,&n ,&m ,&k) ;num = 0 ;memset(head ,-1 ,sizeof(head)) ;for(int i = 1 ;i < n ; ++i){scanf("%d%d" ,&u ,&v) ;addedge(u ,v) ;}dep[1] = siz[0] = 0 ;dfs_find(1 ,1) ;idx = 1 ;dfs_time(1 ,1) ;for(int i = 0 ;i < k ; ++i)build(i ,1 ,1 ,n) ;printf("Case#%d:\n" ,cse++) ;for(int i = 0 ;i < m ; ++i){scanf("%d%d" ,&pos ,&x) ;if(pos == 1){scanf("%d" ,&val) ;for(int j = 0 ;j < k ; ++j){dt = (dep[x] - 1)%k ;  // 根所在的集合cnt = (dt + j)%k ;section(cnt ,1 ,ti[x] ,to[x] ,(j+1)*val) ;}}elseprintf("%d\n" ,Query((dep[x]-1)%k ,1 ,ti[x])) ;}}return 0 ;
}

  

这篇关于FZOJ 2176 easy problem ( 树链剖分 )的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1014330

相关文章

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

uva 10025 The ? 1 ? 2 ? ... ? n = k problem(数学)

题意是    ?  1  ?  2  ?  ...  ?  n = k 式子中给k,? 处可以填 + 也可以填 - ,问最小满足条件的n。 e.g k = 12  - 1 + 2 + 3 + 4 + 5 + 6 - 7 = 12 with n = 7。 先给证明,令 S(n) = 1 + 2 + 3 + 4 + 5 + .... + n 暴搜n,搜出当 S(n) >=

LibSVM学习(六)——easy.py和grid.py的使用

我们在“LibSVM学习(一)”中,讲到libSVM有一个tools文件夹,里面包含有四个python文件,是用来对参数优选的。其中,常用到的是easy.py和grid.py两个文件。其实,网上也有相应的说明,但很不系统,下面结合本人的经验,对使用方法做个说明。        这两个文件都要用python(可以在http://www.python.org上下载到,需要安装)和绘图工具gnup

FZOJ 2166 inversion

算是一道模拟题吧,,,此题就逆序数不用归并排序之类的,而是用双层for循环模拟求解,,因为题意中要求任意两个数交换位置,所以只需在双层for循环中模拟两个数交换即可,在模拟的过程中,不需要准确的求出模拟之后的逆序数,只需要考虑逆序数的变化量为多少,最后,求出变化量最小的,用最初的逆序数求得最后结果。。。 其中模拟两个数交换之后逆序数的变化量解法: 如下: 有   a,b,c,d,e,f,g

11991 - Easy Problem from Rujia Liu?

题意: 输入一串整型数列,再输入两个数k,v,输出第k个v的序号。不存在则输出0,如第一个样例 8 41 3 2 2 4 3 2 11 3 //第1个3,序号为2,输出22 4 //第2个4,不存在,输出03 2 //第3个2,序号为7,输出74 2 思路: struct num {

【开发工具】开发过程中,怎么通过Easy JavaDoc快速生成注释。

文章目录 引言什么是Easy JavaDoc?Easy JavaDoc用来干什么?如何使用Easy JavaDoc?安装Easy JavaDoc配置Easy JavaDoc使用Easy JavaDoc生成注释 Easy JavaDoc与IDEA自带注释的区别IDEA自带注释Easy JavaDoc Easy JavaDoc的优缺点优点缺点 步骤 1:打开设置步骤 2:找到Easy JavaD

HDU 1016 Prime Ring Problem (深搜)

OJ题目 : click here ~~ 大概题意:n个数,形成一个环,使得相邻两个数的和为素数。以1开始,按字典序输出序列。 很简单的深搜。 AC_CODE int n;int visit[22];int num[22];int len;bool Is_prime(int x){for(int i = 2;i*i <= x;i++)if(x%i == 0) return

LVM 'Can’t open /dev/sdb1 exclusively. Mounted filesystem?' Problem

在将几块盘做LVM时,遇到一个之前都没遇到过的问题: root@ubuntu:~# pvcreate /dev/sdc1Can't open /dev/sdc1 exclusively. Mounted filesystem? 首先第一反应就是查看这个分区是否已经在使用了,但是没有。 查看硬盘的一些信息: root@ubuntu:~# cat /proc/partitionsmajo

easy简化封装

//confirm function Confirm(msg, control) {$.messager.confirm('确认', msg, function (r) {if (r) {eval(control.toString().slice(11));}});return false;}//loadfunction Load() {$("<div class=\"datagrid-ma

树链剖分 + 后缀数组 - E. Misha and LCP on Tree

E. Misha and LCP on Tree  Problem's Link   Mean:  给出一棵树,每个结点上有一个字母。每个询问给出两个路径,问这两个路径的串的最长公共前缀。 analyse: 做法:树链剖分+后缀数组. 记录每条链的串,正反都需要标记,组成一个长串. 然后记录每条链对应的串在大串中的位置,对大串求后缀数组,最后询问就是在一些链