[HAOI2006] 均分数据

题目描述：

N个元素分M组，总方差最小

题目分析：

不会啊
虽说大佬们都讲是个膜你退火
但我感觉是个带贪心的随机乱搞！
数据只有 20 我们每次就去搞随机分组。
当然 在初期温度很高的时候 我们需要找一组方差最小的去换
加入退火的可能去找更差解的随机，搞个10000-20000 基本上就可以A了

题目链接：

BZOJ 2428
Luogu 2503

Ac 代码：

#include <ctime>
#include <iostream>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cmath>
const double eps=1e-3;
const double r=0.9;
double minans=1e30,p=0;
const int maxm=200;
int n,m,a[maxm];
int sum[maxm];
int pos[maxm];
void SA(double T)
{
    memset(sum,0,sizeof(sum));
    for(int i=1;i<=n;i++)
     {
        pos[i]=rand()%m+1;
        sum[pos[i]]+=a[i];
     }
    double ans=0;
    for(int i=1;i<=m;i++)
     ans+=(double)(sum[i]-p)*(double)(sum[i]-p);
    while(T>eps)
    {
        int t=rand()%n+1,x=pos[t],y;
        if(T>500)
         y=std::min_element(sum+1,sum+m+1)-sum;
        else y=rand()%m+1;
        if(x==y) continue;
        double tmp=ans;
        ans-=(sum[x]-p)*(sum[x]-p);
        ans-=(sum[y]-p)*(sum[y]-p);
        sum[x]-=a[t],sum[y]+=a[t];
        ans+=(sum[x]-p)*(sum[x]-p);
        ans+=(sum[y]-p)*(sum[y]-p);
        if(ans<tmp||rand()%10000<=T) pos[t]=y;
        else sum[x]+=a[t],sum[y]-=a[t],ans=tmp;
        minans=std::min(minans,ans);
        T*=r;
    }
}
int main()
{
    srand(19990720);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]),p+=a[i];
    p/=(double)m;
    for(int i=1;i<=10000;i++) SA(10000);
    printf("%.2lf\n",std::sqrt(minans/m));
    return 0;
}