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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A087717 Start with x=n, then iterate the map x -> A322982(x) with A322982(x)=2*x-1 if x is noncomposite, otherwise A322982(x) = A032742(x), the largest proper divisor of x. If this iteration leads to a fixed point then a(n) is the value of that fixed point. If the iteration leads to a cycle, a(n) is the smallest value in the cycle. If the iteration never becomes periodic then a(n)=0. 3

%I

%S 1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,19,3,3,3,3,3,3,3,3,3,19,3,3,3,3,

%T 3,3,3,19,19,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,19,19,3,3,3,3,3,3,3,

%U 3,19,3,3,3,3,3,19,19,3,19,3,3,3,3,3,3,3,3,3,3,19,3,3,3,3,3,3,3,19,3,3,3,3,3

%N Start with x=n, then iterate the map x -> A322982(x) with A322982(x)=2*x-1 if x is noncomposite, otherwise A322982(x) = A032742(x), the largest proper divisor of x. If this iteration leads to a fixed point then a(n) is the value of that fixed point. If the iteration leads to a cycle, a(n) is the smallest value in the cycle. If the iteration never becomes periodic then a(n)=0.

%C Conjecture. For n > 1, the iteration given in the definition above always leads to the 3-cycle {3,5,9,3} or the 6-cycle {19,37,73,145,29,57,19}, thus a(n) takes on only the values 3 or 19 for n=2,3,4,.... This has been verified to n=1000000.

%C In range 2..100000 term 3 occurs 77630 times, while 19 occurs 22369 times. - _Antti Karttunen_, Jan 03 2019

%H Antti Karttunen, <a href="/A087717/b087717.txt">Table of n, a(n) for n = 1..16384</a>

%H Antti Karttunen, <a href="/A087717/a087717.txt">Data supplement: n, a(n) computed for n = 1..100000</a>

%t Which[Length@ Union@ #[[-2 ;; -1]] == 1, Last@ #, MemberQ[{3, 5, 9}, Last@ #], 3, MemberQ[{19, 37, 73, 145, 29, 57}, Last@ #], 19, True, 0] & /@ Array[NestWhileList[If[CompositeQ@ #, Divisors[#][[-2]], 2 # - 1] &, #, UnsameQ[##] &, All] &, 106] (* _Michael De Vlieger_, Jan 03 2019 *)

%o (PARI)

%o A322982(n) = if((1==n)||isprime(n),n+n-1,n/vecmin(factor(n)[,1]));

%o A087717(n) = { my(visited = Map(), visited_at_step = Map(), j=0, m=0, t); while(!mapisdefined(visited, n), mapput(visited, n, j); mapput(visited_at_step, j, n); j++; n = A322982(n)); for(k=mapget(visited,n), j-1, t = mapget(visited_at_step,k); if(!m || (t<m), m=t)); (m); }; \\ _Antti Karttunen_, Jan 03 2019

%Y Cf. A322982.

%K nonn

%O 1,2

%A _John W. Layman_, Sep 29 2003

%E Name edited and the term a(1) = 1 prepended by _Antti Karttunen_, Jan 03 2019

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 March 23 05:36 EDT 2019. Contains 321422 sequences. (Running on oeis4.)