SA 60pts~80pts求调

#include <bits/stdc++.h>
using namespace std;
int n,m,a[25],sum[10],q[25],v[10];
double ans;
const double bt=1000,et=1e-15,lt=0.9996,tim=0.99;
double work()
{
	double pjs,mans;
	for (int i=1; i<=m;i++)
	{
		pjs+=v[i];
	}
	pjs=pjs/double(m);
	for (int i=1;i <=m; i++)
	{
		mans+=(pjs-double(v[i]))*(pjs-double(v[i]));
	}
	mans=mans/double(m);
	return mans;
}
int main()
{
	srand(time(0));
	srand(rand());
	cin>>n>>m;
	for (int i=1; i<=n ;i++)
	{
		cin>>a[i];
	}
    for (int i=1;i <=n; i++)
    {
    	int minn=INT_MAX;
    	for (int j=1; j<=m; j++)
    	{
    		if (v[j]<minn)
    		{
    			minn=v[j];
    			q[i]=j;
			}
		}
		v[q[i]]+=a[i];
	}
	/*for (int i=1;i<=n; i++)
	{
		cout<<q[i]<<" "<<i<<endl;
	}*/
	ans=work();
	while (clock()/(1.0*CLOCKS_PER_SEC)<=tim)
	{
		double t=bt;
		//int cnt=0;
		while (t>et)
		{
			//cout<<++cnt<<" "<<1<<endl;
			cout<<"";
			int x,y;
			do
			{
				x=rand()%n+1;
				y=rand()%m+1;
			}while (q[x]==y);
			int kl=q[x];
			v[q[x]]-=a[x];
			v[y]+=a[x];
			q[x]=y;
			double nans=work();
			double de=nans-ans;
			if (de<0)
			{
				ans=nans;
			}
			else if (double(rand())/(RAND_MAX*1.0)>exp(double(-de)/bt))
			{
				q[x]=kl;
				v[q[x]]+=a[x];
				v[y]-=a[x];
			}
			t*=lt;
			//cout<<"t"<<" "<<t<<endl;
		}
	}
	cout<<fixed<<setprecision(2)<<sqrt(ans);
}