绍兴一中模拟赛3.19——白驹过隙

博客围绕类仙人掌展开,定义类仙人掌后,目标是求解类仙人掌上的最大独立集大小,其中n<=50000,m<=100000,还给出代码相关注意事项,如f[p2][p1][0][1]等于f[p1][p2][1][0] 。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Description

定义类仙人掌为:

求类仙人掌上的最大独立集大小,(n&lt;=50000,m&lt;=100000)(n&lt;=50000,m&lt;=100000)(n<=50000,m<=100000)

Solution


Code

注意f[p2][p1][0][1]f[p2][p1][0][1]f[p2][p1][0][1]等于f[p1][p2][1][0]f[p1][p2][1][0]f[p1][p2][1][0],而不是f[p1][p2][0][1]f[p1][p2][0][1]f[p1][p2][0][1]

#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(register int i=(a);i<=(b);++i)
#define dep(i,a,b) for(register int i=(b);i>=(a);--i)
inline char gc(){
	static char buf[100000],*p1=buf,*p2=buf;
	return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int rd(){
	int x=0,fl=1;char ch=gc();
	for(;ch<48||ch>57;ch=gc())if(ch=='-')fl=-1;
	for(;48<=ch&&ch<=57;ch=gc())x=(x<<3)+(x<<1)+(ch^48);
	return x*fl;
}
const int N=100002;
int n,m,x,y,p1,p2,ans;
queue<int>q;
map<int,bool>p[N];
map<int,int>f[2][2][N];
set<int>to[N];
bool vis[N];
int main(){
	n=rd(),m=rd();
	rep(i,1,m){
		x=rd(),y=rd();
		to[x].insert(y),to[y].insert(x);
		p[x][y]=p[y][x]=1;
		f[0][1][x][y]=f[1][0][x][y]=f[0][1][y][x]=f[1][0][y][x]=1;
		f[1][1][x][y]=f[1][1][y][x]=-1e9;
	}
	rep(i,1,n)
		if (to[i].size()==2) q.push(i);
	while (!q.empty()){
		x=q.front(),q.pop();
		if (vis[x] || to[x].size()!=2) continue;
		p1=*to[x].begin(),p2=*(++to[x].begin());
		f[0][0][p1][p2]+=max(f[0][0][x][p1]+f[0][0][x][p2],f[1][0][x][p1]+f[1][0][x][p2]-1);
		f[0][1][p1][p2]+=max(f[0][0][x][p1]+f[0][1][x][p2],f[1][0][x][p1]+f[1][1][x][p2]-1)-(f[0][1][p1][p2]>0);
		f[1][0][p1][p2]+=max(f[0][1][x][p1]+f[0][0][x][p2],f[1][1][x][p1]+f[1][0][x][p2]-1)-(f[1][0][p1][p2]>0);
		if (!p[p1][p2])
		ans=max(ans,f[1][1][p1][p2]+=max(f[0][1][x][p1]+f[0][1][x][p2],f[1][1][x][p1]+f[1][1][x][p2]-1)-(f[1][1][p1][p2]>0)*2);
		ans=max(ans,max(f[0][0][p1][p2],max(f[0][1][p1][p2],f[1][0][p1][p2])));
		f[0][0][p2][p1]=f[0][0][p1][p2],f[1][0][p2][p1]=f[0][1][p1][p2],f[0][1][p2][p1]=f[1][0][p1][p2],f[1][1][p2][p1]=f[1][1][p1][p2];
		vis[x]=1;
		to[p1].erase(x),to[p2].erase(x);
		to[p1].insert(p2),to[p2].insert(p1);
		if (to[p1].size()==2) q.push(p1);
		if (to[p2].size()==2) q.push(p2);
	}
	printf("%d",ans);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值