#include <set>
#include <map>
#include <queue>
#include <ctime>
#include <cmath>
#include <cstdio>
#include <bitset>
#include <vector>
#include <cstring>
#include <cassert>
#include <iostream>
#include <algorithm>
using namespace std;
typedef double db;
typedef long long ll;
typedef long double ld;
typedef vector<int> vi;
typedef pair<int,int> pa;
typedef unsigned int uint;
typedef unsigned long long ull;
#define w1 first
#define ls (x<<1)
#define w2 second
#define ins insert
#define rs (x<<1|1) 
#define mp make_pair
#define pb push_back
#define mid ((l+r)>>1)
#define sqr(x) ((x)*(x))
#define cle(x) ((x).clear())
#define lowbit(x) ((x)&(-x))
#define SZ(x) (int((x).size()))
#define All(x) (x).begin(),(x).end()
#define ms(x,y) memset(x,y,sizeof (x))
#define rep(i,a,b) for(int (i)=(a);(i)<=(b);(i)++)
#define rep2(i,a,b) for(int (i)=(a);(i)<(b);(i)++)
#define per(i,a,b) for(int (i)=(a);(i)>=(b);(i)--)
#define Rep(p,x) for(int (p)=head[(x)];(p);(p)=nxt[(p)]) 
#define Rep2(p,x) for(int (p)=cur[(x)];(p);(p)=nxt[(p)])
template<class T>inline void read(T&num){
	num=0;T f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9')num=num*10+ch-'0',ch=getchar();
	num*=f;
}
inline int getgcd(int x,int y){if(!x)return y;return getgcd(y%x,x);}
inline int power(int x,int k,int p){int res=1;for(;k;k>>=1,x=(ll)x*x%p)if(k&1)res=(ll)res*x%p;return res;}
const double pi=acos(-1);	
inline void judge(){
	freopen("input.txt","r",stdin);
} 
//********************************head*************************************
const int maxn=1e2+5,mod=1e9+7;
int n;
int a[maxn],f[maxn][maxn][maxn];
void upd(int &x,int y){
	x=(x+y)%mod;
}
int main(){
	read(n);
	rep(i,1,2*n-1)read(a[i]);
	sort(a+1,a+2*n-1+1);
	a[0]=a[1];a[2*n]=a[2*n-1];
	f[n][a[n]!=a[n-1]][a[n]!=a[n+1]]=1;
	per(i,n,2)rep(l,0,2*n-1)rep(r,0,2*n-1){
		int fl=(a[i-1]!=a[i-2]),fr=(a[2*n-i+1]!=a[2*n-i+2]);
		rep(k,0,l-1)upd(f[i-1][k+fl][r+fr+1],f[i][l][r]);
		upd(f[i-1][l+fl][r+fr],f[i][l][r]);
		rep(k,0,r-1)upd(f[i-1][l+fl+1][k+fr],f[i][l][r]);
	}
	int ans=0;
	rep(l,0,2*n-1)rep(r,0,2*n-1)upd(ans,f[1][l][r]);
	cout<<ans<<endl;
	return 0;
}