bzoj4428 均分数据

Description

已知N个正整数:A1、A2、……、An 。今要将它们分成M组,使得各组数据的数值和最平均,即各组的均方差最小。均方差公式如下:

,其中σ为均方差,是各组数据和的平均值,xi为第i组数据的数值和。

Input

第一行是两个整数,表示N,M的值(N是整数个数,M是要分成的组数)

第二行有N个整数,表示A1、A2、……、An。整数的范围是1--50。

(同一行的整数间用空格分开)

Output

这一行只包含一个数,表示最小均方差的值(保留小数点后两位数字)。
首先任取一个初始的分组方案,进行两种随机调整操作并模拟退火确定是否进行:

1.随机在一个非空组中选出一个数放入另一组中

2.随机选两个不同组中的数并交换位置

#include<cstdio>
#include<cmath>
#include<cstdlib>
int n,m;
int v[24];
int s[24][24],sp[24],sum[24];
double Ans=1e20;
double calc(){
    double x=0,ans=0;
    for(int i=0;i<m;i++)x+=sum[i];
    x/=m;
    for(int i=0;i<m;i++)ans+=(sum[i]-x)*(sum[i]-x);
    if(ans<Ans)Ans=ans;
    return ans;
}
void run(double k){
    double v=calc();
    double l=1024;
    int a,b,c,c2,d;
    while(l>1){
        a=rand()%m,b=rand()%m;
        if(a==b||!sp[a])continue;
        c=rand()%sp[a],c2,d=s[a][c];
        sum[a]-=d;sum[b]+=d;
        double v1=calc();
        if(v1<v||(rand()&1023)<l){
            s[b][sp[b]++]=d;
            s[a][c]=s[a][--sp[a]];
            v=v1;
        }else sum[a]+=d,sum[b]-=d;

        a=rand()%m,b=rand()%m;
        if(a==b||!sp[a]||!sp[b])continue;
        c=rand()%sp[a],c2=rand()%sp[b],d=s[a][c]-s[b][c2];
        sum[a]-=d;sum[b]+=d;
        v1=calc();
        if(v1<v||(rand()&1023)<l){
            int x=s[b][c2];
            s[b][c2]=s[a][c];
            s[a][c]=x;
            v=v1;
        }else sum[a]+=d,sum[b]-=d;

        l*=k;
    }
}
int main(){
    srand(1844677);
    scanf("%d%d",&n,&m);
    for(int i=0;i<n;i++)scanf("%d",v+i);
    for(int i=0;i<n;i++)s[i%m][sp[i%m]++]=v[i],sum[i%m]+=v[i];
    for(int i=0;i<1000;i++)run(0.95);
    for(int i=0;i<2000;i++)run(0.92);
    for(int i=0;i<2000;i++)run(0.9);
    for(int i=0;i<1000;i++)run(0.8);
    printf("%.2lf",sqrt(Ans/m));
    return 0;
}
时间: 2024-10-14 05:05:44

bzoj4428 均分数据的相关文章

BZOJ2428 均分数据

2428: [HAOI2006]均分数据 Time Limit: 5 Sec  Memory Limit: 128 MB Description 已知N个正整数:A1.A2.…….An .今要将它们分成M组,使得各组数据的数值和最平均,即各组的均方差最小.均方差公式如下: ,其中σ为均方差,是各组数据和的平均值,xi为第i组数据的数值和. Input 第一行是两个整数,表示N,M的值(N是整数个数,M是要分成的组数) 第二行有N个整数,表示A1.A2.…….An.整数的范围是1--50. (同一

[HAOI2006][BZOJ2428] 均分数据

2428: [HAOI2006]均分数据 Time Limit: 5 Sec  Memory Limit: 128 MBSubmit: 801  Solved: 231[Submit][Status][Discuss] Description 已知N个正整数:A1.A2.…….An .今要将它们分成M组,使得各组数据的数值和最平均,即各组的均方差最小.均方差公式如下: ,其中σ为均方差,是各组数据和的平均值,xi为第i组数据的数值和. Input 第一行是两个整数,表示N,M的值(N是整数个数,

2428: [HAOI2006]均分数据

模拟退火.一种十分玄学的随机算法,网上可以查到比较详细的资料. 先随机地把数分成m组,每次随机地选择一个数,一开始直接选最小的一组,后来就随机一组,把这个数换到该组看看答案能不能变小,如果变小则换,如果没有变小,按模拟退火的玄学方式判断一下,也要交换. srand(time(0))在bzoj会RE,不知为何. 非常玄学的正确性,没拍多少数据就会出现差别,然而OJ上可以过的. //Twenty #include<cstdio> #include<cstdlib> #include&l

洛谷P2503 [HAOI2006]均分数据(模拟退火)

题目描述 已知N个正整数:A1.A2.…….An .今要将它们分成M组,使得各组数据的数值和最平均,即各组的均方差最小.均方差公式如下: 输入输出格式 输入格式: 输入文件data.in包括: 第一行是两个整数,表示N,M的值(N是整数个数,M是要分成的组数) 第二行有N个整数,表示A1.A2.…….An.整数的范围是1--50. (同一行的整数间用空格分开) 输出格式: 输出文件data.out包括一行,这一行只包含一个数,表示最小均方差的值(保留小数点后两位数字). 输入输出样例 输入样例#

bzoj 2428 均分数据

题目大意: 已知N个正整数 将它们分成M组,使得各组数据的数值和最平均,即各组的均方差最小 求最小均方差 思路: 模拟退火 1 #include<iostream> 2 #include<cstdio> 3 #include<cmath> 4 #include<cstdlib> 5 #include<cstring> 6 #include<algorithm> 7 #include<queue> 8 #include<

P2503 [HAOI2006]均分数据

题目描述 已知N个正整数:A1.A2.…….An .今要将它们分成M组,使得各组数据的数值和最平均,即各组的均方差最小.均方差公式如下: 输入输出格式 输入格式: 输入文件data.in包括: 第一行是两个整数,表示N,M的值(N是整数个数,M是要分成的组数) 第二行有N个整数,表示A1.A2.…….An.整数的范围是1--50. (同一行的整数间用空格分开) 输出格式: 输出文件data.out包括一行,这一行只包含一个数,表示最小均方差的值(保留小数点后两位数字). 输入输出样例 输入样例#

[模拟退火] HAOI2006 均分数据

题目描述 已知N个正整数:A1.A2.--.An .今要将它们分成M组,使得各组数据的数值和最平均,即各组的均方差最小.均方差公式如下: 输入输出格式 输入格式: 输入文件data.in包括: 第一行是两个整数,表示N,M的值(N是整数个数,M是要分成的组数) 第二行有N个整数,表示A1.A2.--.An.整数的范围是1--50. (同一行的整数间用空格分开) 输出格式: 输出文件data.out包括一行,这一行只包含一个数,表示最小均方差的值(保留小数点后两位数字). 输入输出样例 输入样例#

[HAOI2006] 均分数据 - dp,随机化

将 \(n\) 个数分成 \(m\) 组,使得各组数据的数值和最平均,即各组的均方差最小.\(n\leq 20,m \leq 6\) Solution 考虑如果是按顺序分段,那么设 \(f[i]\) 表示前 \(i\) 个数积攒的最小方差,dp即可 用 std::random_shuffle 随机若干顺序,可以证明在很大概率上能随到答案 当然内部也可以贪心,每次尽量使得每个集合的总和相等 #include <bits/stdc++.h> using namespace std; const i

BZOJ2428 [HAOI2006]均分数据

正确方法是模拟退火? 根据wulala的讲法,我们多random_shuffle几发,对序列贪心一下就好了 反正都是随机化的做法,能过就行2333 1 /************************************************************** 2 Problem: 2428 3 User: rausen 4 Language: C++ 5 Result: Accepted 6 Time:3356 ms 7 Memory:816 kb 8 ***********