链接
https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1824
题解
每次计算平均数,我都对所有的数除以
2
2
2,我也可以过程中不除以
2
2
2,每次只求和,最后再除以
2
n
−
1
2^{n-1}
2n−1
然后就可以发现第
i
i
i个数(下表从
0
0
0开始)对答案的贡献是这个数乘以
C
n
−
1
i
C_{n-1}^i
Cn−1i
也就是说
a
n
s
=
∑
i
=
0
n
−
1
C
n
−
1
i
2
n
−
1
a
i
ans=\sum_{i=0}^{n-1}{C_{n-1}^i\over2^{n-1}}a_i
ans=i=0∑n−12n−1Cn−1iai
这个系数我不可能直接算出来,那样铁定爆
d
o
u
b
l
e
double
double
显然,
C
n
m
+
1
=
C
n
m
×
n
−
m
m
+
1
C_{n}^{m+1}=C_{n}^m\times {n-m\over m+1}
Cnm+1=Cnm×m+1n−m
我从左往右扫,我可以先不把所有的
2
2
2除完,如果我发现除以
2
2
2的过程中,这个数已经很小以至于对答案造不成影响,那我就不必继续除以
2
2
2了,直接把这个没算完的系数作为这个位置的系数,因为分子是递增的,所以分母也是递增的,除以
2
2
2的过程是均摊
O
(
n
)
O(n)
O(n)的。
右半部分可以直接利用组合数的对称性求解。
其它解法
百度之后,发现这个题还可以取对数做
乘除法直接变成对数的加减运算,数学函数真是厉害了
代码
//数学题
#include <bits/stdc++.h>
#define maxn 50010
#define eps 1e-15
using namespace std;
double C, a[maxn], ans;
int N;
void init()
{
int i;
scanf("%d",&N);
for(i=0;i<N;i++)scanf("%lf",a+i);
}
void work()
{
int i, res=N-1;
C=1, ans=0;
for(i=0;i<=N-1-i;i++)
{
while(res and C>eps)C/=2.0, res--;
if(i<N-1-i)ans+=C*(a[i]+a[N-1-i]);
else ans+=C*a[i];
C*=1.0*(N-1-i)/(i+1);
}
printf("%.3lf\n",ans);
}
int main()
{
int kase=0, T;
scanf("%d",&T);
while(T--)
{
printf("Case #%d: ",++kase);
init();
work();
}
return 0;
}