cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 11-13 of 13 results.

A332893 a(1) = 1, a(2n) = n, a(2n+1) = A332819(2n+1).

Original entry on oeis.org

1, 1, 2, 2, 3, 3, 5, 4, 4, 5, 13, 6, 7, 7, 6, 8, 11, 9, 17, 10, 10, 11, 29, 12, 9, 13, 8, 14, 19, 15, 37, 16, 26, 17, 15, 18, 23, 19, 14, 20, 31, 21, 41, 22, 12, 23, 53, 24, 25, 25, 22, 26, 43, 27, 39, 28, 34, 29, 61, 30, 47, 31, 20, 32, 21, 33, 73, 34, 58, 35, 89, 36, 59, 37, 18, 38, 65, 39, 97, 40, 16, 41, 101, 42, 33, 43, 38, 44, 67, 45, 35
Offset: 1

Views

Author

Antti Karttunen, Mar 01 2020

Keywords

Comments

For any node n >= 2 in binary trees like A332815, a(n) gives the parent node of n.

Crossrefs

Cf. also A252463.

Programs

Formula

a(1) = 1, after which a(n) = n/2 for even n, and a(n) = A332819(n) for odd n.

A332899 a(1) = 0, and for n > 2, a(n) = a(A332893(n)) + A000035(n).

Original entry on oeis.org

0, 1, 2, 1, 3, 2, 4, 1, 2, 3, 6, 2, 5, 4, 3, 1, 7, 2, 8, 3, 4, 6, 10, 2, 3, 5, 2, 4, 9, 3, 12, 1, 6, 7, 4, 2, 11, 8, 5, 3, 13, 4, 14, 6, 3, 10, 16, 2, 4, 3, 7, 5, 15, 2, 6, 4, 8, 9, 18, 3, 17, 12, 4, 1, 5, 6, 20, 7, 10, 4, 22, 2, 19, 11, 3, 8, 6, 5, 24, 3, 2, 13, 26, 4, 7, 14, 9, 6, 21, 3, 5, 10, 12, 16, 8, 2, 23, 4, 6, 3, 25, 7, 28, 5, 4
Offset: 1

Views

Author

Antti Karttunen, Mar 04 2020

Keywords

Comments

a(n) tells how many odd numbers are encountered when map x -> A332893(x) is used to traverse from n to 1, the root of the binary tree A332815. This count includes both the starting n itself if it is odd, but excludes 1 where the iteration ends.
a(n) also gives the index of the largest prime factor (A061395) in A332808(n), which is the inverse permutation of A108548 (see also A108546).

Crossrefs

Cf. A000079 (after its initial term, gives the positions of 1's).

Programs

Formula

a(1) = 0, and for n > 1, a(n) = a(A332893(n)) + A000035(n).
a(n) = A000120(A332811(n)).
a(n) = A061395(A332808(n)).
a(n) = A332897(n) + A332898(n).
a(n) <= A332894(n).
For all n > 1, a(n) = 1 + A080791(A332816(n)).

A340323 Multiplicative with a(p^e) = (p + 1) * (p - 1)^(e - 1).

Original entry on oeis.org

1, 3, 4, 3, 6, 12, 8, 3, 8, 18, 12, 12, 14, 24, 24, 3, 18, 24, 20, 18, 32, 36, 24, 12, 24, 42, 16, 24, 30, 72, 32, 3, 48, 54, 48, 24, 38, 60, 56, 18, 42, 96, 44, 36, 48, 72, 48, 12, 48, 72, 72, 42, 54, 48, 72, 24, 80, 90, 60, 72, 62, 96, 64, 3, 84, 144, 68, 54
Offset: 1

Views

Author

Keywords

Comments

Starting with any integer and repeatedly applying the map x -> a(x) reaches the fixed point 12 or the loop {3, 4}.

Examples

			a(2^s) = 3 for all s>0.
		

Crossrefs

Programs

  • Maple
    f:= proc(n) local  t;
      mul((t[1]+1)*(t[1]-1)^(t[2]-1),t=ifactors(n)[2])
    end proc:
    map(f, [$1..100]); # Robert Israel, Jan 07 2021
  • Mathematica
    fa[n_]:=fa[n]=FactorInteger[n];
    phi[1]=1; phi[p_, s_]:= (p + 1)*( p - 1)^(s - 1)
    phi[n_]:=Product[phi[fa[n][[i, 1]], fa[n][[i, 2]]], {i,Length[fa[n]]}];
    Array[phi, 245]
  • PARI
    A340323(n) = if(1==n,n,my(f=factor(n)); prod(i=1,#f~,(f[i,1]+1)*((f[i,1]-1)^(f[i,2]-1)))); \\ Antti Karttunen, Jan 06 2021

Formula

a(n) = A167344(n) / A340368(n) = A048250(n) * A326297(n). - Antti Karttunen, Jan 06 2021
Sum_{k=1..n} a(k) ~ c * n^2, where c = (zeta(6)/(2*zeta(2)*zeta(3))) * Product_{p prime} (1 + 2/p^2) = 0.56361239505... . - Amiram Eldar, Nov 12 2022
Previous Showing 11-13 of 13 results.