模拟退火的

题目链接

体验乱调参数而看天意的奇特体验

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll,ll> pii;
const int inf=0x3f3f3f3f;
const int N=1e5+10;
const int mod=1e9+7;
//#define int long long
const ll INF=2e9+10;

mt19937_64 rd(23333);
uniform_real_distribution<double> drd(0.000001,0.99999);

int n,m,c,res;
int a[300];
int b[30][30];
int nx[4]={0,1,0,-1},ny[4]={1,0,-1,0};
int ans[30][30],q;

void init(){
	int pos=1;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			while(a[pos]==0)
				pos++;
			a[pos]--;
			b[i][j]=pos;
			ans[i][j]=b[i][j];
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			for(int k=0;k<4;k++){
				int xx=i+nx[k];
				int yy=j+ny[k];
				if(xx<1||xx>n||yy<1||yy>m)
					continue;
				if(b[xx][yy]!=b[i][j])
					res++;
			}
		}
	}
	res/=2;
	q=res;
}

int cal(int x,int y){
	int re=0;
	for(int k=0;k<4;k++){
		int xx=x+nx[k];
		int yy=y+ny[k];
		if(xx<1||xx>n||yy<1||yy>m)
			continue;
		if(b[xx][yy]!=b[x][y])
			re++;
	}
	return re;
}

void sa(){
	for(int t=1e15;t>=1e-15;t*=0.9999){
		int x1=(rd()%n)+1;
		int y1=(rd()%m)+1;
		int x2=(rd()%n)+1;
		int y2=(rd()%m)+1;
		if(b[x1][y1]==b[x2][y2])
			continue;
		int pre=cal(x1,y1)+cal(x2,y2);
		swap(b[x1][y1],b[x2][y2]);
		int nex=cal(x1,y1)+cal(x2,y2);
		int temp=res+nex-pre;
		int sub=res-temp;
		if(sub>0){
			res=temp;
			if(res<q){
				for(int i=1;i<=n;i++)
				for(int j=1;j<=m;j++)
					ans[i][j]=b[i][j];
				q=res;
			}
		}
		else if(exp(sub/t)<1.0*drd(rd))
			res=temp;
		else
			swap(b[x1][y1],b[x2][y2]);
	}
}

void solve(){
	cin>>n>>m>>c;
	for(int i=1;i<=c;i++)
		cin>>a[i];
	init();
	while(clock()/(1.0*CLOCKS_PER_SEC)<=4.89) sa();
	//cout<<q<<endl;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			if(j>1) cout<<' ';
			cout<<ans[i][j];
		}
		cout<<endl;
	}
}
signed main(){
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	int t=1;
	//cin>>t;
	while(t--){
		solve();
	}
	return 0;
}