|
%I #81 May 05 2023 05:35:13
%S 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,
%T 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,
%U 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
%N Period of decimal representation of 1/n.
%C Appears to be a divisor of A007733*A007736. - _Henry Bottomley_, Dec 20 2001
%C Primes p such that a(p)=p-1 are in A001913. [_Dmitry Kamenetsky_, Nov 13 2008]
%C 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
%D J. H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, pp. 159 etc.
%H Jon E. Schoenfield, <a href="/A007732/b007732.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from T. D. Noe)
%H Project Euler, <a href="https://projecteuler.net/problem=26">Reciprocal cycles: Problem 26</a>
%H <a href="/index/1#1overn">Index entries for sequences related to decimal expansion of 1/n</a>
%F 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.
%F 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
%F a(n)=a(A132740(n)); a(A132741(n))=a(A003592(n))=1. - _Reinhard Zumkeller_, Aug 27 2007
%p A000732 := proc(n)
%p a132740 := 1 ;
%p for pe in ifactors(n)[2] do
%p if not op(1,pe) in {2,5} then
%p a132740 := a132740*op(1,pe)^op(2,pe) ;
%p end if;
%p end do:
%p if a132740 = 1 then
%p 1 ;
%p else
%p numtheory[order](10,a132740) ;
%p end if;
%p end proc:
%p seq(A000732(n),n=1..50) ; # _R. J. Mathar_, May 05 2023
%t Table[r = n/2^IntegerExponent[n, 2]/5^IntegerExponent[n, 5]; MultiplicativeOrder[10, r], {n, 100}] (* _T. D. Noe_, Oct 17 2012 *)
%o (PARI) a(n)=znorder(Mod(10,n/2^valuation(n,2)/5^valuation(n,5))) \\ _Charles R Greathouse IV_, Jan 14 2013
%o (Sage)
%o def a(n):
%o n = ZZ(n)
%o rad = 2**n.valuation(2) * 5**n.valuation(5)
%o return Zmod(n // rad)(10).multiplicative_order()
%o [a(n) for n in range(1, 20)]
%o # _F. Chapoton_, May 03 2020
%o (Python)
%o from sympy import n_order, multiplicity
%o def A007732(n): return n_order(10,n//2**multiplicity(2,n)//5**multiplicity(5,n)) # _Chai Wah Wu_, Feb 07 2022
%Y Cf. A121341, A066799, A121090, A001913, A084680.
%K nonn,base,easy,nice
%O 1,7
%A _N. J. A. Sloane_, Hal Sampson [ hals(AT)easynet.com ]
%E More terms from _James A. Sellers_, Feb 05 2000
|
