cp's OEIS Frontend

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.

Showing 1-5 of 5 results.

A040017 Prime 3 followed by unique period primes (the period r of 1/p is not shared with any other prime) of the form A019328(r)/gcd(A019328(r),r) in order (periods r are given in A051627).

Original entry on oeis.org

3, 11, 37, 101, 9091, 9901, 333667, 909091, 99990001, 999999000001, 9999999900000001, 909090909090909091, 1111111111111111111, 11111111111111111111111, 900900900900990990990991, 909090909090909090909090909091
Offset: 1

Views

Author

Keywords

Comments

Prime p=3 is the only known example of a unique period prime such that A019328(r)/gcd(A019328(r),r) = p^k with k > 1 (cf. A323748). It is plausible to assume that no other such prime exists. Under this (unproved) assumption, the current sequence lists all unique period primes in order and represents a sorted version of A007615. - Max Alekseyev, Oct 14 2022

Examples

			The decimal expansion of 1/101 is 0.00990099..., having a period of 4 and it is the only prime with that period.
		

References

  • J.-P. Delahaye, Merveilleux nombres premiers ("Amazing primes"), p. 324, Pour la Science Paris 2000.

Crossrefs

Programs

  • Mathematica
    lst = {}; Do[c = Cyclotomic[n, 10]; q = c/GCD[c, n]; If[PrimeQ[q], AppendTo[lst, q]], {n, 62}]; Prepend[Sort[lst], 3] (* Arkadiusz Wesolowski, May 13 2012 *)

Formula

For n >= 2, a(n) = A019328(r) / gcd(A019328(r), r), where r = A051627(n). - Max Alekseyev, Oct 14 2022

Extensions

Missing term a(45) inserted in b-file at the suggestion of Eric Chen by Max Alekseyev, Oct 13 2022
Edited by Max Alekseyev, Oct 14 2022

A036275 The periodic part of the decimal expansion of 1/n. Any initial 0's are to be placed at end of cycle.

Original entry on oeis.org

0, 0, 3, 0, 0, 6, 142857, 0, 1, 0, 90, 3, 769230, 714285, 6, 0, 5882352941176470, 5, 526315789473684210, 0, 476190, 45, 4347826086956521739130, 6, 0, 384615, 370, 571428, 3448275862068965517241379310, 3, 322580645161290, 0, 30, 2941176470588235, 285714, 7
Offset: 1

Views

Author

Keywords

Comments

a(n) = 0 iff n = 2^i*5^j (A003592). - Jon Perry, Nov 19 2014
a(n) = n iff n = 3 or 6 (see De Koninck & Mercier reference). - Bernard Schott, Dec 02 2020

Examples

			1/28 = .03571428571428571428571428571428571428571... and digit-cycle is 571428, so a(28)=571428.
		

References

  • Jean-Marie De Koninck & Armel Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problème 347 pp. 50 and 205, Ellipses, Paris, 2004.

Crossrefs

Programs

  • Maple
    isCycl := proc(n) local ifa,i ; if n <= 2 then RETURN(false) ; fi ; ifa := ifactors(n)[2] ; for i from 1 to nops(ifa) do if op(1,op(i,ifa)) <> 2 and op(1,op(i,ifa)) <> 5 then RETURN(true) ; fi ; od ; RETURN(false) ; end: A036275 := proc(n) local ifa,sh,lpow,mpow,r ; if not isCycl(n) then RETURN(0) ; else lpow:=1 ; while true do for mpow from lpow-1 to 0 by -1 do if (10^lpow-10^mpow) mod n =0 then r := (10^lpow-10^mpow)/n ; r := r mod (10^(lpow-mpow)-1) ; while r*10 < 10^(lpow-mpow) do r := 10*r ; od ; RETURN(r) ; fi ; od ; lpow := lpow+1 ; od ; fi ; end: for n from 1 to 60 do printf("%d %d ",n,A036275(n)) ; od ; # R. J. Mathar, Oct 19 2006
  • Mathematica
    fc[n_]:=Block[{q=RealDigits[1/n][[1,-1]]},If[IntegerQ[q],0,While[First[q]==0,q=RotateLeft[q]];FromDigits[q]]];
    Table[fc[n],{n,36}] (* Ray Chandler, Nov 19 2014, corrected Jun 27 2017 *)
    Table[FromDigits[FindTransientRepeat[RealDigits[1/n,10,120][[1]],3] [[2]]],{n,40}] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Mar 12 2019 *)

Extensions

Corrected and extended by N. J. A. Sloane
Corrected a(92), a(208), a(248), a(328), a(352) and a(488) which missed a trailing zero (see the table). - Philippe Guglielmetti, Jun 20 2017

A112505 Number of primitive prime factors of 10^n-1.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 2, 2, 1, 1, 2, 1, 3, 1, 2, 2, 2, 2, 1, 2, 3, 3, 1, 1, 3, 2, 2, 3, 5, 3, 3, 5, 2, 3, 3, 1, 3, 1, 1, 2, 4, 3, 4, 3, 2, 4, 2, 1, 2, 3, 4, 2, 4, 2, 4, 2, 3, 2, 2, 3, 7, 1, 5, 4, 2, 2, 3, 3, 3, 2, 2, 3, 3, 3, 3, 2, 4, 4, 6, 2, 5, 2, 3, 2, 3, 3, 3, 2, 5, 3, 7, 3, 1, 3, 5, 4, 3, 2, 4, 4
Offset: 1

Views

Author

T. D. Noe, Sep 08 2005

Keywords

Comments

Also the number of primes whose reciprocal is a repeating decimal of length n. The number of numbers in each row of table A046107.
By Zsigmondy's theorem, a(n) >= 1. When a(n)=1, the corresponding prime is called a unique prime (see A007498, A040017 and A051627).

Crossrefs

Cf. A007138 (smallest primitive prime factor of 10^n-1), A102347 (number of distinct prime factors of 10^n-1), A046107.

Programs

  • Mathematica
    pp={}; Table[f=Transpose[FactorInteger[10^n-1]][[1]]; p=Complement[f, pp]; pp=Union[pp, p]; Length[p], {n, 66}]

Extensions

Terms to a(276) in b-file from T. D. Noe, Jun 01 2010
a(277)-a(322) in b-file from Ray Chandler, May 01 2017
a(323)-a(352) in b-file from Max Alekseyev, Apr 28 2022

A247071 Numbers n such that 2^n-1 has only one primitive prime factor, sorted according to the magnitude of the corresponding prime.

Original entry on oeis.org

2, 4, 3, 10, 12, 8, 18, 5, 20, 14, 9, 7, 15, 24, 16, 30, 21, 22, 26, 42, 13, 34, 40, 32, 54, 17, 38, 27, 19, 33, 46, 56, 90, 78, 62, 31, 80
Offset: 1

Views

Author

Eric Chen, Nov 16 2014

Keywords

Comments

Periods associated with A144755 in base 2. The binary analog of A051627.

Examples

			2^12 - 1 = 4095 = 3 * 3 * 5 * 7 * 13, but none of 3, 5, 7 is a primitive prime factor, so the only primitive prime factor of 2^12 - 1 is 13.
		

Crossrefs

Programs

  • Mathematica
    nmax = 65536; primesPeriods = Reap[Do[p = Cyclotomic[n, 2]/GCD[n, Cyclotomic[n, 2]]; If[PrimeQ[p], Print[n]; Sow[{p, n}]], {n, 1, nmax}]][[2, 1]]; Sort[primesPeriods][[All, 2]]

Formula

a(n) = A002326((A144755(n+1)-1)/2). - Max Alekseyev, Feb 11 2024

Extensions

Sequence trimmed to the established terms of A144755 by Max Alekseyev, Feb 11 2024

A072848 Largest prime factor of 10^(6*n) + 1.

Original entry on oeis.org

9901, 99990001, 999999000001, 9999999900000001, 39526741, 3199044596370769, 4458192223320340849, 75118313082913, 59779577156334533866654838281, 100009999999899989999000000010001, 2361000305507449, 111994624258035614290513943330720125433979169
Offset: 1

Views

Author

Rick L. Shepherd, Jul 25 2002

Keywords

Comments

According to the link, there are only 18 "unique primes" below 10^50. The first four terms above are each unique primes, of periods 12, 24, 36 and 48, respectively, according to Caldwell and the cross-referenced sequences. These are precisely the only unique primes (less than 10^50 at least) with this type of digit pattern: m 9's, m-1 0's and 1, in that order. (Also a(10) is a unique prime of period 120.)

Examples

			10^(6*4)+1 = 17 * 5882353 * 9999999900000001, so a(4) = 9999999900000001, the largest prime factor.
		

Crossrefs

Cf. A040017 (unique period primes), A051627 (associated periods).

Programs

  • PARI
    for(n=1,12,v=factor(10^(6*n)+1); print1(v[matsize(v)[1],1],","))

Formula

a(n) = A003021(6n) = A006530(A062397(6n)). - Ray Chandler, May 11 2017
Showing 1-5 of 5 results.