#include <iostream>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <queue>
#include <vector>
#include <algorithm>
#include <map>

int N;
long long M;
std::vector<long long > A(100000);
std::vector<long long > kumisuu(200000, 0); // A_i>=x を満たすiの数
std::vector<long long > souwa(200000, 0); // A_i>=x を満たすiについてA_iの総和

long long get_possible_count(int value){ // 値がvalue以上となる組み合わせ総数を求める
  long long count = 0;
  for (int i=0; i<N; i++){
    count += kumisuu[std::max(value-A[i], (long long)0)];
  }
  return count;
}

int main(){
  scanf("%d",&N);
  scanf("%lld",&M);

  for (int i=0; i<N; i++){
    scanf("%lld",&A[i]);
    kumisuu[A[i]]++;
  }
  sort(A.begin(), A.end(), std::greater<long long>());

  for (int i=0; i<N; i++){
    souwa[i+1] = souwa[i] + A[i];
  }

  /*  累積和を求める */
  for (int i=200000; i>0; i--){
    kumisuu[i-1] += kumisuu[i];
  }

  int left = 0;
  int right = 200001;
  int mid;
  long long ans_temp = 0;
  std::map<long long, long long> mp;
  while(right-left>1){
    //printf("============\n");
    //printf("left=%d, right=%d\n", left, right);
    mid = (left+right)/2;
    ans_temp=0;
    for (int i=0; i<N; i++){
      ans_temp += souwa[kumisuu[std::max(mid-A[i],(long long)0)]];
      ans_temp += A[i]*kumisuu[std::max(mid-A[i],(long long)0)];
      //printf("ans += %lld + %lld\n", souwa[kumisuu[std::max(mid-A[i],(long long)0)]], A[i]*kumisuu[std::max(mid-A[i],(long long)0)]);
    }
    if (get_possible_count(mid)>=M){ // 
      left = mid;
    } else {
      right = mid;
    }
    mp[mid] = get_possible_count(mid);
  }
  //std::cout<<"asn_temp="<<ans_temp << std::endl;
  //printf("left=%d, right=%d\n", left, right);

  /*int target;
  long long target_count;
  if (get_possible_count(right)>=M){
    target = right;
    target_count= get_possible_count(right);
  } else {
    target = left;
    target_count= get_possible_count(left);
    }*/
  //printf("target=%d\n", target);
  /*long long ans;
  for (int i=0; i<N; i++){
    ans += souwa[kumisuu[std::max(target-A[i],(long long)0)]];
    //ans += souwa[std::max(target-A[i],(long long)0)];
    ans += A[i]*kumisuu[std::max(target-A[i],(long long)0)];
    //printf("souwa adding = %lld + %lld\n", souwa[kumisuu[std::max(target-A[i],(long long)0)]], A[i]*kumisuu[std::max(target-A[i],(long long)0)]);
    //printf("souwa adding + A[i] adding = %lld + %d = %lld\n", souwa[std::max(target-A[i],0)], A[i]*kumisuu[std::max(target-A[i],0)], souwa[std::max(target-A[i],0)]+A[i]*kumisuu[std::max(target-A[i],0)]);
    //printf("adding souwa[%d-%d] = souwa[%d] = %lld\n",target,A[i],std::max(target-A[i],0),souwa[std::max(target-A[i],0)]);
    }*/
    /*
  ans -= (long long)target*(target_count-M);
  printf("%lld\n", ans);
  */
  //printf("mp[left]=mp[%d]=%d\n", left, mp[left]);
  //  printf("mp[right]=mp[%d]=%d\n", right, mp[right]);
  if (left == mid){
    /*long long ans;
      for (int i=0; i<N; i++){
      ans += souwa[kumisuu[std::max(left-A[i],(long long)0)]];
      ans += A[i]*kumisuu[std::max(left-A[i],(long long)0)];
      }*/
    //std::cout << ans_temp-abs(M-get_possible_count(left))*left << std::endl;
    std::cout << ans_temp-abs(M-mp[left])*left << std::endl;
  } else if (right==mid){
    /*long long ans;
    for (int i=0; i<N; i++){
      ans += souwa[kumisuu[std::max(right-A[i],(long long)0)]];
      ans += A[i]*kumisuu[std::max(right-A[i],(long long)0)];
      }*/
    //std::cout << ans_temp+abs(M-get_possible_count(right))*left << std::endl;
    std::cout << ans_temp+abs(M-mp[right])*left << std::endl;
  }
  //printf("left=%d, right=%d\n", left, right);
  //printf("left count=%lld, right count=%lld\n", get_possible_count(left), get_possible_count(right));
  /*if(left==mid) {
  } else if(right==mid) {
    }*/
  return 0;
}


Measured memory_bytes: 4224