Problem:
Score : 300 points

Problem Statement

There are 3N participants in AtCoder Group Contest.
The strength of the i-th participant is represented by an integer a_i.
They will form N teams, each consisting of three participants.
No participant may belong to multiple teams.

The strength of a team is defined as the second largest strength among its members.
For example, a team of participants of strength 1, 5, 2 has a strength 2, and a team of three participants of strength 3, 2, 3 has a strength 3.

Find the maximum possible sum of the strengths of N teams.

Constraints

1 ≤ N ≤ 10^5

1 ≤ a_i ≤ 10^{9}

a_i are integers.

Input

Input is given from Standard Input in the following format:

N
a_1 a_2 ... a_{3N}

Output

Print the answer.

Sample Input 1

2
5 2 8 5 1 5

Sample Output 1

10

The following is one formation of teams that maximizes the sum of the strengths of teams:

Team 1: consists of the first, fourth and fifth participants.

Team 2: consists of the second, third and sixth participants.

Sample Input 2

10
1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000 1000000000

Sample Output 2

10000000000

The sum of the strengths can be quite large.

Solution in C++14 (GCC 5.4.1) (judge verdict: Runtime Error)

#include <bits/stdc++.h>
using namespace std;

const int RLEN=1<<18|1;
inline char nc() {
	static char ibuf[RLEN],*ib,*ob;
	(ib==ob) && (ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
	return (ib==ob) ? -1 : *ib++;
}
inline int rd() {
	char ch=nc(); int i=0,f=1;
	while(!isdigit(ch)) {if(ch=='-')f=-1; ch=nc();}
	while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=nc();}
	return i*f;
}

const int N=2e5+50, mod=1e9+7;
inline int mul(int x,int y) {return (long long)x*y%mod;}
inline int power(int a,int b,int rs=1) {for(;b;b>>=1,a=mul(a,a)) if(b&1) rs=mul(rs,a); return rs;}

int n,X,Y,cov[N];
vector <int> c[N];

struct combin {
	int fac[N],ifac[N];
	combin() {
		fac[0]=1;
		for(int i=1;i<N;i++) fac[i]=mul(fac[i-1],i);
		ifac[N-1]=power(fac[N-1],mod-2);
		for(int i=N-2;~i;i--) ifac[i]=mul(ifac[i+1],i+1);
	}
} C;

int main() {
	n=rd(), X=rd(), Y=rd();
	for(int i=1;i<=n;i++) {
		int col=rd(), w=rd();
		c[col].push_back(w);
	}
	for(int i=1;i<=n;i++) sort(c[i].begin(),c[i].end());
	int mn=2e9, id=23333, sec=2e9;
	for(int i=1;i<=n;i++) if(c[i].size() && c[i].front()<mn) mn=c[i].front(), id=i;
	for(int i=1;i<=n;i++) if(c[i].size() && (i==id || c[i].front()+mn<=Y)) {cov[i]=1; sec=min(sec,(i!=id) ? c[i].front() : sec);}
	int prd=1, sum=0;
	for(int i=1;i<=n;i++) {
		if(cov[i]) {
			int j=1;
			while(j<c[i].size() && (c[i][j]+c[i].front()<=X || (i!=id && c[i][j]+mn<=Y) || (i==id && c[i][j]+sec<=Y))) ++j;
			prd=mul(prd,C.ifac[j]); sum=sum+j;
		}
	}
	cout<<mul(C.fac[sum],prd)<<'\n';
}