Problem: Radewoosh is playing a computer game. There are n levels, numbered 1 through n. Levels are divided into k regions (groups). Each region contains some positive number of consecutive levels.

The game repeats the the following process:

  1. If all regions are beaten then the game ends immediately. Otherwise, the system finds the first region with at least one non-beaten level. Let X denote this region.
  2. The system creates an empty bag for tokens. Each token will represent one level and there may be many tokens representing the same level.
    * For each already beaten level i in the region X, the system adds ti tokens to the bag (tokens representing the i-th level). 
    * Let j denote the first non-beaten level in the region X. The system adds tj tokens to the bag. 
  3. Finally, the system takes a uniformly random token from the bag and a player starts the level represented by the token. A player spends one hour and beats the level, even if he has already beaten it in the past. 



Given n, k and values t1, t2, ..., tn, your task is to split levels into regions. Each level must belong to exactly one region, and each region must contain non-empty consecutive set of levels. What is the minimum possible expected number of hours required to finish the game?

Input

The first line of the input contains two integers n and k (1 ≤ n ≤ 200 000, 1 ≤ k ≤ min(50, n)) — the number of levels and the number of regions, respectively.

The second line contains n integers t1, t2, ..., tn (1 ≤ ti ≤ 100 000).

Output

Print one real number — the minimum possible expected value of the number of hours spent to finish the game if levels are distributed between regions in the optimal way. Your answer will be considered correct if its absolute or relative error does not exceed 10 - 4.

Namely: let's assume that your answer is a, and the answer of the jury is b. The checker program will consider your answer correct if <image>.

Examples

Input

4 2
100 3 5 7


Output

5.7428571429


Input

6 2
1 2 4 8 16 32


Output

8.5000000000

Note

In the first sample, we are supposed to split 4 levels into 2 regions. It's optimal to create the first region with only one level (it must be the first level). Then, the second region must contain other three levels.

In the second sample, it's optimal to split levels into two regions with 3 levels each.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <typename A, typename B>
string to_string(pair<A, B> p);
template <typename A, typename B, typename C>
string to_string(tuple<A, B, C> p);
template <typename A, typename B, typename C, typename D>
string to_string(tuple<A, B, C, D> p);
string to_string(const string& s) { return '"' + s + '"'; }
string to_string(const char* s) { return to_string((string)s); }
string to_string(bool b) { return (b ? "true" : "false"); }
string to_string(vector<bool> v) {
  bool first = true;
  string res = "{";
  for (int i = 0; i < static_cast<int>(v.size()); i++) {
    if (!first) {
      res += ", ";
    }
    first = false;
    res += to_string(v[i]);
  }
  res += "}";
  return res;
}
template <size_t N>
string to_string(bitset<N> v) {
  string res = "";
  for (size_t i = 0; i < N; i++) {
    res += static_cast<char>('0' + v[i]);
  }
  return res;
}
template <typename A>
string to_string(A v) {
  bool first = true;
  string res = "{";
  for (const auto& x : v) {
    if (!first) {
      res += ", ";
    }
    first = false;
    res += to_string(x);
  }
  res += "}";
  return res;
}
template <typename A, typename B>
string to_string(pair<A, B> p) {
  return "(" + to_string(p.first) + ", " + to_string(p.second) + ")";
}
template <typename A, typename B, typename C>
string to_string(tuple<A, B, C> p) {
  return "(" + to_string(get<0>(p)) + ", " + to_string(get<1>(p)) + ", " +
         to_string(get<2>(p)) + ")";
}
template <typename A, typename B, typename C, typename D>
string to_string(tuple<A, B, C, D> p) {
  return "(" + to_string(get<0>(p)) + ", " + to_string(get<1>(p)) + ", " +
         to_string(get<2>(p)) + ", " + to_string(get<3>(p)) + ")";
}
void debug_out() { cerr << endl; }
template <typename Head, typename... Tail>
void debug_out(Head H, Tail... T) {
  cerr << " " << to_string(H);
  debug_out(T...);
}
const long long mod = 1e9 + 7;
int sum(int a, int b) { return (a + b) % mod; }
int sub(int a, int b) { return (a + mod - b) % mod; }
int mul(int a, int b) { return (1ll * a * b) % mod; }
int power(int a, int b) {
  int res = 1;
  while (b) {
    if (b & 1) res = mul(res, a);
    b >>= 1;
    a = mul(a, a);
  }
  return res;
}
class ConvexHullTrick {
  struct CHTPoint {
    double x, y;
    int idx;
    double lim;
  };
  std::vector<CHTPoint> pilha;
  inline double get_intersection(CHTPoint a, CHTPoint b);
  bool ccw(CHTPoint p0, CHTPoint p1, CHTPoint p2);

 public:
  void add_line(double a, double b, int idx);
  std::pair<int, double> get_maximum(double x);
};
inline double ConvexHullTrick::get_intersection(CHTPoint a, CHTPoint b) {
  double denom = (double(b.x) - a.x);
  double num = (double(b.y) - a.y);
  return -num / denom;
}
bool ConvexHullTrick::ccw(CHTPoint p0, CHTPoint p1, CHTPoint p2) {
  return ((p1.y - p0.y) * (p2.x - p0.x) > (p2.y - p0.y) * (p1.x - p0.x));
}
void ConvexHullTrick::add_line(double a, double b, int idx) {
  CHTPoint novo = {a, b, idx, 0};
  int tam = pilha.size();
  while (tam >= 2 && !ccw(pilha[tam - 2], pilha[tam - 1], novo)) {
    pilha.pop_back();
    tam--;
  }
  while (tam >= 1 && fabs(pilha[tam - 1].x - a) < 1e-8) {
    pilha.pop_back();
    tam--;
  }
  pilha.push_back(novo);
  if (tam >= 1)
    pilha[tam - 1].lim = get_intersection(pilha[tam - 1], pilha[tam]);
}
std::pair<int, double> ConvexHullTrick::get_maximum(double x) {
  int st = 0, ed = pilha.size() - 1;
  while (st < ed) {
    int mid = (st + ed) / 2;
    if (pilha[mid].lim < x)
      st = mid + 1;
    else
      ed = mid;
  }
  return std::make_pair(pilha[st].idx, pilha[st].x * x + pilha[st].y);
}
int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  int N, K;
  cin >> N >> K;
  vector<int> A(N + 1);
  for (int i = 1; i < N + 1; i++) cin >> A[i];
  vector<double> sum(N + 1), pre(N + 1), rev(N + 1);
  for (int i = 1; i < N + 1; i++) {
    sum[i] = A[i] + sum[i - 1];
    rev[i] = (1.0 / double(A[i])) + rev[i - 1];
    pre[i] = pre[i - 1] + (sum[i] / double(A[i]));
  }
  vector<ConvexHullTrick> CH(K + 1);
  vector<vector<double>> dp(N + 1, vector<double>(K + 1, 0));
  CH[0].add_line(0, 0, 0);
  for (int j = 1; j < N + 1; j++) {
    for (int k = 1; k < min(j, K) + 1; k++) {
      pair<int, double> act = CH[k - 1].get_maximum(rev[j]);
      double add = -act.second;
      dp[j][k] = pre[j] + add;
      double m = -sum[j];
      double b = dp[j][k] - pre[j] + rev[j] * sum[j];
      CH[k].add_line(-m, -b, 0);
    }
  }
  cout << fixed << setprecision(10) << dp[N][K] << '\n';
  return 0;
}