login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A097248 a(n) = the eventual stable point reached when iterating k -> A097246(k), starting from k = n. 10
1, 2, 3, 3, 5, 6, 7, 6, 5, 10, 11, 5, 13, 14, 15, 5, 17, 10, 19, 15, 21, 22, 23, 10, 7, 26, 15, 21, 29, 30, 31, 10, 33, 34, 35, 15, 37, 38, 39, 30, 41, 42, 43, 33, 7, 46, 47, 15, 11, 14, 51, 39, 53, 30, 55, 42, 57, 58, 59, 7, 61, 62, 35, 15, 65, 66, 67, 51, 69, 70, 71, 30, 73, 74, 21 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

a(n) = r(n,m) with m such that r(n,m)=r(n,m+1), where r(n,k) = A097246(r(n,k-1)), r(n,0)=n. (The original definition.)

A097248(n) = r(n,a(n)).

From Antti Karttunen, Nov 15 2016: (Start)

The above remark could be interpreted to mean that A097249(n) <= a(n).

All terms are squarefree, and the squarefree numbers are the fixed points.

These are also fixed points eventually reached when iterating A277886.

(End)

LINKS

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

FORMULA

a(A005117(n)) = A005117(n).

From Antti Karttunen, Nov 15 2016: (Start)

If A008683(n) <> 0 [when n is squarefree], a(n) = n, otherwise a(n) = a(A097246(n)).

If A277885(n) = 0, a(n) = n, otherwise a(n) = a(A277886(n)).

A007913(a(n)) = a(n).

a(A007913(n)) = A007913(n).

A048675(a(n)) = A048675(n).

a(A260443(n)) = A019565(n).

(End)

MATHEMATICA

Table[FixedPoint[Times @@ Map[#1^#2 & @@ # &, Partition[#, 2, 2] &@ Flatten[FactorInteger[#] /. {p_, e_} /; e >= 2 :> {If[OddQ@ e, {p, 1}, {1, 1}], {NextPrime@ p, Floor[e/2]}}]] &, n], {n, 75}] (* Michael De Vlieger, Mar 18 2017 *)

PROG

(PARI)

A097246(n) = { my(f=factor(n)); prod(i=1, #f~, (nextprime(f[i, 1]+1)^(f[i, 2]\2))*((f[i, 1])^(f[i, 2]%2))); };

A097248(n) = { my(k=A097246(n)); while(k<>n, n = k; k = A097246(k)); k; };

\\ Antti Karttunen, Mar 18 2017

(Scheme, with memoization-macro definec)

;; Two implementations:

(definec (A097248 n) (if (not (zero? (A008683 n))) n (A097248 (A097246 n))))

(definec (A097248 n) (if (zero? (A277885 n)) n (A097248 (A277886 n))))

;; Antti Karttunen, Nov 15 2016

(Python)

from sympy import factorint, nextprime

from operator import mul

def a097246(n):

    f=factorint(n)

    return 1 if n==1 else reduce(mul, [(nextprime(i)**int(f[i]/2))*(i**(f[i]%2)) for i in f])

def a(n):

    k=a097246(n)

    while k!=n:

        n=k

        k=a097246(k)

    return k # Indranil Ghosh, May 15 2017

CROSSREFS

Cf. A007913, A005117, A008683, A019565, A048675, A097246, A097247, A097249, A260443, A277885, A277886, A283475, A283478.

Sequence in context: A085310 A055653 A155918 * A097247 A097246 A277886

Adjacent sequences:  A097245 A097246 A097247 * A097249 A097250 A097251

KEYWORD

nonn

AUTHOR

Reinhard Zumkeller, Aug 03 2004

EXTENSIONS

Name changed and the original definition moved to the Comments section by Antti Karttunen, Nov 15 2016

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 22 09:01 EST 2019. Contains 319362 sequences. (Running on oeis4.)