OFFSET
1,1
COMMENTS
Is the sequence infinite?
a(11) > 5 * 10^13. - Hiroaki Yamanouchi, Aug 24 2015
PROG
(Python)
TOP = 5000*5000
a = [0]*TOP
for y in range(2, 5000):
z = TOP//y + 1
for x in range(y, z):
n = x*y + x + y
if n>=TOP: break
a[n] += x+y
for n in range(1, TOP):
if a[n]==a[n-1]+1 and a[n]>0: print(n, end=', ')
(PARI) a259124(n)=sum(y=2, sqrtint(n+1)-1, my(x=(n-y)/(y+1)); if(denominator(x)==1, x+y));
isok(n)= a259124(n) == a259124(n-1) + 1; \\ Michel Marcus, Aug 05 2015
CROSSREFS
KEYWORD
nonn,more,changed
AUTHOR
Alex Ratushnyak, Jul 15 2015
EXTENSIONS
a(8)-a(9) from Hiroaki Yamanouchi, Aug 24 2015
a(10) from Hiroaki Yamanouchi, Aug 25 2015
STATUS
approved