#include <cassert>
#include <algorithm>
#include <array>
#include <bitset>
#include <cctype>
#include <cstdint>
#include <cmath>
#include <complex>
#include <chrono>
#include <deque>
#include <functional>
#include <iomanip>
#include <iostream>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <unordered_map>
#include <vector>
#include <random>
#include <memory>
#include <utility>
#include <limits>
#include "limits.h"
 
#define rep(i, a, b) for (long long (i) = (a); (i) < (b); (i)++)
#define all(i) i.begin(), i.end()
#define debug(...) std::cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__)

void debug_out(){std::cerr<<std::endl;}

template<typename Head,typename... Tail>
void debug_out(Head h,Tail... t){
  std::cerr<<" "<<h;
  debug_out(t...);
}

template <typename T1, typename T2>
std::ostream& operator<<(std::ostream& os, std::pair<T1, T2> pa) {
  return os << pa.first << " " << pa.second;
}
 
template <typename T>
std::ostream& operator<<(std::ostream& os, std::vector<T> vec) {
  for (int i = 0; i < vec.size(); i++)os << vec[i] << (i + 1 == vec.size() ? "" : " ");
  return os;
}
 
template<typename T1,typename T2>
inline bool chmax(T1& a,T2 b){return a<b && (a=b,true);}
 
template<typename T1,typename T2>
inline bool chmin(T1& a,T2 b){return a>b && (a=b,true);}
 
long long pow_mod(long long a, long long b, long long mod=-1) {
  if(b==0)return 1;
  if ((a == 0)||(mod!=-1&&(a+mod)%mod==0)) {
    return 0;
  }
  long long x = 1;
  while (b > 0) {
    if (b & 1) {
      x = (mod!=-1)?(x * a) % mod:x*a;
    }
    a = (mod!=-1)?(a * a) % mod:a*a;
    b >>= 1;
  }
  return x;
}

// const long long MOD = 998244353;
const long long MOD = 1e9 + 7;

using ll = long long;
using P=std::pair<long long,long long>;

ll calc2(ll n,ll k,std::vector<ll> a){

  using namespace std;
   	const ll INF = (ll)(1e18)+100;
    sort(a.begin(), a.end());
    ll count[3] = {0};  //0:- 1:0 2:+   負、0、正それぞれの個数
    rep(i, 0,n){
        if(a[i] < 0)    count[0]++;
        else if(a[i] > 0)   count[2]++;
        else            count[1]++;
    }
    ll cdic[3] = {0};  //0:-  1:0  2:+ となるペアの個数、累積をとる
    cdic[0] = count[0]*count[2];
    cdic[1] = cdic[0] + count[1]*(count[0]+count[1]+count[2]) - count[1]*(count[1]+1)/2;
    cdic[2] = cdic[1] + count[0]*(count[0]-1)/2 + count[2]*(count[2]-1)/2;
    ll ans = 0;
    if(k <= cdic[0]){
        //ansが負
        vector<ll>pos;
        vector<ll>neg;
        rep(i, 0,n){
            if(a[i] < 0)        neg.push_back(a[i]);
            else if(a[i] > 0)   pos.push_back(a[i]);
        }
        ll ok = -INF, ng = 0;
        while(ok + 1 < ng){
            ll x = (ok+ng)/2;   //x未満がk未満
            ll ck = 0;
            rep(i, 0,pos.size()){ //正を固定して負を二分探索
                ll ok2 = -1, ng2 = neg.size();
                while(ok2 + 1 < ng2){
                    ll t = (ok2+ng2)/2;
                    if(pos[i]*neg[t] < x)   ok2 = t;
                    else                    ng2 = t;
                }
                ck += ng2;
            }
            if(ck < k)  ok = x;
            else        ng = x;
        }
        ans = ok;
    }else if(k <= cdic[1]){
        //ansが0
    }else{
        //ansが正
        k -= cdic[1];
        vector<ll>pos;
        vector<ll>neg;
        rep(i,0,n){
            if(a[i] < 0)        neg.push_back(a[i]);
            else if(a[i] > 0)   pos.push_back(a[i]);
        }
        rep(i,0, neg.size())  neg[i] = 0 - neg[i];
        sort(neg.begin(), neg.end());
        ll ok = 0, ng = INF;
        while(ok + 1 < ng){
            ll x = (ok+ng)/2;   //x未満がk未満
            ll ck = 0;
            //pos
            rep(i, 0,pos.size()){
                ll ok2 = -1, ng2 = pos.size();
                while(ok2 + 1 < ng2){
                    ll t = (ok2+ng2)/2;
                    if(pos[i]*pos[t] < x)   ok2 = t;
                    else                    ng2 = t;
                }
                ck += ng2;
                if(pos[i]*pos[i] < x)   ck--;
            }
            //neg
            rep(i, 0,neg.size()){
                ll ok2 = -1, ng2 = neg.size();
                while(ok2 + 1 < ng2){
                    ll t = (ok2+ng2)/2;
                    if(neg[i]*neg[t] < x)   ok2 = t;
                    else                    ng2 = t;
                }
                ck += ng2;
                if(neg[i]*neg[i] < x)    ck--;
            }
            if(ck/2 < k)  ok = x;
            else        ng = x;
        }
        ans = ok;
    }
    return ans;
}

int main(){
	ll n,k;
    std::cin>>n>>k;
    std::vector<ll> a(n);
    rep(i,0,n)std::cin>>a[i];
    
    std::cout<<calc2(n,k,a)<<"\n";
}


Measured memory_bytes: 5620