a(n+1) is the smallest nonsquarefree number m such that Diophantine equation S(x,y) = (x+y) + (x-y) + (x*y) + (x/y) = m has exactly n solutions, for n >= 0 (A353282); example: a(4) = 36 and 36 is the smallest number m such that equation S(x,y) = m has exactly 3 solutions: (9,1), (8,2), (5,5). - Bernard Schott, Apr 13 2022
a(n) is the square of the smallest integer having exactly n divisors (see formula with proof). - Bernard Schott, Oct 01 2022
Proof: Suppose a(n) = Product p_i^(2*e_i), where the p_i are primes. Then the n square divisors are all of the form d = Product p_i^(2*k_i) with 0 <= k_i <= e_i. As a(n) = Product (p_i^e_i)^2 = (Product (p_i^e_i))^2, we get that sqrt(a(n)) = Product (p_i^e_i). This is the prime decomposition of sqrt(a(n)). As there is a bijection between prime factors p_i^(2*k_i) and (p_i^k_i), there is also bijection between square divisors of a(n) and divisors of sqrt(a(n)). We conclude that sqrt(a(n)) is the smallest integer that has exactly n divisors. (End)
PROG
(PARI) a(n) = my(k=1); while(sumdiv(k, d, issquare(d)) != n, k++); k; \\ Michel Marcus, Jul 15 2019