bzoj2428 [HAOI2006]均分数据

原创于 2017-04-18 14:01:26 发布 · 400 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

本文通过一道具体题目介绍了模拟退火算法的基本原理及其应用。详细解释了如何通过模拟物质从高温到低温的状态变化来求解优化问题，并给出了完整的C++实现代码。

一开始想了一个用优先队列的贪心，只能在洛谷过三个点（我还是too young too simple）。。。所以在网上看了一下

模拟退火

模拟退火之前在网上也看到过这个算法吧，但是一直都没有学，觉得难就没学，其实仔细看一下还是比较简单的

这题就是模拟退火的一个裸题~~，先上代码，详细的后面再补~~

模拟退火就是模拟物理上物体从高温到低温的一个过程，具体的网上也很多，兔子的那个比喻就很恰当

高温的时候物体内部的粒子比较混乱，所以它所在的地方可能性很多（大概是这个意思吧。。），那么在程序里面高温的时候状态转移过去的概率就比较大

反之，低温的时候物体内部粒子倾向于稳定，那么这个时候转移的概率比较小，除非当前状态比之前更优，那么直接模拟了

//Decision's template
#include<cstdio>
#include<cstring>
#include<iostream>
#include<cstdlib>
#include<vector>
#include<queue>
#include<stack>
#include<algorithm>
#include<string>
#include<cmath>
#include<map>
#include<set>
using namespace std;

#define DP_maxn 16
#define maxn 100000
#define INF 10000007
#define MOD 1000000007
#define mst(s,k) memset(s,k,sizeof(s))

typedef long long ll;

struct Edge{
    int from,to,dist;
    Edge(int u,int v,int d):from(u),to(v),dist(d){}
};

/*-------------------------------template End--------------------------------*/

int a[30],n,k,pos[30];
double minans = (double)INF,ans,ave = 0.0,sum[30];

double sqr(double a) {return a*a;}

void  Fire()
{
	mst(sum,0);
	double r= 0.9,T = 10000.00;
	ans = 0.0;
	for(int i = 1;i<=n;i++)
	{
		pos[i] = rand()%k+1;
	}
	for(int i = 1;i<=n;i++) sum[pos[i]] += a[i];
	for(int i = 1;i<=k;i++)
	{
		ans+=(sum[i]-ave)*(sum[i]-ave);
	}
	//cout<<ans<<endl;
	while(T>0.1)
	{
		T*=r;
		int t = rand()%n+1;
		int x = pos[t],y;
		if(T>500) y = min_element(sum+1,sum+k+1)-sum;
		else y = rand()%k+1; 
		if(x==y) continue;
		double tmp = ans;
		ans -= (sum[x]-ave)*(sum[x]-ave);
		ans -= (sum[y]-ave)*(sum[y]-ave);
		sum[x]-=a[t];sum[y]+=a[t];
		ans += (sum[x]-ave)*(sum[x]-ave);
		ans += (sum[y]-ave)*(sum[y]-ave);
		if(ans<=tmp||rand()%10000<=T){
			pos[t] = y;
		} 
		else{
			ans = tmp;
			sum[x]+=a[t];
			sum[y]-=a[t];
		}
	}
	minans = min(ans,minans);
	//cout<<minans<<endl;
}

int main()
{
	//freopen("std.in","r",stdin);
	//freopen("std.out","w",stdout);
	mst(sum,0);
	cin>>n>>k;
	for(int i = 1;i<=n;i++)
	{
		cin>>a[i];
		ave+=a[i];
		swap(a[i],a[rand()%i+1]);
	}
	ave /= k;
	for(int i = 1;i<=10000;i++) Fire();
	printf("%.2f\n",sqrt(minans/k));
	return 0;
}