A281188 -id:A281188 - cp's OEIS frontend

A033950 Refactorable numbers: number of divisors of k divides k. Also known as tau numbers.

1, 2, 8, 9, 12, 18, 24, 36, 40, 56, 60, 72, 80, 84, 88, 96, 104, 108, 128, 132, 136, 152, 156, 180, 184, 204, 225, 228, 232, 240, 248, 252, 276, 288, 296, 328, 344, 348, 360, 372, 376, 384, 396, 424, 441, 444, 448, 450, 468, 472, 480, 488, 492, 504, 516, 536

Offset: 1

Simon Colton (simonco(AT)cs.york.ac.uk)

Kennedy and Cooper show that this sequence has density zero.
Spiro showed more precisely that the number of refactorable numbers less than x is asymptotic to (x/sqrt(log x))(log(log x))^(-1+o(1)). - David Eppstein, Aug 25 2014
Numbers k such that the equation gcd(k,x) = tau(k) has solutions. - Benoit Cloitre, Jun 10 2002
Refactorable numbers are the fixed points of A009230. - Labos Elemer, Nov 18 2002
Let ref(n) denote the characteristic function of the refactorable numbers. Then ref(n) = 1 + floor(n/d(n)) - ceiling(n/d(n)), where d(n) is the number of divisors of n. - Wesley Ivan Hurt, Jan 09 2013, Feb 15 2013
An odd number with an even number of divisors cannot be in the sequence by definition. Therefore all odd terms are squares (A000290). - Ivan N. Ianakiev, Aug 25 2013
A054008(k) = k mod A000005(k). - Reinhard Zumkeller, Sep 17 2014
The only squarefree terms are 1 and 2: if x is a squarefree number that is a product of n distinct primes, its number of divisors is 2^n, so x is refactorable if it contains 2^n as a factor, but that makes it nonsquarefree unless n = 0, 1, hence x = 1, 2. - Waldemar Puszkarz, Jun 10 2016
Every positive integer k occurs as tau(m) for some m in the sequence. If the factorization of k is Product p_i^e_i, then Product p_i^(p_i^e_i-1) has the specified property. For k prime, this is the only such number. - Franklin T. Adams-Watters, Jan 14 2017
Zelinsky (2002) proved that for any j > 0 and for sufficiently large m the number of terms not exceeding m is > j*pi(m), where pi(m) = A000720(m). - Amiram Eldar, Feb 20 2021
Numbers m such that the ratio (number of non-divisors of m)/(number of divisors of m) = A049820(m)/A000005(m) is an integer. - Michel Lagneau, Apr 04 2025

Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004, Section B12, pp. 102-103.
New Scientist, Sep 05 1998, p. 17, para. 3.

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Kushagr Ahuja, Patrick Lei and Dylan Pentland, Tau ideals in number fields, PROMYS 2017.
Alan Bundy, Simon Colton and Toby Walsh, HR - A system for Machine Discovery in Finite Algebras, ECAI 1998.
Simon Colton, Refactorable Numbers - A Machine Invention, J. Integer Sequences, Vol. 2 (1999), Article 99.1.2.
Simon Colton, HR - Automatic Theory Formation in Pure Mathematics.
Robert E. Kennedy and Curtis N. Cooper, Tau numbers, natural density and Hardy and Wright's Theorem 437, International Journal of Mathematics and Mathematical Sciences, Vol. 13, No. 2 (1990), pp. 383-386.
Claudia Spiro, How often is the number of divisors of n a divisor of n?, J. Number Theory, Vol. 21, No. 1 (1985), pp. 81-100.
Joshua Zelinsky, Tau Numbers: A Partial Proof of a Conjecture and Other Results , Journal of Integer Sequences, Vol. 5 (2002), Article 02.2.8.

Cf. A000005, A000720, A039819, A036762, A051278, A051279, A051280, A036763.
Cf. A235353 (subsequence).
Cf. A054008, A281188.

a033950 n = a033950_list !! (n-1)
a033950_list = [x | x <- [1..], x `mod` a000005 x == 0]
-- Reinhard Zumkeller, Dec 28 2011

[ n: n in [1..540] | n mod #Divisors(n) eq 0 ]; // Klaus Brockhaus, Apr 29 2009

with(numtheory):
A033950 := proc(n)
    option remember:
    local k:
    if n=1 then
        return 1:
    else
        for k from procname(n-1)+1 do
            if type(k/tau(k), integer) then
                return k:
            end if:
        end do:
    end if:
end proc:
seq(A033950(n), n=1..56); # Nathaniel Johnston, May 04 2011

Do[If[IntegerQ[n/DivisorSigma[0, n]], Print[n]], {n, 1, 1000}]
Select[ Range[559], Mod[ #, DivisorSigma[0, # ]] == 0 &]
Select[Range[550], Divisible[ #, DivisorSigma[0, # ]]&] (* Waldemar Puszkarz, Jun 10 2016 *)

isA033950(n)=n%numdiv(n)==0 \\ Charles R Greathouse IV, Jun 10 2011

from sympy import divisor_count
print([n for n in range(1, 1001) if not n % divisor_count(n)]) # Indranil Ghosh, May 03 2017

More terms from Erich Friedman

A335489 Number of strict permutations of the prime indices of n.

Original entry on oeis.org

1, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 0, 1, 2, 2, 0, 1, 0, 1, 0, 2, 2, 1, 0, 0, 2, 0, 0, 1, 6, 1, 0, 2, 2, 2, 0, 1, 2, 2, 0, 1, 6, 1, 0, 0, 2, 1, 0, 0, 0, 2, 0, 1, 0, 2, 0, 2, 2, 1, 0, 1, 2, 0, 0, 2, 6, 1, 0, 2, 6, 1, 0, 1, 2, 0, 0, 2, 6, 1, 0, 0, 2, 1, 0, 2, 2, 2

Offset: 1

Gus Wiseman, Jun 19 2020

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Also the number of (1,1)-avoiding permutations of the prime indices of n.

Wikipedia, Permutation pattern
Gus Wiseman, Sequences counting and ranking compositions by the patterns they match or avoid.

Positions of first appearances are A002110 with 2 replaced by 4.
Permutations of prime indices are counted by A008480.
The contiguous version is A335451.
Anti-run permutations of prime indices are counted by A335452.
(1,1,1)-avoiding permutations of prime indices are counted by A335511.
Cf. A056239, A056986, A106356, A112798, A238279, A281188, A333221, A335456, A335460, A335462, A335465.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Select[Permutations[primeMS[n]],!MatchQ[#,{_,x_,_,x_,_}]&]],{n,100}]

If n is squarefree, a(n) = A001221(n)!; otherwise a(n) = 0.

a(n != 4) = A281188(n); a(4) = 0.

A335520 Number of (1,2,3)-matching permutations of the prime indices of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 3, 0, 0, 0

Offset: 1

Gus Wiseman, Jun 19 2020

nonn

We define a pattern to be a finite sequence covering an initial interval of positive integers. Patterns are counted by A000670 and ranked by A333217. A sequence S is said to match a pattern P if there is a not necessarily contiguous subsequence of S whose parts have the same relative order as P. For example, (3,1,1,3) matches (1,1,2), (2,1,1), and (2,1,2), but avoids (1,2,1), (1,2,2), and (2,2,1).

			The a(n) permutations for n = 30, 60, 120, 210, 180, 480:
  (123)  (1123)  (11123)  (1234)  (11223)  (1111123)
         (1213)  (11213)  (1243)  (11232)  (1111213)
         (1231)  (11231)  (1324)  (12123)  (1111231)
                 (12113)  (1342)  (12132)  (1112113)
                 (12131)  (1423)  (12213)  (1112131)
                 (12311)  (2134)  (12231)  (1112311)
                          (2314)  (12312)  (1121113)
                          (2341)  (12321)  (1121131)
                          (3124)  (21123)  (1121311)
                          (4123)  (21213)  (1123111)
                                  (21231)  (1211113)
                                           (1211131)
                                           (1211311)
                                           (1213111)
                                           (1231111)

Wikipedia, Permutation pattern
Gus Wiseman, Sequences counting and ranking compositions by the patterns they match or avoid.

Positions of nonzero terms are A000977.
These permutations are ranked by A335479.
These compositions are counted by A335514.
Patterns matching this pattern are counted by A335515.
The complement A335521 is the avoiding version.
Permutations of prime indices are counted by A008480.
Patterns are counted by A000670 and ranked by A333217.
Anti-run permutations of prime indices are counted by A335452.
Cf. A056239, A056986, A112798, A238279, A281188, A333221, A333755, A335456, A335460, A335462, A335463.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Select[Permutations[primeMS[n]],MatchQ[#,{_,x_,_,y_,_,z_,_}/;x

For n > 0, a(n) + A335521(n) = A008480(n).

A039819 Number of divisors of n-th refactorable number (A033950(n)).

Original entry on oeis.org

1, 2, 4, 3, 6, 6, 8, 9, 8, 8, 12, 12, 10, 12, 8, 12, 8, 12, 8, 12, 8, 8, 12, 18, 8, 12, 9, 12, 8, 20, 8, 18, 12, 18, 8, 8, 8, 12, 24, 12, 8, 16, 18, 8, 9, 12, 14, 18, 18, 8, 24, 8, 12, 24, 12, 8, 20, 12, 8, 8, 24, 18, 5, 8, 12, 16, 8, 24, 18, 12, 8, 30, 12, 8, 24, 12, 8, 8, 18, 12, 8, 24, 8

Offset: 1

N. J. A. Sloane

nonn

A number n is refactorable if the number of divisors of n divides n.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A033950, A000005, A036762, A281188.

v:=[ n: n in [1..900] | n mod NumberOfDivisors(n) eq 0 ]; [NumberOfDivisors(v[i]): i in [1..#v]]; // Marius A. Burtea, Jul 02 2019

fQ[n_] := Mod[n, DivisorSigma[0, n]] == 0; DivisorSigma[0, # ] & /@ Select[ Range[1000], fQ[ # ] &] (* Robert G. Wilson v *)

A033950(n)/a(n) = A036762(n).

a(n) = A000005(A033950(n)). - Omar E. Pol, Jan 17 2017

More terms from Robert G. Wilson v, Oct 29 2005

Minor edits by Franklin T. Adams-Watters, Jan 17 2017

A335511 Number of (1,1,1)-avoiding permutations of the prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 0, 1, 2, 1, 3, 1, 2, 2, 0, 1, 3, 1, 3, 2, 2, 1, 0, 1, 2, 0, 3, 1, 6, 1, 0, 2, 2, 2, 6, 1, 2, 2, 0, 1, 6, 1, 3, 3, 2, 1, 0, 1, 3, 2, 3, 1, 0, 2, 0, 2, 2, 1, 12, 1, 2, 3, 0, 2, 6, 1, 3, 2, 6, 1, 0, 1, 2, 3, 3, 2, 6, 1, 0, 0, 2, 1, 12, 2, 2

Offset: 1

Gus Wiseman, Jun 19 2020

nonn

We define a pattern to be a finite sequence covering an initial interval of positive integers. Patterns are counted by A000670 and ranked by A333217. A sequence S is said to match a pattern P if there is a not necessarily contiguous subsequence of S whose parts have the same relative order as P. For example, (3,1,1,3) matches (1,1,2), (2,1,1), and (2,1,2), but avoids (1,2,1), (1,2,2), and (2,2,1).

Wikipedia, Permutation pattern
Gus Wiseman, Sequences counting and ranking compositions by the patterns they match or avoid.

Patterns avoiding this pattern are counted by A080599.
These compositions are counted by A232432.
The (1,1)-avoiding version is A335451.
The complement A335510 is the matching version.
These permutations are ranked by A335513.
Patterns are counted by A000670 and ranked by A333217.
Permutations of prime indices are counted by A008480.
Anti-run permutations of prime indices are counted by A335452.
Cf. A056239, A056986, A112798, A238279, A281188, A333221, A333755, A335456, A335460, A335462, A335463.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Select[Permutations[primeMS[n]],!MatchQ[#,{_,x_,_,x_,_,x_,_}]&]],{n,100}]

If n is cubefree, a(n) = A008480(n), otherwise a(n) = 0.

A335510 Number of (1,1,1)-matching permutations of the prime indices of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 4, 0, 4, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 10, 0, 0, 0, 0, 0, 0, 0, 5, 1, 0, 0, 0, 0, 0, 0

Offset: 1

Gus Wiseman, Jun 19 2020

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

We define a pattern to be a finite sequence covering an initial interval of positive integers. Patterns are counted by A000670 and ranked by A333217. A sequence S is said to match a pattern P if there is a not necessarily contiguous subsequence of S whose parts have the same relative order as P. For example, (3,1,1,3) matches (1,1,2), (2,1,1), and (2,1,2), but avoids (1,2,1), (1,2,2), and (2,2,1).

Wikipedia, Permutation pattern
Gus Wiseman, Sequences counting and ranking compositions by the patterns they match or avoid.

Patterns matching this pattern are counted by A335508.
These compositions are counted by A335455.
The (1,1)-matching version is A335487.
The complement A335511 is the avoiding version.
These permutations are ranked by A335512.
Permutations of prime indices are counted by A008480.
Patterns are counted by A000670 and ranked by A333217.
Anti-run permutations of prime indices are counted by A335452.
Cf. A056239, A056986, A112798, A238279, A281188, A333221, A333755, A335456, A335460, A335462, A335463.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Select[Permutations[primeMS[n]],MatchQ[#,{_,x_,_,x_,_,x_,_}]&]],{n,0,100}]

If n is cubefree, a(n) = 0; otherwise a(n) = A008480(n).

A335521 Number of (1,2,3)-avoiding permutations of the prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 2, 2, 1, 1, 3, 1, 3, 2, 2, 1, 4, 1, 2, 1, 3, 1, 5, 1, 1, 2, 2, 2, 6, 1, 2, 2, 4, 1, 5, 1, 3, 3, 2, 1, 5, 1, 3, 2, 3, 1, 4, 2, 4, 2, 2, 1, 9, 1, 2, 3, 1, 2, 5, 1, 3, 2, 5, 1, 10, 1, 2, 3, 3, 2, 5, 1, 5, 1, 2, 1, 9, 2, 2, 2

Offset: 1

Gus Wiseman, Jun 19 2020

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

We define a pattern to be a finite sequence covering an initial interval of positive integers. Patterns are counted by A000670 and ranked by A333217. A sequence S is said to match a pattern P if there is a not necessarily contiguous subsequence of S whose parts have the same relative order as P. For example, (3,1,1,3) matches (1,1,2), (2,1,1), and (2,1,2), but avoids (1,2,1), (1,2,2), and (2,2,1).

			The a(n) permutations for n = 1, 6, 12, 24, 30, 36, 60, 72, 120:
  ()  (12)  (112)  (1112)  (132)  (1122)  (1132)  (11122)  (11132)
      (21)  (121)  (1121)  (213)  (1212)  (1312)  (11212)  (11312)
            (211)  (1211)  (231)  (1221)  (1321)  (11221)  (11321)
                   (2111)  (312)  (2112)  (2113)  (12112)  (13112)
                           (321)  (2121)  (2131)  (12121)  (13121)
                                  (2211)  (2311)  (12211)  (13211)
                                          (3112)  (21112)  (21113)
                                          (3121)  (21121)  (21131)
                                          (3211)  (21211)  (21311)
                                                  (22111)  (23111)
                                                           (31112)
                                                           (31121)
                                                           (31211)
                                                           (32111)

Wikipedia, Permutation pattern
Gus Wiseman, Sequences counting and ranking compositions by the patterns they match or avoid.

These compositions are counted by A102726.
Patterns avoiding this pattern are counted by A226316.
The complement A335520 is the matching version.
Permutations of prime indices are counted by A008480.
Patterns are counted by A000670 and ranked by A333217.
Anti-run permutations of prime indices are counted by A335452.
Cf. A056239, A056986, A112798, A238279, A281188, A333221, A333755, A335456, A335460, A335462, A335463.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Select[Permutations[primeMS[n]],!MatchQ[#,{_,x_,_,y_,_,z_,_}/;x

For n > 0, a(n) + A335520(n) = A008480(n).

A335487 Number of (1,1)-matching permutations of the prime indices of n.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 3, 0, 0, 0, 1, 0, 3, 0, 3, 0, 0, 0, 4, 1, 0, 1, 3, 0, 0, 0, 1, 0, 0, 0, 6, 0, 0, 0, 4, 0, 0, 0, 3, 3, 0, 0, 5, 1, 3, 0, 3, 0, 4, 0, 4, 0, 0, 0, 12, 0, 0, 3, 1, 0, 0, 0, 3, 0, 0, 0, 10, 0, 0, 3, 3, 0, 0, 0, 5, 1, 0, 0, 12, 0, 0

Offset: 1

Gus Wiseman, Jun 14 2020

nonn

Depends only on sorted prime signature (A118914).
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
We define a pattern to be a finite sequence covering an initial interval of positive integers. Patterns are counted by A000670. A sequence S is said to match a pattern P if there is a not necessarily contiguous subsequence of S whose parts have the same relative order as P. For example, (3,1,1,3) matches (1,1,2), (2,1,1), and (2,1,2), but avoids (1,2,1), (1,2,2), and (2,2,1).

			The a(n) permutations for n = 4, 12, 24, 48, 36, 72, 60:
  (11)  (112)  (1112)  (11112)  (1122)  (11122)  (1123)
        (121)  (1121)  (11121)  (1212)  (11212)  (1132)
        (211)  (1211)  (11211)  (1221)  (11221)  (1213)
               (2111)  (12111)  (2112)  (12112)  (1231)
                       (21111)  (2121)  (12121)  (1312)
                                (2211)  (12211)  (1321)
                                        (21112)  (2113)
                                        (21121)  (2131)
                                        (21211)  (2311)
                                        (22111)  (3112)
                                                 (3121)
                                                 (3211)

Wikipedia, Permutation pattern
Gus Wiseman, Sequences counting and ranking compositions by the patterns they match or avoid.

Positions of zeros are A005117 (squarefree numbers).
The case where the match must be contiguous is A333175.
The avoiding version is A335489.
The (1,1,1)-matching case is A335510.
Patterns are counted by A000670.
Permutations of prime indices are counted by A008480.
(1,1)-matching patterns are counted by A019472.
(1,1)-matching compositions are counted by A261982.
STC-numbers of permutations of prime indices are A333221.
Patterns matched by standard compositions are counted by A335454.
Dimensions of downsets of standard compositions are A335465.
(1,1)-matching compositions are ranked by A335488.
Cf. A000961, A056239, A056986, A112798, A181796, A335451, A335452, A335460, A335462.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Select[Permutations[primeMS[n]],!UnsameQ@@#&]],{n,100}]

a(n) = 0 if n is squarefree, otherwise a(n) = A008480(n).

a(n) = A008480(n) - A281188(n) for n != 4.

A341781 Refactorable numbers (or tau numbers, A033950) k such that k/tau(k) is even, where tau(k) = A000005(k).

Original entry on oeis.org

8, 12, 36, 72, 80, 96, 128, 180, 240, 252, 288, 384, 396, 448, 468, 480, 560, 612, 640, 672, 684, 720, 828, 864, 880, 896, 972, 1040, 1044, 1056, 1116, 1152, 1200, 1248, 1332, 1344, 1360, 1408, 1440, 1476, 1520, 1548, 1620, 1632, 1664, 1680, 1692, 1800, 1824

Offset: 1

Amiram Eldar, Feb 19 2021

nonn

Zelinsky (2002) called these numbers p-generators. He proved that these are the tau numbers k such that for any prime p, if p does not divide k then p*k is also a tau number. He used these numbers to prove that the number of tau numbers not exceeding m is > pi(m)/2 for all m > 7.42*10^13, where pi(m) = A000720(m).

			8 is a term since 8/tau(8) = 8/4 = 2 is even.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Joshua Zelinsky, Tau Numbers: A Partial Proof of a Conjecture and Other Results, Journal of Integer Sequences, Vol. 5 (2002), Article 02.2.8.

Cf. A000005, A000720, A033950, A036762, A039819, A281188.

q[n_] := Divisible[n, (d = DivisorSigma[0, n])] && EvenQ[n/d]; Select[Range[2000], q]

isok(k) = my(q=k/numdiv(k)); (denominator(q)==1) && ((q%2) == 0); \\ Michel Marcus, Feb 20 2021

cp's OEIS Frontend

A033950 Refactorable numbers: number of divisors of k divides k. Also known as tau numbers.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Extensions

A335489 Number of strict permutations of the prime indices of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A335520 Number of (1,2,3)-matching permutations of the prime indices of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A039819 Number of divisors of n-th refactorable number (A033950(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

Formula

Extensions

A335511 Number of (1,1,1)-avoiding permutations of the prime indices of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A335510 Number of (1,1,1)-matching permutations of the prime indices of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A335521 Number of (1,2,3)-avoiding permutations of the prime indices of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica