import sys
import numpy as np

input = sys.stdin.buffer.readline

N, M = map(int, input().split())
S = np.array(tuple(tuple(map(int, input().split())) for _ in range(N)), dtype=bool)

dp = [0] * (1 << M)
alive = np.zeros((1 << M, N), dtype=bool)
dp[0] = 1
alive[0] = True

for i in range(M):
    for b in range(1 << i):
        alive[b | (1 << i)] = alive[b] & S[:, i]
alive_count = np.count_nonzero(alive, axis=1)

for b in range(1 << M):
    for i in range(M):
        if (b >> i) & 1:
            continue
        c = b | (1 << i)
        if alive_count[c] >= alive_count[b] - alive_count[c]:
            dp[c] += dp[b]
print(dp[-1])