Bivariate Lagrange interpolation at the Padua points: the ideal theory approach