模拟退火。
一种十分玄学的随机算法,网上可以查到比较详细的资料。
先随机地把数分成m组,每次随机地选择一个数,一开始直接选最小的一组,后来就随机一组,把这个数换到该组看看答案能不能变小,如果变小则换,如果没有变小,按模拟退火的玄学方式判断一下,也要交换。
srand(time(0))在bzoj会RE,不知为何。
非常玄学的正确性,没拍多少数据就会出现差别,然而OJ上可以过的。
//Twenty
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
#include<ctime>
#define Ti 1e4
const int maxn=;
using namespace std;
int n,from[maxn],m;
double tp,ave,ans,p[maxn],sum[maxn];
double pf(double x) {return x*x;}
double solve() {
double res=;
for(int i=;i<=m;i++) sum[i]=;
for(int i=;i<=n;i++) {
from[i]=rand()%m+;
sum[from[i]]+=p[i];
}
for(int i=;i<=m;i++)
res+=pf(sum[i]-ave);
double T=Ti;
while(T>0.1) {
int x=rand()%n+,y;
y=T>?min_element(sum+,sum+m+)-sum:rand()%m+;
tp=res-pf(sum[y]-ave)+pf(sum[y]+p[x]-ave)-pf(sum[from[x]]-ave)+pf(sum[from[x]]-p[x]-ave);
if(tp<res||exp((res-tp)*Ti/T)>(double)rand()/RAND_MAX) {
res=tp;
sum[y]+=p[x];
sum[from[x]]-=p[x];
from[x]=y;
}
T*=0.9;
}
return res;
} int main()
{
srand();
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++) {
scanf("%lf",&p[i]);
ave+=p[i];
}
ave/=m;
ans=solve();
for(int i=;i<=Ti;i++)
ans=min(ans,solve());
printf("%.2lf\n",sqrt(ans/m));
return ;
}
其实不懂能不能随便放别人博客,侵删