zoukankan      html  css  js  c++  java
  • 【省选模拟】—Cactus(圆方树+dfs序)

    传送门

    割点也建一个方点要好做一点
    建出圆方树后就差不多了吧。。。
    处理一下lcalca如果是方点会对上面那个圆点有贡献
    开个桶记一下就完了

    #include<bits/stdc++.h>
    using namespace std;
    const int RLEN=1<<18|1;
    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;
    }
    const int N=1e6+5;
    const int mod=998244353;
    inline int add(int a,int b){
    	return a+b>=mod?a+b-mod:a+b;
    }
    inline int dec(int a,int b){
    	return a<b?a-b+mod:a-b;
    }
    inline void ad(int &a,int b){
    	a=a+b,a<mod?0:a-=mod;
    }
    struct Graph{
    	int adj[N],nxt[N<<1],to[N<<1],cnt;
    	inline void addedge(int u,int v){
    		nxt[++cnt]=adj[u],adj[u]=cnt,to[cnt]=v;
    		nxt[++cnt]=adj[v],adj[v]=cnt,to[cnt]=u;
    	}
    }T,G;
    int low[N],dfn[N],tot,bel,ori[N];
    stack<int> stk;
    void tarjan(int u){
    	low[u]=dfn[u]=++tot;stk.push(u);
    	for(int e=G.adj[u];e;e=G.nxt[e]){
    		int v=G.to[e];
    		if(!dfn[v]){
    			tarjan(v),low[u]=min(low[u],low[v]);
    			if(low[v]>=dfn[u]){
    				int tmp;bel++;
    				T.addedge(bel,u);
    				do{
    					tmp=stk.top(),stk.pop();
    					T.addedge(bel,tmp);
    				}while(tmp!=v);
    			}
    		}
    		else low[u]=min(low[u],dfn[v]);
    	}
    }
    int n,m,q;
    namespace SLPF{
    	int pos[N],out[N],fa[N],top[N],dep[N],siz[N],son[N],tot,val[N];
    	void dfs1(int u){
    		siz[u]=1;
    		for(int e=T.adj[u];e;e=T.nxt[e]){
    			int v=T.to[e];
    			if(v==fa[u])continue;
    			fa[v]=u,dep[v]=dep[u]+1;
    			dfs1(v),siz[u]+=siz[v];
    			if(siz[v]>=siz[son[u]])son[u]=v;
    		}
    	}
    	void dfs2(int u,int tp){
    		top[u]=tp;if(u>n)pos[u]=++tot;else pos[u]=pos[fa[u]];
    		if(son[u])dfs2(son[u],tp);
    		for(int e=T.adj[u];e;e=T.nxt[e]){
    			int v=T.to[e];
    			if(v==fa[u]||v==son[u])continue;
    			dfs2(v,v);
    		}
    		out[u]=tot;
    	}
    	inline int Lca(int u,int v){
    		while(top[u]!=top[v]){
    			if(dep[top[u]]<dep[top[v]])swap(u,v);
    			u=fa[top[u]];
    		}
    		return dep[u]>dep[v]?v:u;
    	}
    	int tr[N];
    	#define lowbit(x) (x&(-x))
    	inline void update(int p,int k){
    		for(;p&&p<=tot;p+=lowbit(p))ad(tr[p],k);
    	}
    	inline int qry(int p,int res=0){
    		for(;p;p-=lowbit(p))ad(res,tr[p]);return res;
    	}
    	inline int query(int l,int r){
    		return dec(qry(r),qry(l-1));
    	}
    	#undef lowbit
    	inline void update(int u,int v,int k){
    		int lca=Lca(u,v);
    		update(pos[u],k),update(pos[v],k);
    		update(pos[lca],mod-k);if(fa[lca])update(pos[fa[lca]],mod-k);
    		if(lca<=n)ad(val[lca],k);
    		else ad(val[fa[lca]],k);
    	}
    	inline int query(int u){
    		return add(val[u],query(pos[fa[u]],out[fa[u]]));
    	}
    }
    signed main(){
    	bel=n=read(),m=read(),q=read();
    	for(int i=1;i<=m;i++){
    		int u=read(),v=read();
    		G.addedge(u,v);
    	}
    	tarjan(1);
    	SLPF::dfs1(bel),SLPF::dfs2(bel,bel);
    	while(q--){
    		int op=read();
    		if(op==0){
    			int u=read(),v=read(),k=read();
    			SLPF::update(u,v,k);
    		}
    		else {
    			int u=read();
    			cout<<SLPF::query(u)<<'
    ';
    		}
    	}
    }
    
  • 相关阅读:
    梅森数 洛谷p1045
    高精度算法
    并查集(算法)
    各种背包模板
    爬虫作业
    numpy、scipy、pandas、matplotlib的读书报告:
    体育竞技分析羽毛球
    聊斋相关的分词,出现次数最高的20个
    在终端输出如下信息,(两种以上方法)
    7段数码管绘制
  • 原文地址:https://www.cnblogs.com/stargazer-cyk/p/11145578.html
Copyright © 2011-2022 走看看