题目链接

模拟退火:

模拟退火!每次随机一个位置加给sum[]最小的组。

参数真特么玄学啊。。气的不想调了(其实就是想刷刷最优解)

如果用DP去算好像更准。。

//832kb	428ms
#include <cmath>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <assert.h>
#include <algorithm>
#define gc() getchar()
#define D (0.997)//big enough!
#define eps (1e-3)
#define Rand(x) (rand()%x+1)
const int N=22; int n,K,A[N],sum[7],bel[N];
double Ans=1e14,Aver,Tmp; inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline double Squ(double x){
return x*x;
}
inline void Move(int fr,int to,int tar){
sum[fr]-=A[tar], sum[bel[tar]=to]+=A[tar];
}
//inline void Update(double &x){
// x=0.0; for(int i=1; i<=K; ++i) x+=Squ(sum[i]);
//} void SA()
{
memset(sum,0,sizeof sum);
for(int i=1; i<=n; ++i) sum[bel[i]=Rand(K)]+=A[i];
double ans=0,nxt;
for(int i=1; i<=K; ++i) ans+=Squ(sum[i]);
Ans=std::min(Ans,ans+Tmp);
for(double T=1e7/*big enough*/; T>eps; T*=D)//Times:6898
{
int tar=Rand(n), fr=bel[tar], to=std::min_element(sum+1,sum+1+K)-sum;
if(fr==to) continue;
nxt=ans-Squ(sum[fr])-Squ(sum[to]);
Move(fr,to,tar);
nxt+=Squ(sum[fr])+Squ(sum[to]);
if(nxt<ans || (exp((ans-nxt)/T)*RAND_MAX>rand())) ans=nxt;//!
else Move(to,fr,tar);
Ans=std::min(Ans,ans+Tmp);//这个放里头!
}
} int main()
{
n=read(),K=read();
int sum=0;
for(int i=1; i<=n; ++i) sum+=(A[i]=read());
Aver=1.0*sum/K, Tmp=K*Aver*Aver-2.0*sum*Aver; // std::random_shuffle(A+1,A+1+n);//这东西。。有时好有时坏
for(int i=1; i<=30; ++i) SA();
printf("%.2lf",sqrt(Ans/(double)K)); return 0;
}

简单好写(错误率高)的裸随机化贪心:

使每组(sum_i-Average)尽量平均(也就是使Σ(sum_i^2)最小)。数据范围这么小,而且只保留两位。。

不连续分组很难办,但是random_shuffle()一下连续分组很多次就可以达到伪不连续分组的效果了。。

具体,我们可以随便分啊按照某种策略来分,比如依次分给当前sum最小的组。

被随机数种子一直卡一个点是怎样的体验。。我特么不设了。

//832kb	2448ms
#include <cmath>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
const int N=21; int n,K,A[N],sum[7];
double Ans=1e15,Aver,Tmp; inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
void Work()
{
std::random_shuffle(A+1,A+1+n);
memset(sum,0,sizeof sum);
for(int i=1; i<=n; ++i)
sum[std::min_element(sum+1,sum+1+K)-sum]+=A[i];
double res=0;
for(int i=1; i<=K; ++i) res+=1.0*sum[i]*sum[i];//Aver可以提出来。。
Ans=std::min(Ans,res+Tmp);
} int main()
{
n=read(),K=read();
int sum=0;
for(int i=1; i<=n; ++i) sum+=(A[i]=read());
Aver=1.0*sum/K, Tmp=K*Aver*Aver-2.0*sum*Aver;
for(int i=1; i<=300000; ++i) Work();
printf("%.2lf",sqrt(Ans/(double)K)); return 0;
}

BZOJ.2428.[HAOI2006]均分数据(随机化贪心/模拟退火)的更多相关文章

  1. bzoj 2428: [HAOI2006]均分数据 随机化

    2428: [HAOI2006]均分数据 Time Limit: 20 Sec  Memory Limit: 256 MB 题目连接 http://www.lydsy.com/JudgeOnline/ ...

  2. 洛谷 P2503 [HAOI2006]均分数据 随机化贪心

    洛谷P2503 [HAOI2006]均分数据(随机化贪心) 现在来看这个题就是水题,但模拟赛时想了1个小时贪心,推了一堆结论,最后发现贪心做 不了, 又想了半个小时dp 发现dp好像也做不了,在随机化 ...

  3. bzoj 2428: [HAOI2006]均分数据 && bzoj 3680 : 吊打XXX 模拟退火

    每次把元素随便扔随机一个初始解,退火时每次随机拿一个元素扔到随机一个集合里,当温度高时因为状态不稳定扔到那个元素和最小的里边. 如果新解优,更新ans. 把原式拆一下,就可以用int存了. bzoj ...

  4. bzoj 2428: [HAOI2006]均分数据

    #include<cstdio> #include<iostream> #include<cstdlib> #include<ctime> #inclu ...

  5. 2428: [HAOI2006]均分数据

    模拟退火.一种十分玄学的随机算法,网上可以查到比较详细的资料. 先随机地把数分成m组,每次随机地选择一个数,一开始直接选最小的一组,后来就随机一组,把这个数换到该组看看答案能不能变小,如果变小则换,如 ...

  6. bzoj2428: [HAOI2006]均分数据

    模拟退火.挺好理解的.然后res打成了ans一直WA一直WA...!!!一定要注意嗷嗷嗷一定要注意嗷嗷嗷一定要注意嗷嗷嗷. 然后我就一直卡一直卡...发现最少1800次的时候就可以出解了.然后我就去调 ...

  7. bzoj2428 [HAOI2006]均分数据 模拟退火

    [HAOI2006]均分数据 Time Limit: 5 Sec  Memory Limit: 128 MBSubmit: 3434  Solved: 1091[Submit][Status][Dis ...

  8. P2503 [HAOI2006]均分数据

    P2503 [HAOI2006]均分数据 模拟退火+dp (不得不说,我今天欧气爆棚) 随机出1个数列,然后跑一遍dp统计 #include<iostream> #include<c ...

  9. BZOJ 2428 JZYZOJ1533 : [HAOI2006]均分数据 模拟退火 随机化

    http://www.lydsy.com/JudgeOnline/problem.php?id=2428 http://172.20.6.3/Problem_Show.asp?id=1533 http ...

随机推荐

  1. UNIX网络编程 第3章 套接字编程简介

    套接字结构类型和相关的格式转换函数

  2. 32 Profiling Go Programs 分析go语言项目

    Profiling Go Programs  分析go语言项目 24 June 2011 At Scala Days 2011, Robert Hundt presented a paper titl ...

  3. 关于 poScreenCenter 与 poDesktopCenter

    主要是窗体水平方向与垂直方向的的 居中问题,由于水平方向 没有什么,所以不探讨.而垂直方向由于底部有个工具栏,工具栏自身有个高度,所以垂直方向的居中问题,需要探讨下. 结论: poScreenCent ...

  4. IntelliJ IDEA + Tomcat 部署问题

    首先要了解下 tomcat的 几种部署方式(大致分为静态部署和动态部署),可以百度,博客:http://qsfwy.iteye.com/blog/466461 IntelliJ IDEA 下部署项目的 ...

  5. Oracle学习笔记:实现select top N的方法

    由于Oracle不支持select top N语句,所以在Oracle中需要利用order by和rownum的组合来实现select top N的查询. rownum是记录表中数据编号的一个隐藏字段 ...

  6. Effective STL 学习笔记:19 ~ 20

    Effective STL 学习笔记:19 ~ 20 */--> div.org-src-container { font-size: 85%; font-family: monospace; ...

  7. XML文件解析-SaxReader

    一.前言 java解析xml文件有几种方式,这里介绍一下用SaxReader来解析Xml的方法. 二.准备 如果用SaxReader的话,需要引入jar包dom4j, 版本的话官网下载一个就好,这里用 ...

  8. 关于ARM指令那些你必须知道的东西

    1.32位ARM指令每一位都有其作用,具体如下: 低12为第二操作数, 12~15位为目的寄存器, 16~19位为第一操作数, 20~27就是操作码, 28~31就是条件域. 2.多寄存器load和s ...

  9. Pytest里,mark装饰器的使用,双引号,没引号,这种差别很重要

    按最新版的pytest测试框架. 如果只是单一的mark,不要加任何引号. 如果是要作and ,not之类的先把,一定要是双引号! 这个要记清楚,好像和以前版本的书上介绍的不一样,切记! import ...

  10. ssh远程执行nohup命令不退出

    https://blog.csdn.net/oneinmore/article/details/50073443