【SDOI2018】战略游戏(圆方树)(虚树)(树上差分)

传送门


题解:

圆方树的性质不再赘述,两个圆点路径上的圆点就是原图中两点之间的必经点。

那么我们要统计的就是询问的点构成的虚树上有多少个圆点。

并不需要把虚树建出来,直接利用dfs序树上差分即可。

跑得还挺快的。


代码:

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

namespace IO{
	inline char gc(){
		static cs int Rlen=1<<22|1;static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	template<typename T>T get(){
		char c;T num;while(!isdigit(c=gc()));num=c^48;
		while(isdigit(c=gc()))num=((num+(num<<2))<<1)+(c^48);
		return num;
	}inline int gi(){return get<int>();}
}
using namespace IO;

using std::cerr;
using std::cout;

cs int N=1e5+7,M=2e5+7;

int n,m,ext;

namespace T{//Tree
	cs int N=2e5+7;
	int el[N],nxt[N];
	void adde(int u,int v){
		nxt[v]=el[u];el[u]=v;
	}
	int dfn[N],dfc;
	int siz[N],son[N],top[N];
	int fa[N],d[N],dis[N];
	void dfs1(int u,int p){
		fa[u]=p,d[u]=d[p]+1,dis[u]=dis[p]+(u<=n);
		for(int re v=el[u];v;v=nxt[v]){
			dfs1(v,u);siz[u]+=siz[v];
			if(siz[v]>siz[son[u]])son[u]=v;
		}++siz[u];
	}
	void dfs2(int u,int tp){
		top[u]=tp;dfn[u]=++dfc;
		if(son[u])dfs2(son[u],tp);
		for(int re v=el[u];v;v=nxt[v])
		if(v!=son[u])dfs2(v,v);
	}
	int LCA(int u,int v){
		while(top[u]!=top[v])d[top[u]]<d[top[v]]?v=fa[top[v]]:u=fa[top[u]];
		return d[u]<d[v]?u:v;
	}
	void build(){dfc=0;dfs1(1,0);dfs2(1,1);}

	void solve(){
		static int nd[N];int ans=0,s=gi();
		for(int re i=1;i<=s;++i)nd[i]=gi();
		std::sort(nd+1,nd+s+1,[](int u,int v){return dfn[u]<dfn[v];});
		for(int re i=1;i<=s;++i)ans+=dis[nd[i]];
		for(int re i=1;i<s;++i)ans-=dis[LCA(nd[i],nd[i+1])];
		cout<<ans-dis[fa[LCA(nd[1],nd[s])]]-s<<"\n";
	} 
	void init(){
		memset(siz,0,(n+1)<<3);
		memset(son,0,(n+1)<<3);
		memset(el,0,(n+1)<<3);
		memset(dfn,0,(n+1)<<3);
	}
}

namespace G{//Graph
	int el[N],nxt[M<<1],to[M<<1],ec;
	void adde(int u,int v){
		nxt[++ec]=el[u],el[u]=ec,to[ec]=v;
		nxt[++ec]=el[v],el[v]=ec,to[ec]=u;
	}
	int st[N],tp;
	int dfn[N],low[N],dfc;
	void tarjan(int u){
		dfn[u]=low[u]=++dfc;st[++tp]=u;
		for(int re e=el[u],v;e;e=nxt[e])
		if(!dfn[v=to[e]]){
			tarjan(v);low[u]=std::min(low[u],low[v]);
			if(low[v]>=dfn[u]){
				T::adde(u,++ext);
				do{T::adde(ext,st[tp--]);}while(st[tp+1]!=v);
			}
		}else low[u]=std::min(low[u],dfn[v]);
	}
	void init(){
		memset(el,0,(n+1)<<2);
		memset(dfn,0,(n+1)<<2);
		memset(low,0,(n+1)<<2);
		ec=dfc=tp=0;
	}
}

void work(){
	n=gi(),m=gi();G::init(),T::init();
	for(int re i=1;i<=m;++i)G::adde(gi(),gi());
	ext=n;G::tarjan(1);T::build();
	int Q=gi();while(Q--)T::solve();
}

void file(){
#ifdef zxyoi
	freopen("game.in","r",stdin);
#endif
}
signed main(){file();int T=gi();while(T--)work();return 0;}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值