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-4 of 4 results.

A073904 Smallest multiple k*n of n having n divisors.

Original entry on oeis.org

1, 2, 9, 8, 625, 12, 117649, 24, 36, 80, 25937424601, 60, 23298085122481, 448, 2025, 384, 48661191875666868481, 180, 104127350297911241532841, 240, 35721, 11264, 907846434775996175406740561329, 360, 10000, 53248, 26244, 1344
Offset: 1

Views

Author

Amarnath Murthy, Aug 18 2002

Keywords

Comments

Smallest refactorable number, m, such that m=k*n has n divisors. - Robert G. Wilson v, Oct 31 2005

Examples

			Smallest multiple a(n)=k*n; a(1)=1*1, a(2)=1*2, a(3)=3*3, a(4)=2*4, a(5)=125*5, a(6)=2*6, ... having d(k*n)=n divisors; d(1)=1, d(2)=2, d(3^2)=3, d(2^3)=4, d(5^4)=5, d(2^2*3)=3*2=6, ...
		

Crossrefs

Cf. A033950 (refactorable numbers, also known as tau numbers).
Cf. A110821 (SuperRefactorable numbers).

Programs

  • Mathematica
    f[n_] := Block[{k = 1}, If[ PrimeQ[n], n^(n - 1), While[d = DivisorSigma[0, k*n]; d != n, k++ ]; k*n]]; Table[ f[n], {n, 28}] (* Robert G. Wilson v *)

Formula

If p is a prime then a(p) = p^(p-1). If n = p^2 then a(n) = 2^(p-1)*p^(p-1).
a(p^r) = (2*3*5*...*p_r)^(p-1) for r < p <= p_r. a(p^r) = (2*3*...*p_(r-1))^(p-1)*p^(p-1) for p > p_r. Else a(p^r) = ...? for r >= p. Problem a(2^r) = ...? Cf. A005179(p^n)=(2*3*...*p_n)^(p-1) for p_n < 2^p. - Thomas Ordowski, Aug 20 2005
a(p^r) = (2*3...*p_(r-1)*p)^(p-1) for p > p_r; else a(p^r) = (2*3...*p...*p_m)^(p-1)*p^(p^k-p) for p <= p_r and p_m < 2^p, where m=r-k+1 for smallest k such that p^k > r, so k=floor(log(r)/log(p))+1 and p > log(p_m)/log(2). Examples: If k=1 then a(p^r) = (2*3*...*p_r)^(p-1) for r < p <= p_r. If p=2 then a(2^r) = (2*3*...*p_m)*2^(2^k-2) for r < 5. For instance, let r=4 so k=3, m=2 and a(2^4)=384. - Thomas Ordowski, Aug 22 2005
If p is a prime and n=p^r then a(p^r) = (s_1*s_2*...*s_r)^(p-1) where (s_r) is a permutation of the (ascending sequence) numbers of the form q^(p^j) for every prime q and j>=0; permutation such that s_(p^j)=p^(p^j) and shifted remainder. For example, if p=3 then (s_r): 3, 2, 3^3, 5, 7, 2^3, 11, 13, 3^9, 17, 19, ... so a(3^r) = (3*2*27*5*...*s_r)^2. - Thomas Ordowski, Aug 29 2005
If n=2^r then a(2^r) is the product of the first r members of the A109429 sequence. - Thomas Ordowski, Aug 29 2005
a(n) = n * A076931(n). - Thomas Ordowski, Oct 07 2005
a(4) = 8; a(2*prime(n)) = A299795(n), for n>1. - Bernard Schott, Nov 06 2022

Extensions

a(12) corrected by Thomas Ordowski, Aug 18 2005
Further corrections from Thomas Ordowski, Oct 07 2005
a(21), a(27) & a(28) from Robert G. Wilson v, Oct 31 2005

A262983 Terms of A005179 divisible by their indices in order of appearance in A005179.

Original entry on oeis.org

1, 2, 12, 24, 36, 60, 180, 240, 360, 720, 1260, 1680, 3600, 6720, 5040, 10080, 32400, 15120, 20160, 25200, 60480, 55440, 810000, 100800, 181440, 110880, 226800, 221760, 277200, 907200, 665280, 1587600, 720720, 5670000, 1108800, 3548160, 1995840, 1441440, 2494800, 6350400
Offset: 1

Views

Author

Vladimir Letsko, Oct 06 2015

Keywords

Comments

A005179(n) is in this sequence iff it is divisible by n. Thus this is a subsequence of A005179 indexed by A262981.
Also this sequence is the intersection of A033950 and A005179. Hence this sequence has density zero. - Vladimir Letsko, Dec 16 2016
It seems that this sequence is a subsequence of A262981.
This sequence is not in ascending order as terms of A005179 divisible by their number of divisors do not occur in ascending order. For terms sorted in ascending order see A110821. - David A. Corneth, Dec 10 2021

Examples

			12 is a term since it is the smallest positive integer having exactly 6 divisors and divisible by 6.
		

Crossrefs

Programs

  • Mathematica
    Take[#, 33] &@ DeleteCases[#, 0] &@ Function[s, ReplacePart[#, Flatten@ Map[{# -> Function[k, k Boole[Divisible[k, #]]]@ Lookup[s, #]} &, Keys@ s]] &@ ConstantArray[0, Max@ Keys@ s]]@ Map[First, KeySort@ PositionIndex@ Table[DivisorSigma[0, n], {n, 10^7}]] (* Michael De Vlieger, Dec 11 2016, Version 10 *)

Formula

a(n) = A005179(A262981(n)).
A000005(a(n)) = A262981(n).

Extensions

Name clarified by David A. Corneth, Dec 10 2021

A359964 Refactorable numbers (A033950) having more divisors than all smaller refactorable numbers.

Original entry on oeis.org

1, 2, 8, 12, 24, 36, 60, 180, 240, 360, 720, 1260, 1680, 3360, 5040, 10080, 15120, 20160, 25200, 30240, 55440, 100800, 110880, 221760, 277200, 443520, 665280, 720720, 1108800, 1441440, 2494800, 2882880, 3603600, 5765760, 8648640, 12972960, 14414400, 25945920, 28828800
Offset: 1

Views

Author

Amiram Eldar, Jan 20 2023

Keywords

Comments

The corresponding numbers of divisors are 1, 2, 4, 6, 8, 9, 12, 18, 20, 24, ... .
This sequence if infinite since there are refactorable numbers with arbitrarily large number of divisors. E.g., for any prime p, p^(p-1) is a refactorable number with p divisors.

Crossrefs

Subsequence of A033950.
Similar sequences: A002182, A335317, A356078, A359963.

Programs

  • Mathematica
    seq[nmax_] := Module[{s = {}, dm = 0, d}, Do[d = DivisorSigma[0, n]; If[d > dm && Divisible[n, d], dm = d; AppendTo[s, n]], {n, 1, nmax}]; s]; seq[10^6]
  • PARI
    lista(nmax) = {my(dm = 0, d); for(n = 1, nmax, d = numdiv(n); if(d > dm && n%d == 0, dm = d; print1(n, ", "))); }

A174457 Infinitely refactorable numbers: numbers k such that each iteration under the map x -> A000005(x) produces a divisor of k.

Original entry on oeis.org

1, 2, 12, 24, 36, 60, 72, 84, 96, 108, 132, 156, 180, 204, 228, 240, 252, 276, 288, 348, 360, 372, 396, 444, 468, 480, 492, 504, 516, 564, 600, 612, 636, 640, 672, 684, 708, 720, 732, 792, 804, 828, 852, 864, 876, 936, 948, 972, 996, 1044, 1056, 1068, 1116, 1152
Offset: 1

Views

Author

Matthew Vandermast, Dec 04 2010

Keywords

Comments

In other words, let d^1(n) = A000005(n) and, for all positive integers k, let d^(k+1)(n) = A000005(d^k(n)). Sequence lists numbers n with the property that every such value of d^k(n) divides n.
A141586 is a subsequence. Is A110821 a subsequence?
Not a subsequence of A141551: 504 is the smallest term in this sequence not member of A141551.
a(n) is even for all n, since for any n >= 2, d^k(n) = 2 for some k. Proof: {d^k(n)} is a nonincreasing sequence of k, so it must stablize at a fixed point of the map x -> A000005(x), namely x = 1 or 2. But d^k(n) = 1 for some k implies that n = 1. - Jianing Song, Apr 20 2022

Examples

			9 has 3 divisors, and 9 is a multiple of 3. But 3 has 2 divisors, and 9 is not a multiple of 2. Hence, 9 does not belong to this sequence.
36 has 9 divisors, 9 has 3 divisors, 3 has 2 divisors, and 9, 3, and 2 are all divisors of 36. (Since 2 has 2 divisors, all further steps produce a value of 2.) Hence, 36 belongs to this sequence.
		

Crossrefs

Cf. A036459 (number of steps of the map), A000005 (d(n): number of divisors).
Cf. A010553 (d(d(n))), A036450 (d^3(n)), A036452 (d^4(n)), A036453 (d^5(n)).
Subsequence of A033950 (refactorable numbers: d(n) | n) and A141113 (d(d(n))| n).

Programs

  • PARI
    is_A174457(n, d=n)=!until(d<3, n%(d=numdiv(d)) && return) \\ M. F. Hasler, Dec 05 2010, updated PARI syntax Apr 16 2022

Extensions

Edited by M. F. Hasler, Apr 16 2022
Showing 1-4 of 4 results.