1 #include<algorithm>
 2 #include<iostream>
 3 #include<cstring>
 4 #include<cstdio>
 5 using namespace std;
 6 inline int read(){
 7     int sum(0);
 8     char ch(getchar());
 9     for(;ch<'0'||ch>'9';ch=getchar());
10     for(;ch>='0'&&ch<='9';sum=sum*10+(ch^48),ch=getchar());
11     return sum;
12 }
13 typedef long long L;
14 const L mod(1000000007);
15 L n,k;
16 L n1;
17 L cnt,luc[1031];
18 inline void dfs(int dep,int x){
19     if(x)
20         luc[++cnt]=x;
21     if(dep==10)
22         return;
23     dfs(dep+1,x*10+4);
24     dfs(dep+1,x*10+7);
25 }
26 L sum[1031];
27 L f[100001],c[100001];
28 inline L pw(L x,L p){
29     L ret(1);
30     while(p){
31         if(p&1)
32             ret=(ret*x)%mod;
33         x=(x*x)%mod;
34         p>>=1;
35     }
36     return ret;
37 }
38 int main(){
39     dfs(1,0);
40     sort(luc+1,luc+cnt+1);
41     n=read(),k=read();
42     n1=n;
43     for(int i=1;i<=n;i++){
44         int a(read());
45         int pos(lower_bound(luc+1,luc+cnt+1,a)-luc);
46         if(luc[pos]==a){
47             sum[pos]++;
48             if(sum[pos]==2)
49                 n1-=2;
50             if(sum[pos]>2)
51                 n1--;
52         }
53     }
54     c[0]=f[0]=1;
55     for(int i=1;i<=cnt;i++)
56         if(sum[i]>=2)
57             for(int j=i;j>0;j--)
58                 f[j]=(f[j]+f[j-1]*sum[i])%mod;
59     for(int i=1;i<=n1;i++)
60         c[i]=c[i-1]*(n1-i+1)%mod*pw(i,mod-2)%mod;
61     L ans(0);
62     for(int i=0;i<=k;i++)
63         ans=(ans+c[i]*f[k-i]%mod)%mod;
64     printf("%lld",ans);
65 }