usaco 4.2 Cowcycles

 搜索题目，只需要两个优化即可，其他的效果并不明显
1）因为递归搜索出的数列是有序的，所以将分子升序，分母降序排列后，再按顺序相除，得到的数列将“较为有序”，采用插入排

序比较好，因为c++的sort函数在数组小的时候会自动insert_sort()，所以还是可以用sort的。这个优化是这道题目的关键！
2）计算方差可以采用 D(x) = E(x*x) - E(x)*E(x)

这道题目花了很长时间，总结如下：
1）对于搜索最终更新答案的函数要尽量优化，对较为有序的数列排序，速度很快
2）当遇到像编译器错误的问题时，一般是数组大小没开够，或者开错了，这次就是把ans[2][N]写成了ans[N][2]导致错误

 

1. #include <iostream>
2. #include <algorithm>
3. #include <map>
4. #include <vector>
5. using namespace std;
7. /*
8. PROG: cowcycle
9. LANG: C++
10. ID: heben991
11. */
13. const int N = 60;
15. int f[N], r[N], F, R, f1, f2, r1, r2;
16. int ans[2][N];
17. double a[N], mini=1e15;
20. void update()
21. {
22.     double mean=0, vari=0, t;
23.     int i, j, k=0, p;
25.     for(i = 1; i <= F; ++i)
26.     for(j = R; j >= 1; --j)
27.     {
28.         t = (double)f[i]/r[j];
29.         a[++k] = t;
30.     }
31.     sort(a+1,a+k+1);
33.     --k;
34.     for(i = 1; i <= k; ++i)
35.     {
36.         a[i]=a[i+1]-a[i];
37.         mean += a[i];
38.         vari += a[i]*a[i];
39.     }
40.     mean /= k;
41.     vari = vari/k - mean*mean;
43.     if(vari < mini)
44.     {
45.         mini = vari;
47.         memcpy(ans[0], f+1, F*sizeof(int));
48.         memcpy(ans[1], r+1, R*sizeof(int));
52.     }
53. }
55. void dfs_rear(int k, int low)
56. {
57.     if(k > R)
58.     {
59.         int i, j;
60.         if( f[F]*r[R] < 3*f[1]*r[1] ) return;
62.         update();
63.     }
64.     else
65.     {
66.         for(r[k] = low; r[k]-R+k <= r2; ++r[k])
67.         {
68.             dfs_rear(k+1,r[k]+1);
69.         }
70.     }
71. }
73. void dfs_front(int k, int low)
74. {
75.     if(k > F)
76.     {
77.         dfs_rear(1, r1);
78.     }
79.     else
80.     {
81.         for(f[k] = low; f[k]-F+k <= f2; ++f[k])
82.         {
83.             dfs_front(k+1, f[k]+1);
84.         }
85.     }
86. }
88. int main()
89. {
90.     int i, j, k, t, h;
92.     freopen("cowcycle.in", "r", stdin);
93.     freopen("cowcycle.out","w",stdout);
95.     scanf("%d%d%d%d%d%d", &F,&R,&f1,&f2,&r1,&r2);
97.     dfs_front(1, f1);
99.     printf("%d", ans[0][0]);
100.     for(i = 1; i < F; ++i) printf(" %d", ans[0][i]);
101.     puts("");
102.     printf("%d", ans[1][0]);
103.     for(i = 1; i < R; ++i) printf(" %d", ans[1][i]);
104.     puts("");
105.     return 0;
106. }