uo什么意思_2018五大学科集训队

uo什么意思_2018五大学科集训队uoj#418. 【集训队作业2018】三角形(线段树合并)

大家好,又见面了,我是你们的朋友全栈君。

传送门

好迷啊……膜一下ljz

考虑每个操作,如果把操作按先后顺序放到序列上的话,操作一就是把\(w_i\)的石子放到某个节点,那么就是在序列末端加入\(w_i\),然后根据贪心肯定要把它所有儿子的石子拿走,也就是要减去\(\sum w_{son}\)

那么每个点的答案就是序列的最大前缀

因为父亲节点的操作一要在儿子之后进行,很麻烦,那么可以每次在自己这里把\(w_i\)减掉,到父亲的时候再加回去

\((x,y)\)为一个二元组,\(x\)表示当前位置的最大前缀和,\(y\)表示最小后缀和,然后定义一个运算\((ax,by)=(x+max(0,a+y),b+min(0,a+y)\),大概能看出是个什么东西,注意这个运算不满足交换律

然后考虑一下这些二元组在序列中的顺序,如果\(x+y<0\),那么肯定得放前面,因为可以让之后的前缀和减小。

那么当\(x+y<0\)时,按\(x\)排序,这样能使前缀和不断减小。如果\(x+y>0\),按\(y\)排序就好了

维护二元组的话,用线段树合并,如果当前二元组的\(x+y<0\)且有比它最大前缀大的二元组,那么就已经是最优的了否则将在它前面的二元组与它合并,合并完后加进线段树里就好了。注意合并完之后节点没了要记得清空

//minamoto
#include<bits/stdc++.h>
#define R register
#define int long long
#define loli 200000000000000
#define fp(i,a,b) for(R int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(R int i=a,I=b-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){
    R int res,f=1;R char ch;
    while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
    for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
    return res*f;
}
char sr[1<<21],z[20];int C=-1,Z=0;
inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;}
void print(R int x){
    if(C>1<<20)Ot();if(x<0)sr[++C]='-',x=-x;
    while(z[++Z]=x%10+48,x/=10);
    while(sr[++C]=z[Z],--Z);sr[++C]=' ';
}
const int N=4e5+5;
struct eg{int v,nx;}e[N<<1];int head[N],tot;
inline void add(R int u,R int v){e[++tot]={v,head[u]},head[u]=tot;}
struct node{
    int mx,mn;
    node(){}
    node(R int mx,R int mn):mx(mx),mn(mn){}
    inline node operator +(const node &b){return node(mx+max(0ll,mn+b.mx),b.mn+min(0ll,mn+b.mx));}
}tr[N<<5];
int w[N],rt[N],ans[N],ls[N<<5],rs[N<<5];
int n,cnt,fa;
void ins(int &p,int l,int r,node x){
    if(!p)p=++cnt;if(l==r)return (void)(tr[p]=tr[p]+x);
    int mid=(l+r)>>1;
    x.mx<=mid?ins(ls[p],l,mid,x):ins(rs[p],mid+1,r,x);
    tr[p]=tr[ls[p]]+tr[rs[p]];
}
int merge(int x,int y,int l,int r){
    if(!x||!y)return x|y;if(l==r)return tr[x]=tr[x]+tr[y],x;
    int mid=(l+r)>>1;
    ls[x]=merge(ls[x],ls[y],l,mid),rs[x]=merge(rs[x],rs[y],mid+1,r);
    tr[x]=tr[ls[x]]+tr[rs[x]];
    return x;
}
void update(int &p,int l,int r,node &x,bool &flag){
    if(!p||flag)return;
    if(l==r){
        if(x.mx+x.mn<0&&x.mx<l)return (void)(flag=1);
        x=x+tr[p],p=0;return;
    }
    int mid=(l+r)>>1;
    update(ls[p],l,mid,x,flag),update(rs[p],mid+1,r,x,flag);
    if(ls[p]||rs[p])tr[p]=tr[ls[p]]+tr[rs[p]];
    else p=0;
}
void dfs(int u){
    node res=node(0,-w[u]);
    go(u){
        res.mx+=w[v],dfs(v);
        rt[u]=merge(rt[u],rt[v],0,loli);
    }
    node qaq=res;bool flag=0;
    qaq.mx+=w[u],ans[u]=(qaq+tr[rt[u]]).mx;
    update(rt[u],0,loli,res,flag),ins(rt[u],0,loli,res);
}
signed main(){
//  freopen("testdata.in","r",stdin);
    read(),n=read();
    fp(i,2,n)fa=read(),add(fa,i);
    fp(i,1,n)w[i]=read();
    dfs(1);
    fp(i,1,n)print(ans[i]);
    return Ot(),0;
}

转载于:https://www.cnblogs.com/bztMinamoto/p/10231616.html

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/107104.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • JS数组扁平化_扁平化js

    JS数组扁平化_扁平化jsJS 扁平化 (flatten) 数组

    2022年4月21日
    38
  • POJ 1322 Chocolate

    POJ 1322 Chocolate

    2021年12月1日
    43
  • pycharm怎么导入数据库_pycharm导入python

    pycharm怎么导入数据库_pycharm导入python记录一下自己最近一段时间可能用不上的代码,免得以后再找起来麻烦。pycharm连接数据库首先在找到Database,选择连接自己下载的数据库软件。填写完自己的数据库密码后,如果正确的话,就会出现下面这张图。点下“刷新”按钮后即可连接数据库。测试是否连接成功。到这里基本完成了外部数据导入步骤。2.代码importpymysqlimportmath#用来操作数据库的类classclimb_test_DataInput(object):#…

    2022年8月27日
    0
  • Mysql锁详解(行锁、表锁、意向锁、Gap锁、插入意向锁)

    Mysql锁详解(行锁、表锁、意向锁、Gap锁、插入意向锁)锁:对“某种范围”的数据上“某种锁”1.“某种范围”:行、表2.“某种锁”2.1共享锁SharedLocks(S锁)1、兼容性:加了S锁的记录,允许其他事务再加S锁,不允许其他事务再加X锁2、加锁方式:select…lockinsharemode2.2排他锁ExclusiveLocks(X锁)1、兼容性:加了X锁的记录,不允许其他事务再加S锁或者X锁2、加锁方式…

    2022年5月9日
    69
  • 傻瓜攻略(一)——MATLAB主成分分析(PCA)代码及结果分析实例「建议收藏」

    傻瓜攻略(一)——MATLAB主成分分析(PCA)代码及结果分析实例「建议收藏」主成分分析主成分分析法(PCA)是一种高效处理多维数据的多元统计分析方法,将主成分分析用于多指标(变量)的综合评价较为普遍。该方法的基本思想是运用较少的变量去解释原始数据中的大部分变异,通过对原始数据相关矩阵内部结构关系的分析和计算,产生一系列互不相关的新变量。根据需要从中选取比原始变量个数少的几个新变量,这些新的变量就是所谓的主成分,它们能够充分解释原始数据的变化。因此,主成分分析法本质上是…

    2022年5月29日
    42
  • VWware15安装windows XP纯净版虚拟机「建议收藏」

    VWware15安装windows XP纯净版虚拟机「建议收藏」需要准备的工具:VMware15或其他版本winXP.iso(需要的留言)安装步骤:1.创建虚拟机选择“典型”选择“安装光盘映像文件”winXP.iso选择操作系统类型选择安装位置及虚拟机名称选择磁盘大小,最少30G选择处理器数量和内存,根据个人电脑情况而定,就是虚拟机运行速度的问题完成后启动虚拟机等待进入winPE系统选择打开“分区工具…

    2022年4月19日
    502

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号