A051626 -id:A051626 - cp's OEIS frontend

A003592 Numbers of the form 2^i*5^j with i, j >= 0.

1, 2, 4, 5, 8, 10, 16, 20, 25, 32, 40, 50, 64, 80, 100, 125, 128, 160, 200, 250, 256, 320, 400, 500, 512, 625, 640, 800, 1000, 1024, 1250, 1280, 1600, 2000, 2048, 2500, 2560, 3125, 3200, 4000, 4096, 5000, 5120, 6250, 6400, 8000, 8192, 10000, 10240, 12500, 12800

Offset: 1

N. J. A. Sloane

These are the natural numbers whose reciprocals are terminating decimals. - David Wasserman, Feb 26 2002
A132726(a(n), k) = 0 for k <= a(n); A051626(a(n)) = 0; A132740(a(n)) = 1; A132741(a(n)) = a(n). - Reinhard Zumkeller, Aug 27 2007
Where record values greater than 1 occur in A165706: A165707(n) = A165706(a(n)). - Reinhard Zumkeller, Sep 26 2009
Also numbers that are divisible by neither 10k - 7, 10k - 3, 10k - 1 nor 10k + 1, for all k > 0. - Robert G. Wilson v, Oct 26 2010
A204455(5*a(n)) = 5, and only for these numbers. - Wolfdieter Lang, Feb 04 2012
Since p = 2 and q = 5 are coprime, sum_{n >= 1} 1/a(n) = sum_{i >= 0} sum_{j >= 0} 1/p^i * 1/q^j = sum_{i >= 0} 1/p^i q/(q - 1) = p*q/((p-1)*(q-1)) = 2*5/(1*4) = 2.5. - Franklin T. Adams-Watters, Jul 07 2014
Conjecture: Each positive integer n not among 1, 4 and 12 can be written as a sum of finitely many numbers of the form 2^a*5^b + 1 (a,b >= 0) with no one dividing another. This has been verified for n <= 3700. - Zhi-Wei Sun, Apr 18 2023
1,2 and 4,5 are the only consecutive terms in the sequence. - Robin Jones, May 03 2025

Albert H. Beiler, Recreations in the theory of numbers, New York, Dover, (2nd ed.) 1966. See p. 73.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Vaclav Kotesovec, Graph - the asymptotic ratio (200000 terms)
Math Stack Exchange, Are (1,2) and (4,5) the only two consecutive pairs in A003592, the integers of the form 2^i5^j?
Eric Weisstein's World of Mathematics, Regular Number
Eric Weisstein's World of Mathematics, Decimal Expansion

Complement of A085837. Cf. A094958, A022333 (list of j), A022332 (list of i).
Cf. A003586, A003591, A003593, A003594, A003595, A257997.
Cf. A164768 (difference between consecutive terms)

Filtered([1..10000],n->PowerMod(10,n,n)=0); # Muniru A Asiru, Mar 19 2019

import Data.Set (singleton, deleteFindMin, insert)
a003592 n = a003592_list !! (n-1)
a003592_list = f $ singleton 1 where
   f s = y : f (insert (2 * y) $ insert (5 * y) s')
               where (y, s') = deleteFindMin s
-- Reinhard Zumkeller, May 16 2015

[n: n in [1..10000] | PrimeDivisors(n) subset [2,5]]; // Bruno Berselli, Sep 24 2012

isA003592 := proc(n)
      if n = 1 then
        true;
    else
        return (numtheory[factorset](n) minus {2,5} = {} );
    end if;
end proc:
A003592 := proc(n)
     option remember;
     if n = 1 then
        1;
    else
        for a from procname(n-1)+1 do
            if isA003592(a) then
                return a;
            end if;
        end do:
    end if;
end proc: # R. J. Mathar, Jul 16 2012

twoFiveableQ[n_] := PowerMod[10, n, n] == 0; Select[Range@ 10000, twoFiveableQ] (* Robert G. Wilson v, Jan 12 2012 *)
twoFiveableQ[n_] := Union[ MemberQ[{1, 3, 7, 9}, # ] & /@ Union@ Mod[ Rest@ Divisors@ n, 10]] == {False}; twoFiveableQ[1] = True; Select[Range@ 10000, twoFiveableQ] (* Robert G. Wilson v, Oct 26 2010 *)
maxExpo = 14; Sort@ Flatten@ Table[2^i * 5^j, {i, 0, maxExpo}, {j, 0, Log[5, 2^(maxExpo - i)]}] (* Or *)
Union@ Flatten@ NestList[{2#, 4#, 5#} &, 1, 7] (* Robert G. Wilson v, Apr 16 2011 *)

list(lim)=my(v=List(),N);for(n=0,log(lim+.5)\log(5),N=5^n;while(N<=lim,listput(v,N);N<<=1));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jun 28 2011

# A003592.py
from heapq import heappush, heappop
def A003592():
    pq = [1]
    seen = set(pq)
    while True:
        value = heappop(pq)
        yield value
        seen.remove(value)
        for x in 2*value, 5*value:
            if x not in seen:
                heappush(pq, x)
                seen.add(x)
sequence = A003592()
A003592_list = [next(sequence) for _ in range(100)]

from sympy import integer_log
def A003592(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum((x//5**i).bit_length() for i in range(integer_log(x,5)[0]+1))
    return bisection(f,n,n) # Chai Wah Wu, Feb 24 2025

def isA003592(n) :
    return not any(d != 2 and d != 5 for d in prime_divisors(n))
@CachedFunction
def A003592(n) :
    if n == 1 : return 1
    k = A003592(n-1) + 1
    while not isA003592(k) : k += 1
    return k
[A003592(n) for n in (1..48)]  # Peter Luschny, Jul 20 2012

The characteristic function of this sequence is given by Sum_{n >= 1} x^a(n) = Sum_{n >= 1} mu(10*n)*x^n/(1 - x^n), where mu(n) is the Möbius function A008683. Cf. with the formula of Hanna in A051037. - Peter Bala, Mar 18 2019
a(n) ~ exp(sqrt(2*log(2)*log(5)*n)) / sqrt(10). - Vaclav Kotesovec, Sep 22 2020
a(n) = 2^A022332(n) * 5^A022333(n). - R. J. Mathar, Jul 06 2025

Incomplete Python program removed by David Radcliffe, Jun 27 2016

A002371 Period of decimal expansion of 1/(n-th prime) (0 by convention for the primes 2 and 5).

Original entry on oeis.org

0, 1, 0, 6, 2, 6, 16, 18, 22, 28, 15, 3, 5, 21, 46, 13, 58, 60, 33, 35, 8, 13, 41, 44, 96, 4, 34, 53, 108, 112, 42, 130, 8, 46, 148, 75, 78, 81, 166, 43, 178, 180, 95, 192, 98, 99, 30, 222, 113, 228, 232, 7, 30, 50, 256, 262, 268, 5, 69, 28, 141, 146, 153, 155, 312, 79, 110

Offset: 1

N. J. A. Sloane

a(n) is the minimum solution x of modular equation 10^x == 1 (mod p), where p = prime(n). - Carmine Suriano, Oct 10 2012
a(n) = smallest m such that 111...11 (m 1's) is divisible by the n-th prime, or 0 if no such m exists (with the exception that a(2) = 3 instead of 1). E.g., the 5th prime, 11, divides 11, so a(5) = 2. - N. J. A. Sloane, Oct 03 2013 [Comment corrected by Derek Orr, Jun 14 2014]
Numbers n such that A071126(n) = A000040(n) - 1. - Hugo Pfoertner, Mar 18 2003
Except for n = 1 and 3, a(n) divides A006093(n). - Robert Israel, Jul 15 2016

			A002371(11) = 15 because the 11th prime is 31, and 1/31 = 0.03225806451612903225806451612903225806452... has period 15. - _Richard F. Lyon_, Mar 29 2022

Albert H. Beiler, Recreations in the Theory of Numbers, 2nd ed. New York: Dover, 1966, pages 65, 309. ISBN 0-486-21096-0.
John H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, 1996, p. 162. ISBN 978-0-387-97993-9.
D. H. Lehmer, Guide to Tables in the Theory of Numbers. Bulletin No. 105, National Research Council, Washington, DC, 1941, p. 15.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Alois P. Heinz, Table of n, a(n) for n = 1..20000 (first 1000 terms from T. D. Noe)
C. K. Caldwell, The Prime Glossary, Period of a prime
Matt Parker and Brady Haran, The Reciprocals of Primes, Numberphile video (2022).
William Shanks, On the number of figures in the period of the reciprocal of every prime number below 20,000, Proc. Royal Soc. London, 22 (1874), 200-210. See also on JSTOR.
Eric Weisstein's World of Mathematics, Decimal Expansion
Index entries for sequences related to decimal expansion of 1/n

See A048595 for another version. Cf. A006883, A007732, A051626, A071126, A000040, A002275, A097443.
Cf. A001913 (full repetend primes), A060257 (1/prime(n) has period prime(n) - 1).
In other bases: A014664, A062117, A211241, A211242, A211243, A211244, A211245, A002371.

seq(subs(FAIL=0,numtheory:-order(10, ithprime(n))),n=1..100); # Robert Israel, Jul 15 2016

Table[ Length[ RealDigits[1 / Prime[n]] [[1, 1]]], {n, 1, 70}]
Table[If[IntegerQ[#], #, 0] &[MultiplicativeOrder[10, Prime[n]]], {n, 1, 70}] (* Jan Mangaldan, Jul 07 2020 *)

a(n)=if(n<4,n==2,znorder(Mod(10, prime(n))))

from sympy import prime, n_order
def A002371(n): return 0 if n == 1 or n == 3 else n_order(10,prime(n)) # Chai Wah Wu, Feb 07 2022

From Alexander Adamchuk, Jan 28 2007: (Start)
a(A000720(p)) = p - 1 for primes p in A001913.
a(A060257(n)) = prime(A060257(n)) - 1. (End)

More terms from Arlin Anderson (starship1(AT)gmail.com)

Edited by Charles R Greathouse IV, Mar 24 2010

A001913 Full reptend primes: primes with primitive root 10.

Original entry on oeis.org

7, 17, 19, 23, 29, 47, 59, 61, 97, 109, 113, 131, 149, 167, 179, 181, 193, 223, 229, 233, 257, 263, 269, 313, 337, 367, 379, 383, 389, 419, 433, 461, 487, 491, 499, 503, 509, 541, 571, 577, 593, 619, 647, 659, 701, 709, 727, 743, 811, 821, 823, 857, 863, 887, 937, 941, 953, 971, 977, 983

Offset: 1

N. J. A. Sloane

Primes p such that the decimal expansion of 1/p has period p-1, which is the greatest period possible for any integer.
Primes p such that the corresponding entry in A002371 is p-1.
Pieter Moree writes (Oct 20 2004): Assuming the Generalized Riemann Hypothesis it can be shown that the density of primes p such that a prescribed integer g has order (p-1)/t, with t fixed exists and, moreover, it can be computed. This density will be a rational number times the so-called Artin constant. For 2 and 10 the density of primitive roots is A, the Artin constant itself.
R. K. Guy writes (Oct 20 2004): MR 2004j:11141 speaks of the unearthing by Lenstra & Stevenhagen of correspondence concerning the density of this sequence between the Lehmers & Artin.
Also called long period primes, long primes or maximal period primes.
The base-10 cyclic numbers A180340, (b^(p-1) - 1) / p, with b = 10, are obtained from the full reptend primes p. - Daniel Forgues, Dec 17 2012
The number of terms < 10^n: A086018(n). - Robert G. Wilson v, Aug 18 2014

			7 is in the sequence because 1/7 = 0.142857142857... and the period = 7-1 = 6.

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 864.
Albert H. Beiler, Recreations in the Theory of Numbers, 2nd ed. New York: Dover, 1966, pages 65, 309.
John H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, p. 161.
C. F. Gauss, Disquisitiones Arithmeticae, Yale, 1965; see p. 380.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 115.
M. Kraitchik, Recherches sur la Théorie des Nombres. Gauthiers-Villars, Paris, Vol. 1, 1924, Vol. 2, 1929, see Vol. 1, p. 61.
H. Rademacher and O. Toeplitz, Von Zahlen und Figuren (Springer 1930, reprinted 1968), Ch. 19, 'Die periodischen Dezimalbrüche'.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Robert Israel, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
B. Chanco, Full Reptend Prime
Sebastian M. Cioabă and Werner Linde, A Bridge to Advanced Mathematics: from Natural to Complex Numbers, Amer. Math. Soc. (2023) Vol. 58, see page 186.
L. J. Goldstein, Density questions in algebraic number theory, Amer. Math. Monthly, 78 (1971), 342-349.
Pieter Moree, Artin's primitive root conjecture - a survey
Katsuya Mori, On a Certain Inverse Problem for Carousel Numbers, INTEGERS 20 (2020), #A77.
OEIS Wiki, Full reptend primes
Matt Parker and Brady Haran, The Reciprocals of Primes, Numberphile video (2022)
Eric Weisstein's World of Mathematics, Cyclic Number.
Eric Weisstein's World of Mathematics, Decimal Expansion.
Eric Weisstein's World of Mathematics, Full Reptend Prime.
D. Williams, Primitive Roots (Check) [Dead link]
Chai Wah Wu, Pigeonholes and repunits, Amer. Math. Monthly, 121 (2014), 529-533.
Index entries for primes by primitive root
Index entries for sequences related to decimal expansion of 1/n

Apart from initial term, identical to A006883.
Other definitions of cyclic numbers: A003277, A001914, A180340.
Cf. A005596, A001122, A048296, A051626.

A001913 := proc(n) local st, period:
st := ithprime(n):
period := numtheory[order](10,st):
if (st-1 = period) then
   RETURN(st):
fi: end:  seq(A001913(n), n=1..200); # Jani Melik, Feb 25 2011

pr=10; Select[Prime[Range[200]], MultiplicativeOrder[pr, # ] == #-1 &]
(* Second program: *)
Join[{7},Select[Prime[Range[300]],PrimitiveRoot[#,10]==10&]] (* Harvey P. Dale, Feb 01 2018 *)

forprime(p=7,1e3,if(znorder(Mod(10,p))+1==p,print1(p", "))) \\ Charles R Greathouse IV, Feb 27 2011

is(n)=Mod(10,n)^(n\2)==-1 && isprime(n) && znorder(Mod(10,n))+1==n \\ Charles R Greathouse IV, Oct 24 2013

from itertools import count, islice
from sympy import nextprime, n_order
def A001913_gen(startvalue=1): # generator of terms >= startvalue
    p = max(startvalue-1,1)
    while (p:=nextprime(p)):
        if p!=2 and p!=5 and n_order(10,p)==p-1:
            yield p
A001913_list = list(islice(A001913_gen(),20)) # Chai Wah Wu, Mar 03 2025

A007732 Period of decimal representation of 1/n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 6, 1, 1, 1, 2, 1, 6, 6, 1, 1, 16, 1, 18, 1, 6, 2, 22, 1, 1, 6, 3, 6, 28, 1, 15, 1, 2, 16, 6, 1, 3, 18, 6, 1, 5, 6, 21, 2, 1, 22, 46, 1, 42, 1, 16, 6, 13, 3, 2, 6, 18, 28, 58, 1, 60, 15, 6, 1, 6, 2, 33, 16, 22, 6, 35, 1, 8, 3, 1, 18, 6, 6, 13, 1, 9, 5, 41, 6, 16, 21, 28, 2, 44, 1

Offset: 1

N. J. A. Sloane, Hal Sampson [ hals(AT)easynet.com ]

Appears to be a divisor of A007733*A007736. - Henry Bottomley, Dec 20 2001
Primes p such that a(p) = p-1 are in A001913. - Dmitry Kamenetsky, Nov 13 2008
When 1/n has a finite decimal expansion (namely, when n = 2^a*5^b), a(n) = 1 while A051626(n) = 0. - M. F. Hasler, Dec 14 2015
a(n.n) >= a(n) where n.n is A020338(n). - Davide Rotondo, Jun 13 2024

J. H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, pp. 159 etc.

Jon E. Schoenfield, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Project Euler, Reciprocal cycles: Problem 26
Index entries for sequences related to decimal expansion of 1/n

Cf. A121341, A066799, A121090, A001913, A084680.

A007732 := proc(n)
    a132740 := 1 ;
    for pe in ifactors(n)[2] do
        if not op(1,pe) in {2,5} then
            a132740 := a132740*op(1,pe)^op(2,pe) ;
        end if;
    end do:
    if a132740 = 1 then
        1 ;
    else
        numtheory[order](10,a132740) ;
    end if;
end proc:
seq(A007732(n),n=1..50) ; # R. J. Mathar, May 05 2023

Table[r = n/2^IntegerExponent[n, 2]/5^IntegerExponent[n, 5]; MultiplicativeOrder[10, r], {n, 100}] (* T. D. Noe, Oct 17 2012 *)

a(n)=znorder(Mod(10,n/2^valuation(n,2)/5^valuation(n,5))) \\ Charles R Greathouse IV, Jan 14 2013

from sympy import n_order, multiplicity
def A007732(n): return n_order(10,n//2**multiplicity(2,n)//5**multiplicity(5,n)) # Chai Wah Wu, Feb 07 2022

def a(n):
    n = ZZ(n)
    rad = 2**n.valuation(2) * 5**n.valuation(5)
    return Zmod(n // rad)(10).multiplicative_order()
[a(n) for n in range(1, 20)]
# F. Chapoton, May 03 2020

Note that if n=r*s where r is a power of 2 and s is odd then a(n)=a(s). Also if n=r*s where r is a power of 5 and s is not divisible by 5 then a(n) = a(s). So we just need a(n) for n not divisible by 2 or 5. This is the smallest number m such that n divides 10^m - 1; m is a divisor of phi(n), where phi = A000010.
phi(n) = n-1 only if n is prime and since a(n) divides phi(n), a(n) can only equal n-1 if n is prime. - Scott Hemphill (hemphill(AT)alumni.caltech.edu), Nov 23 2006
a(n)=a(A132740(n)); a(A132741(n))=a(A003592(n))=1. - Reinhard Zumkeller, Aug 27 2007

More terms from James Sellers, Feb 05 2000

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

Floor van Lamoen

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

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

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.

Philippe Guglielmetti, Table of n, a(n) for n = 1..1001 (first 500 terms from T. D. Noe)
Index entries for sequences related to decimal expansion of 1/n

Cf. A007732, A051626, A002371, A048595, A006883, A007498, A007615, A040017, A051627.
See also A060282, A060283, A060251.
A051628 is length of preamble.

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

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 *)

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

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

Jud McCranie

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

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

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

Robert G. Wilson v, Table of n, a(n) for n = 1..47
Chris Caldwell, The Prime Glossary, Unique prime.
C. K. Caldwell, "Top Twenty" page, Unique.
Chris K. Caldwell and Harvey Dubner, Unique-Period Primes, J. Recreational Math., 29:1 (1998) 43-48.
Ernest G. Hibbs, Component Interactions of the Prime Numbers, Ph. D. Thesis, Capitol Technology Univ. (2022), see p. 33.
Eric Weisstein's World of Mathematics, Unique Prime.
Wikipedia, Unique prime.
Index entries for sequences related to decimal expansion of 1/n

Cf. A007615, A007498, A002371, A048595, A006883, A007732, A051626, A051627, A323748.

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 *)

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

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

A132740 Largest divisor of n coprime to 10.

Original entry on oeis.org

1, 1, 3, 1, 1, 3, 7, 1, 9, 1, 11, 3, 13, 7, 3, 1, 17, 9, 19, 1, 21, 11, 23, 3, 1, 13, 27, 7, 29, 3, 31, 1, 33, 17, 7, 9, 37, 19, 39, 1, 41, 21, 43, 11, 9, 23, 47, 3, 49, 1, 51, 13, 53, 27, 11, 7, 57, 29, 59, 3, 61, 31, 63, 1, 13, 33, 67, 17, 69, 7, 71, 9, 73, 37, 3, 19, 77, 39, 79, 1, 81

Offset: 1

Reinhard Zumkeller, Aug 27 2007

Or: n with all factors of 2 and 5 removed. - M. F. Hasler, Apr 25 2017

			a(1050) = a(2*3*5*5*7) = 3*7 = 21.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A000265, A003592, A069105, A070021, A070022, A070023, A132739, A132741.

a132740 = a132739 . a000265  -- Reinhard Zumkeller, Apr 08 2011

A132740 := proc(n) n/A132741(n) ; end proc: # R. J. Mathar, Sep 06 2011

a[n_] := FixedPoint[ Quotient[#, GCD[#, 10]]& , n]; Table[a[n], {n, 1, 81}] (* Jean-François Alcover, Sep 06 2011, after Vladimir Joseph Stephan Orlovsky *)
Table[SelectFirst[Reverse[Divisors[n]],CoprimeQ[#,10]&],{n,90}] (* Uses the SelectFirst function from Mathematica version 10. - Harvey P. Dale, Mar 22 2015 *)
a[n_] := n / Times @@ ({2, 5}^IntegerExponent[n, {2, 5}]); Array[a, 100] (* Amiram Eldar, Jun 12 2022 *)

a(n)=n/5^valuation(n,5)>>valuation(n,2) \\ Charles R Greathouse IV, Sep 06 2011

a(n) = A000265(A132739(n)) = A132739(A000265(n)) = n / A132741(n);
A051626(a(n)) = A051626(n); A007732(a(n)) = A007732(n);
a(A003592(n)) = 1.
Multiplicative with a(2^e) = 1, a(5^e) = 1 and a(p^e) = p^e for p = 3 and p >= 7.
Dirichlet g.f. zeta(s-1)*(2^s-2)*(5^s-5)/((2^s-1)*(5^s-1)). - R. J. Mathar, Sep 06 2011
Sum_{k=1..n} a(k) ~ (5/18) * n^2. - Amiram Eldar, Nov 28 2022

Edited by M. F. Hasler, Apr 25 2017

A007498 Unique period lengths of primes mentioned in A007615.

Original entry on oeis.org

1, 2, 3, 4, 9, 10, 12, 14, 19, 23, 24, 36, 38, 39, 48, 62, 93, 106, 120, 134, 150, 196, 294, 317, 320, 385, 586, 597, 654, 738, 945, 1031, 1172, 1282, 1404, 1426, 1452, 1521, 1752, 1812, 1836, 1844, 1862, 2134, 2232, 2264, 2667, 3750, 3903, 3927, 4274, 4354

Offset: 1

N. J. A. Sloane, Robert G. Wilson v

Let {Zs(m, 10, 1)} be the Zsigmondy numbers for a = 10, b = 1: Zs(m, 10, 1) is the greatest divisor of 10^m - 1^m that is coprime to 10^r - 1^r for all positive integers r < m. Then this sequence gives m such that Zs(m, 10, 1) is a prime power (e.g., Zs(1, 10, 1) = 9 = 3^2, Zs(2, 10, 1) = 11, Zs(3, 10, 1) = 37, Zs(4, 10, 1) = 101). It is very likely that Zs(m, 10, 1) is prime if m > 1 is in this sequence (note that the Mathematica and PARI programs below are based on this assumption). - Jianing Song, Aug 12 2020

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Samuel Yates, Period Lengths of Exactly One or Two Prime Numbers, J. Rec. Math., 18 (1985), 22-24.

Ray Chandler, Table of n, a(n) for n = 1..106
Chris K. Caldwell, Unique (period) primes and the factorization of cyclotomic polynomials minus one, Mathematica Japonica, 26 (1997), 189-195.
Chris K. Caldwell & H. Dubner, Unique-Period Primes, Table 2 in Journal of Recreational Mathematics 29(1) 46 1998.
Robert G. Wilson v, Notes, n.d.
Index entries for sequences related to decimal expansion of 1/n

Cf. A007615, A002371, A048595, A006883, A007732, A051626, A046108.

Cf. A161508 (unique period lengths in base 2).

lst={1}; Do[p=Cyclotomic[n, 10]/GCD[n, Cyclotomic[n, 10]]; If[PrimeQ[p], AppendTo[lst, n]], {n, 3000}]; lst (* T. D. Noe, Sep 08 2005 *)

isok(n) = if (n==1, 1, my(p = polcyclo(n, 10)); isprime(p/gcd(p, n))); \\ Michel Marcus, Jun 20 2018

More terms from T. D. Noe, Sep 08 2005

a(48)-a(52) from Ray Chandler, Jul 09 2008

A132741 Largest divisor of n having the form 2^i*5^j.

Original entry on oeis.org

1, 2, 1, 4, 5, 2, 1, 8, 1, 10, 1, 4, 1, 2, 5, 16, 1, 2, 1, 20, 1, 2, 1, 8, 25, 2, 1, 4, 1, 10, 1, 32, 1, 2, 5, 4, 1, 2, 1, 40, 1, 2, 1, 4, 5, 2, 1, 16, 1, 50, 1, 4, 1, 2, 5, 8, 1, 2, 1, 20, 1, 2, 1, 64, 5, 2, 1, 4, 1, 10, 1, 8, 1, 2, 25, 4, 1, 2, 1, 80, 1, 2, 1, 4, 5, 2, 1, 8, 1, 10, 1, 4, 1, 2, 5, 32, 1, 2

Offset: 1

Reinhard Zumkeller, Aug 27 2007

The range of this sequence, { a(n); n>=0 }, is equal to A003592. - M. F. Hasler, Dec 28 2015

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A003592, A006519, A007732, A007814, A051626, A060904, A112765, A132740.
Cf. A001620, A082633.
Cf. A379003 (ordinal transform), A379004 (rgs-transform).
Cf. also A355582.

a132741 = f 2 1 where
   f p y x | r == 0    = f p (y * p) x'
           | otherwise = if p == 2 then f 5 y x else y
           where (x', r) = divMod x p
-- Reinhard Zumkeller, Nov 19 2015

A132741 := proc(n) local f,a; f := ifactors(n)[2] ; a := 1; for f in ifactors(n)[2] do if op(1,f) =2 then a := a*2^op(2,f) ; elif op(1,f) =5 then a := a*5^op(2,f) ; end if; end do;a; end proc: # R. J. Mathar, Sep 06 2011

a[n_] := SelectFirst[Reverse[Divisors[n]], MatchQ[FactorInteger[#], {{1, 1}} | {{2, }} | {{5, }} | {{2, }, {5, }}]&]; Array[a, 100] (* Jean-François Alcover, Feb 02 2018 *)
a[n_] := Times @@ ({2, 5}^IntegerExponent[n, {2, 5}]); Array[a, 100] (* Amiram Eldar, Jun 12 2022 *)

A132741(n)=5^valuation(n,5)<M. F. Hasler, Dec 28 2015

a(n) = n / A132740(n).
a(A003592(n)) = A003592(n).
A051626(a(n)) = 0.
A007732(a(n)) = 1.
From R. J. Mathar, Sep 06 2011: (Start)
Multiplicative with a(2^e)=2^e, a(5^e)=5^e and a(p^e)=1 for p=3 or p>=7.
Dirichlet g.f. zeta(s)*(2^s-1)*(5^s-1)/((2^s-2)*(5^s-5)). (End)
a(n) = A006519(n)*A060904(n) = 2^A007814(n)*5^A112765(n). - M. F. Hasler, Dec 28 2015
Sum_{k=1..n} a(k) ~ n*(12*log(n)^2 + (24*gamma + 36*log(2) - 24)*log(n) + 24 - 24*gamma - 36*log(2) + 36*gamma*log(2) + 2*log(2)^2 - 18*log(5) + 18*gamma*log(5) + 27*log(2)*log(5) + 2*log(5)^2 + 18*log(5)*log(n) - 24*gamma_1)/(60*log(2)*log(5)), where gamma is Euler's constant (A001620) and gamma_1 is the first Stieltjes constant (A082633). - Amiram Eldar, Jan 26 2023

A007615 Primes with unique period length (the periods are given in A007498).

Original entry on oeis.org

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

Offset: 1

N. J. A. Sloane, Robert G. Wilson v, Mira Bernstein

Additional terms are Phi(n,10)/gcd(n,Phi(n,10)) for the n in A007498, where Phi(n,10) is the n-th cyclotomic polynomial evaluated at 10.

			3 is the only prime p such that decimal expansion of 1/p has (nontrivial) period exactly 1.

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Samuel Yates, Period Lengths of Exactly One or Two Prime Numbers, J. Rec. Math., 18 (1985), 22-24.

Max Alekseyev, Table of n, a(n) for n = 1..98 (terms 1..25 from T. D. Noe; terms 26..31 from Ray Chandler)
C. K. Caldwell, The Prime Glossary, unique prime
Makoto Kamada, Factorizations of Phi_n(10)
Index entries for sequences related to decimal expansion of 1/n

Cf. A007498, A040017, A002371, A048595, A006883, A007732, A051626, A061075, A006530, A019328.

nmax = 50; periods = Reap[ Do[ p = Cyclotomic[n, 10] / GCD[n, Cyclotomic[n, 10]]; If[ PrimeQ[p], Sow[n]], {n, 1, nmax}]][[2, 1]]; Cyclotomic[#, 10] / GCD[#, Cyclotomic[#, 10]]& /@ periods // Prepend[#, 3]& (* Jean-François Alcover, Mar 28 2013 *)

a(n) = A061075(A007498(n)). - Max Alekseyev, Oct 16 2010

a(n) = A006530(A019328(A007498(n))). - Ray Chandler, May 10 2017

cp's OEIS Frontend

A003592 Numbers of the form 2^i*5^j with i, j >= 0.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

Python

Python

Sage

Formula

Extensions

A002371 Period of decimal expansion of 1/(n-th prime) (0 by convention for the primes 2 and 5).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A001913 Full reptend primes: primes with primitive root 10.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

A007732 Period of decimal representation of 1/n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Sage

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica