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.

Previous Showing 11-20 of 53 results. Next

A316428 Heinz numbers of integer partitions such that every part is divisible by the number of parts.

Original entry on oeis.org

1, 2, 3, 5, 7, 9, 11, 13, 17, 19, 21, 23, 29, 31, 37, 39, 41, 43, 47, 49, 53, 57, 59, 61, 67, 71, 73, 79, 83, 87, 89, 91, 97, 101, 103, 107, 109, 111, 113, 125, 127, 129, 131, 133, 137, 139, 149, 151, 157, 159, 163, 167, 169, 173, 179, 181, 183, 191, 193, 197
Offset: 1

Views

Author

Gus Wiseman, Jul 02 2018

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			93499 is the Heinz number of (12,8,8,4) and belongs to the sequence because each part is divisible by 4.
Sequence of partitions such that every part is divisible by the number of parts begins (1), (2), (3), (4), (2,2), (5), (6), (7), (8), (4,2), (9).
		

Crossrefs

Programs

  • Mathematica
    Select[Range[200],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>Divisible[PrimePi[p],PrimeOmega[#]]]&]

A009998 Triangle in which j-th entry in i-th row is (j+1)^(i-j).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 8, 9, 4, 1, 1, 16, 27, 16, 5, 1, 1, 32, 81, 64, 25, 6, 1, 1, 64, 243, 256, 125, 36, 7, 1, 1, 128, 729, 1024, 625, 216, 49, 8, 1, 1, 256, 2187, 4096, 3125, 1296, 343, 64, 9, 1, 1, 512, 6561, 16384, 15625, 7776, 2401, 512, 81, 10, 1
Offset: 0

Views

Author

Keywords

Comments

Read as a square array this is the Hilbert transform of triangle A123125 (see A145905 for the definition of this term). For example, the fourth row of A123125 is (0,1,4,1) and the expansion (x + 4*x^2 + x^3)/(1-x)^4 = x + 8*x^2 + 27*x^3 + 64*x^4 + ... generates the entries in the fourth row of this array read as a square. - Peter Bala, Oct 28 2008

Examples

			Triangle begins:
  1;
  1,  1;
  1,  2,  1;
  1,  4,  3,  1;
  1,  8,  9,  4,  1;
  1, 16, 27, 16,  5,  1;
  1, 32, 81, 64, 25,  6,  1;
  ...
From _Gus Wiseman_, May 01 2021: (Start)
The rows of the triangle are obtained by reading antidiagonals upward in the following table of A(k,n) = n^k, with offset k = 0, n = 1:
         n=1:     n=2:     n=3:     n=4:     n=5:     n=6:
   k=0:   1        1        1        1        1        1
   k=1:   1        2        3        4        5        6
   k=2:   1        4        9       16       25       36
   k=3:   1        8       27       64      125      216
   k=4:   1       16       81      256      625     1296
   k=5:   1       32      243     1024     3125     7776
   k=6:   1       64      729     4096    15625    46656
   k=7:   1      128     2187    16384    78125   279936
   k=8:   1      256     6561    65536   390625  1679616
   k=9:   1      512    19683   262144  1953125 10077696
  k=10:   1     1024    59049  1048576  9765625 60466176
(End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 24.

Crossrefs

Row sums give A026898.
Column n = 2 of the array is A000079.
Column n = 3 of the array is A000244.
Row k = 2 of the array is A000290.
Row k = 3 of the array is A000578.
Diagonal n = k of the array is A000312.
Diagonal n = k + 1 of the array is A000169.
Diagonal n = k + 2 of the array is A000272.
The transpose of the array is A009999.
The numbers of divisors of the entries are A343656 (row sums: A343657).
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.

Programs

  • Haskell
    a009998 n k = (k + 1) ^ (n - k)
    a009998_row n = a009998_tabl !! n
    a009998_tabl = map reverse a009999_tabl
    -- Reinhard Zumkeller, Feb 02 2014
    
  • Maple
    E := (n,x) -> `if`(n=0,1,x*(1-x)*diff(E(n-1,x),x)+E(n-1,x)*(1+(n-1)*x));
    G := (n,x) -> E(n,x)/(1-x)^(n+1);
    A009998 := (n,k) -> coeff(series(G(n-k,x),x,18),x,k);
    seq(print(seq(A009998(n,k),k=0..n)),n=0..6);
    # Peter Luschny, Aug 02 2010
  • Mathematica
    Flatten[Table[(j+1)^(i-j),{i,0,20},{j,0,i}]] (* Harvey P. Dale, Dec 25 2012 *)
  • PARI
    T(i,j)=(j+1)^(i-j) \\ Charles R Greathouse IV, Feb 06 2017

Formula

T(n,n) = 1; T(n,k) = (k+1)*T(n-1,k) for k=0..n-1. - Reinhard Zumkeller, Feb 02 2014
T(n,m) = (m+1)*Sum_{k=0..n-m}((n+1)^(k-1)*(n-m)^(n-m-k)*(-1)^(n-m-k)*binomial(n-m-1,k-1)). - Vladimir Kruchinin, Sep 12 2015

Extensions

a(62) corrected to 512 by T. D. Noe, Dec 20 2007

A077592 Table by antidiagonals of tau_k(n), the k-th Piltz function (see A007425), or n-th term of the sequence resulting from applying the inverse Möbius transform (k-1) times to the all-ones sequence.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 6, 2, 1, 1, 6, 5, 10, 3, 4, 1, 1, 7, 6, 15, 4, 9, 2, 1, 1, 8, 7, 21, 5, 16, 3, 4, 1, 1, 9, 8, 28, 6, 25, 4, 10, 3, 1, 1, 10, 9, 36, 7, 36, 5, 20, 6, 4, 1, 1, 11, 10, 45, 8, 49, 6, 35, 10, 9, 2, 1, 1, 12, 11, 55, 9, 64, 7, 56, 15, 16, 3, 6, 1
Offset: 1

Views

Author

Henry Bottomley, Nov 08 2002

Keywords

Comments

As an array with offset n=0, k=1, also the number of length n chains of divisors of k. - Gus Wiseman, Aug 04 2022

Examples

			T(6,3) = 9 because we have: 1*1*6, 1*2*3, 1*3*2, 1*6*1, 2*1*3, 2*3*1, 3*1*2, 3*2*1, 6*1*1. - _Geoffrey Critzer_, Feb 16 2015
From _Gus Wiseman_, May 03 2021: (Start)
Array begins:
       k=1 k=2 k=3 k=4 k=5 k=6 k=7 k=8
  n=0:  1   1   1   1   1   1   1   1
  n=1:  1   2   2   3   2   4   2   4
  n=2:  1   3   3   6   3   9   3  10
  n=3:  1   4   4  10   4  16   4  20
  n=4:  1   5   5  15   5  25   5  35
  n=5:  1   6   6  21   6  36   6  56
  n=6:  1   7   7  28   7  49   7  84
  n=7:  1   8   8  36   8  64   8 120
  n=8:  1   9   9  45   9  81   9 165
The triangular form T(n,k) = A(n-k,k) gives the number of length n - k chains of divisors of k. It begins:
  1
  1  1
  1  2  1
  1  3  2  1
  1  4  3  3  1
  1  5  4  6  2  1
  1  6  5 10  3  4  1
  1  7  6 15  4  9  2  1
  1  8  7 21  5 16  3  4  1
  1  9  8 28  6 25  4 10  3  1
  1 10  9 36  7 36  5 20  6  4  1
  1 11 10 45  8 49  6 35 10  9  2  1
(End)
		

Crossrefs

Columns include (with multiplicity and some offsets) A000012, A000027, A000027, A000217, A000027, A000290, A000027, A000292, A000217, A000290, A000027, A002411, A000027, A000290, A000290, A000332 etc.
Cf. A077593.
Row n = 2 of the array is A007425.
Row n = 3 of the array is A007426.
Row n = 4 of the array is A061200.
The diagonal n = k of the array (central column of the triangle) is A163767.
The transpose of the array is A334997.
Diagonal n = k of the array is A343939.
Antidiagonal sums of the array (or row sums of the triangle) are A343940.
A067824(n) counts strict chains of divisors starting with n.
A074206(n) counts strict chains of divisors from n to 1.
A146291(n,k) counts divisors of n with k prime factors (with multiplicity).
A251683(n,k) counts strict length k + 1 chains of divisors from n to 1.
A253249(n) counts nonempty chains of divisors of n.
A334996(n,k) counts strict length k chains of divisors from n to 1.
A337255(n,k) counts strict length k chains of divisors starting with n.

Programs

  • Maple
    with(numtheory):
    A:= proc(n,k) option remember; `if`(k=1, 1,
          add(A(d, k-1), d=divisors(n)))
        end:
    seq(seq(A(n, 1+d-n), n=1..d), d=1..14);  # Alois P. Heinz, Feb 25 2015
  • Mathematica
    tau[n_, 1] = 1; tau[n_, k_] := tau[n, k] = Plus @@ (tau[ #, k - 1] & /@ Divisors[n]); Table[tau[n - k + 1, k], {n, 14}, {k, n, 1, -1}] // Flatten (* Robert G. Wilson v *)
    tau[1, k_] := 1; tau[n_, k_] := Times @@ (Binomial[Last[#] + k - 1, k - 1] & /@ FactorInteger[n]); Table[tau[k, n - k + 1], {n, 1, 13}, {k, 1, n}] // Flatten (* Amiram Eldar, Sep 13 2020 *)
    Table[Length[Select[Tuples[Divisors[k],n-k],And@@Divisible@@@Partition[#,2,1]&]],{n,12},{k,1,n}] (* TRIANGLE, Gus Wiseman, May 03 2021 *)
    Table[Length[Select[Tuples[Divisors[k],n-1],And@@Divisible@@@Partition[#,2,1]&]],{n,6},{k,6}] (* ARRAY, Gus Wiseman, May 03 2021 *)

Formula

If n = Product_i p_i^e_i, then T(n,k) = Product_i C(k+e_i-1, e_i). T(n,k) = Sum_d{d|n} T(n-1,d) = A077593(n,k) - A077593(n-1,k).
Columns are multiplicative.
Dirichlet g.f. for column k: Zeta(s)^k. - Geoffrey Critzer, Feb 16 2015
A(n,k) = A334997(k,n). - Gus Wiseman, Aug 04 2022

Extensions

Typo in formula fixed by Geoffrey Critzer, Feb 16 2015

A340609 Numbers whose number of prime factors (A001222) is divisible by their greatest prime index (A061395).

Original entry on oeis.org

2, 4, 6, 8, 9, 16, 20, 24, 30, 32, 36, 45, 50, 54, 56, 64, 75, 81, 84, 96, 125, 126, 128, 140, 144, 160, 176, 189, 196, 210, 216, 240, 256, 264, 294, 315, 324, 350, 360, 384, 396, 400, 416, 440, 441, 486, 490, 512, 525, 540, 576, 594, 600, 616, 624, 660, 686
Offset: 1

Views

Author

Gus Wiseman, Jan 27 2021

Keywords

Comments

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.
If n is a term, then so is n^k for k > 1. - Robert Israel, Feb 08 2021

Examples

			The sequence of terms together with their prime indices begins:
      2: {1}             64: {1,1,1,1,1,1}      216: {1,1,1,2,2,2}
      4: {1,1}           75: {2,3,3}            240: {1,1,1,1,2,3}
      6: {1,2}           81: {2,2,2,2}          256: {1,1,1,1,1,1,1,1}
      8: {1,1,1}         84: {1,1,2,4}          264: {1,1,1,2,5}
      9: {2,2}           96: {1,1,1,1,1,2}      294: {1,2,4,4}
     16: {1,1,1,1}      125: {3,3,3}            315: {2,2,3,4}
     20: {1,1,3}        126: {1,2,2,4}          324: {1,1,2,2,2,2}
     24: {1,1,1,2}      128: {1,1,1,1,1,1,1}    350: {1,3,3,4}
     30: {1,2,3}        140: {1,1,3,4}          360: {1,1,1,2,2,3}
     32: {1,1,1,1,1}    144: {1,1,1,1,2,2}      384: {1,1,1,1,1,1,1,2}
     36: {1,1,2,2}      160: {1,1,1,1,1,3}      396: {1,1,2,2,5}
     45: {2,2,3}        176: {1,1,1,1,5}        400: {1,1,1,1,3,3}
     50: {1,3,3}        189: {2,2,2,4}          416: {1,1,1,1,1,6}
     54: {1,2,2,2}      196: {1,1,4,4}          440: {1,1,1,3,5}
     56: {1,1,1,4}      210: {1,2,3,4}          441: {2,2,4,4}
		

Crossrefs

Note: Heinz numbers are given in parentheses below.
The case of equality is A047993 (A106529).
These are the Heinz numbers of certain partitions counted by A168659.
The reciprocal version is A340610, with strict case A340828 (A340856).
If all parts (not just the greatest) are divisors we get A340693 (A340606).
A001222 counts prime factors.
A006141 counts partitions whose length equals their minimum (A324522).
A056239 adds up prime indices.
A061395 selects the maximum prime index.
A067538 counts partitions whose length divides their sum (A316413).
A067538 counts partitions whose maximum divides their sum (A326836).
A112798 lists the prime indices of each positive integer.
A200750 counts partitions with length coprime to maximum (A340608).

Programs

  • Maple
    filter:= proc(n) local F,m,g,t;
      F:= ifactors(n)[2];
      m:= add(t[2],t=F);
      g:= numtheory:-pi(max(seq(t[1],t=F)));
      m mod g = 0;
    end proc:
    seelect(filter, [$2..1000]); # Robert Israel, Feb 08 2021
  • Mathematica
    Select[Range[2,100],Divisible[PrimeOmega[#],PrimePi[FactorInteger[#][[-1,1]]]]&]

Formula

A061395(a(n)) divides A001222(a(n)).

A343656 Array read by antidiagonals where A(n,k) is the number of divisors of n^k.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 5, 2, 1, 1, 6, 5, 7, 3, 4, 1, 1, 7, 6, 9, 4, 9, 2, 1, 1, 8, 7, 11, 5, 16, 3, 4, 1, 1, 9, 8, 13, 6, 25, 4, 7, 3, 1, 1, 10, 9, 15, 7, 36, 5, 10, 5, 4, 1, 1, 11, 10, 17, 8, 49, 6, 13, 7, 9, 2, 1, 1, 12, 11, 19, 9, 64, 7, 16, 9, 16, 3, 6, 1
Offset: 1

Views

Author

Gus Wiseman, Apr 28 2021

Keywords

Comments

First differs from A343658 at A(4,2) = 5, A343658(4,2) = 6.
As a triangle, T(n,k) = number of divisors of k^(n-k).

Examples

			Array begins:
       k=0 k=1 k=2 k=3 k=4 k=5 k=6 k=7
  n=1:  1   1   1   1   1   1   1   1
  n=2:  1   2   3   4   5   6   7   8
  n=3:  1   2   3   4   5   6   7   8
  n=4:  1   3   5   7   9  11  13  15
  n=5:  1   2   3   4   5   6   7   8
  n=6:  1   4   9  16  25  36  49  64
  n=7:  1   2   3   4   5   6   7   8
  n=8:  1   4   7  10  13  16  19  22
  n=9:  1   3   5   7   9  11  13  15
Triangle begins:
  1
  1  1
  1  2  1
  1  3  2  1
  1  4  3  3  1
  1  5  4  5  2  1
  1  6  5  7  3  4  1
  1  7  6  9  4  9  2  1
  1  8  7 11  5 16  3  4  1
  1  9  8 13  6 25  4  7  3  1
  1 10  9 15  7 36  5 10  5  4  1
  1 11 10 17  8 49  6 13  7  9  2  1
  1 12 11 19  9 64  7 16  9 16  3  6  1
  1 13 12 21 10 81  8 19 11 25  4 15  2  1
For example, row n = 8 counts the following divisors:
  1  64  243  256  125  36  7  1
     32  81   128  25   18  1
     16  27   64   5    12
     8   9    32   1    9
     4   3    16        6
     2   1    8         4
     1        4         3
              2         2
              1         1
		

Crossrefs

Columns k=1..9 of the array give A000005, A048691, A048785, A344327, A344328, A344329, A343526, A344335, A344336.
Row n = 6 of the array is A000290.
Diagonal n = k of the array is A062319.
Array antidiagonal sums (row sums of the triangle) are A343657.
Dominated by A343658.
A000312 = n^n.
A007318 counts k-sets of elements of {1..n}.
A009998(n,k) = n^k (as an array, offset 1).
A059481 counts k-multisets of elements of {1..n}.

Programs

  • Mathematica
    Table[DivisorSigma[0,k^(n-k)],{n,10},{k,n}]
  • PARI
    A(n, k) = numdiv(n^k); \\ Seiichi Manyama, May 15 2021

Formula

A(n,k) = A000005(A009998(n,k)), where A009998(n,k) = n^k is the interpretation as an array.
A(n,k) = Sum_{d|n} k^omega(d). - Seiichi Manyama, May 15 2021

A298426 Regular triangle where T(n,k) is number of k-ary rooted trees with n nodes.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 2, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 3, 0, 1, 0, 0, 0, 1, 0, 1, 0, 2, 0, 0, 0, 0, 0, 1, 0, 1, 6, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 11, 4, 2, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 23, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1
Offset: 1

Views

Author

Gus Wiseman, Jan 19 2018

Keywords

Comments

Row sums are A298422.

Examples

			Triangle begins:
1
0  1
0  1  1
0  1  0  1
0  1  1  0  1
0  1  0  0  0  1
0  1  2  1  0  0  1
0  1  0  0  0  0  0  1
0  1  3  0  1  0  0  0  1
0  1  0  2  0  0  0  0  0  1
0  1  6  0  0  1  0  0  0  0  1
0  1  0  0  0  0  0  0  0  0  0  1
0  1  11 4  2  0  1  0  0  0  0  0  1
0  1  0  0  0  0  0  0  0  0  0  0  0  1
0  1  23 0  0  0  0  1  0  0  0  0  0  0  1
0  1  0  8  0  2  0  0  0  0  0  0  0  0  0  1
		

Crossrefs

Programs

  • Mathematica
    nn=16;
    arut[n_,k_]:=If[n===1,{{}},Join@@Function[c,Union[Sort/@Tuples[arut[#,k]&/@c]]]/@Select[IntegerPartitions[n-1],Length[#]===k&]]
    Table[arut[n,k]//Length,{n,nn},{k,0,n-1}]

A340692 Number of integer partitions of n of odd rank.

Original entry on oeis.org

0, 0, 2, 0, 4, 2, 8, 4, 14, 12, 26, 22, 44, 44, 76, 78, 126, 138, 206, 228, 330, 378, 524, 602, 814, 950, 1252, 1466, 1900, 2238, 2854, 3362, 4236, 5006, 6232, 7356, 9078, 10720, 13118, 15470, 18800, 22152, 26744, 31456, 37772, 44368, 53002, 62134, 73894
Offset: 0

Views

Author

Gus Wiseman, Jan 29 2021

Keywords

Comments

The Dyson rank of a nonempty partition is its maximum part minus its length. The rank of an empty partition is undefined.

Examples

			The a(0) = 0 through a(9) = 12 partitions (empty columns indicated by dots):
  .  .  (2)   .  (4)     (32)   (6)       (52)     (8)         (54)
        (11)     (31)    (221)  (33)      (421)    (53)        (72)
                 (211)          (51)      (3211)   (71)        (432)
                 (1111)         (222)     (22111)  (422)       (441)
                                (411)              (431)       (621)
                                (3111)             (611)       (3222)
                                (21111)            (3221)      (3321)
                                (111111)           (3311)      (5211)
                                                   (5111)      (22221)
                                                   (22211)     (42111)
                                                   (41111)     (321111)
                                                   (311111)    (2211111)
                                                   (2111111)
                                                   (11111111)
		

Crossrefs

Note: A-numbers of Heinz-number sequences are in parentheses below.
The case of length/maximum instead of rank is A027193 (A026424/A244991).
The case of odd positive rank is A101707 is (A340604).
The strict case is A117193.
The even version is A340601 (A340602).
The Heinz numbers of these partitions are (A340603).
A072233 counts partitions by sum and length.
A168659 counts partitions whose length is divisible by maximum.
A200750 counts partitions whose length and maximum are relatively prime.
- Rank -
A047993 counts partitions of rank 0 (A106529).
A063995/A105806 count partitions by Dyson rank.
A064173 counts partitions of positive/negative rank (A340787/A340788).
A064174 counts partitions of nonpositive/nonnegative rank (A324521/A324562).
A101198 counts partitions of rank 1 (A325233).
A101708 counts partitions of even positive rank (A340605).
A257541 gives the rank of the partition with Heinz number n.
A324520 counts partitions with rank equal to least part (A324519).
- Odd -
A000009 counts partitions into odd parts (A066208).
A026804 counts partitions whose least part is odd.
A058695 counts partitions of odd numbers (A300063).
A067659 counts strict partitions of odd length (A030059).
A160786 counts odd-length partitions of odd numbers (A300272).
A339890 counts factorizations of odd length.
A340385 counts partitions of odd length and maximum (A340386).

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],OddQ[Max[#]-Length[#]]&]],{n,0,30}]

Formula

Having odd rank is preserved under conjugation, and self-conjugate partitions cannot have odd rank, so a(n) = 2*A101707(n) for n > 0.

A340852 Numbers that can be factored in such a way that every factor is a divisor of the number of factors.

Original entry on oeis.org

1, 4, 16, 27, 32, 64, 96, 128, 144, 192, 216, 256, 288, 324, 432, 486, 512, 576, 648, 729, 864, 972, 1024, 1296, 1458, 1728, 1944, 2048, 2560, 2592, 2916, 3125, 3888, 4096, 5120, 5184, 5832, 6144, 6400, 7776, 8192, 9216, 11664, 12288, 12800, 13824, 15552
Offset: 1

Views

Author

Gus Wiseman, Feb 04 2021

Keywords

Comments

Also numbers that can be factored in such a way that the length is divisible by the least common multiple.

Examples

			The sequence of terms together with their prime indices begins:
    1: {}
    4: {1,1}
   16: {1,1,1,1}
   27: {2,2,2}
   32: {1,1,1,1,1}
   64: {1,1,1,1,1,1}
   96: {1,1,1,1,1,2}
  128: {1,1,1,1,1,1,1}
  144: {1,1,1,1,2,2}
  192: {1,1,1,1,1,1,2}
  216: {1,1,1,2,2,2}
  256: {1,1,1,1,1,1,1,1}
  288: {1,1,1,1,1,2,2}
  324: {1,1,2,2,2,2}
  432: {1,1,1,1,2,2,2}
For example, 24576 has three suitable factorizations:
  (2*2*2*2*2*2*2*2*2*2*2*12)
  (2*2*2*2*2*2*2*2*2*2*4*6)
  (2*2*2*2*2*2*2*2*2*3*4*4)
so is in the sequence.
		

Crossrefs

Partitions of this type are counted by A340693 (A340606).
These factorizations are counted by A340851.
The reciprocal version is A340853.
A143773 counts partitions whose parts are multiples of the number of parts.
A320911 can be factored into squarefree semiprimes.
A340597 have an alt-balanced factorization.
A340656 lack a twice-balanced factorization, complement A340657.
- Factorizations -
A001055 counts factorizations, with strict case A045778.
A316439 counts factorizations by product and length.
A339846 counts factorizations of even length.
A339890 counts factorizations of odd length.
A340101 counts factorizations into odd factors, odd-length case A340102.
A340653 counts balanced factorizations.
A340831/A340832 count factorizations with odd maximum/minimum.
A340785 counts factorizations into even numbers, even-length case A340786.
A340854 cannot be factored with odd least factor, complement A340855.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Select[Range[1000],Select[facs[#],And@@IntegerQ/@(Length[#]/#)&]!={}&]

A343658 Array read by antidiagonals where A(n,k) is the number of ways to choose a multiset of k divisors of n.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 6, 2, 1, 1, 6, 5, 10, 3, 4, 1, 1, 7, 6, 15, 4, 10, 2, 1, 1, 8, 7, 21, 5, 20, 3, 4, 1, 1, 9, 8, 28, 6, 35, 4, 10, 3, 1, 1, 10, 9, 36, 7, 56, 5, 20, 6, 4, 1, 1, 11, 10, 45, 8, 84, 6, 35, 10, 10, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Apr 29 2021

Keywords

Comments

First differs from A343656 at A(4,2) = 6, A343656(4,2) = 5.
As a triangle, T(n,k) = number of ways to choose a multiset of n - k divisors of k.

Examples

			Array begins:
       k=0 k=1 k=2 k=3 k=4 k=5 k=6 k=7 k=8
  n=1:  1   1   1   1   1   1   1   1   1
  n=2:  1   2   3   4   5   6   7   8   9
  n=3:  1   2   3   4   5   6   7   8   9
  n=4:  1   3   6  10  15  21  28  36  45
  n=5:  1   2   3   4   5   6   7   8   9
  n=6:  1   4  10  20  35  56  84 120 165
  n=7:  1   2   3   4   5   6   7   8   9
  n=8:  1   4  10  20  35  56  84 120 165
  n=9:  1   3   6  10  15  21  28  36  45
Triangle begins:
   1
   1   1
   1   2   1
   1   3   2   1
   1   4   3   3   1
   1   5   4   6   2   1
   1   6   5  10   3   4   1
   1   7   6  15   4  10   2   1
   1   8   7  21   5  20   3   4   1
   1   9   8  28   6  35   4  10   3   1
   1  10   9  36   7  56   5  20   6   4   1
   1  11  10  45   8  84   6  35  10  10   2   1
For example, row n = 6 counts the following multisets:
  {1,1,1,1,1}  {1,1,1,1}  {1,1,1}  {1,1}  {1}  {}
               {1,1,1,2}  {1,1,3}  {1,2}  {5}
               {1,1,2,2}  {1,3,3}  {1,4}
               {1,2,2,2}  {3,3,3}  {2,2}
               {2,2,2,2}           {2,4}
                                   {4,4}
Note that for n = 6, k = 4 in the triangle, the two multisets {1,4} and {2,2} represent the same divisor 4, so they are only counted once under A343656(4,2) = 5.
		

Crossrefs

Row k = 1 of the array is A000005.
Column n = 4 of the array is A000217.
Column n = 6 of the array is A000292.
Row k = 2 of the array is A184389.
The distinct products of these multisets are counted by A343656.
Antidiagonal sums of the array (or row sums of the triangle) are A343661.
A000312 = n^n.
A009998(n,k) = n^k (as an array, offset 1).
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.

Programs

  • Mathematica
    multchoo[n_,k_]:=Binomial[n+k-1,k];
    Table[multchoo[DivisorSigma[0,k],n-k],{n,10},{k,n}]
  • PARI
    A(n,k) = binomial(numdiv(n) + k - 1, k)
    { for(n=1, 9, for(k=0, 8, print1(A(n,k), ", ")); print ) } \\ Andrew Howroyd, Jan 11 2024

Formula

A(n,k) = ((A000005(n), k)) = A007318(A000005(n) + k - 1, k).
T(n,k) = ((A000005(k), n - k)) = A007318(A000005(k) + n - k - 1, n - k).

A353399 Numbers whose product of prime exponents equals the product of prime shadows of its prime indices.

Original entry on oeis.org

1, 2, 12, 20, 36, 44, 56, 68, 100, 124, 164, 184, 208, 236, 240, 268, 332, 436, 464, 484, 508, 528, 608, 628, 688, 716, 720, 752, 764, 776, 816, 844, 880, 964, 1108, 1132, 1156, 1168, 1200, 1264, 1296, 1324, 1344, 1360, 1412, 1468, 1488, 1584, 1604, 1616, 1724
Offset: 1

Views

Author

Gus Wiseman, May 17 2022

Keywords

Comments

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 the prime shadow A181819(n) to be the product of primes indexed by the exponents in the prime factorization of n. For example, 90 = prime(1)*prime(2)^2*prime(3) has prime shadow prime(1)*prime(2)*prime(1) = 12.

Examples

			The terms together with their prime indices begin:
     1: {}
     2: {1}
    12: {1,1,2}
    20: {1,1,3}
    36: {1,1,2,2}
    44: {1,1,5}
    56: {1,1,1,4}
    68: {1,1,7}
   100: {1,1,3,3}
   124: {1,1,11}
   164: {1,1,13}
   184: {1,1,1,9}
   208: {1,1,1,1,6}
   236: {1,1,17}
   240: {1,1,1,1,2,3}
		

Crossrefs

Product of prime indices is A003963, counted by A339095.
The LHS (product of exponents) is A005361, counted by A266477.
The RHS (product of shadows) is A353394, first appearances A353397.
A related comparison is A353395, counted by A353396.
The partitions are counted by A353398.
Taking indices instead of exponents on the LHS gives A353503.
A001222 counts prime factors with multiplicity, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914.
A130091 lists numbers with distinct prime exponents, counted by A098859.
A181819 gives prime shadow, with an inverse A181821.
A325131 lists numbers relatively prime to their prime shadow.
Numbers divisible by their prime shadow:
- counted by A325702
- listed by A325755
- co-recursive version A325756
- nonprime recursive version A353389
- recursive version A353393
- recursive version counted by A353426

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    red[n_]:=If[n==1,1,Times@@Prime/@Last/@FactorInteger[n]];
    Select[Range[100],Times@@red/@primeMS[#]==Times@@Last/@FactorInteger[#]&]

Formula

A005361(a(n)) = A353394(a(n)).
Previous Showing 11-20 of 53 results. Next