A379004 -id:A379004 - cp's OEIS frontend

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

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

A379005 Lexicographically earliest infinite sequence such that a(i) = a(j) => v_2(i) = v_2(j), v_3(i) = v_3(j) and v_5(i) = v_5(j), for all i, j, where v_2 (A007814), v_3 (A007949) and v_5 (A112765) give the 2-, 3- and 5-adic valuations of n respectively.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 1, 7, 8, 9, 1, 10, 1, 2, 11, 12, 1, 13, 1, 14, 3, 2, 1, 15, 16, 2, 17, 4, 1, 18, 1, 19, 3, 2, 5, 20, 1, 2, 3, 21, 1, 6, 1, 4, 22, 2, 1, 23, 1, 24, 3, 4, 1, 25, 5, 7, 3, 2, 1, 26, 1, 2, 8, 27, 5, 6, 1, 4, 3, 9, 1, 28, 1, 2, 29, 4, 1, 6, 1, 30, 31, 2, 1, 10, 5, 2, 3, 7, 1, 32, 1, 4, 3, 2, 5, 33, 1, 2, 8, 34, 1, 6, 1, 7, 11

Offset: 1

Antti Karttunen, Dec 15 2024

nonn

Restricted growth sequence transform of A355582.
For all i, j:
  A379001(i) = A379001(j) => a(i) = a(j),
  a(i) = a(j) => A322026(i) = A322026(j),
  a(i) = a(j) => A379004(i) = A379004(j).

Antti Karttunen, Table of n, a(n) for n = 1..100000

Cf. A007814, A007949, A112765, A355582, A379006 (ordinal transform).

Cf. A322026, A379001, A379004.

up_to = 100000;
rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
v379005 = rgs_transform(vector(up_to, n, [valuation(n,2), valuation(n,3), valuation(n,5)]));
A379005(n) = v379005[n];

A379003 Ordinal transform of A132741, where A132741 is the largest divisor of n having the form 2^i*5^j. a(0) = 0 by convention.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 2, 3, 1, 4, 1, 5, 2, 6, 3, 2, 1, 7, 4, 8, 1, 9, 5, 10, 2, 1, 6, 11, 3, 12, 2, 13, 1, 14, 7, 3, 4, 15, 8, 16, 1, 17, 9, 18, 5, 4, 10, 19, 2, 20, 1, 21, 6, 22, 11, 5, 3, 23, 12, 24, 2, 25, 13, 26, 1, 6, 14, 27, 7, 28, 3, 29, 4, 30, 15, 2, 8, 31, 16, 32, 1, 33, 17, 34, 9, 7, 18, 35, 5, 36, 4, 37, 10

Offset: 0

Antti Karttunen, Dec 15 2024

nonn

Ordinal transform of the ordered pair [A007814(n), A112765(n)].

This sequence and A379004 are ordinal transforms of each other (if the initial 0 is discarded).

Antti Karttunen, Table of n, a(n) for n = 0..20000

Cf. A007814, A112765, A132741, A379004 (ordinal transform of this sequence after the initial 0).

Cf. also A126760, A379006.

up_to = 20000;
ordinal_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), pt); for(i=1, length(invec), if(mapisdefined(om,invec[i]), pt = mapget(om, invec[i]), pt = 0); outvec[i] = (1+pt); mapput(om,invec[i],(1+pt))); outvec; };
v379003 = ordinal_transform(vector(up_to, n, [valuation(n,2), valuation(n,5)]));
A379003(n) = if(!n,n,v379003[n]);

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A379005 Lexicographically earliest infinite sequence such that a(i) = a(j) => v_2(i) = v_2(j), v_3(i) = v_3(j) and v_5(i) = v_5(j), for all i, j, where v_2 (A007814), v_3 (A007949) and v_5 (A112765) give the 2-, 3- and 5-adic valuations of n respectively.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A379003 Ordinal transform of A132741, where A132741 is the largest divisor of n having the form 2^i*5^j. a(0) = 0 by convention.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI