zoukankan      html  css  js  c++  java
  • 联考20200521 T3 cac



    分析:
    胡乱分析一波就是两点经过的所有环上点权值加
    于是可以想到圆方树+树链剖分维护
    对于每个环和每个桥都建立一个方点
    路径上圆方点全部加权
    查询一个点的权值可以将他相邻方点的权值和减去他自身的权值(如果不是端点,一个点被经过,他相邻的两个方点加权2w,他自身也加权w,计算时用2w-w=w)
    区间加单点查树状数组维护
    这个做法的瓶颈在于查询单点时的复杂度时(O(Dlogn))(D为度数)
    对于查询点,他的父亲和重儿子暴力查询,其余轻儿子在加权时直接加到查询点上,由于树链剖分一条路径的轻儿子时(O(logn))级别的,可以接受

    #include<cstdio>
    #include<cstring>
    #include<cmath>
    #include<algorithm>
    #include<vector>
    #include<iostream>
    #include<map>
    #include<string>
    
    #define maxn 600005
    #define MOD 998244353
    
    using namespace std;
    
    inline long long getint()
    {
    	long long num=0,flag=1;char c;
    	while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    	while(c>='0'&&c<='9')num=num*10+c-48,c=getchar();
    	return num*flag;
    }
    
    int n,m,q;
    int fir[maxn],nxt[maxn<<1],to[maxn<<1],cnt;
    vector<int>G[maxn];
    int sz[maxn],tp[maxn],fa[maxn],dpt[maxn],son[maxn];
    int pos[maxn],cur;
    int dfn[maxn],low[maxn],S[maxn],tim,top,Tot;
    int t[maxn];
    
    inline int add(int x,int y){return x+y<MOD?x+y:x+y-MOD;}
    inline void newnode(int u,int v)
    {to[++cnt]=v,nxt[cnt]=fir[u],fir[u]=cnt;}
    
    inline void update(int x,int num)
    {for(int i=x;i<=Tot;i+=i&(-i))t[i]=add(t[i],num);}
    inline int getans(int x)
    {int num=0;for(int i=x;i;i-=i&(-i))num=add(num,t[i]);return num;}
    
    inline void add(int u,int v,int w)
    {
    	while(tp[u]!=tp[v])
    	{
    		if(dpt[tp[u]]<dpt[tp[v]])swap(u,v);
    		update(pos[tp[u]],w),update(pos[u]+1,MOD-w);
    		if(tp[u]>n)update(pos[fa[tp[u]]],MOD-w),update(pos[fa[tp[u]]]+1,w);
    		u=fa[tp[u]];
    	}
    	if(dpt[u]>dpt[v])swap(u,v);
    	update(pos[u],w),update(pos[v]+1,MOD-w);
    	if(u>n&&tp[u]==u&&fa[u])update(pos[fa[u]],MOD-w),update(pos[fa[u]]+1,w);
    }
    
    inline void Tarjan(int u)
    {
    	dfn[u]=low[u]=++tim;S[++top]=u;
    	for(int i=0;i<G[u].size();i++)
    	{
    		int v=G[u][i];
    		if(!dfn[v])
    		{
    			Tarjan(v);low[u]=min(low[u],low[v]);
    			if(low[v]>=dfn[u])
    			{
    				newnode(++Tot,u),newnode(u,Tot);int x;
    				do{x=S[top--],newnode(Tot,x),newnode(x,Tot);}while(v!=x);
    			}
    		}
    		else low[u]=min(low[u],dfn[v]);
    	}
    }
    
    inline void dfs1(int u)
    {
    	sz[u]=1;
    	for(int i=fir[u];i;i=nxt[i])if(to[i]!=fa[u])
    	{
    		fa[to[i]]=u,dpt[to[i]]=dpt[u]+1;
    		dfs1(to[i]),sz[u]+=sz[to[i]];
    		if(sz[to[i]]>sz[son[u]])son[u]=to[i];
    	}
    }
    
    inline void dfs2(int u,int ac)
    {
    	tp[u]=ac,pos[u]=++cur;
    	if(son[u])dfs2(son[u],ac);
    	for(int i=fir[u];i;i=nxt[i])if(to[i]!=fa[u]&&to[i]!=son[u])dfs2(to[i],to[i]);
    }
    
    inline void solve()
    {
    	for(int i=1;i<=m;i++)
    	{
    		int u=getint(),v=getint();
    		G[u].push_back(v),G[v].push_back(u);
    	}
    	Tarjan(1),dfs1(1),dfs2(1,1);
    	while(q--)
    	{
    		int op=getint();
    		if(op)
    		{
    			int u=getint();
    			int ans=MOD-getans(pos[u]);
    			if(son[u])ans=add(ans,getans(pos[son[u]]));
    			if(fa[u])ans=add(ans,getans(pos[fa[u]]));
    			printf("%d
    ",ans);
    		}
    		else
    		{
    			int u=getint(),v=getint(),w=getint();
    			if(u==v)update(pos[u],MOD-w),update(pos[u]+1,w);
    			else add(u,v,w),update(pos[u],MOD-w),update(pos[u]+1,w),update(pos[v],MOD-w),update(pos[v]+1,w);
    		}
    	}
    }
    
    int main()
    {
    	n=getint(),m=getint(),q=getint();Tot=n;
    	solve();
    }
    

  • 相关阅读:
    SPOJ 8222 NSUBSTR
    bzoj千题计划284:bzoj2882: 工艺
    bzoj千题计划283:bzoj4516: [Sdoi2016]生成魔咒(后缀数组)
    bzoj千题计划282:bzoj4517: [Sdoi2016]排列计数
    bzoj千题计划281:bzoj4558: [JLoi2016]方
    bzoj千题计划280:bzoj4592: [Shoi2015]脑洞治疗仪
    bzoj千题计划279:bzoj4591: [Shoi2015]超能粒子炮·改
    bzoj千题计划278:bzoj4590: [Shoi2015]自动刷题机
    bzoj千题计划277:bzoj4513: [Sdoi2016]储能表
    bzoj千题计划276:bzoj4515: [Sdoi2016]游戏
  • 原文地址:https://www.cnblogs.com/Darknesses/p/12931518.html
Copyright © 2011-2022 走看看