// #pragma GCC optimize("Ofast")
// // #pragma GCC optimize("unroll-loops")
// // #pragma GCC push_options
// // #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
// #include<bits/stdc++.h>
// #include <xmmintrin.h>
// #include <immintrin.h>
// using namespace::std;
// __attribute__((constructor))void init(){cin.tie(0);ios::sync_with_stdio(false);cout<<fixed<<setprecision(15);}
// #include<ext/pb_ds/assoc_container.hpp>
// #include<ext/pb_ds/tree_policy.hpp>
// #include <ext/pb_ds/priority_queue.hpp>
// #include<ext/pb_ds/tag_and_trait.hpp>
// // #include <boost/multiprecision/cpp_dec_float.hpp>
// // #include <boost/multiprecision/cpp_int.hpp>
// // namespace mp = boost::multiprecision;
// // typedef mp::number<mp::cpp_dec_float<0>> cdouble;
// // typedef mp::cpp_int cint;
// template<typename T>using pbds=__gnu_pbds::tree<T,__gnu_pbds::null_type,less<T>,__gnu_pbds::rb_tree_tag,__gnu_pbds::tree_order_statistics_node_update>;
// template<typename T>using pbds_map=__gnu_pbds::tree<T,T,less<T>,__gnu_pbds::rb_tree_tag,__gnu_pbds::tree_order_statistics_node_update>;
// template<typename T,typename E>using hash_map=__gnu_pbds::gp_hash_table<T,E>;
// template<typename T>using pqueue =__gnu_pbds::priority_queue<T, greater<T>,__gnu_pbds::rc_binomial_heap_tag>;
// typedef long long lint;
// #define INF (1LL<<60)
// #define IINF (1<<30)
// #define EPS (1e-10)
// #define MOD 1000000007LL
// //#define MOD 998244353LL
// typedef vector<lint> vec;
// typedef vector<vector<lint>> mat;
// typedef vector<vector<vector<lint>>> mat3;
// typedef vector<string> svec;
// typedef vector<vector<string>> smat;
// template<typename T>inline void numout(T t){bool f=0;for(auto i:t){cout<<(f?" ":"")<<i<INF/2?i:"INF";f=1;}cout<<endl;}
// template<typename T>inline void numout2(T t){for(auto i:t)numout(i);}
// template<typename T>inline void output(T t){bool f=0;for(auto i:t){cout<<(f?" ":"")<<i;f=1;}cout<<endl;}
// template<typename T>inline void output2(T t){for(auto i:t)output(i);}
// template<typename T>inline void _output(T t){bool f=0;for(lint i=0;i<t.size();i++){cout<<f?"":" "<<t[i];f=1;}cout<<endl;}
// template<typename T>inline void _output2(T t){for(lint i=0;i<t.size();i++)output(t[i]);}
// #define rep(i,n) for(lint i=0;i<lint(n);++i)
// #define repi(i,a,b) for(lint i=lint(a);i<(lint)(b);++i)
// #define rrep(i,n) for(lint i=lint(n)-1;i>=0;--i)
// #define rrepi(i,a,b) for(lint i=lint(b)-1;i>=lint(a);--i)
// #define irep(i) for(lint i=0;;++i)
// #define all(n) begin(n),end(n)
// #define dist(a,b,c,d) sqrt(pow(a-c,2)+pow(b-d,2))
// inline lint gcd(lint A,lint B){return B?gcd(B,A%B):A;}
// inline lint lcm(lint A,lint B){return A/gcd(A,B)*B;}
// // inline cint cgcd(cint A,cint B){return B?cgcd(B,A%B):A;}
// // inline cint clcm(cint A,cint B){return A/cgcd(A,B)*B;}
// //inline bool chmin(auto& s,const auto& t){bool res=s>t;s=min(s,t);return res;}
// //inline bool chmax(auto& s,const auto& t){bool res=s<t;s=max(s,t);return res;}
// const vector<lint> dx={1,0,-1,0,1,1,-1,-1};
// const vector<lint> dy={0,1,0,-1,1,-1,1,-1};
// #define SUM(v) accumulate(all(v),0LL)
// //auto call=[&](auto f,auto... args){return f(f,args...);};

// lint bs(lint mn,lint mx,function<bool(lint)> func){
//     mn--;
//     mx++;
// 	while(mx-mn>1){
// 		lint mid=(mn+mx)/2;
// 		if(!func(mid))mx=mid;
// 		else mn=mid;
// 	}
//     return mn;
// }
 
 
// int main(){
//     lint n,k;
//     cin>>n>>k;
//     vector<lint> a(n);
//     rep(i,n)cin>>a[i];
//     sort(all(a));
//     vec b=a;
//     reverse(all(b));
//     auto f=[&](lint x){
//         lint cnt=0;
//         rep(i,n){
//             if(a[i]*a[i]<=x)cnt--;
//             if(a[i]<0){
//                 cnt+=bs(0,n-1,[&](lint y){return b[y]*a[i]<=x;})+1;
//             }else if(a[i]>0){
//                 cnt+=bs(0,n-1,[&](lint y){return a[y]*a[i]<=x;})+1;
//             }else{
//                 if(x>=0)cnt+=n;
//             }
//         }
//         return cnt/2<k;
//     };
//     cout<<bs(-INF,INF,f)+1<<endl;
// }

#include<iostream>
#include<string>
#include<algorithm>
#include<cmath>
#include<functional>
#include<vector>

using namespace std;

typedef long long ll;

ll N, K;
ll a[200000 + 100];

int my_b_search(ll key, ll val, bool mode) {
	//mode -> trueならば昇順での二分探索 falseなら降順での二分探索
	ll ok = mode ? -1 : N;
	ll ng = mode ? N : -1;
	while (abs(ok - ng) > 1) {
		ll mid = (ok + ng) / 2;
		if (val * a[mid] <= key)ok = mid;
		else ng = mid;
	}
	return ok;
}

bool f(ll key) {
	ll cnt = 0;
	for (int i = 0; i < N; i++) {
		if (a[i] >= 0) {
			cnt += my_b_search(key, a[i], true) + 1;
		}
		else {
			cnt += N - my_b_search(key, a[i], false);
		}
		if (key >= a[i] * a[i])cnt--;
	}
	return cnt >= 2 * K;
}

int main() {
	cin >> N >> K;
	for (int i = 0; i < N; i++)cin >> a[i];

	sort(a, a + N);
	
	ll ng = -1e18 - 100, ok = 1e18 + 100;
	while (abs(ok - ng) > 1) {
		ll mid = (ok + ng) / 2;
		if (f(mid))ok = mid;
		else ng = mid;
	}
	cout << ok << endl;
	return 0;
}


Measured memory_bytes: 1792