A280391 - OEIS

%I #18 Jul 29 2020 12:12:45

%S 1,12,25,57,81,141,169,259,297,413,441,621,625,825,873,1079,1089,1403,

%T 1369,1739,1729,2021,2025,2507,2433,2859,2905,3301,3249,4029,3721,

%U 4509,4305,4793,4989,5551,5329,6027,6025,6807,6561,7917,7225,8357,8121,8677,8649,9843,9481,10889

%N Number of 2 X 2 matrices with all elements in {0,...,n} with permanent = determinant * n.

%C All the values except a(1) are odd.

%C From _Robert Israel_, Jan 02 2017: (Start)

%C Number of solutions to (n+1)*x*y = (n-1)*z*w for x,y,z,w in [0..n].

%C a(n) >= (2n+1)^2, with equality if n+1 is an odd prime. (End)

%H Robert Israel and Indranil Ghosh, <a href="/A280391/b280391.txt">Table of n, a(n) for n = 0..1600</a> (n = 0..200 from Indranil Ghosh)

%p g:= proc(r,n) if r = 0 then 2*n+1 else nops(select(t -> t <= n and r <= t*n, numtheory:-divisors(r))) fi end proc:

%p f:= proc(n) local c;

%p     if n::even then (2*n+1)^2 + add(g((n+1)*c,n)*g((n-1)*c,n), c=1..n-1)

%p     else (2*n+1)^2 + add(g((n+1)/2*c,n) * g((n-1)/2*c,n), c=1..2*n-1)

%p     fi

%p end proc:

%p map(f, [$0..100]); # _Robert Israel_, Jan 02 2017

%t g[r_, n_] := If[r == 0, 2n + 1, Length[Select[Divisors[r], # <= n && r <= # n&]]];

%t f[n_] := If[EvenQ[n], (2n + 1)^2 + Sum[g[(n + 1)c, n] g[(n - 1)c, n], {c, 1, n - 1}], (2n + 1)^2 + Sum[g[(n + 1)/2 c, n] g[(n - 1)/2 c, n], {c, 1, 2n - 1}]];

%t f /@ Range[0, 100] (* _Jean-François Alcover_, Jul 29 2020, after _Robert Israel_ *)

%o (Python)

%o def t(n):

%o     s=0

%o     for a in range(0,n+1):

%o         for b in range(0,n+1):

%o             for c in range(0,n+1):

%o                 for d in range(0,n+1):

%o                     if (a*d-b*c)*n==(a*d+b*c):

%o                         s+=1

%o     return s

%o for i in range(0,201):

%o     print str(i)+" "+str(t(i))

%Y Cf. A280321 (Number of 2 X 2 matrices with all elements in {0,..,n} with permanent*n = determinant).

%Y Cf. A015237 (Number of 2 X 2 matrices having all elements in {0..n} with determinant = permanent).

%Y Cf. A016754 (Number of 2 X 2 matrices having all elements in {0..n} with determinant =2* permanent).

%Y Cf. A280364 (Number of 2 X 2 matrices having all elements in {0..n} with determinant^n = permanent).

%K nonn

%O 0,2

%A _Indranil Ghosh_, Jan 02 2017