The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A357588 The compositional inverse of n -> n^[isprime(n)], where [b] is the Iverson bracket of b. 17
 1, -2, 5, -11, 6, 146, -1295, 7712, -36937, 141514, -357676, -322973, 12078666, -102218510, 623243991, -3041134727, 11440387382, -23657862864, -95377084665, 1570488584608, -12255377466362, 72288056416374, -340793435817068, 1186234942871544, -1525020468715715 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,2 LINKS Table of n, a(n) for n=1..25. MAPLE # REVERT from N. J. A. Sloane's 'Transforms' (see the footer of the page). REVERT([seq(if isprime(k) then k else 1 fi, k = 1..25)]); # Alternative: CompInv := proc(len, seqfun) local n, k, m, g, M, A; A := [seq(seqfun(i), i=1..len)]; M := Matrix(len+1, shape=triangular[lower]); M[1, 1] := 1; for m from 2 to len + 1 do M[m, m] := M[m - 1, m - 1]/A[1]; for k from m-1 by -1 to 2 do M[m, k] := M[m-1, k-1] - add(A[i+1]*M[m, k+i], i=1..m-k)/A[1] od od; seq(M[k, 2], k=2..len + 1) end: CompInv(25, n -> if isprime(n) then n else 1 fi); CROSSREFS Cf. A089026, A357368. Sequence in context: A363554 A127011 A170868 * A276736 A069162 A300677 Adjacent sequences: A357585 A357586 A357587 * A357589 A357590 A357591 KEYWORD sign AUTHOR _Peter Luschny_, Oct 04 2022 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified March 3 17:49 EST 2024. Contains 370512 sequences. (Running on oeis4.)