La méthode que j'ai utilisé n'est probablement pas la plus simple :
Soient i appartenant à 〚0,n〛, X et Y deux parties de E telles que card(X∩Y)=i.
Soit j appartenant à 〚0,n-i〛tel que card(X)=i+j.
Il y a j parmi n-i manières de choisir ces j éléments.
Lorsque les j éléments de X-Y sont déterminés, on a une liberté de k éléments de Y qui ne sont pas dans X (entre 0 et n-i-j donc) parmi n-i-j. On les somme tous pour k allant 0 à n-i-j et on le multiplie par la quantité précédente (j parmi n-i).
En sommant alors cette nouvelle quantité pour j allant de 0 à n-i, on obtient le nombre de manière de choisir les j éléments de X et les k éléments de Y qui ne sont pas dans X∩Y.
Il suffit alors de multiplier ceci par i parmi n pour avoir le nombre de possibilité d'avoir une intersection à i éléments. En remultipliant par i, on obtient la somme des cardinaux des intersections à i éléments, et il suffit de sommer cette quantité pour i allant de 0 à n.
Ensuite, c'est plein de binômes de Newton qui apparaissent, et ça se simplifie... mais j'avoue que l'explication est pas évidente à l'écrit 