Problem: Every true king during his life must conquer the world, hold the Codeforces world finals, win pink panda in the shooting gallery and travel all over his kingdom.

King Copa has already done the first three things. Now he just needs to travel all over the kingdom. The kingdom is an infinite plane with Cartesian coordinate system on it. Every city is a point on this plane. There are n cities in the kingdom at points with coordinates (x1, 0), (x2, 0), ..., (xn, 0), and there is one city at point (xn + 1, yn + 1). 

King starts his journey in the city number k. Your task is to find such route for the king, which visits all cities (in any order) and has minimum possible length. It is allowed to visit a city twice. The king can end his journey in any city. Between any pair of cities there is a direct road with length equal to the distance between the corresponding points. No two cities may be located at the same point.

Input

The first line contains two integers n and k (1 ≤ n ≤ 105, 1 ≤ k ≤ n + 1) — amount of cities and index of the starting city. The second line contains n + 1 numbers xi. The third line contains yn + 1. All coordinates are integers and do not exceed 106 by absolute value. No two cities coincide.

Output

Output the minimum possible length of the journey. Your answer must have relative or absolute error less than 10 - 6.

Examples

Input

3 1
0 1 2 1
1


Output

3.41421356237309490000

Input

3 1
1 0 2 1
1


Output

3.82842712474619030000

Input

4 5
0 5 -1 -5 2
3


Output

14.24264068711928400000

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int MAXN = 100000;
double x[MAXN + 10], y;
int n, k;
double getans(int c) {
  double ret = 1e18;
  if (c == 1) {
    ret =
        ((ret) < (x[n] - x[1] +
                  (sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) + ((y) * (y)))))
             ? (ret)
             : (x[n] - x[1] +
                (sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) + ((y) * (y))))));
    for (int i = 1; i < n; i++)
      ret = ((ret) < (x[i] - x[1] +
                      (sqrt(((x[i] - x[n + 1]) * (x[i] - x[n + 1])) +
                            ((y) * (y)))) +
                      x[n] - x[i + 1] +
                      (((sqrt(((x[i + 1] - x[n + 1]) * (x[i + 1] - x[n + 1])) +
                              ((y) * (y))))) <
                               ((sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) +
                                      ((y) * (y)))))
                           ? ((sqrt(((x[i + 1] - x[n + 1]) *
                                     (x[i + 1] - x[n + 1])) +
                                    ((y) * (y)))))
                           : ((sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) +
                                    ((y) * (y)))))))
                 ? (ret)
                 : (x[i] - x[1] +
                    (sqrt(((x[i] - x[n + 1]) * (x[i] - x[n + 1])) +
                          ((y) * (y)))) +
                    x[n] - x[i + 1] +
                    (((sqrt(((x[i + 1] - x[n + 1]) * (x[i + 1] - x[n + 1])) +
                            ((y) * (y))))) <
                             ((sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) +
                                    ((y) * (y)))))
                         ? ((sqrt(
                               ((x[i + 1] - x[n + 1]) * (x[i + 1] - x[n + 1])) +
                               ((y) * (y)))))
                         : ((sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) +
                                  ((y) * (y))))))));
  } else
    ret =
        ((ret) < (x[n] - x[c] +
                  (sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) +
                        ((y) * (y)))) +
                  x[c - 1] - x[1] +
                  (((sqrt(((x[c - 1] - x[n + 1]) * (x[c - 1] - x[n + 1])) +
                          ((y) * (y))))) <
                           ((sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) +
                                  ((y) * (y)))))
                       ? ((sqrt(
                             ((x[c - 1] - x[n + 1]) * (x[c - 1] - x[n + 1])) +
                             ((y) * (y)))))
                       : ((sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) +
                                ((y) * (y)))))))
             ? (ret)
             : (x[n] - x[c] +
                (sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) + ((y) * (y)))) +
                x[c - 1] - x[1] +
                (((sqrt(((x[c - 1] - x[n + 1]) * (x[c - 1] - x[n + 1])) +
                        ((y) * (y))))) <
                         ((sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) +
                                ((y) * (y)))))
                     ? ((sqrt(((x[c - 1] - x[n + 1]) * (x[c - 1] - x[n + 1])) +
                              ((y) * (y)))))
                     : ((sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) +
                              ((y) * (y))))))));
  if (c == n) {
    ret =
        ((ret) < (x[n] - x[1] +
                  (sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) + ((y) * (y)))))
             ? (ret)
             : (x[n] - x[1] +
                (sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) + ((y) * (y))))));
    for (int i = 2; i <= n; i++)
      ret = ((ret) < (x[n] - x[i] +
                      (sqrt(((x[i] - x[n + 1]) * (x[i] - x[n + 1])) +
                            ((y) * (y)))) +
                      x[i - 1] - x[1] +
                      (((sqrt(((x[i - 1] - x[n + 1]) * (x[i - 1] - x[n + 1])) +
                              ((y) * (y))))) <
                               ((sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) +
                                      ((y) * (y)))))
                           ? ((sqrt(((x[i - 1] - x[n + 1]) *
                                     (x[i - 1] - x[n + 1])) +
                                    ((y) * (y)))))
                           : ((sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) +
                                    ((y) * (y)))))))
                 ? (ret)
                 : (x[n] - x[i] +
                    (sqrt(((x[i] - x[n + 1]) * (x[i] - x[n + 1])) +
                          ((y) * (y)))) +
                    x[i - 1] - x[1] +
                    (((sqrt(((x[i - 1] - x[n + 1]) * (x[i - 1] - x[n + 1])) +
                            ((y) * (y))))) <
                             ((sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) +
                                    ((y) * (y)))))
                         ? ((sqrt(
                               ((x[i - 1] - x[n + 1]) * (x[i - 1] - x[n + 1])) +
                               ((y) * (y)))))
                         : ((sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) +
                                  ((y) * (y))))))));
  } else
    ret =
        ((ret) < (x[c] - x[1] +
                  (sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) +
                        ((y) * (y)))) +
                  x[n] - x[c + 1] +
                  (((sqrt(((x[c + 1] - x[n + 1]) * (x[c + 1] - x[n + 1])) +
                          ((y) * (y))))) <
                           ((sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) +
                                  ((y) * (y)))))
                       ? ((sqrt(
                             ((x[c + 1] - x[n + 1]) * (x[c + 1] - x[n + 1])) +
                             ((y) * (y)))))
                       : ((sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) +
                                ((y) * (y)))))))
             ? (ret)
             : (x[c] - x[1] +
                (sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) + ((y) * (y)))) +
                x[n] - x[c + 1] +
                (((sqrt(((x[c + 1] - x[n + 1]) * (x[c + 1] - x[n + 1])) +
                        ((y) * (y))))) <
                         ((sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) +
                                ((y) * (y)))))
                     ? ((sqrt(((x[c + 1] - x[n + 1]) * (x[c + 1] - x[n + 1])) +
                              ((y) * (y)))))
                     : ((sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) +
                              ((y) * (y))))))));
  return ret;
}
void init() {
  scanf("%d%d", &n, &k);
  for (int i = 1; i <= n + 1; i++) scanf("%lf", &x[i]);
  scanf("%lf", &y);
  double ans = 1e18;
  if (k == n + 1) {
    sort(x + 1, x + n + 1);
    ans =
        (((sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) + ((y) * (y))))) <
                 ((sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) + ((y) * (y)))))
             ? ((sqrt(((x[1] - x[n + 1]) * (x[1] - x[n + 1])) + ((y) * (y)))))
             : ((sqrt(((x[n] - x[n + 1]) * (x[n] - x[n + 1])) +
                      ((y) * (y)))))) +
        x[n] - x[1];
  } else {
    int st = x[k];
    sort(x + 1, x + n + 1);
    for (int i = 1; i <= n; i++)
      if (x[i] == st) {
        k = i;
        break;
      }
    for (int i = 1; i <= n; i++) {
      double now = abs(x[i] - x[k]) + getans(i);
      ans = ((ans) < (now) ? (ans) : (now));
    }
  }
  printf("%.10lf\n", ans);
  return;
}
int main() {
  init();
  return 0;
}