模拟退火60pt求调
查看原帖
模拟退火60pt求调
852827
Mr_yeh楼主2025/8/29 15:49
#include<bits/stdc++.h>
using namespace std;
const double dlt=0.995;
const int N=30;
int n,m,x,y,mp[N][N],a[N],as[N],ans=0x3f;
int calm(){
	int sum=0;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(mp[i][j]&&a[i]!=a[j]) sum++;
		}
	}
	return sum/2;
}
void SA(){
	double bt=3000,et=1e-15;
	while(bt>et){
		int sr1,sr2;
		while(1){
			sr1=rand()%n+1,sr2=rand()%n+1;
			if(a[sr1]!=a[sr2]) break;
		}
		swap(a[sr1],a[sr2]);
		int res=calm();
		int diff=res-ans;
		if(diff<0){
			ans=res;
			int tot=0;
			for(int i=1;i<=n;i++){
				if(a[i]) as[++tot]=i;
			}
		}else if(exp(-1.0*diff/bt)*RAND_MAX<rand()){
			swap(a[sr1],a[sr2]);
		}
		bt*=dlt;
	}
}
int main(){
	srand(time(0));
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	cin>>n>>m;
	for(int i=1;i<=n/2;i++) a[i]=1;
	while(m--){
		cin>>x>>y;
		mp[x][y]=mp[y][x]=1;
	}
	while((double)clock()/CLOCKS_PER_SEC<0.95){
		SA();
	}
	for(int i=1;i<=n/2;i++) cout<<as[i]<<" ";
	cout<<'\n';
	return 0;
}
2025/8/29 15:49
加载中...