如何優雅地生成仙人掌圖


用途

如果某個無向連通圖的任意一條邊至多只出現在一條簡單回路里,我們就稱這張圖為仙人掌圖。

所謂簡單回路就是指在圖上不重復經過任何一個頂點的回路。

在某些情況下,我們會需要生成仙人掌圖來檢驗代碼的正確性。

隨機連邊的話效率太低,而且生成的圖也可能不合法。

看上去似乎不大好實現,但實際上有一種比較優秀的做法:

對於原樹進行隨機鏈剖分,隨機幾個點向鏈頂連邊即可。

代碼

#include<bits/stdc++.h>
using namespace std;
#define rg register
const int maxn=2e5+5;
int h[maxn],tot=1,x[maxn],y[maxn];
struct asd{
	int to,nxt;
}b[maxn<<1];
void ad(rg int aa,rg int bb){
	b[tot].to=bb;
	b[tot].nxt=h[aa];
	h[aa]=tot++;
}
int tp[maxn];
bool vis[maxn];
void dfs(rg int now,rg int lat){
	rg int jud=0;
	for(rg int i=h[now];i!=-1;i=b[i].nxt){
		rg int u=b[i].to;
		if(u==lat) continue;
		if(!jud) tp[u]=tp[now],jud=1;
		else tp[u]=u;
		dfs(u,now);
	}
}
int main(){
	freopen("/dev/urandom","r",stdin);
	srand(getchar()*getchar()*getchar()*time(0));
	freopen("data.in","w",stdout);
	memset(h,-1,sizeof(h));
	int n=rand()%10+2,m=n-1;
	rg int aa,bb;
	tp[1]=1;
	for(rg int i=2;i<=n;i++){
		aa=rand()%(i-1)+1,bb=i;
		x[i-1]=aa,y[i-1]=bb;
		ad(aa,bb),ad(bb,aa);
	}
	dfs(1,0);
	rg int cnt=rand()%10+1;
	for(rg int i=1;i<=cnt;i++){
		aa=rand()%n+1;
		bb=tp[aa];
		if(aa!=bb && !vis[bb]){
			x[++m]=aa,y[m]=bb;
			vis[bb]=1;
		}
	}
	printf("%d %d\n",n,m);
	for(rg int i=1;i<=m;i++){
		printf("%d %d\n",x[i],y[i]);
	}
	return 0;
}


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM