#include <bits/stdc++.h>

#define INF ll_MAX
#define UINF Ull_MAX
#define MOD 1000000007
#define REP(i,n) for(int i=0;i<n;i++)
#define all(x) (x).begin(),(x).end()

typedef long long ll;
typedef long double ld;
typedef std::pair<ll,ll> pii;

using namespace std;

// Cooley–Tukey FFT algorithm O(N log N)
vector<complex<ld>> fft(vector<complex<ld>> a, bool inverse = false) {
    ll n = a.size();
    ll h = 0; // h = log_2(n)
    for (ll i = 0; 1 << i < n; i++) h++;
    // バタフライ演算用の配置入れ替え
    for (ll i = 0; i < n; i++) {
        ll j = 0;
        for (ll k = 0; k < h; k++) j |= (i >> k & 1) << (h - 1 - k);
        if (i < j) swap(a[i], a[j]);
    }
    // バタフライ演算
    for (ll b = 1; b < n; b *= 2) {
        // 第 log_2(b) + 1 段
        // ブロックサイズ = b * 2
        for (ll j = 0; j < b; j++) {
            // ブロック内 j 個目
            // 重み w = (1 の原始 2b 乗根の j 乗)
            complex<ld> w =
                polar(1.0, (2 * M_PI) / (2 * b) * j * (inverse ? 1 : -1));
            for (ll k = 0; k < n; k += b * 2) {
                // k を先頭とするブロック
                complex<ld> s = a[j + k];         // 前
                complex<ld> t = a[j + k + b] * w; // 後
                a[j + k] = s + t;                     // 前の更新
                a[j + k + b] = s - t;                 // 後の更新
            }
        }
    }
    // 逆変換時にサイズで割る調整
    if (inverse)
        for (ll i = 0; i < n; i++) a[i] /= n;
    return a;
}
// Cooley–Tukey FFT algorithm O(N log N)
vector<complex<ld>> fft(vector<ld> a, bool inverse = false) {
    vector<complex<ld>> a_complex(a.size());
    for (ll i = 0; i < a.size(); i++) a_complex[i] = complex<ld>(a[i], 0);
    return fft(a_complex, inverse);
}
 
// FFT による畳み込み O(N log N)
vector<ld> convolve(vector<ld> a, vector<ld> b) {
    ll s = a.size() + b.size() - 1; // 畳み込み結果のサイズ
    ll t = 1; // FFT に使う配列のサイズ（2 の累乗）
    while (t < s) t *= 2;
    a.resize(t); // FFT するためにリサイズ
    b.resize(t); // FFT するためにリサイズ
    vector<complex<ld>> A = fft(a);
    vector<complex<ld>> B = fft(b);
    for (ll i = 0; i < t; i++) {
        A[i] *= B[i]; // 畳み込み結果の FFT 結果を得る
    }
    A = fft(A, true); // IFFT で畳み込み結果を得る
    a.resize(s);      // 畳み込み結果を入れるためにリサイズ
    for (ll i = 0; i < s; i++) a[i] = A[i].real(); // 実部が答え
    return a;
}

int main() {
    ll max_num = 100100;
    ll N,M;
    cin >> N >> M;
    vector<ld> a(max_num,0);
    a[0] = 0; // 0 円のメニューは存在しない
    for (ll i = 0; i < N; i++) {
        ll aa = 0;
        cin >> aa;
        a[aa]+=1.0; 
    }
 
    // 畳み込みで答えを計算
    vector<ld> ret = convolve(a, a);
 
    //ld → ll
    vector<ll> ret_ll(max_num*2,0);
    for(ll i = 0;i < ret_ll.size();i++){
        ret_ll[i] = (ll)(ret[i] + 0.5);
    }

    ll ans = 0;
    ll count = 0;

    for (ll i = ret_ll.size() - 1; i > 0 && count < M; i--) {
        if(ret_ll[i]){
            // for(ll j = 0;j < ret_ll[i] && count < M;j++){
            //     ans += i;
            //     count++;
            // }
            if(M >= count + (ll)ret_ll[i]){
                ans += (ll)ret_ll[i] * i;
                count += (ll)ret_ll[i];
            }else{
                ans += (ll)(M - count) * i ;
                break;
            }
        }
    }

    cout << ans << endl;

    return 0;
}

Measured memory_bytes: 38712