[BZOJ2428][HAOI2006]均分数据（模拟退火）

=== ===

这里放传送门

=== ===

题解

感觉这题才是正儿八经的模拟退火啊。。。先给每个元素随机一个位置，然后对于每个元素随机调整它所在的组，调整完了以后计算当前方差，如果当前方差更优那么就直接更改，否则就有一定的概率接受更改，并且这个概率是随着时间的推移逐渐变小的。

需要注意的问题是首先它这个东西似乎只做一遍无法得到最优解，所以要做好多遍。。还有一个问题就是如果直接无脑退火会有点GG，因为温度比较高的时候不确定性比较大，接受较差解的概率比较大，那如果无脑随机的话就可能做很多不必要的移动，浪费很多机会。所以当温度比较大的时候直接把当前这个元素扔到和最小的那一组就可以了，当温度较低的时候就可以随机一个位置来进行微调。

代码

请自行忽略奇怪的随机种子。。ATP当时本来想随便按一个日期的但是手一滑就按成了这样子结果这个数还莫名好用。。

#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define random(x)(rand()%(x)+1)
using namespace std;
const double eps=1e-3;
int n,m,a[50],pos[50];
double Min,sum[50],aver,now,temp;
double calc(double x){return x*x;}
int getMin(){
    int pos=1;
    for (int i=2;i<=m;i++)
      if (sum[i]<sum[pos]) pos=i;
    return pos;
}
void work(){
    int lp,np,loc;
    double tmp;
    aver=now=0;
    memset(sum,0,sizeof(sum));
    for (int i=1;i<=n;i++) pos[i]=random(m);
    for (int i=1;i<=n;i++){
        sum[pos[i]]+=a[i];
        aver+=a[i];
    }
    aver/=m;temp=10000;
    for (int i=1;i<=m;i++) now+=calc(sum[i]-aver);
    Min=min(Min,now);
    while (temp>eps){
        loc=random(n);lp=pos[loc];
        if (temp>50) np=getMin();
        else np=random(m);
        if (lp==np) continue;
        tmp=now;
        tmp-=calc(sum[lp]-aver);
        tmp-=calc(sum[np]-aver);
        tmp+=calc(sum[lp]-a[loc]-aver);
        tmp+=calc(sum[np]+a[loc]-aver);
        if (tmp<=now||temp>random(10000)){
            pos[loc]=np;now=tmp;
            sum[lp]-=a[loc];sum[np]+=a[loc];
        }
        Min=min(Min,now);
        temp*=0.9;
    }
    Min=min(Min,now);
}
int main()
{
    scanf("%d%d",&n,&m);
    Min=0x7fffffff;srand(20020987);
    for (int i=1;i<=n;i++) scanf("%d",&a[i]);
    for (int i=1;i<=12000;i++) work();
    printf("%.2lf\n",sqrt(Min/(double)m));
    return 0;
}