OFFSET
1,9
COMMENTS
That is a generalization of a problem proposed by French site Diophante in link.
Some results:
x and y satisfy (x+1)*(y+1) = 2^n + 1.
x and y are both even, so 2 <= x <= y < 2^n.
There is only one case such that x = y, it is for n = 3 with x = y = 2 (Examples).
a(n) = 0 iff 2^n+1 is Fermat prime (A019434), hence iff n = 1, 2, 4, 8, 16.
a(n) = 1 iff 2^n+1 is semiprime (n is in A092559).
LINKS
Diophante, A4921 - Zéro, une, mille et plus (in French).
FORMULA
For n<>3, the number of positive pairs solution (x,y) is a(n) = (tau(2^n+1) - 2)/2.
For n = 3, there is only one pair solution and a(3) = (tau(2^3+1) - 1)/2 = 1, with (x, y) = (2, 2).
a(n) = ceiling((tau(2^n+1) - 2)/2) = ceiling((A046798(n)-2)/2) is the general formula.
EXAMPLE
For n = 3, only (x=y=2) satisfy s = 2+2 = 4, p = 2*2 = 4 and s+p = 8 = 2^3, hence a(3) = 1.
For n = 6, only (x=4, y=12) satisfy s = 4+12 = 16, p = 4*12 = 48 and s+p = 64 = 2^6 hence a(6) = 1.
For n = 9, (2,170), (8,56), (18,26) are the 3 solutions, with 172+340=512=2^9, 64+448=512, 44+468=512, hence a(9) = 3.
For n = 10, (4, 204) and (24, 40) are the 2 solutions, with 208+816=1024=2^10 and 64+960=1024, hence a(10) = 2.
MAPLE
with(numtheory):
M := seq(ceil((tau(2^n+1)-2)/2), n=1..100);
MATHEMATICA
a[3] = 1; a[n_] := DivisorSigma[0, 2^n + 1]/2 - 1; Array[a, 80] (* Amiram Eldar, Oct 05 2021 *)
PROG
(PARI) a(n) = ceil((numdiv(2^n+1) - 2)/2); \\ Michel Marcus, Oct 11 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Bernard Schott, Oct 05 2021
STATUS
approved