OFFSET
1,2
COMMENTS
f is the composite (P o D), where P(x,y)=((x+y)^2+3*x+y)/2 is the Cantor polynomial and D(x,y)=((x+1)*(y-1),(x-1)*(y+1)) is a divisor plot built so as to fit the first quadrant. This sequence can be viewed as an irregular table where the length of row n is A000005(n), the number of divisors of n.
LINKS
Luc Rousseau, Illustration
FORMULA
EXAMPLE
1: 0
. .
2: 3 . 5
. . . .
3: 10 . . . 14
. . . . . .
4: 21 . . 24 . . 27
. . . . . . . .
5: 36 . . . . . . . 44
. . . . . . . . . .
6: 55 . . . 59 . 61 . . . 65
. . . . . . . . . . . .
7: 78 . . . . . . . . . . . 90
. . . . . . . . . . . . . .
8: 105 . . . . 110 . . . 114 . . . . 119
. . . . . . . . . . . . . . . .
9: 136 . . . . . . . 144 . . . . . . . 152
. . . . . . . . . . . . . . . . . .
(...)
This sequence is what remains when one removes the dots or "unoccupied integers" from the above schema, result of the superposition of a divisor plot on a triangle of numbered points.
PROG
(PARI)
f(x, y)=2*x*y*(x*y-1)-x+y
for(n=1, 20, fordiv(n, d, print1(f(n/d, d), ", ")))
CROSSREFS
KEYWORD
nonn,tabf
AUTHOR
Luc Rousseau, Jun 25 2019
STATUS
approved