SP6779 GSS7 - Can you answer these queries VII-程序员宅基地

技术标签: 数据结构与算法  

纯数据结构题,没有思维难度。直接用线段树求最大子段和的方法完成树上路径的合并。注意链上合并顺序要符合序列的前后顺序。

#include <cstdio>
#include <cstring>
#define cd w<<1
const int S=200030;
int n,Q,h[S],v[S],nx[S],d[S],e[S],s[S],eg=1,fa[S],t[S],f[S],_=0,rk[S],a[S],dep[S];
inline int ma(int a,int b){return a>b?a:b;}
inline int mi(int a,int b){return a>b?b:a;}
inline void swap(int &a,int &b){a^=b^=a^=b;}
inline void egadd(int uu,int vv)
{
    nx[++eg]=h[uu];h[uu]=eg;
    v[eg]=vv;
}
struct node
{
    int a,ls,rs,s,tag;
    
};
inline node merge(node b,node c)
{
    node a;
    a.ls=ma(b.ls,b.s+c.ls);
    a.rs=ma(c.rs,c.s+b.rs);
    a.s=b.s+c.s;
    a.a=ma(b.a,c.a);
    a.a=ma(ma(a.ls,a.rs),ma(b.rs+c.ls,a.a));
    a.a=ma(a.a,0);
    a.ls=ma(a.ls,0);
    a.rs=ma(a.rs,0);
    return a;
}   
const int inf=(1<<30)-1;
struct Sgtr
{
    node sg[S<<2];
    inline int pos(int l,int r){return (l+r)|(l!=r);}
    void build(int l,int r)
    {
        int w=pos(l,r);
        // printf("Build (%d,%d)=%d\n",l,r,w);
        if (l==r)
        {
            int x=a[rk[l]];
            // printf("--build (%d %d)\n",rk[l],x);
            if (x>0)
                sg[w]=(node){x,x,x,x,inf};
            else
                sg[w]=(node){0,0,0,x,inf};
            return;
        }
        int mid=(l+r)>>1,lc=pos(l,mid),rc=pos(mid+1,r);
        build(l,mid);
        build(mid+1,r);
        sg[w]=merge(sg[lc],sg[rc]);
        sg[w].tag=inf;
    }
    void init()
    {
        memset(sg,0,sizeof(sg));
        build(1,n);
    }
    inline void pushdown(int w,int l,int mid,int r)
    {
        int lc=pos(l,mid),rc=pos(mid+1,r),x=sg[w].tag;
        if (x!=inf)
        {
            sg[lc].tag=sg[rc].tag=x;
            if (x>=0)
            {
                sg[lc]=(node){(mid-l+1)*x,(mid-l+1)*x,(mid-l+1)*x,(mid-l+1)*x,x};
                sg[rc]=(node){(r-mid)*x,(r-mid)*x,(r-mid)*x,(r-mid)*x,x};
            }
            else
            {
                sg[lc]=(node){0,0,0,(mid-l+1)*x,x};
                sg[rc]=(node){0,0,0,(r-mid)*x,x};
            }
        }
        sg[w].tag=inf;
    }
    void modify(int l,int r,int ll,int rr,int x)
    {
        int w=pos(l,r);
        // printf("MODIFY %d %d %d %d %d\n",l,r,ll,rr,x);
        if (ll<=l && r<=rr)
        {
            if (x>=0)
            {
                sg[w].a=sg[w].ls=sg[w].rs=sg[w].s=(r-l+1)*x;
                sg[w].tag=x;
            }
            else
            {
                sg[w].a=sg[w].ls=sg[w].rs=0;
                sg[w].s=(r-l+1)*x;
                sg[w].tag=x;
            }
            return;
        }
        int mid=(l+r)>>1,lc=pos(l,mid),rc=pos(mid+1,r);
        pushdown(w,l,mid,r);
        if (ll<=mid) modify(l,mid,ll,rr,x);
        if (rr>mid) modify(mid+1,r,ll,rr,x);
        sg[w]=merge(sg[lc],sg[rc]);
        sg[w].tag=inf;
    }
    node calcu(int l,int r,int ll,int rr)
    {
        int w=pos(l,r);
        // printf("(%d %d)=%d  %d %d\n",l,r,w,ll,rr);
        if (ll<=l && r<=rr)
        {
            return sg[w];
        }
        int mid=(l+r)>>1;
        pushdown(w,l,mid,r);
        if (ll>mid) return calcu(mid+1,r,ll,rr);
        if (rr<=mid) return calcu(l,mid,ll,rr);
        node a,b,c;
        b=calcu(l,mid,ll,rr);
        c=calcu(mid+1,r,ll,rr);
        a=merge(b,c);
        return a;
    }
}sg;
void dfs_1(int x)
{
    s[x]=1;
    for (int i=h[x];i;i=nx[i])
        if (v[i]!=fa[x])
        {
            fa[v[i]]=x;
            dep[v[i]]=dep[x]+1;
            dfs_1(v[i]);
            s[x]+=s[v[i]];
            if (s[v[i]]>s[e[x]])
                e[x]=v[i];
        }
}
void dfs_2(int x,int top)
{
    f[x]=++_;
    rk[_]=x;
    t[x]=top;
    if (!e[x]) return;
    dfs_2(e[x],top);
    for (int i=h[x];i;i=nx[i])
        if (v[i]!=fa[x] && v[i]!=e[x])
            dfs_2(v[i],v[i]);
}
void _modify(int x,int y,int o)
{
    while (t[x]!=t[y])
    {
        if (dep[t[x]]<dep[t[y]])
            swap(x,y);
        sg.modify(1,n,f[t[x]],f[x],o);
        // printf("modify %d~%d\n",t[x],x);
        x=fa[t[x]];
    }
    if (dep[x]>dep[y])
        swap(x,y);
    // printf("modify %d~%d\n",x,y);
    sg.modify(1,n,f[x],f[y],o);
}
int _calcu(int x,int y)
{
    node hida,migi,ans;
    hida=migi=(node){0,0,0,0,inf};
    while (t[x]!=t[y])
    {
        if (dep[t[x]]>dep[t[y]])
        {
            // printf("Calcu: x:%d~%d   value=%d \n",f[t[x]],f[x],hida.a);
            // printf("o=%d ",o.a);
            hida=merge(sg.calcu(1,n,f[t[x]],f[x]),hida);
            // printf(" v2=%d\n",hida.a);
            x=fa[t[x]];
        }
        else
        {
            migi=merge(sg.calcu(1,n,f[t[y]],f[y]),migi);
            // printf("Calcu: y:%d~%d\n",t[y],y);
            y=fa[t[y]];
        }
    }
    if (dep[x]<dep[y])
    {
        migi=merge(sg.calcu(1,n,f[x],f[y]),migi);
        y=x;
    }
    else
    {
        hida=merge(sg.calcu(1,n,f[y],f[x]),hida);
        x=y;
    }
    // printf("Lca=%d\n",x);
    swap(hida.ls,hida.rs);
    ans=merge(hida,migi);
    return ans.a;
}
int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
        scanf("%d",a+i);
    for (int i=1,uu,vv;i<n;i++)
    {
        scanf("%d%d",&uu,&vv);
        egadd(uu,vv);
        egadd(vv,uu);
    }
    dfs_1(1);
    dfs_2(1,1);
    sg.init();
    // for (int i=1;i<=n;i++)
        // printf("(f=%d,%d,%d) ",f[i],a[i],sg.sg[sg.pos(f[i],f[i])].a);puts("");
    for (scanf("%d",&Q);Q--;)
    {
        int o,x,y,z;
        scanf("%d%d%d",&o,&x,&y);
        if (o==1)
            printf("%d\n",_calcu(x,y));
        else
        {
            scanf("%d",&z);
            _modify(x,y,z);
        }
        // puts("---------------");sg.travel(1,n);
    }
}

转载于:https://www.cnblogs.com/Algebra-hy/p/11145852.html

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weixin_30732825/article/details/97463883

智能推荐

React学习记录-程序员宅基地

文章浏览阅读936次,点赞22次,收藏26次。React核心基础

Linux查磁盘大小命令,linux系统查看磁盘空间的命令是什么-程序员宅基地

文章浏览阅读2k次。linux系统查看磁盘空间的命令是【df -hl】,该命令可以查看磁盘剩余空间大小。如果要查看每个根路径的分区大小,可以使用【df -h】命令。df命令以磁盘分区为单位查看文件系统。本文操作环境:red hat enterprise linux 6.1系统、thinkpad t480电脑。(学习视频分享:linux视频教程)Linux 查看磁盘空间可以使用 df 和 du 命令。df命令df 以磁..._df -hl

Office & delphi_range[char(96 + acolumn) + inttostr(65536)].end[xl-程序员宅基地

文章浏览阅读923次。uses ComObj;var ExcelApp: OleVariant;implementationprocedure TForm1.Button1Click(Sender: TObject);const // SheetType xlChart = -4109; xlWorksheet = -4167; // WBATemplate xlWBATWorksheet = -4167_range[char(96 + acolumn) + inttostr(65536)].end[xlup]

若依 quartz 定时任务中 service mapper无法注入解决办法_ruoyi-quartz无法引入ruoyi-admin的service-程序员宅基地

文章浏览阅读2.3k次。上图为任务代码,在任务具体执行的方法中使用,一定要写在方法内使用SpringContextUtil.getBean()方法实例化Spring service类下边是ruoyi-quartz模块中util/SpringContextUtil.java(已改写)import org.springframework.beans.BeansException;import org.springframework.context.ApplicationContext;import org.s..._ruoyi-quartz无法引入ruoyi-admin的service

CentOS7配置yum源-程序员宅基地

文章浏览阅读2w次,点赞10次,收藏77次。yum,全称“Yellow dog Updater, Modified”,是一个专门为了解决包的依赖关系而存在的软件包管理器。可以这么说,yum 是改进型的 RPM 软件管理器,它很好的解决了 RPM 所面临的软件包依赖问题。yum 在服务器端存有所有的 RPM 包,并将各个包之间的依赖关系记录在文件中,当管理员使用 yum 安装 RPM 包时,yum 会先从服务器端下载包的依赖性文件,通过分析此文件从服务器端一次性下载所有相关的 RPM 包并进行安装。_centos7配置yum源

智能科学毕设分享(算法) 基于深度学习的抽烟行为检测算法实现(源码分享)-程序员宅基地

文章浏览阅读828次,点赞21次,收藏8次。今天学长向大家分享一个毕业设计项目毕业设计 基于深度学习的抽烟行为检测算法实现(源码分享)毕业设计 深度学习的抽烟行为检测算法实现通过目前应用比较广泛的 Web 开发平台,将模型训练完成的算法模型部署,部署于 Web 平台。并且利用目前流行的前后端技术在该平台进行整合实现运营车辆驾驶员吸烟行为检测系统,方便用户使用。本系统是一种运营车辆驾驶员吸烟行为检测系统,为了降低误检率,对驾驶员视频中的吸烟烟雾和香烟目标分别进行检测,若同时检测到则判定该驾驶员存在吸烟行为。进行流程化处理,以满足用户的需要。

随便推点

STM32单片机示例:多个定时器同步触发启动_stm32 定时器同步-程序员宅基地

文章浏览阅读3.7k次,点赞3次,收藏14次。多个定时器同步触发启动是一种比较实用的功能,这里将对此做个示例说明。_stm32 定时器同步

android launcher分析和修改10,Android Launcher分析和修改9——Launcher启动APP流程(转载)...-程序员宅基地

文章浏览阅读348次。出处 : http://www.cnblogs.com/mythou/p/3187881.html本来想分析AppsCustomizePagedView类,不过今天突然接到一个临时任务。客户反馈说机器界面的图标很难点击启动程序,经常点击了没有反应,Boss说要优先解决这问题。没办法,只能看看是怎么回事。今天分析一下Launcher启动APP的过程。从用户点击到程序启动的流程,下面针对WorkSpa..._回调bubbletextview

Ubuntu 12 最快的两个源 个人感觉 163与cn99最快 ubuntu安装源下包过慢_un.12.cc-程序员宅基地

文章浏览阅读6.2k次。Ubuntu 12 最快的两个源 个人感觉 163与cn99最快 ubuntu下包过慢 1、首先备份Ubuntu 12.04源列表 sudo cp /etc/apt/sources.list /etc/apt/sources.list.backup (备份下当前的源列表,有备无患嘛) 2、修改更新源 sudo gedit /etc/apt/sources.list (打开Ubuntu 12_un.12.cc

vue动态路由(权限设置)_vue动态路由权限-程序员宅基地

文章浏览阅读5.8k次,点赞6次,收藏86次。1.思路(1)动态添加路由肯定用的是addRouter,在哪用?(2)vuex当中获取到菜单,怎样展示到界面2.不管其他先试一下addRouter找到router/index.js文件,内容如下,这是我自己先配置的登录路由现在先不管请求到的菜单是什么样,先写一个固定的菜单通过addRouter添加添加以前注意:addRoutes()添加的是数组在export defult router的上一行图中17行写下以下代码var addRoute=[ { path:"/", name:"_vue动态路由权限

JSTL 之变量赋值标签-程序员宅基地

文章浏览阅读8.9k次。 关键词: JSTL 之变量赋值标签 /* * Author Yachun Miao * Created 11-Dec-06 */关于JSP核心库的set标签赋值变量,有两种方式: 1.日期" />2. 有种需求要把ApplicationResources_zh_CN.prope

VGA带音频转HDMI转换芯片|VGA转HDMI 转换器方案|VGA转HDMI1.4转换器芯片介绍_vga转hdmi带音频转换器,转接头拆解-程序员宅基地

文章浏览阅读3.1k次,点赞3次,收藏2次。1.1ZY5621概述ZY5621是VGA音频到HDMI转换器芯片,它符合HDMI1.4 DV1.0规范。ZY5621也是一款先进的高速转换器,集成了MCU和VGA EDID芯片。它还包含VGA输入指示和仅音频到HDMI功能。进一步降低系统制造成本,简化系统板上的布线。ZY5621方案设计简单,且可以完美还原输入端口的信号,此方案设计广泛应用于投影仪、教育多媒体、视频会议、视频展台、工业级主板显示、手持便携设备、转换盒、转换线材等产品设计上面。1.2 ZY5621 特性内置MCU嵌入式VGA_vga转hdmi带音频转换器,转接头拆解

推荐文章

热门文章

相关标签