739B - Alyona and a tree倍增＋差分

最新推荐文章于 2020-12-10 23:36:11 发布
最新推荐文章于 2020-12-10 23:36:11 发布
阅读量213 收藏
点赞数
分类专栏： 树上差分
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_34927456/article/details/79476009
版权

发现cf文本读确实变快了，但是用了快读树剖还是过不了2e5，好像常数大了4？

后来才想起来这是离线的树链改值，差分就行了

跪了的树剖：

#include<bits/stdc++.h>  
//#pragma comment(linker, "/STACK:1024000000,1024000000")   
#include<stdio.h>  
#include<algorithm>  
#include<queue>  
#include<string.h>  
#include<iostream>  
#include<math.h>                    
#include<stack>
#include<set>  
#include<map>  
#include<vector>  
#include<iomanip> 
#include<bitset>
using namespace std;         //

#define ll long long  
#define ull unsigned long long
#define pb push_back  
#define FOR(a) for(int i=1;i<=a;i++) 
#define sqr(a) (a)*(a)
#define dis(a,b) sqrt(sqr(a.x-b.x)+sqr(a.y-b.y))
ll qp(ll a,ll b,ll mod){
	ll t=1;while(b){if(b&1)t=t*a%mod;b>>=1;a=a*a%mod;}return t;
}
struct DOT{ll x;ll y;};
//inline void read(int &x){int k=0;char f=1;char c=getchar();for(;!isdigit(c);c=getchar())if(c=='-')f=-1;for(;isdigit(c);c=getchar())k=k*10+c-'0';x=k*f;} 
const int dx[4]={0,0,-1,1};
const int dy[4]={1,-1,0,0};
const int inf=0x3f3f3f3f; 
const ll Linf=0x3f3f3f3f3f3f3f3f;
const ll mod=1e9+7;;

namespace fastIO{    
    #define BUF_SIZE 100000    
    #define OUT_SIZE 100000    
    #define ll long long    
    //fread->read    
    bool IOerror=0;    
    inline char nc(){    
        static char buf[BUF_SIZE],*p1=buf+BUF_SIZE,*pend=buf+BUF_SIZE;    
        if (p1==pend){    
            p1=buf; pend=buf+fread(buf,1,BUF_SIZE,stdin);    
        if (pend==p1){IOerror=1;return -1;}    
        //{printf("IO error!\n");system("pause");for (;;);exit(0);}    
    }    
    return *p1++;    
}    
inline bool blank(char ch){return ch==' '||ch=='\n'||ch=='\r'||ch=='\t';}    
inline void read(int &x){    
    bool sign=0; char ch=nc(); x=0;    
    for (;blank(ch);ch=nc());    
    if (IOerror)return;    
    if (ch=='-')sign=1,ch=nc();    
    for (;ch>='0'&&ch<='9';ch=nc())x=x*10+ch-'0';    
    if (sign)x=-x;    
}    
inline void read(ll &x){    
    bool sign=0; char ch=nc(); x=0;    
    for (;blank(ch);ch=nc());    
    if (IOerror)return;    
    if (ch=='-')sign=1,ch=nc();    
    for (;ch>='0'&&ch<='9';ch=nc())x=x*10+ch-'0';    
    if (sign)x=-x;    
}    
inline void read(double &x){    
    bool sign=0; char ch=nc(); x=0;    
    for (;blank(ch);ch=nc());    
    if (IOerror)return;    
    if (ch=='-')sign=1,ch=nc();    
    for (;ch>='0'&&ch<='9';ch=nc())x=x*10+ch-'0';    
    if (ch=='.'){    
        double tmp=1; ch=nc();    
        for (;ch>='0'&&ch<='9';ch=nc())tmp/=10.0,x+=tmp*(ch-'0');    
    }    
    if (sign)x=-x;    
}    
inline void read(char *s){    
    char ch=nc();    
    for (;blank(ch);ch=nc());    
    if (IOerror)return;    
    for (;!blank(ch)&&!IOerror;ch=nc())*s++=ch;    
    *s=0;    
}    
inline void read(char &c){    
    for (c=nc();blank(c);c=nc());    
    if (IOerror){c=-1;return;}    
}   
#undef OUT_SIZE    
#undef BUF_SIZE    
}; using namespace fastIO;  





const int maxn=2e5+9;  

int n;
 
struct EDGE{
    int v;int w;int next;
}G[maxn<<2];
int head[maxn],tot;

inline void addedge(int u,int v,int w){
    ++tot;G[tot].v=v;G[tot].next=head[u];G[tot].w=w;head[u]=tot;
}
 
int top[maxn];
int pre[maxn];
int dep[maxn];
int num[maxn];	//子树尺寸
int p[maxn];    //v的对应位置
int out[maxn];  //退出时间戳
int fp[maxn];   //访问序列
int son[maxn];  //重儿子
int pos;		//计时器

inline void init(){
    memset(head,-1,sizeof head);tot=0;
	memset(son,-1,sizeof son);
}

inline void dfs1(int u,int fa,int d){
    dep[u]=d;
    pre[u]=fa;
    num[u]=1;
    for(int i=head[u];~i;i=G[i].next){
        int v=G[i].v;
        if(v==fa)continue;
        dfs1(v,u,d+G[i].w);
        num[u]+=num[v];
        if(son[u]==-1||num[v]>num[son[u]])son[u]=v;
    }
}
inline void getpos(int u,int sp){
    top[u]=sp;
    p[u]=out[u]=++pos;

	//cout<<u<<"ww"<<p[u]<<endl;

    fp[p[u]]=u;
    if(son[u]==-1)return;
    getpos(son[u],sp);
    for(int i=head[u];~i;i=G[i].next){
        int v=G[i].v;
        if(v!=son[u]&&v!=pre[u])getpos(v,v);
    }
    out[u]=pos;
}

int a[maxn];

struct NODE{
	int l,r;
	ll sum;
	ll add;
}ST[maxn<<2];
inline void pushup(int rt){
	ST[rt].sum=ST[rt<<1].sum+ST[rt<<1|1].sum;
}
inline void pushdown(int rt){
	if(ST[rt].add){
		ST[rt<<1].sum+=ST[rt].add*(ST[rt<<1].r-ST[rt<<1].l+1);
		ST[rt<<1|1].sum+=ST[rt].add*(ST[rt<<1|1].r-ST[rt<<1|1].l+1);
		ST[rt<<1].add+=ST[rt].add;
		ST[rt<<1|1].add+=ST[rt].add;
		ST[rt].add=0;
	}
}
inline void build(int l,int r,int rt){
	if(l==r){ST[rt].l=ST[rt].r=l;ST[rt].sum=0;return;}
	ST[rt].l=l;ST[rt].r=r;
	int m=l+r>>1;
	build(l,m,rt<<1);build(m+1,r,rt<<1|1);pushup(rt);
}
inline void update(int a,int b,int c,int l,int r,int rt){
	if(a<=l&&b>=r){
		ST[rt].sum+=(r-l+1)*c;
		ST[rt].add+=c;
		return;
	}
	int m=l+r>>1;
	pushdown(rt);
	if(a<=m)update(a,b,c,l,m,rt<<1);
	if(b> m)update(a,b,c,m+1,r,rt<<1|1);
	pushup(rt);
}
inline ll query(int x,int l,int r,int rt){
	if(l==r){return ST[rt].sum;}
	int m=l+r>>1;
	pushdown(rt);
	if(x<=m)return query(x,l,m,rt<<1);
	else return query(x,m+1,r,rt<<1|1);
}

inline void solve(int u,int v){

    while(top[u]!=top[v]){

		if(dep[top[u]]<dep[top[v]])swap(u,v);
        update(p[top[u]],p[u],1,1,n,1);
		u=pre[top[u]];
    }

	if(dep[u]<dep[v])swap(u,v);
	
	//cout<<p[u]<<"ww"<<p[v]<<endl;
    update(p[v],p[u],1,1,n,1);
}

int fa[maxn][50];
ll falen[maxn][50];

int main(){
	//scanf("%d",&n);
	read(n);
	init();
	for(int i=1;i<=n;i++){
		//scanf("%d",&a[i]);
		read(a[i]);
	}
	for(int i=2,x,y;i<=n;i++){
		//scanf("%d%d",&x,&y);
		read(x);read(y);
		addedge(x,i,y);
		fa[i][0]=x;
		falen[i][0]=y;
	}
	dfs1(1,1,0);
	getpos(1,1);
	build(1,n,1);

	for(int i=0;i<=32;i++){
		fa[1][i]=0;
		falen[0][i]=Linf;
		falen[1][i]=Linf;
	}
	for(int i=1;i<=32;i++){
		for(int j=1;j<=n;j++){
			fa[j][i]=fa[fa[j][i-1]][i-1];
			if(falen[fa[j][i-1]][i-1]==Linf){
				falen[j][i]=Linf;continue;
			}
			falen[j][i]=falen[fa[j][i-1]][i-1]+falen[j][i-1];
		}
	}

	for(int i=2;i<=n;i++){
		int now=i;
		for(int j=32;j>=0;j--){
			if(falen[now][j] <= a[i]){
				a[i]-=falen[now][j];
				now=fa[now][j];
			}
			//cout<<i<<"ww"<<now<<"ww"<<a[i]<<endl;
		}

		if(i==now)continue;
		
		solve(fa[i][0],now);	
	}
	for(int i=1;i<=n;i++)printf("%lld ",query(p[i],1,n,1));
	
}

Ａ了的差分，都用的倍增

#include<bits/stdc++.h>  
//#pragma comment(linker, "/STACK:1024000000,1024000000")   
#include<stdio.h>  
#include<algorithm>  
#include<queue>  
#include<string.h>  
#include<iostream>  
#include<math.h>                    
#include<stack>
#include<set>  
#include<map>  
#include<vector>  
#include<iomanip> 
#include<bitset>
using namespace std;         //

#define ll long long  
#define ull unsigned long long
#define pb push_back  
#define FOR(a) for(int i=1;i<=a;i++) 
#define sqr(a) (a)*(a)
#define dis(a,b) sqrt(sqr(a.x-b.x)+sqr(a.y-b.y))
ll qp(ll a,ll b,ll mod){
	ll t=1;while(b){if(b&1)t=t*a%mod;b>>=1;a=a*a%mod;}return t;
}
struct DOT{ll x;ll y;};
inline void read(int &x){int k=0;char f=1;char c=getchar();for(;!isdigit(c);c=getchar())if(c=='-')f=-1;for(;isdigit(c);c=getchar())k=k*10+c-'0';x=k*f;} 
const int dx[4]={0,0,-1,1};
const int dy[4]={1,-1,0,0};
const int inf=0x3f3f3f3f; 
const ll Linf=0x3f3f3f3f3f3f3f3f;
const ll mod=1e9+7;;

const int maxn=2e5+9;  

int n;
 
struct EDGE{
    int v;int w;int next;
}G[maxn<<2];
int head[maxn],tot;

inline void addedge(int u,int v,int w){
    ++tot;G[tot].v=v;G[tot].next=head[u];G[tot].w=w;head[u]=tot;
}
 
/*
int top[maxn];
int pre[maxn];
int dep[maxn];
int num[maxn];	//子树尺寸
int p[maxn];    //v的对应位置
int out[maxn];  //退出时间戳
int fp[maxn];   //访问序列
int son[maxn];  //重儿子
int pos;		//计时器
*/
inline void init(){
    memset(head,-1,sizeof head);tot=0;
	//memset(son,-1,sizeof son);
}
/*
inline void dfs1(int u,int fa,int d){
    dep[u]=d;
    pre[u]=fa;
    num[u]=1;
    for(int i=head[u];~i;i=G[i].next){
        int v=G[i].v;
        if(v==fa)continue;
        dfs1(v,u,d+G[i].w);
        num[u]+=num[v];
        if(son[u]==-1||num[v]>num[son[u]])son[u]=v;
    }
}
inline void getpos(int u,int sp){
    top[u]=sp;
    p[u]=out[u]=++pos;

    fp[p[u]]=u;
    if(son[u]==-1)return;
    getpos(son[u],sp);
    for(int i=head[u];~i;i=G[i].next){
        int v=G[i].v;
        if(v!=son[u]&&v!=pre[u])getpos(v,v);
    }
    out[u]=pos;
}
*/
ll add[maxn];
int a[maxn];

/*
struct NODE{
	int l,r;
	ll sum;
	ll add;
}ST[maxn<<2];
inline void pushup(int rt){
	ST[rt].sum=ST[rt<<1].sum+ST[rt<<1|1].sum;
}
inline void pushdown(int rt){
	if(ST[rt].add){
		ST[rt<<1].sum+=ST[rt].add*(ST[rt<<1].r-ST[rt<<1].l+1);
		ST[rt<<1|1].sum+=ST[rt].add*(ST[rt<<1|1].r-ST[rt<<1|1].l+1);
		ST[rt<<1].add+=ST[rt].add;
		ST[rt<<1|1].add+=ST[rt].add;
		ST[rt].add=0;
	}
}
inline void build(int l,int r,int rt){
	if(l==r){ST[rt].l=ST[rt].r=l;ST[rt].sum=0;return;}
	ST[rt].l=l;ST[rt].r=r;
	int m=l+r>>1;
	build(l,m,rt<<1);build(m+1,r,rt<<1|1);pushup(rt);
}
inline void update(int a,int b,int c,int l,int r,int rt){
	if(a<=l&&b>=r){
		ST[rt].sum+=(r-l+1)*c;
		ST[rt].add+=c;
		return;
	}
	int m=l+r>>1;
	pushdown(rt);
	if(a<=m)update(a,b,c,l,m,rt<<1);
	if(b> m)update(a,b,c,m+1,r,rt<<1|1);
	pushup(rt);
}
inline ll query(int x,int l,int r,int rt){
	if(l==r){return ST[rt].sum;}
	int m=l+r>>1;
	pushdown(rt);
	if(x<=m)return query(x,l,m,rt<<1);
	else return query(x,m+1,r,rt<<1|1);
}

inline void solve(int u,int v){

    while(top[u]!=top[v]){

		if(dep[top[u]]<dep[top[v]])swap(u,v);
        update(p[top[u]],p[u],1,1,n,1);
		u=pre[top[u]];
    }

	if(dep[u]<dep[v])swap(u,v);
	
	//cout<<p[u]<<"ww"<<p[v]<<endl;
    update(p[v],p[u],1,1,n,1);
}
*/

int fa[maxn][50];
ll falen[maxn][50];

void dfs(int u){

	for(int i=head[u];~i;i=G[i].next){
		dfs(G[i].v);
		add[u]+=add[G[i].v];
	}
}

int main(){
	//scanf("%d",&n);
	read(n);
	init();
	for(int i=1;i<=n;i++){
		//scanf("%d",&a[i]);
		read(a[i]);
	}
	for(int i=2,x,y;i<=n;i++){
		//scanf("%d%d",&x,&y);
		read(x);read(y);
		addedge(x,i,y);
		fa[i][0]=x;
		falen[i][0]=y;
	}
	//dfs1(1,1,0);
	//getpos(1,1);
	//build(1,n,1);

	for(int i=0;i<=32;i++){
		fa[1][i]=0;
		falen[0][i]=Linf;
		falen[1][i]=Linf;
	}
	for(int i=1;i<=32;i++){
		for(int j=1;j<=n;j++){
			fa[j][i]=fa[fa[j][i-1]][i-1];
			if(falen[fa[j][i-1]][i-1]==Linf){
				falen[j][i]=Linf;continue;
			}
			falen[j][i]=falen[fa[j][i-1]][i-1]+falen[j][i-1];
		}
	}

	for(int i=2;i<=n;i++){
		int now=i;
		for(int j=32;j>=0;j--){
			if(falen[now][j] <= a[i]){
				a[i]-=falen[now][j];
				now=fa[now][j];
			}
			//cout<<i<<"ww"<<now<<"ww"<<a[i]<<endl;
		}

		if(i==now)continue;
		
		//solve(fa[i][0],now);	
		add[fa[i][0]]++;add[fa[now][0]]--;
	}
	dfs(1);
	for(int i=1;i<=n;i++)printf("%lld ",add[i]);
	//for(int i=1;i<=n;i++)printf("%lld ",query(p[i],1,n,1));
	
}


Drenight
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
vue-20210119_alyona75
02-18
КурсVue.js 。 Использованиезадачника Инструкцияпоработенепосредственносзадачникомподробноописанав00 。...
codeforces 739B B. Alyona and a tree
ZT在努力
12-01 1471
这个题目想了好久都没有想出好的方法，最终看了别人的思路解决的。写一篇博客纪念一下做这个题的想法和学到的一些新东西。读这个题的时候以为会是dp，仔细想了想，又不太符合dp的特征。然后，就按照数据结构的想法思考，想了好久，还是没什么思路，而且用队列，从叶子节点模拟了一发，以为可以过。提交的时候才想起来，肯定会TLE呀，结果果断TLE，想想也是醉了。最后想着想着，觉得时间复杂度怎么也应该时n*lgn才能
B. Alyona and a tree
weixin_30402085的博客
08-20 110
二分+dfs序+前缀和+瞎几把yy #include "a.h" int n; aLL a; struct edge { int to, nx, w; } e[N << 1]; arr fi; int ce = 1; aLL sum, ans; aLL vec, up, down, res; ll tot = 1; ll rs = ...
Codeforces 739B Alyona and a tree （树上差分+二分）
K键盘里的青春K
04-19 1083
D. Alyona and a tree time limit per test 2 seconds memory limit per test 256 megabytes input standard input output standard output Alyona has a tree with n vertices. The
CF - 739B (Alyona and a tree)
Joker & Liar的博客
03-10 156
题意：给定一颗树(1为根节点)，每个点有个权值，每条边也有边权，现在给出所有的权值和边权，定义一个点u控制另一个点v的条件： ① u是 v的祖先节点； ② u 到 v 上的边权之和 sum 不大于 v 的点权 a[v] ，即 sum<=a[v]； 求每个点控制的节点的数量； 分析：先考虑对于每个子节点 v ，找到它最大的祖先 u 满足边权之和 sum<=a[v]，那么...
[CF739B] Alyona and a tree 详解
banmi4671的博客
03-12 278
B. Alyona and a tree Alyona has a tree with n vertices. The root of the tree is the vertex 1. In each vertex Alyona wrote an positive integer, in the vertex i she wrote ai. Moreover, the girl wro...
string2wav:字符串到 wav。 “text2wave”应用程序的模拟，基于 Microsoft SAPI（语音 API），通过选择的语音从文件或字符串中读出文本
06-09
最常用于俄语： Alyona（阿卡佩拉集团） Katerina 和 Milena (ScanSoft - RealSpeak) 对于工作，您必须安装 SAPI（Windows 7 及更高版本默认包括 SAPI）在 Windows 7、Windows 2008 和 SAPI 5.4 上进行测试用法 ： ...
Codeforces739B. Alyona and a tree（树+二分）
KeyboardPianist的博客
12-18 1439
B. Alyona and a tree time limit per test 2 seconds memory limit per test 256 megabytes input standard input output standard output Alyona has a tree with n vertices. The root of the tree is t
【题解】CF739B Alyona and a tree
鹭天
12-05 161
题目描述： [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-dM1R7HPn-1607157484643)(C:\Users\Stuedent\AppData\Roaming\Typora\typora-user-images\image-20201127211106195.png)] Solution 一、O(n2)O(n^2)O(n2)暴力 枚举每一个点，扫一遍她的子树，合法的点累加进答案。 二、O(n2)O(n^2)O(n2)暴力 我们换一个角度来打暴力 考虑每一个点对
Alyona and a tree 树上差分 +倍增
tlyzxc的博客
12-10 101
传送门 分析 这道题跟我前段时间写的一道数位DP的思路有点一致，但是居然没想出 来首先我们确定一下，任何一个平衡数的支点，都仅存在一个支点，所以，我们可以去枚举每一个支点，计算每一个支点下，平衡数的数量，然后相加即可 最后只需要注意一下，如果每一位都取0，那么任何一位都可以作为支点，但合法的数字只有一个0，所以最后答案需要特殊处理一下 代码 #include <iostream> #include <cstdio> #include <cmath> #include &l
Codeforces 740D Alyona and a tree 二分+树上差分
Masamiiiiii
04-13 632
点击打开链接 题意:n个结点,带权的树,v控制u当且仅当 u为v的子孙&& dist[u,v] 设d[i]为i->rt的距离,dist[u,v]=d[u]-d[v],根据u来更新答案,控制u的v必须满足 ,d[v]>=d[u]-a[u]  d[i]随着深度递增,若在dfs过程中保存u的祖先的d值,可以二分找到第一个满足的pos,则ans[pos~fa[u]]都应该+1 怎么判断祖先
C. Alyona and the Tree
ao9795的博客
08-28 102
http://codeforces.com/contest/682/problem/C dfs，没什么好说的，直接看代码吧 1 import java.util.ArrayList; 2 import java.util.Scanner; 3 4 5 public class Main { 6 static long[] a =...
#Codeforces 381 [div2] D. Alyona and a tree 【树链剖分 + 优先队列】
算法本子
11-24 1253
题目： D. Alyona and a tree time limit per test  2 seconds memory limit per test  256 megabytes input  standard input output  standard output Alyona has atree with n vertices. The root of the tree
Codeforces Round #381 (Div. 1) B Alyona and a tree 树状数组
来自Fop_zz的辣鸡题解或教程
06-19 396
B. Alyona and a tree 题意：给定一棵树，树上的边有权值为val[i],点有权值为a[i],定义dist(a,b)为a到b的路径上的边权的和 定义v控制u，当且仅当v是u的祖先且dist(u,v) 第一反应树形dp嘛，可行，不过感觉略麻烦。。。而且二分的还要判断一些奇怪的东西（好吧懒得打 因为u是v的后代嘛。所以dist(u,v)=dist(1,u)-di
CF739B Alyona and a tree 树上差分+倍增
liangzihao1的博客
08-28 193
题目描述： Alyona有一棵有 nnn 个节点的树。这棵树的根节点是 111。在每个节点里，Alyona写了一个正整数，在节点 iii 她写了正整数 aia_iai​。另外，她在这棵树上的每条边上写了一个正整数（不同边上可能有不同的数）。 让我们定义 dist(v,u)dist(v,u)dist(v,u) 作为从 vvv 到 uuu 的简单路径上的边权和。 当且仅当 uuu 在 vvv 的子树中并且 dist(v,u)≤audist(v,u)≤a_udist(v,u)≤au​，顶点 vvv 控制顶点 u(
CF 739B - Alyona and a tree
Green_Days
11-27 416
题目大意：一棵树有n(n≤2∗105)n(n\leq 2*10^5)个点，编号为1−n1-n，固定11号节点为根。每个点有一个值aia_i，树上的每一条边的长度为wiw_i。对每一个点定义一个支配值cic_i，cic_i为ii点为根的子树中的满足到ii的距离不超过aja_j的点j的个数。 求每个点的cic_i。题解：比较容易想到，对于每一个点ii，二分它的祖先节点，找到第一个到ii的距离大于aia_
Codeforces Round #381 (Div. 2) D. Alyona and a tree
survivorone的博客
11-25 804
D. Alyona and a tree time limit per test 2 seconds memory limit per test 256 megabytes Alyona has a tree with n vertices. The root of the tree is the vertex 1. In each vertex Alyona wr
codeforces 739b Alyona and a tree
rua~
11-27 882
题意给定一棵树，每个节点有一个值a(u)，每条边有一个权值w，定义节点u控制节点v当且仅当dis(u,v) <= a(v)。要求每个节点控制的点数。 链接思路首先求出每个节点到根节点的前缀边权和pre[u],那么dis(u,v) = dis[v] - dis[u] -> dis[v] - dis[u] <= a[v] -> dis[v] - a[v] <= dis[u].问题转化为了求树上节点的子

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值