|
%I #26 Apr 03 2023 10:36:09
%S 1,2,3,4,10,12,9,14,24,36,48,38,19,23,39,62,120,150,106,93,134,294,
%T 196,320,654,738,385,586,317,597,1404,945,1452,1836,1752,1172,1812,
%U 1282,1426,2232,1862,1844,1521,2134,3750,1031,2264,2667,4354,3927,4274,6522,3903,6022,6682,6135,9550,5877
%N Periods associated with A040017.
%C The numbers in A007498 sorted according to the magnitude of the corresponding prime. - _T. D. Noe_, Sep 08 2005
%H Ray Chandler, <a href="/A051627/b051627.txt">Table of n, a(n) for n = 1..99</a>
%H C. K. Caldwell, <a href="https://t5k.org/glossary/page.php?sort=UniquePrime">Unique Primes</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/UniquePrime.html">Unique Prime</a>
%H <a href="/index/1#1overn">Index entries for sequences related to decimal expansion of 1/n</a>
%F a(n) = A002371(A000720(A040017(n))). - _Max Alekseyev_, Oct 14 2022
%e The decimal expansion of 1/101 is 0.00990099..., having a period of 4 and it is the only prime with that period.
%t nmax = 10000; primesPeriods = Reap[Do[p = Cyclotomic[n, 10]/GCD[n, Cyclotomic[n, 10]]; If[PrimeQ[p], Print[n]; Sow[{p, n}]], {n, 1, nmax}]][[2, 1]]; Sort[primesPeriods][[All, 2]] // Prepend[#, 1]& // Take[#, 58]& (* _Jean-François Alcover_, Mar 29 2013 *)
%K nonn,nice,base
%O 1,2
%A _N. J. A. Sloane_
%E More terms from _Jud McCranie_
%E More terms from _T. D. Noe_, Sep 08 2005
%E Corrected a(45)=3750 and extended by _Ray Chandler_, Oct 13 2008
|
