zoukankan      html  css  js  c++  java
  • 【CSP-S 2019模拟】T1 —看门人(长链剖分+线段树)

    传送门


    签到水题

    长剖维护一下,每次在链顶把链底向上一个个合并轻儿子
    用一颗线段树维护单点修改和区间询问即可

    考试的时候看错题了…
    以为只找向下一条路径
    还写了个线段树合并对拍了2h2h
    我是个瓜皮

    #include<bits/stdc++.h>
    using namespace std;
    const int RLEN=1<<22|5;
    inline char gc(){
    	static char ibuf[RLEN],*ib,*ob;
    	(ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
    	return (ib==ob)?EOF:*ib++;
    }
    inline int read(){
    	char ch=gc();
    	int res=0,f=1;
    	while(!isdigit(ch))f^=ch=='-',ch=gc();
    	while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
    	return f?res:-res;
    }
    #define re register
    #define cs const
    #define ll long long
    #define pb push_back
    #define poly vector<int>
    #define pii pair<int,int>
    #define fi first
    #define se second
    inline void chemx(ll &a,ll b){a<b?a=b:0;}
    inline void chemn(int &a,int b){a>b?a=b:0;}
    cs int mod=998244353;
    inline int Add(int &a,int b){return (a+=b)>=mod?a-=mod:0;}
    inline int mul(int a,int b){return 1ll*a*b%mod;}
    cs int N=1000005;
    int n,m,pw[N];
    ll ans[N];
    vector<pii> e[N];
    vector<int> chain[N];
    int dep[N],mxdep[N],son[N],L[N],R[N],fa[N],top[N];
    ll dis[N];
    ll *f[N],tmp[N<<2],*id;
    void dfs1(int u){
    	for(int i=0;i<e[u].size();i++){
    		int v=e[u][i].fi;
    		dep[v]=dep[u]+1,dis[v]=dis[u]+e[u][i].se;
    		dfs1(v);
    		if(mxdep[v]>mxdep[son[u]])son[u]=v;
    	}
    	mxdep[u]=mxdep[son[u]]+1;
    }
    void dfs2(int u,int tp){
    	top[u]=tp;
    	if(son[u])dfs2(son[u],tp);
    	for(int i=0;i<e[u].size();i++){
    		int v=e[u][i].fi;
    		if(v==son[u])continue;
    		dfs2(v,v);
    	}
    	chain[top[u]].pb(u);
    }
    namespace Seg{
    	ll mx[N<<2],a[N];
    	int n;
    	#define lc (u<<1)
    	#define rc ((u<<1)|1)
    	#define mid ((l+r)>>1)
    	inline void pushup(int u){
    		mx[u]=max(mx[lc],mx[rc]);
    	}
    	void build(int u,int l,int r){
    		if(l==r){mx[u]=a[l];return;}
    		build(lc,l,mid),build(rc,mid+1,r);
    		pushup(u);
    	}
    	inline void init(int len,ll *f){
    		n=len-1;
    		for(int i=0;i<=n;i++)a[i]=f[i];
    		build(1,0,n);
    	}
    	void insert(int u,int l,int r,int p,ll k){
    		if(l==r){chemx(mx[u],k);return;}
    		if(p<=mid)insert(lc,l,mid,p,k);
    		else insert(rc,mid+1,r,p,k);
    		pushup(u);
    	}
    	inline void update(int p,ll k){
    		insert(1,0,n,p,k);
    	}
    	ll query(int u,int l,int r,int st,int des){
    		if(st>r||l>des)return -1;
    		if(st<=l&&r<=des)return mx[u];
    		ll res=-1;
    		if(st<=mid)chemx(res,query(lc,l,mid,st,des));
    		if(mid<des)chemx(res,query(rc,mid+1,r,st,des));
    		return res;
    	}
    	inline ll ask(int l,int r){
    		
    		ll res=query(1,0,n,l,r);
    		return res;
    	}
    	#undef lc 
    	#undef rc
    	#undef mid
    }
    inline ll Addchain(int u,int tt,int v){
    	ll res=-1;re int del=dep[v]-dep[u];
    	for(int i=0;i<mxdep[v];i++){
    		int l=dep[tt]-dep[u]+L[tt]-i-1,r=dep[tt]-dep[u]+R[tt]-i-1;
    		ll now=-1;
    		if(r>=dep[tt]-dep[u])now=Seg::ask(l,r);
    		if(now!=-1)chemx(res,f[v][i]+now);
    	}
    	for(re int i=0;i<mxdep[v];i++){
    		if(f[v][i]>f[u][del+i])f[u][del+i]=f[v][i],Seg::update(del+i,f[u][del+i]);
    	}
    	if(res!=-1)res-=2*dis[tt];
    	return res;
    }
    inline ll ask(int u,int v){
    	int del=dep[v]-dep[u];
    	ll res=Seg::ask(del+L[v],del+R[v]);
    	if(res!=-1)res-=dis[v];
    	return res;
    }
    inline void calc_ans(int u){
    	Seg::init(mxdep[u],f[u]);
    	for(re int i=0,v,j,t;i<chain[u].size();i++){
    		v=chain[u][i];
    		for(j=0;j<e[v].size();j++){
    			t=e[v][j].fi;
    			if(t==son[v])continue;
    			chemx(ans[v],Addchain(u,v,t));
    		}
    		chemx(ans[v],ask(u,v));
    	}
    }
    void dp(int u){
    	f[u][0]=dis[u];
    	if(son[u]){
    		f[son[u]]=f[u]+1,dp(son[u]);
    	}
    	for(re int i=0,v;i<e[u].size();i++){
    		v=e[u][i].fi;
    		if(v==son[u])continue;
    		f[v]=id,id+=mxdep[v];
    		dp(v);
    	}
    	if(u==top[u])calc_ans(u);
    }
    signed main(){
    	int size=100<<20;
        __asm__ ("movq %0,%%rsp
    "::"r"((char*)malloc(size)+size));//提交用这个 
    	n=read();memset(ans,-1,sizeof(ans));
    	for(int i=1;i<=n;i++)L[i]=read(),R[i]=read();
    	for(int i=2,w;i<=n;i++)
    		fa[i]=read(),w=read(),e[fa[i]].pb(pii(i,w));
    	dep[1]=pw[0]=1;
    	for(int i=1;i<=n;i++)pw[i]=mul(pw[i-1],23333);
    	dfs1(1),dfs2(1,1);
    	id=tmp,f[1]=id,id+=mxdep[1];
    	dp(1);int res=0;
    	for(int i=1;i<=n;i++)Add(res,mul(pw[n-i],(ans[i]%mod+mod)%mod));
    	cout<<res;exit(0);
    }
    
  • 相关阅读:
    SAP 关于ALV 双击 打开其它事务码
    记一次PGSQL数据库恢复过程
    SAP 供应商主数据维护(字段示例)
    SAP 高仿 CKM3 物料成本价格分析报表SQL
    SAP scc4更改设置 (不允许更改资源库和跨客户端定)
    SAP 主数据 查重常用SQL
    SAP S4/HANA 连接SQL SERVER 数据库
    SAP VBKD业务表的情况处理
    SAP 物料主数据维护说明
    SAP CO成本控制 相关表
  • 原文地址:https://www.cnblogs.com/stargazer-cyk/p/12328528.html
Copyright © 2011-2022 走看看