bitarray.jl

function solve(n, m, ts, xss, ss)
    mat = BitArray(undef, n, m)
    for i in 1:n
        for j in 1:ts[i]
            mat[i, xss[i][j]] = 1
        end
    end
    println(mat)
    s = BitArray(ss)
    c = 1
    for i in 1:m
        for j in c:n
            if mat[j, i] == 1
                if s[i] == 1
                    s .⊻= mat[j, :]
                end
                mat[j, :], mat[c, :] = mat[c, :], mat[j, :]
                c += 1
                break
            end
        end
    end
end

function main()
    n, m, ts, xss, ss = (9, 8, [4, 5, 8, 5, 1, 4, 8, 2, 2], [[1, 2, 3, 4], [1, 5, 6, 7, 8], [1, 2, 3, 4, 5, 6, 7, 8], [1, 2, 3, 5, 7], [4], [1, 5, 6, 8], [1, 2, 3, 4, 5, 6, 7, 8], [2, 4], [1, 8]], Bool[1, 1, 0, 1, 0, 0, 0, 1])
    println(solve(n, m, ts, xss, ss))
end

main()