bzoj 2428: [HAOI2006]均分数据

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cstdlib>
 4 #include<ctime>
 5 #include<cstring>
 6 #include<cmath>
 7 #include<algorithm>
 8 #define M 100
 9 using namespace std;
10 int n,m,a[M],pos[M];
11 double rev,sum[M],ans,mn=1e30,pre;
12 void solve()
13 {
14     memset(sum,0,sizeof(sum));
15     ans=0;
16     for(int i=1;i<=n;i++)
17       pos[i]=rand()%m+1;
18     for(int i=1;i<=n;i++)
19       sum[pos[i]]+=a[i];
20     for(int i=1;i<=m;i++)
21       ans+=(sum[i]-rev)*(sum[i]-rev);
22     double T=10000;
23     for(;T>0.1;)
24       {
25         int t1=rand()%n+1,x=pos[t1],y=rand()%m+1;
26         if(x==y)
27           continue;
28         T*=0.9;
29         pre=ans;
30         ans-=(sum[x]-rev)*(sum[x]-rev);
31         ans-=(sum[y]-rev)*(sum[y]-rev);
32         sum[x]-=a[t1];
33         sum[y]+=a[t1];
34         ans+=(sum[x]-rev)*(sum[x]-rev);
35         ans+=(sum[y]-rev)*(sum[y]-rev);
36         if(ans>pre)
37           {
38             ans=pre;
39             sum[x]+=a[t1];
40             sum[y]-=a[t1];
41           }
42         else
43           pos[t1]=y;
44       }
45     if(ans<mn)
46       mn=ans;
47     return;
48 }
49 int main()
50 {
51     scanf("%d%d",&n,&m);
52     for(int i=1;i<=n;i++)
53     {
54     scanf("%d",&a[i]);
55     rev+=a[i];
56     swap(a[i],a[rand()%i+1]);
57     }
58     rev/=m;
59     for(int i=1;i<=10000;i++)solve();
60     printf("%.2lf",sqrt(mn/m));
61     return 0;
62 }

模拟退火 乱搞

时间: 2024-12-28 11:46:04

bzoj 2428: [HAOI2006]均分数据的相关文章

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包括一行,这一行只包含一个数,表示最小均方差的值(保留小数点后两位数字). 输入输出样例 输入样例#

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 ***********

[BZOJ2503][HAOI2006]均分数据

BZOJ Luogu sol 如果已经确定了一个序列,现要求把这个序列分成m个连续段作为答案,那么就可以用一个显而易见的DP DP显然可以得到当前序列下的最优解. 所以模拟退火瞎JB改一改序列每次DP一下就可以了 据说这题random_shuffle可以AC code #include<cstdio> #include<cstring> #include<algorithm> #include<cmath> #include<ctime> usi

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

[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是整数个数,

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. (同一