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.
%I A080822 #14 Jul 01 2025 19:27:47 %S A080822 23,61,67,137,149,233,547,701,709,887,1583,1723,1777,2239,2459,2477, %T A080822 2693,3067,3167,3329,3881,3947,3967,4073,5333,5449,5839,6217,6709, %U A080822 6833,8221,8273,8501,8527,8641,9227,9341,10193,10729,10771,11939,12197,12347 %N A080822 Primes in A058633. %H A080822 Paolo Xausa, <a href="/A080822/b080822.txt">Table of n, a(n) for n = 1..10000</a> %t A080822 A006577[n_] := Length[NestWhileList[If[OddQ[#], 3*# + 1, #/2] &, n, # > 1 &]] - 1; %t A080822 Select[Accumulate[Array[A006577, 500]], PrimeQ] (* _Paolo Xausa_, Oct 01 2024 *) %o A080822 (PARI) xpcount(n,p) = { ct=0; sr=0; for(x=1,n, p1 = x; while(p1>1, if(p1%2==0,p1/=2; ct++,p1 = p1*p+1; ct++) ); if(isprime(ct),print1(ct" "); sr+=1.0/ct) ); print(); print(sr) } %Y A080822 Cf. A006577, A058633. %K A080822 easy,nonn %O A080822 1,1 %A A080822 _Cino Hilliard_, Mar 26 2003 %E A080822 Edited by Hauke Worpel and _N. J. A. Sloane_, Jun 03 2003