P3038 [USACO11DEC]牧草种植Grass Planting
题目描述
Farmer John has N barren pastures (2 <= N <= 100,000) connected by N-1 bidirectional roads, such that there is exactly one path between any two pastures. Bessie, a cow who loves her grazing time, often complains about how there is no grass on the roads between pastures. Farmer John loves Bessie very much, and today he is finally going to plant grass on the roads. He will do so using a procedure consisting of M steps (1 <= M <= 100,000).
At each step one of two things will happen:
FJ will choose two pastures, and plant a patch of grass along each road in between the two pastures, or,
- Bessie will ask about how many patches of grass on a particular road, and Farmer John must answer her question.
Farmer John is a very poor counter -- help him answer Bessie's questions!
给出一棵n个节点的树,有m个操作,操作为将一条路径上的边权加一或询问某条边的权值。
输入输出格式
输入格式:
* Line 1: Two space-separated integers N and M
* Lines 2..N: Two space-separated integers describing the endpoints of a road.
* Lines N+1..N+M: Line i+1 describes step i. The first character of the line is either P or Q, which describes whether or not FJ is planting grass or simply querying. This is followed by two space-separated integers A_i and B_i (1 <= A_i, B_i <= N) which describe FJ's action or query.
输出格式:
* Lines 1..???: Each line has the answer to a query, appearing in the same order as the queries appear in the input.
输入输出样例
输入样例#1: 复制
4 6
1 4
2 4
3 4
P 2 3
P 1 3
Q 3 4
P 1 4
Q 2 4
Q 1 4
输出样例#1: 复制
2
1
2
//树剖是在点上操作的,这道题是边
//那么怎么把边权转成点权呢?
//根据树的性质可以知道,一个点可以有多个儿子,但是只会有一个爸爸,
//所以我们可以把这个点和它爸爸之间的那条边的边权转移到这个点上来
//用这个点的点权来表示这条边的权值
//因为根节点没有爸爸,所以它不表示任何边权,点权为0
//但是我们怎么样才能不把两个点的公共祖先的权值算进去啊?
//node[fx].s+1? 不行,这是它的重儿子的位置
// 考虑一下,我们在Query或者Modify的时候,都是当x和y同时处于一条链了之后就break
//然后再把这条链加上,最近公共祖先不就是这条链的top嘛!
//所以,我们在while循环外边写node[x].s+1就可以不算上公共祖先了。
//但是也要注意,如果最后是条轻边,我们就要if特判一下,不能让他进线段树查询了
//因为如果是轻边的话,最后的那条链退化成了最近公共祖先这一个点,不能要! #include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std; const int N=1e5+; int n,m;
int w[N];
int head[N],num_edge;
struct Edge
{
int v,nxt;
}edge[N<<];
struct Node
{
int fa,son;
int dep,top;
int size;
int s,t;
}node[N];
struct TREE
{
TREE *lson,*rson;
int l,r,mid,len;
int num,lazy;
}tree[N<<]; typedef TREE* Tree;
Tree Root,now_node=tree; inline int read()
{
char c=getchar();int num=;
for(;!isdigit(c);c=getchar())
if(c=='P') return ;
else if(c=='Q') return ;
for(;isdigit(c);c=getchar())
num=num*+c-'';
return num;
} inline void add_edge(int u,int v)
{
edge[++num_edge].v=v;
edge[num_edge].nxt=head[u];
head[u]=num_edge;
} void dfs1(int u)
{
node[u].size=;
for(int i=head[u],v;i;i=edge[i].nxt)
{
v=edge[i].v;
if(v==node[u].fa)
continue;
node[v].fa=u;
node[v].dep=node[u].dep+;
dfs1(v);
node[u].size+=node[v].size;
if(node[v].size>node[node[u].son].size)
node[u].son=v;
}
} int bound;
void dfs2(int u,int top)
{
node[u].top=top;
node[u].s=++bound;
if(node[u].son)
{
dfs2(node[u].son,top);
for(int i=head[u],v;i;i=edge[i].nxt)
{
v=edge[i].v;
if(v==node[u].son||v==node[u].fa)
continue;
dfs2(v,v);
}
}
node[u].t=bound;
} void build(Tree &root,int l,int r)
{
root=++now_node;
root->l=l,root->r=r,root->mid=l+r>>,root->len=r-l+;
if(l==r)
return;
build(root->lson,l,root->mid);
build(root->rson,root->mid+,r);
} inline void pushdown(Tree root)
{
if(root->lazy)
{
root->lson->lazy+=root->lazy;
root->rson->lazy+=root->lazy;
root->lson->num+=root->lson->len*root->lazy;
root->rson->num+=root->rson->len*root->lazy;
root->lazy=;
}
} void update(Tree root,int l,int r)
{
if(root->l==l&&r==root->r)
{
root->num+=root->len;
root->lazy+=;
return;
}
pushdown(root);
if(r<=root->mid)
update(root->lson,l,r);
else if(l>root->mid)
update(root->rson,l,r);
else
{
update(root->lson,l,root->mid);
update(root->rson,root->mid+,r);
}
root->num=root->lson->num+root->rson->num;
} int query(Tree root,int l,int r)
{
if(root->l==l&&root->r==r)
return root->num;
pushdown(root);
if(r<=root->mid)
return query(root->lson,l,r);
else if(l>root->mid)
return query(root->rson,l,r);
else
return query(root->lson,l,root->mid)+query(root->rson,root->mid+,r);
} inline void Modify(int x,int y)
{
int fx=node[x].top,fy=node[y].top;
while(fx!=fy)
{
if(node[fx].dep>node[fy].dep)
{
update(Root,node[fx].s,node[x].s);
x=node[fx].fa;
fx=node[x].top;
}
else
{
update(Root,node[fy].s,node[y].s);
y=node[fy].fa;
fy=node[y].top;
}
}
if(x!=y)
{
if(node[x].dep>node[y].dep)
update(Root,node[y].s+,node[x].s);
else
update(Root,node[x].s+,node[y].s);
}
} inline int Query(int x,int y)
{
int fx=node[x].top,fy=node[y].top;
int ans=;
while(fx!=fy)
{
if(node[fx].dep>node[fy].dep)
{
ans+=query(Root,node[fx].s,node[x].s);
x=node[fx].fa;
fx=node[x].top;
}
else
{
ans+=query(Root,node[fy].s,node[y].s);
y=node[fy].fa;
fy=node[y].top;
}
}
if(x!=y)
{
if(node[x].dep>node[y].dep)
return ans+query(Root,node[y].s+,node[x].s);
else
return ans+query(Root,node[x].s+,node[y].s);
}
return ans;
} int opt,u,v;
int main()
{
n=read(),m=read();
for(int i=;i<n;++i)
{
u=read(),v=read();
add_edge(u,v);
add_edge(v,u);
}
dfs1();
dfs2(,);
build(Root,,n);
for(int i=;i<=m;++i)
{
opt=read(),u=read(),v=read();
if(opt==)
Modify(u,v);
else
printf("%d\n",Query(u,v));
}
return ;
}
P3038 [USACO11DEC]牧草种植Grass Planting的更多相关文章
- 洛谷P3038 [USACO11DEC]牧草种植Grass Planting
题目描述 Farmer John has N barren pastures (2 <= N <= 100,000) connected by N-1 bidirectional road ...
- 洛谷 P3038 [USACO11DEC]牧草种植Grass Planting
题目描述 Farmer John has N barren pastures (2 <= N <= 100,000) connected by N-1 bidirectional road ...
- 洛谷 P3038 [USACO11DEC]牧草种植Grass Planting(树链剖分)
题解:仍然是无脑树剖,要注意一下边权,然而这种没有初始边权的题目其实和点权也没什么区别了 代码如下: #include<cstdio> #include<vector> #in ...
- AC日记——[USACO11DEC]牧草种植Grass Planting 洛谷 P3038
题目描述 Farmer John has N barren pastures (2 <= N <= 100,000) connected by N-1 bidirectional road ...
- 树链剖分【p3038】[USACO11DEC]牧草种植Grass Planting
表示看不太清. 概括题意 树上维护区间修改与区间和查询. 很明显树剖裸题,切掉,细节处错误T了好久 TAT 代码 #include<cstdio> #include<cstdlib& ...
- [USACO11DEC]牧草种植Grass Planting
图很丑.明显的树链剖分,需要的操作只有区间修改和区间查询.不过这里是边权,我们怎么把它转成点权呢?对于E(u,v),我们选其深度大的节点,把边权扔给它.因为这是树,所以每个点只有一个父亲,所以每个边权 ...
- 【LuoguP3038/[USACO11DEC]牧草种植Grass Planting】树链剖分+树状数组【树状数组的区间修改与区间查询】
模拟题,可以用树链剖分+线段树维护. 但是学了一个厉害的..树状数组的区间修改与区间查询.. 分割线里面的是转载的: ----------------------------------------- ...
- 洛谷P3038 牧草种植Grass Planting
思路: 首先,这道题的翻译是有问题的(起码现在是),查询的时候应该是查询某一条路径的权值,而不是某条边(坑死我了). 与平常树链剖分题目不同的是,这道题目维护的是边权,而不是点权,那怎么办呢?好像有点 ...
- USACO Grass Planting
洛谷 P3038 [USACO11DEC]牧草种植Grass Planting 洛谷传送门 JDOJ 2282: USACO 2011 Dec Gold 3.Grass Planting JDOJ传送 ...
随机推荐
- vue 下拉刷新数据的插件的使用:
1.安装: npm i vue-scroller -S npm install vue-scroller -D 2.在需要加载的页面中引入,或在公共js文件中引入: import VueScrolle ...
- linux 磁盘占用的排查流程
Linux 服务器在使用过程中可能会遇到各种问题,其中之一就是"没有可用空间". 遇到这种情况,就需要进行排查,定位到消耗了磁盘的那个文件夹. 流程如下: 1. df -h df ...
- 使用Jenkins编译打包SpringCloud微服务中的个别目录
意义说明: 使用Jenkins从Gogs拉取SpringCloud微服务,拉取的是整个仓库的内容,分好多个模块文件夹,但是使用maven编译打包的话只编译打包指定的模块文件夹 Gogs Webhook ...
- robot framework 如何处理循环条件下面的变量自增
下面举了一个基础栗子,可以运行的.${num}就是我需要的自增变量.有人也许会问为什么不用${i},不是我不想用,而是我${i}有其他用处,必须另外定义一个变量,需要注意的是定义变量的时候,应该在循环 ...
- 前端编译原理 笔记 -- BISON
前面总结的差不多了,这边记录下,零碎的相关阅读可以备忘的一些知识点 Bsion文档,下面是中文的地址 https://blog.csdn.net/chinamming/article/details ...
- shopxo代码审计
由于工作原因,分析了很多的cms也都写过文章,不过觉得好像没什么骚操作都是网上的基本操作,所以也就没发表在网站上,都保存在本地.最近突然发现自己博客中实战的东西太少了,决定将以前写的一些文章搬过来,由 ...
- php框架——laravel
准备工作:见php语法篇末尾环境部署 一.[创建一个表+插入数据]:手动创建或者执行sql语句 DROP TABLE IF EXISTS `tb_business`; CREATE TABLE IF ...
- Django的结构
一.Django的结构 二.静态文件的配置 STATIC_URL = '/static/' # HTML中使用的静态文件夹前缀 STATICFILES_DIRS = [ os.path.join(BA ...
- mysql(单表查询,多表查询,MySQl创建用户和授权,可视化工具Navicat的使用)
单表查询 语法: 一.单表查询的语法 SELECT 字段1,字段2... FROM 表名 WHERE 条件 GROUP BY field HAVING 筛选 ORDER BY field LIMIT ...
- CentOS6和7启动流程
CentOS6启动流程 https://linux.cn/article-8807-1.html BIOS 开机自检,硬件自检 MBR MBR磁盘分区是一种使用最为广泛的分区结构,它也被称为DOS分区 ...