「SDOI2019」热闹的聚会与尴尬的聚会 ,LOJ3113,图论好题

本文探讨了如何通过寻找最小度数最大的连通子图和无向图的最大独立集来最大化pq值的问题。介绍了两种策略:一是删去当前最小度数的点以找到最大度数的连通子图;二是删去度数最小的点及其周围点,直到满足特定条件,以求得最大独立集。文章还提供了一段C++代码实现。

正题

      Portal

      看看题,发现我们要做的事情就是最大化pq。

      那么我们最大化p和q就好了呀!

      接着你就会发现,你要解决两个事情:

      1.找一个最小度数最大的连通子图

      2.找一个无向图的最大独立集。

      第一个貌似可做,第二个不是NP问题吗?

      对!

      首先说第一个怎么做。

      有一个贪心策略就是删去当前最小度数的点。

      因为删去的如果不是最小度数的点,那么这个最小度数不会增长。

      每次删去一个点和边就可以了,由于度数只会变小,所以这里完全不用可删除堆。

      最大独立集发现有一个限制,只要满足q<=\left \lfloor \frac{n}{p+1} \right \rfloor

      我们又有一个策略,每次删去度数最小的那个点和他周围的点,可以这些点的总和不超过p+1,否则当前就是一个最小度数更大的连通子图。

      那么直接这样取就可以了,这个也不用可删除堆,通过这题可以发现4倍常数可以决定一切。

      

#include<bits/stdc++.h>
using namespace std;

int T,n,m,ans1,ans2;
const int N=10010,M=100010;
struct edge{
	int y,next;
}s[M<<1];
struct node{
	int x,d;
	bool operator < (const node q)const{
		return d!=q.d?d>q.d:x>q.x;
	}
};
priority_queue<node> V;
int first[N],len=0;
vector<int> f;
bool tf[N];
int d[N],now[N];

#define nc() getchar()
inline int _read(){
    char ch=nc();int sum=0;
    while(!(ch>='0'&&ch<='9'))ch=nc();
    while(ch>='0'&&ch<='9')sum=sum*10+ch-48,ch=nc();
    return sum;
}

void ins(int x,int y){
	s[++len]=(edge){y,first[x]};first[x]=len;
	s[++len]=(edge){x,first[y]};first[y]=len;
	d[x]++;d[y]++;
}

void clear(priority_queue<node>&Q){
	static priority_queue<node> empty;
	Q=empty;
}

int main(){
//	freopen("party6.in","r",stdin);
//	freopen("party6.out","w",stdout);
	T=_read();
	while(T--){
		n=_read();m=_read();ans1=ans2=0;len=0;
		for(int i=1;i<=n;i++) first[i]=d[i]=tf[i]=0;
		int x,y;
		for(int i=1;i<=m;i++) x=_read(),y=_read(),ins(x,y);
		f.clear();clear(V);
		for(int i=1;i<=n;i++) now[i]=d[i],V.push((node){i,d[i]});
		int mmax=0,tot=-1,t,ans=0;
		while(!V.empty()){
			node X=V.top();V.pop();
			if(tf[X.x]) continue;
			f.push_back(X.x);
			if(mmax<=X.d) mmax=X.d,t=tot;
			tot++;
			tf[X.x]=true;
			for(int i=first[X.x];i!=0;i=s[i].next) if(!tf[s[i].y])
				now[s[i].y]--,V.push((node){s[i].y,now[s[i].y]});
		}
		memset(tf,false,sizeof(tf)); 
		for(int i=0;i<=t;i++) tf[f[i]]=true;
		printf("%d ",n-(t+1));for(int i=1;i<=n;i++) if(!tf[i]) printf("%d ",i);printf("\n");
		memset(tf,false,sizeof(tf));
		f.clear();clear(V);for(int i=1;i<=n;i++) now[i]=d[i],V.push((node){i,d[i]});
		while(!V.empty()){
			node X=V.top();V.pop();
			if(tf[X.x]) continue;
			f.push_back(X.x);tf[X.x]=true;
			for(int i=first[X.x];i!=0;i=s[i].next) if(!tf[s[i].y]){
				tf[s[i].y]=true;
				for(int j=first[s[i].y];j;j=s[j].next) if(!tf[s[j].y]){
					ans++;now[s[j].y]--;
					V.push((node){s[j].y,now[s[j].y]});
				}
			}
		}
		printf("%d ",f.size());for(int i=0;i<f.size();i++) printf("%d ",f[i]);printf("\n");
	}
}

 

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值