n,k=map(int,input().split())
a=list(map(int,input().split()))
from math import*
from collections import*
mod=10**9+7
d=defaultdict(int)
d[1]=1
for i in a:
  d2=d.copy()
  for j in d2:
    m=gcd(k,i*j)
    d[m]=(d2[j]+d[m])%mod
print(d[k]%mod-(k==1))