OFFSET
1,1
COMMENTS
The sequence gives the index of the first element of a pair of consecutive triangular numbers with the same sum of divisors.
The even terms of this sequence are exactly twice the terms of A002961. - Amiram Eldar, Sep 16 2024
EXAMPLE
Sum of the numbers from 1 to 28 -> 406; from 1 to 29 -> 435;
sigma(406) = sigma(435) = 720.
MAPLE
with(numtheory): P:=proc(q) local a, j, n, t; j:=[]; t:=0;
for n from 1 to q do a:=n*(n+1)/2; if sigma(a)=t then j:=[op(j), n-1]; fi;
t:=sigma(a); od; op(j); end: P(143370);
MATHEMATICA
Select[Range[145000], DivisorSigma[1, #(#+1)/2]==DivisorSigma[1, (#+1)(#+2)/2] &] (* Stefano Spezia, Sep 16 2024 *)
PROG
(Python)
from itertools import count, islice
from sympy import divisor_sigma
def A375819_gen(): # generator of terms
s, a = 1, 1
return (k-1 for k in count(2) if a==(a:=divisor_sigma(s:=s+k)))
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paolo P. Lava, Sep 16 2024
STATUS
approved