// template version 1.6
// {{{ define basic macro
using namespace std;
#include <iostream>
#include <bits/stdc++.h>
// varibable settings
#define int long long
#define INF 1000000000000000000LL
// #define INF 2147483647
#define MOD 1000000007LL
// #define MOD 998244353LL
#define infile "../test/sample-1.in"

#define _overload3(_1,_2,_3,name,...) name
#define _rep(i,n) repi(i,0,n)
#define repi(i,a,b) for(int i=(int)(a);i<(int)(b);++i)
#define rep(...) _overload3(__VA_ARGS__,repi,_rep,)(__VA_ARGS__)
#define _rrep(i,n) rrepi(i,0,n)
#define rrepi(i,a,b) for(int i=(int)(b-1);i>=(int)(a);--i)
#define rrep(...) _overload3(__VA_ARGS__,rrepi,_rrep,)(__VA_ARGS__)
#define all(x) (x).begin(),(x).end()
#define sz(x) ((int)(x).size())
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define uni(x) sort(all(x));x.erase(unique(all(x)),x.end())
#define ten(n) ((int)1e##n)
template<class T, size_t N> size_t sza(const T (&array)[N]) { return N; }
template<class T> inline void chmax(T &a, const T &b) { if(a < b) a = b; }
template<class T> inline void chmin(T &a, const T &b) { if(a > b) a = b; }
template<class T=int> T in() {T x; cin>>x; return (x);}
struct Fast { Fast(){ std::cin.tie(0); ios::sync_with_stdio(false); } } fast;

// dump macro
#ifdef PCM
  #include "dump.hpp"
#else
  #define dump(...) 42
  #define dump_1d(...) 42
  #define dump_2d(...) 42
#endif

typedef long long ll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef long double ld;
typedef pair<int,int> pii;
typedef tuple<int,int,int> iii;
template<typename T> using PQ = priority_queue<T, vector<T>, greater<T>>;
int dx[]={1, -1, 0, 0};
int dy[]={0, 0, 1, -1};
#define fi first
#define se second

// }}}

int N,X;
int S[205];
int fac[205];
int memo[100005];

//BEGIN CUT HERE
// mod lib {{{
#define MAX_P 210
ll fact[MAX_P]; // fact[n]: n!
ll inv[MAX_P];  // inv[n]: n^(-1)
ll finv[MAX_P]; // finv[n]: (n!)^(-1)

void init(ll mod){  // 階乗関連テーブルの初期化 O(n)
  fact[0]=1;
  for(ll i=1;i<MAX_P;i++)
    fact[i]=(fact[i-1]*i)%mod;

  inv[1]=1;
  for(ll i=2;i<MAX_P;i++)
    inv[i]=inv[mod%i]*(mod-mod/i)%mod;

  finv[0]=1;  //　高速(10倍程度)な逆元階乗テーブルの作成。
  for(ll i=1;i<MAX_P;i++)
    finv[i]=finv[i-1]*inv[i]%mod;
}

ll extgcd(ll a,ll b,ll& x,ll& y){ //O(log(max(a,b)))
  ll d=a;
  if(b!=0){
    d=extgcd(b,a%b,y,x);
    y-=(a/b)*x;
  }else{
    x=1;y=0;
  }
  return d;
}
ll mod_inverse(ll a,ll mod){  // ユークリッドの拡張互除法による逆元
  ll x,y;
  extgcd(a,mod,x,y);
  return (mod+x%mod)%mod;
}

ll mod_pow(ll x,ll n,ll mod){ // 二分累乗法による階乗
  ll res=1;
  while(n){
    if(n&1) (res*=x)%=mod;
    (x*=x)%=mod;
    n>>=1;
  }
  return res;
}

ll mod_inverse2(ll a,ll mod){  // フェルマーの小定理による逆元
  return mod_pow(a,mod-2,mod);
}


ll mod_fact(ll n,ll mod,ll& e){ // n!=a*p^eとしたときのa mod pを返す。O(p)
  e=0;
  if(n==0) return 1;
  ll res=mod_fact(n/mod,mod,e);
  e+=n/mod;
  if(n/mod%2!=0)return res*(mod-fact[n%mod]) %mod;
  return res*fact[n%mod]%mod;
}

ll mod_comb(ll n,ll k,ll mod){
  if(n==k||k==0) return 1;
  ll e1,e2,e3;
  ll a1=mod_fact(n,mod,e1),a2=mod_fact(k,mod,e2),a3=mod_fact(n-k,mod,e3);
  if(e1>e2+e3) return 0;
  return a1*mod_inverse(a2*a3%mod,mod)%mod;
}

ll mod_comb2(ll n,ll k,ll mod){
  ll res=1;
  for(ll i=0;i<k;i++){
    res*=(n-i)%mod;
    res%=mod;
    res*=mod_inverse(i+1,mod);
    res%=mod;
  }
  return res;
}

//only for prime mod
ll mod_comb_pm(ll n,ll k,ll mod){  // テーブルがあればO(1)
  if(k<0||k>n) return 0;
  return fact[n]*finv[k]%mod*finv[n-k]%mod;
}

ll montmort(ll n,ll mod){
  ll res=0,inv=1;
  for(ll k=2;k<=n;k++){
    (inv*=mod_inverse(k,mod))%=mod;
    if(k%2) (res+=mod-inv)%=mod;
    else (res+=inv)%=mod;
  }
  for(ll i=1;i<=n;i++)
    (res*=i)%=mod;
  return res;
}

// calculate P(t) from given points in [0,N]
ll LagrangePolynomial(vector<ll> &y,ll t,const ll mod){
  init(mod);
  ll n=y.size()-1;
  ll num=1;
  for(ll i=0;i<=n;i++)
    num=num*((t-i)%mod)%mod;
  ll res=0;
  for(ll i=0;i<=n;i++){
    ll tmp=(y[i]*num%mod)*mod_inverse((t-i)%mod,mod)%mod;
    tmp=tmp*finv[i]%mod;
    tmp=tmp*finv[n-i]%mod;
    if((n-i)&1) tmp=mod-tmp;
    res=(res+tmp)%mod;
  }
  return res;
}
// }}}
//END CUT HERE

int rec(int x, int n){
  dump(x, n);
  if (n==0) return x;
  if (memo[x]!=-1) return memo[x];

  int res=0;
  rep(i, n){
    int m = upper_bound(S, S+n-1, x%S[i])-S;
    int tmp = (rec(x%S[i], m) * fac[n-m-1])%MOD;
    res +=  (tmp * mod_comb2(n-1, m, 1e9+7))%MOD;
    res %= MOD;
  }
  memo[x] = res;
  return res;
}

int solve(){
  cin>>N>>X;
  rep(i, N) {cin>>S[i];}
  sort(S, S+N);
  dump_1d(S, N);

  init(1e9+7);
  fac[0]=1; fac[1]=1;
  rep(i, 2, 205){
    fac[i] = (i*fac[i-1])%MOD;
  }

  rep(i, 100005) memo[i]=-1;

  dump_1d(memo, 6);
  cout << rec(X,N) << endl;

  return 0;
}

signed main() { //{{{
#ifdef INPUT_FROM_FILE
  std::ifstream in(infile);
  std::cin.rdbuf(in.rdbuf());
#endif
  solve();
  return 0;
} //}}}


Measured memory_bytes: 1024