题目描述:
N个元素分M组,总方差最小
题目分析:
不会啊
虽说大佬们都讲是个膜你退火
但我感觉是个带贪心的随机乱搞!
数据只有 20 我们每次就去搞随机分组。
当然 在初期温度很高的时候 我们需要找一组方差最小的去换
加入退火的可能去找更差解的随机,搞个10000-20000 基本上就可以A了
题目链接:
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;
}