#include using namespace std; #define ll long long #define all(x) x.begin(),x.end() #define rep(i,a,b) for(int i=a;i<=b;i++) #define rep_r(i,a,b) for(int i=a;i>=b;i--) #define each(a,x) for (auto& x : a) using pi = pair; using pl = pair; using vi = vector; using vl = vector; #define sz(x) int(x.size()) #define so(x) sort(all(x)) #define so_r(x) sort(all(x),greater()) #define lb lower_bound #define ub upper_bound const char nl = '\n'; int dx[4] = {1,-1,0,0}; int dy[4] = {0,0,1,-1}; int bit_cnt(int x){ return __builtin_popcount(x); } ll bex(ll a, ll b, ll mod = 1e9 + 7){ll res = 1LL; while(b){ if (b&1) res = res * a % mod; a = a * a % mod; b >>= 1;} return res;} template bool chmax(t&a,u b){if(a bool chmin(t&a,u b){if(b fac(int x) { vector f; for (int i = 1; 1LL * i * i <= n; i++){ if (x % i == 0) { f.push_back(i); if (x / i != i) f.push_back(x/i); } } sort(all(f)); return f; } map dp[2005]; void solve(){ cin >> n >> k; rep(i,1,n) cin >> a[i]; dp[0][1] = 1; rep(i,1,n){ ll G = __gcd(a[i], 1LL * k); each(dp[i-1],p){ ll uoc = p.first; ll g = __gcd(G * uoc, 1LL * k); // choose a[i] dp[i][g] += dp[i-1][uoc]; dp[i][g] %= MOD; // not choose a[i] dp[i][uoc] += dp[i-1][uoc]; dp[i][uoc] %= MOD; } } cout << dp[n][k] << nl; } int main(){ ios_base::sync_with_stdio(false); cin.tie(NULL); int t = 1; // cin >> t; while (t--){ solve(); } }