import itertools

def all_relations(n):
    for I in itertools.product((0,1),repeat=n**2):
        yield [I[n*k:n*(k+1)] for k in range(n)]

def is_symmetric(rel):
    n = len(rel)
    return all(rel[i][j] == rel[j][i] for i in range(n) for j in range(n))

print(len(list(filter(is_symmetric, all_relations(3)))))