|
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Is the sequence infinite?
|
|
LINKS
|
|
|
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 str(n)+', ',
(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
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|