A265751 -id:A265751 - cp's OEIS frontend

A266116 The last nonzero term on each row of A265751.

7, 7, 13, 7, 8, 7, 13, 7, 8, 13, 20, 13, 25, 13, 20, 19, 24, 19, 25, 19, 20, 37, 25, 37, 24, 25, 40, 37, 28, 37, 50, 37, 40, 33, 50, 37, 36, 37, 50, 43, 40, 43, 49, 43, 50, 67, 49, 67, 56, 49, 50, 67, 52, 67, 68, 55, 56, 67, 68, 67, 136, 67, 68, 63, 64, 67, 66, 67, 68, 79, 74, 79, 136, 79, 74, 75, 103, 79, 98, 79, 88, 103, 98, 103, 136, 85

Offset: 0

Antti Karttunen, Dec 21 2015

nonn

Starting from j = n, search for a smallest number k such that k - d(k) = j, and if found such a number, replace j with k and repeat the procedure. When eventually such k is no longer found, then the (last such) j must be one of the terms of A045765, and it is set as the value of a(n).

			Starting from n = 21, we get the following chain: 21 -> 23 -> 27 -> 29 -> 31 -> 35 -> 37, with A082284 iterated 6 times before the final term 37 (for which A060990(37) = A082284(37) = 0) is encountered. Thus a(21) = 37.

Antti Karttunen, Table of n, a(n) for n = 0..124340

Cf. A000005, A045765, A060990, A082284, A265751.
Cf. A266110 (gives the number of iterations of A082284 needed before a(n) is found).
Cf. also tree A263267 (and its illustration).

(definec (A266116 n) (cond ((A082284 n) => (lambda (lad) (if (zero? lad) n (A266116 lad))))))
;; Alternatively:
(define (A266116 n) (A265751bi n (A266110 n))) ;; Code for A265751bi given in A265751.

a(n) = A265751(n, A266110(n)).
If A060990(n) = 0, a(n) = n, otherwise a(n) = a(A082284(n)), where A082284(n) = smallest number k such that k - d(k) = n, or 0 if no such number exists, and d(n) = the number of divisors of n (A000005).
Other identities and observations. For all n >= 0:
a(n) >= n.
A060990(a(n)) = 0. [All terms are in A045765.]

A082284 a(n) = smallest number k such that k - tau(k) = n, or 0 if no such number exists, where tau(n) = the number of divisors of n (A000005).

Original entry on oeis.org

1, 3, 6, 5, 8, 7, 9, 0, 0, 11, 14, 13, 18, 0, 20, 17, 24, 19, 22, 0, 0, 23, 25, 27, 0, 0, 32, 29, 0, 31, 34, 35, 40, 0, 38, 37, 0, 0, 44, 41, 0, 43, 46, 0, 50, 47, 49, 51, 56, 0, 0, 53, 0, 57, 58, 0, 0, 59, 62, 61, 72, 65, 68, 0, 0, 67, 0, 0, 0, 71, 74, 73, 84, 77, 0, 0, 81, 79, 82, 0, 88

Offset: 0

Amarnath Murthy, Apr 14 2003

nonn

a(p-2) = p for odd primes p.

Antti Karttunen, Table of n, a(n) for n = 0..124340

Column 1 of A265751.
Cf. A000005, A002182, A002183, A049820, A060990, A261100.
Cf. A262686 (the largest such number), A262511 (positions where these are equal and nonzero).
Cf. A266114 (same sequence sorted into ascending order, with zeros removed).
Cf. A266115 (positive numbers missing from this sequence).
Cf. A266110 (number of iterations before zero is reached), A266116 (final nonzero value reached).
Cf. also tree A263267 and its illustration.

N:= 1000: # to get a(0) .. a(N)
V:= Array(0..N):
for k from 1 to 2*(N+1) do
  v:= k - numtheory:-tau(k);
  if v <= N and V[v] = 0 then V[v]:= k fi
od:
seq(V[n],n=0..N); # Robert Israel, Dec 21 2015

Table[k = 1; While[k - DivisorSigma[0, k] != n && k <= 2 (n + 1), k++]; If[k > 2 (n + 1), 0, k], {n, 0, 80}] (* Michael De Vlieger, Dec 22 2015 *)

allocatemem(123456789);
uplim1 = 2162160 + 320; \\ = A002182(41) + A002183(41).
uplim2 = 2162160;
v082284 = vector(uplim1);
A082284 = n -> if(!n,1,v082284[n]);
for(n=1, uplim1, k = n-numdiv(n); if((0 == A082284(k)), v082284[k] = n));
for(n=0, 124340, write("b082284.txt", n, " ", A082284(n)));
\\ Antti Karttunen, Dec 21 2015

(define (A082284 n) (if (zero? n) 1 (let ((u (+ n (A002183 (+ 2 (A261100 n)))))) (let loop ((k n)) (cond ((= (A049820 k) n) k) ((> k u) 0) (else (loop (+ 1 k))))))))
;; Antti Karttunen, Dec 21 2015

Other identities and observations. For all n >= 0:

a(n) <= A262686(n).

More terms from David Wasserman, Aug 31 2004

A263267 Breadth-first traversal of the tree defined by the edge-relation A049820(child) = parent.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 5, 8, 9, 10, 12, 7, 11, 14, 18, 13, 15, 16, 20, 22, 17, 24, 25, 26, 28, 30, 19, 21, 32, 34, 23, 40, 38, 42, 27, 44, 48, 46, 29, 36, 50, 56, 60, 49, 52, 54, 31, 33, 72, 58, 35, 84, 62, 66, 37, 39, 96, 68, 70, 41, 45, 104, 108, 74, 76, 78, 80, 43, 47, 120, 81, 82, 90, 88, 51, 128, 132, 83, 85, 86, 94, 53, 55, 136, 140, 87, 92, 102

Offset: 0

Antti Karttunen, Nov 27 2015

It is conjectured that the terms of A259934 trace the only infinite path in this tree.

After the root (0), the tree narrows next time to the width of just one node at level A262508(1) = 9236, with vertex 119143.

			Rows 0 - 21 of the table. The lines show the nodes of the tree connected by the edge-relation A049820(child) = parent:
0;
| \
1, 2;
| \  \
3, 4, 6;____
|  |  | \   \
5, 8, 9, 10, 12;
|     |   |   |
7, _ 11, 14, 18;
  /  | \   \   \
13, 15, 16, 20, 22;____
     |  |      / | \   \
    17, 24, 25, 26, 28, 30;
     | \         |      |
    19, 21,     32,     34;
         |       |      | \
        23,     40,    38, 42;____
         |              | \       \
        27,            44, 48,     46;____
         | \            |   | \    |  \   \
        29, 36,        50, 56, 60, 49, 52, 54;
         | \                   |           |
        31, 33,                72,         58;
         |                     |           |  \
        35,                    84,         62, 66;
         | \                   |           |  \
        37, 39,                96,         68, 70;_______
            |  \               |  \           / |  \     \
            41, 45,           104, 108,     74, 76, 78,   80;
            |   |              |                |   |  \    \
            43, 47,           120,             _81, 82, 90, 88;
                |              |  \           / |   |   |
                51,           128, 132,     83, 85, 86, 94;
                 | \            | \          |       |   |
                53, 55        136, 140      87,     92, 102;______
                 |                           | \     |    |  \    \
                57,_                        89, 91, 98, 106,  110, 112;
               / |  \                       /   / \       |     |
             59, 63, 64,                  93, 95, 100,   114,   116;
              |                            |   |          |  \
             61,                          99, 97,       _118, 126;
              |                            |   |       /  |  \
             65,                         101, 105,  121, 122, 124;
(See also _Michael De Vlieger_'s poster in the Links section.)

Antti Karttunen, Table of n, a(n) for n = 0..10425; levels 0 .. 1001 of the tree
Michael De Vlieger, Poster illustrating A259934 and A263267
Index entries for sequences that are permutations of the natural numbers

Inverse permutation: A263268.
Cf. A000005, A049820, A060990, A082284, A155043, A259934, A262508.
Cf. A262507 (number of terms on row/level n), A263260 (total number of terms in levels 0 .. n).
Cf. A264988 (the left edge), this differs from A261089 (the least term on each level) for the first time at level 69.
Cf. A263269 (the right edge).
Cf. A262686 (maximum term on the level n).
Cf. A045765 (the leaves of the tree).
Cf. also permutations A263265 (obtained from this table by sorting each row into ascending order), A263266.
Cf. also arrays A265751 and A263271.
Differs from A263265 for the first time at n=31, where a(31) = 40, while A263265(31) = 38.
Cf. also A088975.

uplim = 125753; \\ = A263260(10001).
checklimit = 1440; \\ Hard limit 1440 good for at least up to A002182(67) = 1102701600 as A002183(67) = 1440.
v263267 = vector(uplim);
A263267 = n -> if(!n,n,v263267[n]);
z = 0; for(n=0, uplim, t = A263267(n); write("b263267.txt", n, " ", t); for(k=t+1, t+checklimit, if((k-numdiv(k)) == t, z++; if(z <= uplim, v263267[z] = k))));

# After David Eppstein's Python-code for A088975.
def A263267():
  '''Breadth-first reading of irregular tree defined by the edge-relation A049820(child) = parent'''
  yield 0
  for x in A263267():
    for k in [x+1 .. 2*(x+1)]:
      if ((k - sloane.A000005(k)) == x): yield k
def take(n,g):
  '''Returns a list composed of the next n elements returned by generator g.'''
  return [next(g) for _ in range(n)]
take(120, A263267())

;; This version creates the list of terms incrementally, using append! function that physically modifies the list at the same time as it is traversed. Otherwise the idea is essentially the same as with Python/Sage-program above:
(define (A263267list_up_to_n_terms_at_least n) (let ((terms-produced (list 0))) (let loop ((startp terms-produced) (endp terms-produced) (k (- n 1))) (cond ((<= k 0) terms-produced) (else (let ((children (children-of-n-in-A049820-tree (car startp)))) (cond ((null? children) (loop (cdr startp) endp k)) (else (begin (append! endp children) (loop (cdr startp) children (- k (length children))))))))))))
(define (children-of-n-in-A049820-tree n) (let loop ((k (A262686 n)) (children (list))) (cond ((<= k n) children) ((= (A049820 k) n) (loop (- k 1) (cons k children))) (else (loop (- k 1) children)))))

A263271 Square array A(row,col): A(row,0) = row and for col >= 1, if A262686(row) is 0, then A(row,col) = 0, otherwise A(row,col) = A(A262686(row),col-1).

Original entry on oeis.org

0, 2, 1, 6, 4, 2, 12, 8, 6, 3, 18, 0, 12, 5, 4, 22, 0, 18, 7, 8, 5, 30, 0, 22, 0, 0, 7, 6, 34, 0, 30, 0, 0, 0, 12, 7, 42, 0, 34, 0, 0, 0, 18, 0, 8, 46, 0, 42, 0, 0, 0, 22, 0, 0, 9, 54, 0, 46, 0, 0, 0, 30, 0, 0, 11, 10, 58, 0, 54, 0, 0, 0, 34, 0, 0, 16, 14, 11

Offset: 0

Antti Karttunen, Nov 29 2015

The square array A(row>=0, col>=0) is read by downwards antidiagonals as: A(0,0), A(0,1), A(1,0), A(0,2), A(1,1), A(2,0), A(0,3), A(1,2), A(2,1), A(3,0), ...

Each row n lists all the nodes in A263267-tree that one encounters when one starts from node n and always chooses the largest possible child of it (A262686), and then the largest possible child of that child, etc, until a leaf-child (one of the terms of A045765) is encountered, after which the rest of the row contains only zeros.

			The top left corner of the array:
   0,  2,  6, 12, 18, 22, 30, 34, 42, 46, 54, 58, 66,  0,  0,  0,  0
   1,  4,  8,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
   2,  6, 12, 18, 22, 30, 34, 42, 46, 54, 58, 66,  0,  0,  0,  0,  0
   3,  5,  7,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
   4,  8,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
   5,  7,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
   6, 12, 18, 22, 30, 34, 42, 46, 54, 58, 66,  0,  0,  0,  0,  0,  0
   7,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
   8,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
   9, 11, 16, 24,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
  10, 14, 20,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
  11, 16, 24,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
  12, 18, 22, 30, 34, 42, 46, 54, 58, 66,  0,  0,  0,  0,  0,  0,  0
  13,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
  14, 20,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
  15, 17, 21, 23, 27, 36,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0
  ...

Antti Karttunen, Table of n, a(n) for n = 0..10439; the first 144 antidiagonals

Column 0: A001477, Column 1: A262686.
Cf. A264971 (number of significant terms on each row, position where the first trailing zero-term occurs).
Cf. A264970.
Cf. also A000005, A045765, A263267.
See also array A265751 constructed in the same way, but obtained by following always the smallest child A082284, instead of the largest child A262686.

(define (A263271 n) (A263271bi (A002262 n) (A025581 n)))
(define (A263271bi row col) (cond ((zero? col) row) ((A262686 row) => (lambda (lad) (if (zero? lad) lad (A263271bi lad (- col 1)))))))
;; An alternative implementation, reflecting the new recurrence:
(define (A263271bi row col) (cond ((zero? col) row) ((and (zero? row) (= 1 col)) 2) ((zero? (A263271bi row (- col 1))) 0) (else (A262686 (A263271bi row (- col 1))))))

A(row,0) = row and for col >= 1, if A262686(row) is 0, then A(row,col) = 0, otherwise A(row,col) = A(A262686(row),col-1).

A(0,0) = 0, A(0,1) = 2; if col = 0, A(row,0) = row; and for col > 0, if A(row,col-1) = 0, then A(row,col) = 0, otherwise A(row,col) = A262686(A(row,col-1)). [Another, perhaps more intuitive recurrence for this array.] - Antti Karttunen, Dec 21 2015

A266111 If A082284(n) = 0, a(n) = 1, otherwise a(n) = 1 + a(A082284(n)), where A082284(n) = smallest number k such that k - d(k) = n, or 0 if no such number exists, and d(n) = the number of divisors of n (A000005).

Original entry on oeis.org

5, 4, 5, 3, 2, 2, 4, 1, 1, 3, 3, 2, 4, 1, 2, 3, 2, 2, 3, 1, 1, 7, 2, 6, 1, 1, 3, 5, 1, 4, 5, 3, 2, 1, 4, 2, 1, 1, 3, 3, 1, 2, 3, 1, 2, 9, 2, 8, 2, 1, 1, 7, 1, 6, 4, 1, 1, 5, 3, 4, 8, 3, 2, 1, 1, 2, 1, 1, 1, 5, 2, 4, 7, 3, 1, 1, 9, 2, 5, 1, 2, 8, 4, 7, 6, 1, 3, 6, 1, 5, 13, 6, 2, 4, 12, 5, 5, 4, 1, 3, 1, 2, 11, 1, 4, 3, 10, 2, 1, 1, 2, 2, 1, 1, 9, 3, 1, 1, 8, 2, 3, 7

Offset: 0

Antti Karttunen, Dec 21 2015

nonn

			Starting from n = 21, we get the following chain: 21 -> 23 -> 27 -> 29 -> 31 -> 35 -> 37, with A082284 iterated 6 times before the final term 37 (for which A060990(37) = A082284(37) = 0) is encountered. Here we count the terms (not steps) in whole chain, thus a(21) = 7.

Antti Karttunen, Table of n, a(n) for n = 0..10000

One more than A266110.
Number of significant terms on row n of A265751 (without its trailing zeros).
Cf. tree A263267 (and its illustration).
Cf. A000005, A060990, A082284.
Cf. also A264971.

If A060990(n) = 0, a(n) = 1, otherwise a(n) = 1 + a(A082284(n)).
Other identities. For all n >= 0:
a(n) = 1 + A266110(n).

cp's OEIS Frontend

A266116 The last nonzero term on each row of A265751.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Scheme

Formula

A082284 a(n) = smallest number k such that k - tau(k) = n, or 0 if no such number exists, where tau(n) = the number of divisors of n (A000005).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Scheme

Formula

Extensions

A263267 Breadth-first traversal of the tree defined by the edge-relation A049820(child) = parent.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Sage

Scheme

A263271 Square array A(row,col): A(row,0) = row and for col >= 1, if A262686(row) is 0, then A(row,col) = 0, otherwise A(row,col) = A(A262686(row),col-1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Scheme

Formula

A266111 If A082284(n) = 0, a(n) = 1, otherwise a(n) = 1 + a(A082284(n)), where A082284(n) = smallest number k such that k - d(k) = n, or 0 if no such number exists, and d(n) = the number of divisors of n (A000005).

Views

Author

Keywords

Examples

Links

Crossrefs

Formula