A288636 -id:A288636 - cp's OEIS frontend

A082090 Length of iteration sequence if function A056239, a pseudo-logarithm is iterated and started at n. Fixed point equals zero for all initial values.

2, 3, 4, 4, 5, 5, 5, 5, 5, 5, 6, 5, 6, 6, 6, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 6, 7, 6, 6, 6, 7, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 6, 7, 6, 6, 6, 6, 6, 7, 6, 7, 6, 7, 6, 7, 7, 6, 6, 6, 6, 7, 6, 6, 7, 7, 6, 6, 7, 6, 6, 7, 6, 6, 7, 7, 6, 7, 6, 7, 6, 6, 6, 7, 6, 7, 6, 6, 7

Offset: 1

Labos Elemer, Apr 09 2003

nonn

From Gus Wiseman, Dec 01 2023: (Start)
Conjecture:
- The position of first appearance of k is n = A007097(k-2).
- The position of last appearance of k is n = A014221(k-2) = 2^^(k-2).
- The number of times k appears is: 1, 1, 2, 8, 435, ...
(End)

			n=127:list={127,31,11,5,3,2,1,0},a[127]=8

Mohammad K. Azarian, On the Fixed Points of a Function and the Fixed Points of its Composite Functions, International Journal of Pure and Applied Mathematics, Vol. 46, No. 1, 2008, pp. 37-44. Mathematical Reviews, MR2433713 (2009c:65129), March 2009. Zentralblatt MATH, Zbl 1160.65015.
Mohammad K. Azarian, Fixed Points of a Quadratic Polynomial, Problem 841, College Mathematics Journal, Vol. 38, No. 1, January 2007, p. 60. Solution published in Vol. 39, No. 1, January 2008, pp. 66-67.

Alois P. Heinz, Table of n, a(n) for n = 1..1000

Cf. A008475, A001414, A082083-A082086, A082091.
A112798 lists prime indices, length A001222, sum A056239.
A304038 lists distinct prime indices, length A001221, sum A066328.
Cf. A000079, A000720, A052410, A181819, A225485, A238747, A288636, A296150.

f:= n-> add (numtheory[pi](i[1])*i[2], i=ifactors(n)[2]):
a:= n-> 1+ `if`(n=1, 1, a(f(n))):
seq (a(n), n=1..120);  # Alois P. Heinz, Aug 09 2012

ffi[x_] := Flatten[FactorInteger[x]] lf[x_] := Length[FactorInteger[x]] ba[x_] := Table[Part[ffi[x], 2*w-1], {w, 1, lf[x]}] ep[x_] := Table[Part[ffi[x], 2*w], {w, 1, lf[x]}] bpi[x_] := Table[PrimePi[Part[ba[x], j]], {j, 1, lf[x]}] api[x_] := Apply[Plus, ep[x]*bpi[x]] Table[Length[FixedPointList[api, w]]-1, {w, 2, 128}]
Table[Length[FixedPointList[Total[PrimePi/@Join@@ ConstantArray@@@FactorInteger[#]]&,n]]-1, {n,100}] (* Gus Wiseman, Dec 01 2023 *)

A294339 Number of ways to write 2^n as a finite power-tower of positive integers greater than one, allowing both left and right nesting of parentheses.

Original entry on oeis.org

1, 2, 2, 5, 2, 6, 2, 12, 5, 6, 2, 19, 2, 6, 6, 32, 2, 19, 2, 19, 6, 6, 2, 56, 5, 6, 12, 19, 2, 26, 2, 79, 6, 6, 6, 71, 2, 6, 6, 56, 2, 26, 2, 19, 19, 6, 2, 169, 5, 19, 6, 19, 2, 56, 6, 56, 6, 6, 2, 101, 2, 6, 19, 203, 6, 26, 2, 19, 6, 26, 2, 237, 2, 6, 19, 19

Offset: 1

Gus Wiseman, Oct 28 2017

nonn

			The a(6) = 6 ways are 64, 8^2, (2^3)^2, 4^3, (2^2)^3, 2^6.

Hans Havermann, Table of n, a(n) for n = 1..10000

Cf. A001055, A001597, A007916, A052409, A052410, A089723, A164336, A277562, A284639, A288636, A294336, A294337, A294338.

f:= proc(n) option remember; local F,t,s,g,a;
  F:= ifactors(n)[2];
  g:= igcd(op(map(t -> t[2],F)));
  t:= 1;
  for s in numtheory:-divisors(g) minus {1} do
    t:= t + procname(mul(a[1]^(a[2]/s),a=F))*procname(s)
  od;
  t
end proc:
seq(f(2^n),n=1..100); # Robert Israel, Dec 01 2017

a[n_]:=1+Sum[a[n^(1/g)]*a[g],{g,Rest[Divisors[GCD@@FactorInteger[n][[All,2]]]]}];
Table[a[2^n],{n,100}]

a(n) = A294338(2^n). - R. J. Mathar, Nov 27 2017

A304481 Turn the power-tower for n upside-down.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 9, 8, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 32, 26, 27, 28, 29, 30, 31, 25, 33, 34, 35, 64, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 128, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 36, 65, 66, 67

Offset: 1

Gus Wiseman, May 13 2018

This is an involution of the positive integers.

The power-tower for n is defined as follows. Let {c(i)} = A007916 denote the sequence of numbers > 1 which are not perfect powers. Every positive integer n has a unique representation as a tower n = c(x_1)^c(x_2)^c(x_3)^...^c(x_k), where the exponents are nested from the right. Then a(n) = c(x_k)^...^c(x_3)^c(x_2)^c(x_1).

			The power tower of 81 is 3^2^2, which turned upside-down is 2^2^3 = 256, so a(81) = 256.

Robert Israel, Table of n, a(n) for n = 1..10000
N. J. A. Sloane, Maple programs for A007916, A278028, A278029, A052409, A089723, A277564

Cf. A052409, A052410, A007916, A089723, A164337, A277562, A277564, A278028, A288636, A289023, A294336, A294337, A304491, A304492, A304495.

f:= proc(n,r) local F,a,y;
     if n = 1 then return 1 fi;
     F:= ifactors(n)[2];
     y:= igcd(seq(t[2],t=F));
     if y = 1 then return n^r fi;
     a:= mul(t[1]^(t[2]/y),t=F);
     procname(y,a^r)
end proc:
seq(f(n,1),n=1..100); # Robert Israel, May 13 2018

tow[n_]:=If[n==1,{},With[{g=GCD@@FactorInteger[n][[All,2]]},If[g===1,{n},Prepend[tow[g],n^(1/g)]]]];
Table[Power@@Reverse[tow[n]],{n,100}]

A304491 Last or deepest exponent in the power-tower for n.

Original entry on oeis.org

1, 2, 3, 2, 5, 6, 7, 3, 2, 10, 11, 12, 13, 14, 15, 2, 17, 18, 19, 20, 21, 22, 23, 24, 2, 26, 3, 28, 29, 30, 31, 5, 33, 34, 35, 2, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 2, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 6, 65, 66, 67, 68, 69

Offset: 1

Gus Wiseman, May 13 2018

nonn

Let {c(i)} = A007916 denote the sequence of numbers > 1 which are not perfect powers. Every positive integer n has a unique representation as a tower n = c(x_1)^c(x_2)^c(x_3)^...^c(x_k), where the exponents are nested from the right. Then a(n) = c(x_k).

			We have 16 = 2^2^2, so a(16) = 2.
We have 64 = 2^6, so a(64) = 6.
We have 81 = 3^2^2, so a(81) = 2.
We have 256 = 2^2^3, so a(256) = 3.

Andrew Howroyd, Table of n, a(n) for n = 1..10000

Cf. A052409, A052410, A007916, A089723, A164337, A277562, A277564, A278028, A288636, A289023, A294336, A294337, A304481, A304492, A304495.

a[n_]:=If[n==1,1,With[{g=GCD@@FactorInteger[n][[All,2]]},If[g==1,n,a[g]]]];
Array[a,100]

a(n)={my(t=n); while(t, n=t; t=ispower(t)); n} \\ Andrew Howroyd, Aug 26 2018

a(n) = A007916(A278028(n, A288636(n))).

A304495 Decapitate the power-tower for n, i.e., remove the last (deepest) exponent.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 2, 3, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 3, 1, 1, 1, 1, 2, 1, 1, 1, 6, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 9, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 10, 1, 1, 1, 1, 1

Offset: 1

Gus Wiseman, May 13 2018

nonn

a(1) = 0 by convention.

Let {c(i)} = A007916 denote the sequence of numbers > 1 which are not perfect powers. Every positive integer n has a unique representation as a tower n = c(x_1)^c(x_2)^c(x_3)^...^c(x_k), where the exponents are nested from the right. Then a(n) = c(x_1)^c(x_2)^c(x_3)^...^c(x_{k-1}).

			We have 64 = 2^6, so a(64) = 2.
We have 216 = 6^3, so a(216) = 6.
We have 256 = 2^2^3, so a(256) = 2^2 = 4.

Antti Karttunen, Table of n, a(n) for n = 1..65537

Cf. A052409, A052410, A001597, A007916, A089723, A164337, A277562, A277564, A278028, A288636, A289023, A294336, A294337, A304481, A304491, A304492.

tow[n_]:=If[n==1,{},With[{g=GCD@@FactorInteger[n][[All,2]]},If[g===1,{n},Prepend[tow[g],n^(1/g)]]]];
Table[If[n==1,0,Power@@Most[tow[n]]],{n,100}]

A304495(n) = if(1==n,0,my(e, r, tow = List([])); while((e = ispower(n,,&r)) > 1, listput(tow, r); n = e;); n = 1; while(length(tow)>0, e = tow[#tow]; listpop(tow); n = e^n;); (n)); \\ Antti Karttunen, Jul 23 2018

a(m) <> 1 if m is a perfect power (A001597). - Michel Marcus, Jul 23 2018

Name edited and more terms from Antti Karttunen, Jul 23 2018

A304492 Position in the sequence of numbers that are not perfect powers (A007916) of the last or deepest exponent in the power-tower for n.

Original entry on oeis.org

1, 2, 3, 2, 4, 5, 6, 3, 2, 7, 8, 9, 10, 11, 12, 2, 13, 14, 15, 16, 17, 18, 19, 20, 2, 21, 3, 22, 23, 24, 25, 4, 26, 27, 28, 2, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 2, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 5, 55, 56, 57, 58, 59, 60

Offset: 1

Gus Wiseman, May 13 2018

nonn

Let {c(i)} = A007916 denote the sequence of numbers > 1 which are not perfect powers. Every positive integer n has a unique representation as a tower n = c(x_1)^c(x_2)^c(x_3)^...^c(x_k), where the exponents are nested from the right. Then a(n) = x_k.

Cf. A052409, A052410, A007916, A089723, A164337, A277562, A277564, A278028, A288636, A289023, A294336, A294337, A304481, A304491, A304495.

nn=100;
a[n_]:=If[n==1,1,With[{g=GCD@@FactorInteger[n][[All,2]]},If[g==1,n,a[g]]]];
rads=Union[Array[a,nn]];
Table[a[n],{n,nn}]/.Table[rads[[i]]->i,{i,Length[rads]}]

a(n) = A278028(n, A288636(n)).

cp's OEIS Frontend

A082090 Length of iteration sequence if function A056239, a pseudo-logarithm is iterated and started at n. Fixed point equals zero for all initial values.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

A294339 Number of ways to write 2^n as a finite power-tower of positive integers greater than one, allowing both left and right nesting of parentheses.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A304481 Turn the power-tower for n upside-down.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A304491 Last or deepest exponent in the power-tower for n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A304495 Decapitate the power-tower for n, i.e., remove the last (deepest) exponent.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A304492 Position in the sequence of numbers that are not perfect powers (A007916) of the last or deepest exponent in the power-tower for n.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

Formula