En fait l'exo repose uniquement sur la multilinéarité du déterminant, rien de plus. Le truc bien dans l'affaire, c'est que XtX a un très petit rang, ce qui va annuler beaucoup de termes dans le développement du dét (à la barbare, comme un binôme de Newton).
Appelons C1,..,Cn les colonnes de A, posons X = (x1,..,xn).
On a : det(A+XtX) = det(C1+x1*X,..,Cn+xn*X) = det(C1,..,Cn) = somme(i=1, n, det(C1,..,xi*X,..,Cn) ) + (termes comprenant au moins deux colonnes de X).
Mais si i et j sont deux indices, dès qu'on a deux colonnes de X dans le dét, on a :
det(C1,..,xi*X,..,xj*X,..,Cn) = 0, car xi*(xj*X) - xj*(xi*X) = 0 !
Il ne reste donc que det(A) + somme(i=1, n, xi*det(C1,..,X,..,Cn) ),
le terme dans la somme n'étant qu'un déterminant de Cramer relativement à X.
Fini. 