a = sps.csc_matrix( (z , (x,y)), shape = (N,N), dtype = int)

count = 0
for i in range(N):
    count += (z == a[:,j]).sum()

a == z