#include using namespace std; template class RangeMinimumQuery { private: int numLeaves; vector nodes; public: RangeMinimumQuery(int n) { numLeaves = 1 << (int)log2(n * 2 - 1); nodes.resize(numLeaves * 2, numeric_limits::max()); } RangeMinimumQuery(const vector &a) { numLeaves = 1 << (int)log2(a.size() * 2 - 1); nodes.resize(numLeaves * 2, numeric_limits::max()); for (int i = 0; i < a.size(); i++) { nodes[i + numLeaves] = a[i]; } for (int i = numLeaves - 1; i >= 1; i--) { nodes[i] = min(nodes[i * 2], nodes[i * 2 + 1]); } } // 0-indexed void setValue(int index, T value) { nodes[index + numLeaves] = value; for (int i = index + numLeaves; i > 1; i >>= 1) { nodes[i >> 1] = min(nodes[i], nodes[i ^ 1]); } } // 0-indexed [l,r) T minimum(int l, int r) { T result = numeric_limits::max(); l += numLeaves; r += numLeaves; for (; l < r; l >>= 1, r >>= 1) { if (l & 1) { result = min(result, nodes[l++]); } if (r & 1) { result = min(result, nodes[--r]); } } return result; } }; int main() { int n, m; cin >> n >> m; const int M = n * 2 - 2; vector c(M * 2); for (int i = 1; i < n - 1; i++) { cin >> c[i]; } for (int i = 0; i < n - 2; i++) { c[n + i] = c[n - 2 - i]; } for (int i = 0; i < M; i++) { c[i + M] = c[i]; } const int T = 10000000 / M; vector E(M * 2), EE(M * 2); for (int ii = 0; ii < T; ii++) { RangeMinimumQuery rmq(M * 2); double sumE = 0; double sumEE = 0; for (int i = 0; i < M; i++) { E[i + M] = E[i]; EE[i + M] = EE[i]; E[i] = 0; EE[i] = 0; rmq.setValue(i + M, EE[i + M] + c[i + M]); } for (int i = 1; i <= m; i++) { sumE += E[M - 1 + i] + c[M - 1 + i]; sumEE += EE[M - 1 + i] + c[M - 1 + i]; } for (int i = M - 1; i >= 0; i--) { if (i == n - 1) { E[i] = 0; EE[i] = 0; } else { double mini = rmq.minimum(i + 1, i + m + 1); E[i] += sumE / m; EE[i] += sumEE / m; E[i] = min(E[i], mini); } rmq.setValue(i, EE[i] + c[i]); sumE += E[i] + c[i]; sumEE += EE[i] + c[i]; sumE -= E[i + m] + c[i + m]; sumEE -= EE[i + m] + c[i + m]; } } printf("%.20f\n", (double)min(E[0], EE[0])); }