#include <bits/stdc++.h>
#pragma GCC optimize("O3")
using namespace std;

#if __cplusplus < 201103L
#error C++11 or higher expected.
#endif
#define CPPV (__cplusplus / 100 % 100) // C++ Version

// define
#define endl "\n"
#ifdef LOCAL
#define lout cerr
#else
ostream devnull(0);
#define lout devnull
#endif
#define dump(x) lout << #x << " = " << (x) << endl;

#define all(a) (a).begin(), (a).end()
#define rall(a) (a).rbegin(), (a).rend()

// rep
// Inspired by:
// https://github.com/kurokoji/.cpp-Template/blob/master/template/template2.cpp
#define GET_MACRO(_1, _2, _3, _4, NAME, ...) NAME
#define rep(...)                                 \
  GET_MACRO(__VA_ARGS__, rep4, rep3, rep2, rep1) \
  (__VA_ARGS__)
#define rep1(n) rep2(_, n)
#define rep2(i, n) rep3(i, 0, n)
#define rep3(i, a, n) rep4(i, a, n, 1)
#define rep4(i, a, n, d) rep5(i, a, n, d, <)
#define rep5(i, a, n, d, o) \
  for (auto i = decltype(n)(a), i##_len = (n); i o i##_len; i += (d))

#define rrep(...)                                    \
  GET_MACRO(__VA_ARGS__, rrep4, rrep3, rrep2, rrep1) \
  (__VA_ARGS__)
#define rrep1(a) rrep2(_, a)
#define rrep2(i, a) rrep3(i, a, 0)
#define rrep3(i, a, n) rrep4(i, a, n, -1)
#define rrep4(i, a, n, d) rep5(i, (a)-1, n, d, >=)

#define repc(...)                                    \
  GET_MACRO(__VA_ARGS__, repc4, repc3, repc2, repc1) \
  (__VA_ARGS__)
#define repc1(n) repc2(_, n)
#define repc2(i, n) repc3(i, 0, n)
#define repc3(i, a, n) repc4(i, a, n, 1)
#define repc4(i, a, n, d) rep5(i, a, n, d, <=)

// alias
using ll = long long;
using lint = long long;

// func
template <class C>
void vcin(C &c)
{
  for (auto &&a : c)
  {
    cin >> a;
  }
}
template <class C, class D>
void vcin(C &c, D &d)
{
  for (auto it1 = c.begin(), it2 = d.begin(); it1 != c.aend() && it2 != d.aend();
       it1++, it2++)
  {
    cin >> *it1 >> *it2;
  }
}
template <class C>
void vcout(const C &c, const string delimiter = " ")
{
  for (auto it = c.begin(); it != c.aend(); it++)
  {
    if (it != c.begin())
      cout << delimiter;
    cout << *it;
  }
}
template <class T = int>
inline T input()
{
  T x;
  cin >> x;
  return (x);
}
template <typename T>
inline void chmax(T &x, const T y)
{
  if (x < y)
    x = y;
}
template <typename T>
inline void chmin(T &x, const T y)
{
  if (x > y)
    x = y;
}

#if CPPV >= 14
// http://beet-aizu.hatenablog.com/entry/2018/04/08/145516
template <typename T>
vector<T> make_v(size_t a) { return vector<T>(a); }
template <typename T, typename... Ts>
auto make_v(size_t a, Ts... ts)
{
  return vector<decltype(make_v<T>(ts...))>(a, make_v<T>(ts...));
}
template <typename T, typename V>
typename enable_if<is_class<T>::value == 0>::type fill_v(T &t, const V &v)
{
  t = v;
}
template <typename T, typename V>
typename enable_if<is_class<T>::value != 0>::type fill_v(T &t, const V &v)
{
  for (auto &e : t)
    fill_v(e, v);
}
#endif

constexpr int INF = 1 << 24;
constexpr lint LINF = 1LL << 48;

inline lint bit(int n)
{
  return 1LL << n;
}

inline bool bitof(lint n, int i)
{
  return (n >> i) & 1;
}

template <typename T>
inline bool between(T a, T x, T b)
{
  return a <= x && x <= b;
}

template <typename T, typename Fx>
T binary_search_func(T ok, T ng, Fx isOK, T eps = 1)
{
  while (abs(ok - ng) > eps)
  {
    T mid = ok + (ng - ok) / 2;
    if (isOK(mid))
      ok = mid;
    else
      ng = mid;
  }

  return ok;
}

void solve();

signed main()
{
  cin.tie(0);
  ios::sync_with_stdio(false);
  cout << fixed << setprecision(15);

  solve();

  return 0;
}

int sgn(lint x) { return (x > 0) - (x < 0); }

void solve()
{
  int n;
  lint k;
  cin >> n >> k;

  vector<int> a(n);
  vcin(a);

  sort(all(a));

  auto f = [&](lint x) {
    dump(x);
    lint cnt = 0;
    rep(i, 1, n)
    {

      auto divfloor = [](lint a, lint b) { return sgn(a) == sgn(b) ? a / b : -((abs(a) + abs(b) - 1) / abs(b)); };
      auto divceil  = [](lint a, lint b) { return sgn(a) != sgn(b) ? a / b :   (abs(a) + abs(b) - 1) / abs(b); };

      auto aend = a.begin() + i;
      if (a[i] > 0)
        cnt += upper_bound(a.begin(), aend, divfloor(x, a[i])) - a.begin();
      if (a[i] < 0)
        cnt += aend - lower_bound(a.begin(), aend, divceil(x, a[i]));
      if (a[i] == 0)
        cnt += (x >= 0) * i;
    }
    return cnt >= k;
  };

  cout << binary_search_func((lint)1e18, -(lint)1e18, f) << endl;
}


Measured memory_bytes: 1024