OFFSET
1,1
COMMENTS
Property of the sequence:
The sequence is of the form a1, a1+a2, a2, a1+a2, a3, a3+a4, a4, a3+a4, a5, a5+a6, a6,... where {a1, a2, a3,...} = A239436 are the members of a pair (m,n) such that sigma(m) = sigma(n) = sigma(n+m), and where the members a1+a2, a3+a4, a5+a6,... are repeated two times.
LINKS
Michel Lagneau, Table of n, a(n) for n = 1..38
EXAMPLE
The pair (1288, 2773) is in the sequence because sigma(1288) = sigma(2773) = 2880 and sigma(2773-1288) = sigma(1485) = 2880.
MATHEMATICA
a[n1_, n2_] := (t = Table[{DivisorSigma[1, n], n}, {n, n1, n2}] // Sort; s = Select[Split[t, #1[[1]] == #2[[1]] &], Length[#] >= 2 &]; f[lst_] := Select[Table[{lst[[i]], lst[[j]]}, {i, 1, Length[lst] - 1}, {j, i + 1, Length[lst]}] // Flatten[#, 1] &, #[[1, 1]] == DivisorSigma[1, #[[1, 2]] - #[[2, 2]]] &]; Select[f /@ s, # != {} &]); Flatten[a[1, 10^5], 2][[All, 2]] (* Program from Jean-François Alcover, adapted for this sequence. See A239436 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Mar 31 2014
STATUS
approved