Gym - 102759E Chemistry——双指针、LCT

E. Chemistry

题目描述

给定一个 $n$ 个点 $m$ 条边的无向图，求有多少个区间 $[l,r]$ 满足标号在 $[l,r]$ 范围内的点的导出子图是一条链。

$n,m\le 250000$。

题解

我们不妨把“链”这个条件拆分一下：

1. 任何一个点的度数≤2；
2. 图中没有环；
3. 点数-边数=1。

容易证明，满足这三个条件的图一定是连通的一条链。前两个条件的优点是：满足单调性，即加入的点和边越多越不容易满足。

当我们固定了区间右端点 $r$ 时，在前两个条件的约束下，合法的 $l$ 一定 $\ge$ 某个边界 $l_{min}$，并且由单调性还可得出 $l_{min}$ 随 $r$ 的增大单调不降，所以我们可以想办法快速判断前两个条件，然后用双指针维护这个 $l_{min}$。

第一个条件用数组简单维护一下即可。第二个条件需要支持动态加边删边和判断连通性，用LCT维护即可。

这个时候我们发现，由于图中已经不存在环，所以点数-边数必定 $\ge 1$，所以我们用一个线段树维护以每个位置为左端点时点数-边数的最小值即可，可以顺便求出最小值的数量。

总复杂度是 $O(n\log n)$，所以不明白为什么要开5秒。

代码

#include<bits/stdc++.h>//JZM yyds!!
#define ll long long
#define uns unsigned
#define IF (it->first)
#define IS (it->second)
#define END putchar('\n')
using namespace std;
const int MAXN=250005;
const ll INF=1e17;
inline ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+(s^48),s=getchar();
	return f?x:-x;
}
int ptf[50],lpt;
inline void print(ll x,char c='\n'){
	if(x<0)putchar('-'),x=-x;
	ptf[lpt=1]=x%10;
	while(x>9)x/=10,ptf[++lpt]=x%10;
	while(lpt)putchar(ptf[lpt--]^48);
	if(c>0)putchar(c);
}
inline ll lowbit(ll x){return x&-x;}

int n,m;

struct spl{
	int fa,h[2],a;bool lz;spl(){}
	spl(int A){fa=h[0]=h[1]=0,a=A,lz=0;}
}t[MAXN];
inline void cover(int x){
	if(!x)return;
	t[x].lz^=1,swap(t[x].h[0],t[x].h[1]);
}
inline void pushd(int x){
	if(t[x].lz)cover(t[x].h[0]),cover(t[x].h[1]),t[x].lz=0;
}
inline void update(int x){
	//do nothing
}
inline bool sd(int x){return x==t[t[x].fa].h[1];}
inline bool isroot(int x){return x!=t[t[x].fa].h[sd(x)];}
inline void lin(int x,int y,bool f){if(x)t[x].h[f]=y;if(y)t[y].fa=x;}
inline void rott(int x){
	if(!t[x].fa||isroot(x))return;
	bool f1=sd(x),f2=sd(t[x].fa);
	int fa=t[x].fa,ff=t[fa].fa,sn=t[x].h[f1^1];
	if(isroot(fa))t[x].fa=ff;else lin(ff,x,f2);
	lin(fa,sn,f1),lin(x,fa,f1^1);
	update(fa),update(x),update(ff);
}
inline void pushtag(int x){
	if(!isroot(x))pushtag(t[x].fa);
	pushd(x);
}
inline void splay(int x){
	pushtag(x);
	while(!isroot(x)){
		if(!isroot(t[x].fa)){
			if(sd(x)==sd(t[x].fa))rott(t[x].fa);
			else rott(x);
		}rott(x);
	}
}
inline void access(int x){
	for(int y=0;x;y=x,x=t[x].fa)splay(x),t[x].h[1]=y,update(x);
}
inline void makeroot(int x){access(x),splay(x),cover(x);}
inline void LINK(int x,int y){
	if(x==y)return;
	makeroot(x),access(y),splay(y);
	if(t[x].fa==0)lin(y,x,1),update(y);
}
inline void CUT(int x,int y){
	makeroot(x),access(x),splay(y);
	if(!t[y].h[0]&&t[y].fa==x)t[y].fa=0;
}
inline bool findedge(int x,int y){
	makeroot(x),access(y),splay(y);
	return (x^y)&&t[x].fa==0;
}

int du[MAXN],n3;
struct itn{
	int a,b;itn(){}
	itn(int A,int B){a=A,b=B;}
	itn operator+(const int&y)const{return itn(a+y,b);}
	void operator+=(const int &y){a+=y;}
};
inline itn MIN(itn a,itn b){
	if(!b.b)return a;
	if(!a.b)return b;
	if(a.a^b.a)return a.a<b.a?a:b;
	else return itn(a.a,a.b+b.b);
}
struct zkw{
	itn f[MAXN<<2];
	int lz[MAXN<<2],p;
	inline void init(int n){
		for(p=1;p<n+2;p<<=1);
		for(int i=p+n+2;i>0;i--)f[i]=itn(0,0),lz[i]=0;
		for(int i=1;i<=n;i++)f[p+i]=itn(0,1);
		for(int i=(p+n)>>1;i;i--)f[i]=MIN(f[i<<1],f[i<<1|1]);
	}
	inline void add(int l,int r,int d){
		if(l>r)return;
		for(l=p+l-1,r=p+r+1;l^1^r;){
			if(~l&1)f[l^1]+=d,lz[l^1]+=d;
			if(r&1)f[r^1]+=d,lz[r^1]+=d;
			l>>=1,r>>=1;
			f[l]=MIN(f[l<<1],f[l<<1|1])+lz[l];
			f[r]=MIN(f[r<<1],f[r<<1|1])+lz[r];
		}
		for(l>>=1;l;l>>=1)
			f[l]=MIN(f[l<<1],f[l<<1|1])+lz[l];
	}
	inline itn query(int l,int r){
		if(l>r)return itn(0,0);
		itn sl=itn(0,0),sr=itn(0,0);
		for(l=p+l-1,r=p+r+1;l^1^r;){
			if(~l&1)sl=MIN(sl,f[l^1]);
			if(r&1)sr=MIN(sr,f[r^1]);
			l>>=1,r>>=1,sl+=lz[l],sr+=lz[r];
		}sl=MIN(sl,sr);
		for(l>>=1;l;l>>=1)sl+=lz[l];
		return sl;
	}
}T;
struct edge{
	int v,to;edge(){}
	edge(int V,int T){v=V,to=T;}
}e[MAXN<<1];
int EN,G[MAXN];
inline void addedge(int u,int v){
	e[++EN]=edge(v,G[u]),G[u]=EN;
	e[++EN]=edge(u,G[v]),G[v]=EN;
}
ll ans;
inline void delp(int x,int r){
	for(int i=G[x];i;i=e[i].to){
		int v=e[i].v;
		if(v<x||v>r)continue;
		du[x]--,du[v]--;
		if(du[x]==2)n3--;
		if(du[v]==2)n3--;
		CUT(x,v);
	}
}
signed main()
{
	n=read(),m=read();
	T.init(n);
	for(int i=1;i<=m;i++)addedge(read(),read());
	for(int x=1,l=1;x<=n;x++){
		T.add(1,x,1);
		for(int i=G[x];i;i=e[i].to){
			int v=e[i].v;
			if(v>=x)continue;
			T.add(1,v,-1);
			if(v>=l){
				du[x]++,du[v]++;
				if(du[x]==3)n3++;
				if(du[v]==3)n3++;
			}
		}
		for(int i=G[x];i;i=e[i].to){
			int v=e[i].v;
			if(v<x&&v>=l){
				while(l<x&&!findedge(x,v))delp(l,x),l++;
				if(v>=l)LINK(x,v);
			}
		}
		while(l<x&&n3>0)delp(l,x),l++;
		itn ad=T.query(l,x);
		if(ad.a==1)ans+=ad.b;
	}
	print(ans);
	return 0;
}