#include <bits/stdc++.h>
using namespace std;
int MOD = 1000000007;
bool lucky(int a) {
  while (a > 0) {
    if (a % 10 != 4 && a % 10 != 7) return false;
    a /= 10;
  }
  return true;
}
long long D[5][1100];
long long F[100047], I[100047];
long long umocni(long long a, long long b) {
  long long res = 1;
  while (b > 0) {
    if (b % 2 == 1) res = (res * a) % MOD;
    a *= a;
    a %= MOD;
    b /= 2;
  }
  return res;
}
void precompute() {
  F[0] = 1;
  for (int i = 1; i < 100042; i++) F[i] = (F[i - 1] * i) % MOD;
  for (int i = 0; i < (100042); i++) I[i] = umocni(F[i], MOD - 2);
}
long long comb(int n, int k) {
  if (k > n) return 0;
  long long res = F[n] * I[k];
  res %= MOD;
  res *= I[n - k];
  return res % MOD;
}
int main() {
  precompute();
  int n, k;
  scanf("%d %d", &n, &k);
  int un = 0;
  map<int, int> M;
  for (int i = 0; i < (n); i++) {
    int x;
    scanf(" %d", &x);
    if (lucky(x))
      M[x]++;
    else
      un++;
  }
  vector<int> A;
  for (map<int, int>::iterator it = M.begin(); it != M.end(); it++)
    A.push_back((*it).second);
  D[0][0] = D[1][0] = 1;
  for (int i = 1; i <= A.size(); i++)
    for (int j = 1; j <= i; j++) {
      D[i % 2][j] = D[(i + 1) % 2][j];
      D[i % 2][j] += D[(i + 1) % 2][j - 1] * (long long)A[i - 1];
      D[i % 2][j] %= MOD;
    }
  long long res = 0;
  for (int i = 0; i < (min((int)A.size() + 1, k + 1)); i++) {
    res += D[A.size() % 2][i] * comb(un, k - i);
    res %= MOD;
  }
  cout << res << endl;
  return 0;
}
