bzoj 2428 均分数据模拟退火

最新推荐文章于 2018-05-02 19:26:00 发布

原创最新推荐文章于 2018-05-02 19:26:00 发布 · 409 阅读

0 ·

CC 4.0 BY-SA版权

bzoj 同时被 2 个专栏收录

61 篇文章

订阅专栏

模拟退火

3 篇文章

订阅专栏

本文介绍了一种使用模拟退火算法解决特定问题的实现方法。通过调整关键参数如初温和降温策略，结合随机数生成来优化解的质量。该算法在外层加入循环进一步提升了解的准确性。

模拟退火

按照自己的思路打了，结果WA，发现退火最关键的就是初温，降温，和修改次数，

这个题还在外层带了一个循环，骚气

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#define inf 0x3f3f3f3f3f3f3f3fll
using namespace std;
double av,ans,minn=inf,sum[30];
int n,m,a[30],be[30];
double Rand(){return (rand()%1000)/1000.0;}
void work(){
	memset(sum,0,sizeof sum);
	ans=0;
	for(int i=1;i<=n;i++)
		be[i]=rand()%m+1,sum[be[i]]+=a[i];
	for(int i=1;i<=m;i++)
		ans+=(sum[i]-av)*(sum[i]-av);
	double t=10000;
	while(t>0.1){
		int x=rand()%n+1,y=be[x],z;
		if(t>500)z=min_element(sum+1,sum+m+1)-sum;
		else z=rand()%m+1;
		double nxt=ans;
		nxt-=(sum[y]-av)*(sum[y]-av);
		nxt-=(sum[z]-av)*(sum[z]-av);
		sum[y]-=a[x]; sum[z]+=a[x];
		nxt+=(sum[y]-av)*(sum[y]-av);
		nxt+=(sum[z]-av)*(sum[z]-av);
		double dE=ans-nxt;
		if(dE>0||exp(dE/t)>Rand())
			ans=nxt,be[x]=z;
		else
			sum[y]+=a[x],sum[z]-=a[x];
		t*=0.9;
	}
	if(ans<minn) minn=ans;
}
int main(){
	srand(20001101);
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
		av+=a[i];
	}av/=m;
	for(int i=1;i<=10000;i++)work();
	printf("%0.2lf\n",sqrt(minn/m));
	return 0;
}