bzoj2428 [HAOI2006]均分数据

来源：互联网发布：什么软件可以升级win10 编辑：程序博客网时间：2024/06/06 02:13

一开始想了一个用优先队列的贪心，只能在洛谷过三个点（我还是too young too simple）。。。所以在网上看了一下

模拟退火

模拟退火之前在网上也看到过这个算法吧，但是一直都没有学，觉得难就没学，其实仔细看一下还是比较简单的

这题就是模拟退火的一个裸题~~，先上代码，详细的后面再补~~

模拟退火就是模拟物理上物体从高温到低温的一个过程，具体的网上也很多，兔子的那个比喻就很恰当

高温的时候物体内部的粒子比较混乱，所以它所在的地方可能性很多（大概是这个意思吧。。），那么在程序里面高温的时候状态转移过去的概率就比较大

反之，低温的时候物体内部粒子倾向于稳定，那么这个时候转移的概率比较小，除非当前状态比之前更优，那么直接模拟了

//Decision's template#include<cstdio>#include<cstring>#include<iostream>#include<cstdlib>#include<vector>#include<queue>#include<stack>#include<algorithm>#include<string>#include<cmath>#include<map>#include<set>using namespace std;#define DP_maxn 16#define maxn 100000#define INF 10000007#define MOD 1000000007#define mst(s,k) memset(s,k,sizeof(s))typedef long long ll;struct Edge{    int from,to,dist;    Edge(int u,int v,int d):from(u),to(v),dist(d){}};/*-------------------------------template End--------------------------------*/int a[30],n,k,pos[30];double minans = (double)INF,ans,ave = 0.0,sum[30];double sqr(double a) {return a*a;}void  Fire(){mst(sum,0);double r= 0.9,T = 10000.00;ans = 0.0;for(int i = 1;i<=n;i++){pos[i] = rand()%k+1;}for(int i = 1;i<=n;i++) sum[pos[i]] += a[i];for(int i = 1;i<=k;i++){ans+=(sum[i]-ave)*(sum[i]-ave);}//cout<<ans<<endl;while(T>0.1){T*=r;int t = rand()%n+1;int x = pos[t],y;if(T>500) y = min_element(sum+1,sum+k+1)-sum;else y = rand()%k+1; if(x==y) continue;double tmp = ans;ans -= (sum[x]-ave)*(sum[x]-ave);ans -= (sum[y]-ave)*(sum[y]-ave);sum[x]-=a[t];sum[y]+=a[t];ans += (sum[x]-ave)*(sum[x]-ave);ans += (sum[y]-ave)*(sum[y]-ave);if(ans<=tmp||rand()%10000<=T){pos[t] = y;} else{ans = tmp;sum[x]+=a[t];sum[y]-=a[t];}}minans = min(ans,minans);//cout<<minans<<endl;}int main(){//freopen("std.in","r",stdin);//freopen("std.out","w",stdout);mst(sum,0);cin>>n>>k;for(int i = 1;i<=n;i++){cin>>a[i];ave+=a[i];swap(a[i],a[rand()%i+1]);}ave /= k;for(int i = 1;i<=10000;i++) Fire();printf("%.2f\n",sqrt(minans/k));return 0;}

0 0