Dunque, io invece facevo una versione leggermente diversa... Chiamiamo "e" il vettore di tutti uni di dimensione opportuna, e I la matrice identità di dimensione opportuna, F2 il campo con due elementi, R i reali (non ho voglia di fare i mathbb...).
Osservazione: ee^T-I è la matrice che ha tutti zeri sulla diagonale e tutti uni fuori.
Lemma 1: in uno spazio vettoriale F2^n di dimensione n pari, ee^T-I ha rango massimo. Dimostrazione: facciamo finta di essere in R^n per un attimo. Allora, gli autovettori di questa matrice sono e (con autovalore n-1) e una qualunque base dello spazio n-1 dimensionale ortogonale a e, ognuno con autovalore -1. Il determinante è il prodotto degli autovalori, quindi viene $ (n-1)(-1)^{n-1} $. Ora, fare il determinante in F2 è come farlo in Z (e quindi in R) e poi ridurlo modulo 2, quindi in F2 il determinante viene $ (n-1)(-1)^{n-1}=1 $ (perché n pari). []
Lemma 2: rango(BC)<=rango(C). Dimostrazione: in tanti modi... per esempio: il rango è la dimensione dell'immagine di una matrice; l'immagine di BC è l'immagine sotto B del sottospazio immagine di C, quindi se Im C ha dimensione al più n, applicare la matrice B potrà mandare a spasso quanto vuole questo sottospazio ma non potrà mai farlo crescere di dimensione. []
Dimostrazione del tuo problema nel caso n dispari: supponiamo per assurdo che ci siano più di n clubs, quindi almeno n+1. Prendiamo A=la matrice nx(n+1) che ha come i-esima colonna il "vettore di incidenza" dell'i-esimo club (cioè: numeriamo gli abitanti e i club, e poniamo a_{ij}=1 sse l'i-esimo abitante sta nel club j). Ora, dalle condizioni dell'ipotesi sappiamo che A^TA = ee^T-I (prova a farti i prodotti riga per colonna uno per uno e viene), dove il LHS e il RHS hanno dimensione n+1 (pari quindi). Qual è il rango del RHS? n+1 (per il lemma 1). Qual è il rango del LHS? Al più n (per il lemma 2, perché il rango di A è al più n).
Dimostrazione del tuo problema nel caso n pari (questa è nuova di questo post, spero che funzioni): supponiamo per assurdo che ci siano n clubs o più, e costruiamo la matrice A di dimensione nxn fatta accostando i vettori di incidenza dei primi n clubs (come sopra). e^TA=0, perché ogni colonna ha un numero pari di uni (dall'ipotesi). Quindi A non ha rango massimo. D'altra parte, come prima abbiamo A^TA=ee^T-I (entrambi i membri di dimensione n pari). Il RHS ha rango massimo per il lemma 1; il LHS non ha rango massimo per il lemma 2, quindi assurdo.
Spero che funzioni tutto...

Se hai altri dubbi o se trovi errori non farti problemi a riscrivere!