#include <bits/stdc++.h>
using namespace std;
#define ldb double

const int N=200050;
const int M=17;
char s[N];
int cnt[1<<M][M+1];
int q[N],m;

ldb ans=0;
ldb dp[1<<M];
bool was[1<<M];
int in[1<<M];
void Solve(int mask,int lvl){
	if(was[mask])return;
	was[mask]=true;
	cnt[0][lvl]=0;
	cnt[mask][lvl]=0;
	int tot=0;
	for(int i=mask;i;i=(i-1)&mask){
		tot+=cnt[i][lvl];
	}
	if(tot==0){
		if(mask&1){
			ans+=dp[mask];
		}
	}else{
		for(int i=(mask-1)&mask;i;i=(i-1)&mask){
			dp[i]+=dp[mask]*(ldb)cnt[i][lvl]/tot;
		}
		for(int i=0;i<m;i++){
			if(mask>>i&1){
				int nxt=mask^(1<<i);
				in[nxt]--;
				if(in[nxt]==0){
					for(int j=mask;j;j=(j-1)&mask){
						cnt[j][lvl+1]=0;
					}
					for(int j=mask;j;j=(j-1)&mask){
						int tmp=j;
						if(j>>i&1)tmp-=1<<i;
						cnt[tmp][lvl+1]+=cnt[j][lvl];
					}
					Solve(nxt,lvl+1);
				}
			}
		}
	}
}
int main(){
	int n;
	scanf("%i %i",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%s",s);
		for(int j=0;j<m;j++){
			if(s[j]=='1'){
				q[i]+=1<<j;
			}
		}
	}
	for(int i=1;i<=n;i++){
		cnt[q[i]][0]++;
	}
	for(int mask=0;mask<(1<<m);mask++){
		for(int i=0;i<m;i++){
			if(mask>>i&1){
				in[mask^(1<<i)]++;
			}
		}
		/*for(int sub=(mask-1)&mask;sub;sub=(sub-1)&mask){
			in[sub]++;
		}*/
	}
	dp[(1<<m)-1]=1;
	Solve((1<<m)-1,0);
	printf("%.12lf\n",ans);
	return 0;
}