BZOJ 2836 魔法树 树链剖分

题目大意:维护一棵有根树,每个节点初始权值为0,支持下列操作:

1.链上+

2.子树求和

。。。。。链剖裸题- -

果然链剖这种东西想要1A实在是不咋现实- -

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 100100
using namespace std;
struct Segtree{
	Segtree *ls,*rs;
	long long val,mark;
	void* operator new (size_t)
	{
		static Segtree mempool[M<<1],*C=mempool;
		return C++;
	}
	void Add(int x,int y,long long val)
	{
		this->val+=val*(y-x+1);
		mark+=val;
	}
	void Push_Down(int x,int y)
	{
		int mid=x+y>>1;
		ls->Add(x,mid,mark);
		rs->Add(mid+1,y,mark);
		mark=0;
	}
	void Build_Tree(int x,int y)
	{
		int mid=x+y>>1;
		if(x==y) return ;
		(ls=new Segtree)->Build_Tree(x,mid);
		(rs=new Segtree)->Build_Tree(mid+1,y);
	}
	void Add(int x,int y,int l,int r,long long val)
	{
		int mid=x+y>>1;
		if(x==l&&y==r)
		{
			Add(x,y,val);
			return ;
		}
		Push_Down(x,y);
		if(r<=mid)
			ls->Add(x,mid,l,r,val);
		else if(l>mid)
			rs->Add(mid+1,y,l,r,val);
		else
			ls->Add(x,mid,l,mid,val),rs->Add(mid+1,y,mid+1,r,val);
		this->val=ls->val+rs->val;
	}
	long long Query(int x,int y,int l,int r)
	{
		int mid=x+y>>1;
		if(x==l&&y==r)
			return val;
		Push_Down(x,y);
		if(r<=mid)
			return ls->Query(x,mid,l,r);
		if(l>mid)
			return rs->Query(mid+1,y,l,r);
		return ls->Query(x,mid,l,mid) + rs->Query(mid+1,y,mid+1,r);
	}
}*tree=new Segtree;
struct abcd{
	int to,next;
}table[M<<1];
int head[M],tot;
int n,m;
int fa[M],son[M],dpt[M],size[M],pos[M],top[M],cnt;
void Add(int x,int y)
{
	table[++tot].to=y;
	table[tot].next=head[x];
	head[x]=tot;
}
void DFS1(int x)
{
	int i;
	dpt[x]=dpt[fa[x]]+1;
	size[x]=1;
	for(i=head[x];i;i=table[i].next)
		if(table[i].to!=fa[x])
		{
			fa[table[i].to]=x;
			DFS1(table[i].to);
			size[x]+=size[table[i].to];
			if(size[table[i].to]>size[son[x]])
				son[x]=table[i].to;
		}
}
void DFS2(int x)
{
	int i;
	if(son[fa[x]]==x)
		top[x]=top[fa[x]];
	else
		top[x]=x;
	pos[x]=++cnt;
	if(son[x]) DFS2(son[x]);
	for(i=head[x];i;i=table[i].next)
		if(table[i].to!=fa[x]&&table[i].to!=son[x])
			DFS2(table[i].to);
}
void Modify(int x,int y,int z)
{
	int fx=top[x],fy=top[y];
	while(fx!=fy)
	{
		if(dpt[fx]<dpt[fy])
			swap(x,y),swap(fx,fy);
		tree->Add(1,n,pos[fx],pos[x],z);
		x=fa[fx];fx=top[x];
	}
	if(dpt[x]<dpt[y])
		swap(x,y);
	tree->Add(1,n,pos[y],pos[x],z);
}
int main()
{
	int i,x,y,z;
	char p[10];
	cin>>n;
	for(i=1;i<n;i++)
	{
		scanf("%d%d",&x,&y);
		++x;++y;
		Add(x,y);Add(y,x);
	}
	DFS1(1);DFS2(1);
	tree->Build_Tree(1,n);
	cin>>m;
	for(i=1;i<=m;i++)
	{
		scanf("%s",p);
		if(p[0]=='A')
		{
			scanf("%d%d%d",&x,&y,&z);
			++x;++y;
			Modify(x,y,z);
		}
		else
		{
			scanf("%d",&x);++x;
			printf("%lld\n",tree->Query(1,n,pos[x],pos[x]+size[x]-1) );
		}
	}
	return 0;
}


  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
题目描述 有一个 $n$ 个点的棋盘,每个点上有一个数字 $a_i$,你需要从 $(1,1)$ 走到 $(n,n)$,每次只能往右或往下走,每个格子只能经过一次,路径上的数字和为 $S$。定义一个点 $(x,y)$ 的权值为 $a_x+a_y$,求所有满足条件的路径中,所有点的权值和的最小值。 输入格式 第一行一个整数 $n$。 接下来 $n$ 行,每行 $n$ 个整数,表示棋盘上每个点的数字。 输出格式 输出一个整数,表示所有满足条件的路径中,所有点的权值和的最小值。 数据范围 $1\leq n\leq 300$ 输入样例 3 1 2 3 4 5 6 7 8 9 输出样例 25 算法1 (形dp) $O(n^3)$ 我们可以先将所有点的权值求出来,然后将其看作是一个有权值的图,问题就转化为了在这个图中求从 $(1,1)$ 到 $(n,n)$ 的所有路径中,所有点的权值和的最小值。 我们可以使用形dp来解决这个问题,具体来说,我们可以将这个图看作是一棵,每个点的父节点是它的前驱或者后继,然后我们从根节点开始,依次向下遍历,对于每个节点,我们可以考虑它的两个儿子,如果它的两个儿子都被遍历过了,那么我们就可以计算出从它的左儿子到它的右儿子的路径中,所有点的权值和的最小值,然后再将这个值加上当前节点的权值,就可以得到从根节点到当前节点的路径中,所有点的权值和的最小值。 时间复杂度 形dp的时间复杂度是 $O(n^3)$。 C++ 代码 算法2 (动态规划) $O(n^3)$ 我们可以使用动态规划来解决这个问题,具体来说,我们可以定义 $f(i,j,s)$ 表示从 $(1,1)$ 到 $(i,j)$ 的所有路径中,所有点的权值和为 $s$ 的最小值,那么我们就可以得到如下的状态转移方程: $$ f(i,j,s)=\min\{f(i-1,j,s-a_{i,j}),f(i,j-1,s-a_{i,j})\} $$ 其中 $a_{i,j}$ 表示点 $(i,j)$ 的权值。 时间复杂度 动态规划的时间复杂度是 $O(n^3)$。 C++ 代码
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值