OFFSET
1,2
COMMENTS
See A243917 for definition of non-twin divisor.
Squares included in the sequence are : 1, 4, 9, 16, 64, 729, ...
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..8692 (a(n) < 5,000,000)
EXAMPLE
The positive divisors of 8 are 1, 2, 4, 8. Of these, 1 and 8 are non-twin divisors. So 8 is in this sequence, which is 1 + 8 = 3^2.
MATHEMATICA
a243984[n_Integer] := Total[Select[Divisors[n], If[And[# <= 2 || Divisible[n, # - 2] == False, Divisible[n, # + 2] == False], True, False] &]]; a243985[n_Integer] := Flatten@Select[Position[Sqrt[a243984 /@ Range[n]], _?IntegerQ], If[Length[#] == 1, True, False] &]; a243985[1500] (* Michael De Vlieger, Aug 17 2014 *)
PROG
(PARI)
A243984(n) = s=0; fordiv(n, d, if(!((d>2 && n%(d-2)==0) || (d<=n-2 && n%(d+2)==0)), s+=d)); s
for(n=1, 200, if(issquare(A243984(n)), print1(n, ", "))) \\ Colin Barker, Jun 29 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Jun 16 2014
EXTENSIONS
Several terms corrected by Colin Barker, Jun 29 2014
STATUS
approved