#include<bits/stdc++.h>
using namespace std;
double dp[1<<17];
int buf[1<<17];
int main(){
    int n,m;
    scanf("%d %d",&n,&m);
    for(int i=0;i<n;i++)
    {
        int sta=0;
        for(int j=0;j<m;j++)
        {
            char ch;
            cin>>ch;
            if(ch=='1')
                sta|=(1<<j);
        }
        buf[sta]++;
    }
    int full=(1<<m)-1;
    dp[full]=1;
    double ans=0;
    int sta=full;
    auto DP=[&](auto self,int id)->void
    {
        if(!id)
        {
            int cnt=0;
            for(int g=((sta-1)&sta);g;g=(g-1)&sta)
                cnt+=buf[g];
            if(!cnt)
            {
                ans+=dp[sta];
                return;
            }
            for(int g=((sta-1)&sta);g;g=(g-1)&sta)
                dp[g]+=dp[sta]*buf[g]/cnt;
            return;
        }
        self(self,id-1);
        sta^=(1<<id);
        for(int g=sta;;g=(g-1)&sta)
        {
            buf[g]+=buf[g|(1<<id)];
            if(!g)
                break;
        }
        self(self,id-1);
        for(int g=sta;;g=(g-1)&sta)
        {
            buf[g]-=buf[g|(1<<id)];
            if(!g)
                break;
        }
        sta^=(1<<id);
    };
    DP(DP,m-1);
    printf("%.12lf\n",ans);
    
    return 0;
}