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-10 of 12 results. Next

A003024 Number of acyclic digraphs (or DAGs) with n labeled nodes.

Original entry on oeis.org

1, 1, 3, 25, 543, 29281, 3781503, 1138779265, 783702329343, 1213442454842881, 4175098976430598143, 31603459396418917607425, 521939651343829405020504063, 18676600744432035186664816926721, 1439428141044398334941790719839535103, 237725265553410354992180218286376719253505
Offset: 0

Views

Author

Keywords

Comments

Also the number of n X n real (0,1)-matrices with all eigenvalues positive. - Conjectured by Eric W. Weisstein, Jul 10 2003 and proved by McKay et al. 2003, 2004
Also the number of n X n real (0,1)-matrices with permanent equal to 1, up to permutation of rows/columns, cf. A089482. - Vladeta Jovovic, Oct 28 2009
Also the number of nilpotent elements in the semigroup of binary relations on [n]. - Geoffrey Critzer, May 26 2022
From Gus Wiseman, Jan 01 2024: (Start)
Also the number of sets of n nonempty subsets of {1..n} such that there is a unique way to choose a different element from each. For example, non-isomorphic representatives of the a(3) = 25 set-systems are:
{{1},{2},{3}}
{{1},{2},{1,3}}
{{1},{2},{1,2,3}}
{{1},{1,2},{1,3}}
{{1},{1,2},{2,3}}
{{1},{1,2},{1,2,3}}
These set-systems have ranks A367908, subset of A367906, for multisets A368101.
The version for no ways is A368600, any length A367903, ranks A367907.
The version for at least one way is A368601, any length A367902.
(End)

Examples

			For n = 2 the three (0,1)-matrices are {{{1, 0}, {0, 1}}, {{1, 0}, {1, 1}}, {{1, 1}, {0, 1}}}.
		

References

  • Archer, K., Gessel, I. M., Graves, C., & Liang, X. (2020). Counting acyclic and strong digraphs by descents. Discrete Mathematics, 343(11), 112041.
  • S. R. Finch, Mathematical Constants, Cambridge, 2003, p. 310.
  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 19, Eq. (1.6.1).
  • R. W. Robinson, Counting labeled acyclic digraphs, pp. 239-273 of F. Harary, editor, New Directions in the Theory of Graphs. Academic Press, NY, 1973.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P Stanley, Enumerative Combinatorics I, 2nd. ed., p. 322.

Crossrefs

Cf. A086510, A081064 (refined by # arcs), A307049 (by # descents).
Cf. A055165, which counts nonsingular {0, 1} matrices and A085656, which counts positive definite {0, 1} matrices.
Cf. A188457, A135079, A137435 (acyclic 3-multidigraphs), A188490.
For a unique sink we have A003025.
The unlabeled version is A003087.
These are the reverse-alternating sums of rows of A046860.
The weakly connected case is A082402.
A reciprocal version is A334282.
Row sums of A361718.

Programs

  • Maple
    p:=evalf(solve(sum((-1)^n*x^n/(n!*2^(n*(n-1)/2)), n=0..infinity) = 0, x), 50); M:=evalf(sum((-1)^(n+1)*p^n/((n-1)!*2^(n*(n-1)/2)), n=1..infinity), 40); # program for evaluation of constants p and M in the asymptotic formula, Vaclav Kotesovec, Dec 09 2013
  • Mathematica
    a[0] = a[1] = 1; a[n_] := a[n] = Sum[ -(-1)^k * Binomial[n, k] * 2^(k*(n-k)) * a[n-k], {k, 1, n}]; Table[a[n], {n, 0, 13}](* Jean-François Alcover, May 21 2012, after PARI *)
    Table[2^(n*(n-1)/2)*n! * SeriesCoefficient[1/Sum[(-1)^k*x^k/k!/2^(k*(k-1)/2),{k,0,n}],{x,0,n}],{n,0,20}] (* Vaclav Kotesovec, May 19 2015 *)
    Table[Length[Select[Subsets[Subsets[Range[n]],{n}],Length[Select[Tuples[#],UnsameQ@@#&]]==1&]],{n,0,5}] (* Gus Wiseman, Jan 01 2024 *)
  • PARI
    a(n)=if(n<1,n==0,sum(k=1,n,-(-1)^k*binomial(n,k)*2^(k*(n-k))*a(n-k)))
    
  • PARI
    {a(n)=polcoeff(1-sum(k=0, n-1, a(k)*x^k/(1+2^k*x+x*O(x^n))^(k+1)), n)} \\ Paul D. Hanna, Oct 17 2009

Formula

a(0) = 1; for n > 0, a(n) = Sum_{k=1..n} (-1)^(k+1)*C(n, k)*2^(k*(n-k))*a(n-k).
1 = Sum_{n>=0} a(n)*exp(-2^n*x)*x^n/n!. - Vladeta Jovovic, Jun 05 2005
a(n) = Sum_{k=1..n} (-1)^(n-k)*A046860(n,k) = Sum_{k=1..n} (-1)^(n-k)*k!*A058843(n,k). - Vladeta Jovovic, Jun 20 2008
1 = Sum_{n=>0} a(n)*x^n/(1 + 2^n*x)^(n+1). - Paul D. Hanna, Oct 17 2009
1 = Sum_{n>=0} a(n)*C(n+m-1,n)*x^n/(1 + 2^n*x)^(n+m) for m>=1. - Paul D. Hanna, Apr 01 2011
log(1+x) = Sum_{n>=1} a(n)*(x^n/n)/(1 + 2^n*x)^n. - Paul D. Hanna, Apr 01 2011
Let E(x) = Sum_{n >= 0} x^n/(n!*2^C(n,2)). Then a generating function for this sequence is 1/E(-x) = Sum_{n >= 0} a(n)*x^n/(n!*2^C(n,2)) = 1 + x + 3*x^2/(2!*2) + 25*x^3/(3!*2^3) + 543*x^4/(4!*2^6) + ... (Stanley). Cf. A188457. - Peter Bala, Apr 01 2013
a(n) ~ n!*2^(n*(n-1)/2)/(M*p^n), where p = 1.488078545599710294656246... is the root of the equation Sum_{n>=0} (-1)^n*p^n/(n!*2^(n*(n-1)/2)) = 0, and M = Sum_{n>=1} (-1)^(n+1)*p^n/((n-1)!*2^(n*(n-1)/2)) = 0.57436237330931147691667... Both references to the article "Acyclic digraphs and eigenvalues of (0,1)-matrices" give the wrong value M=0.474! - Vaclav Kotesovec, Dec 09 2013 [Response from N. J. A. Sloane, Dec 11 2013: The value 0.474 has a typo, it should have been 0.574. The value was taken from Stanley's 1973 paper.]
exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + x + 2*x^2 + 10*x^3 + 146*x^4 + 6010*x^5 + ... appears to have integer coefficients (cf. A188490). - Peter Bala, Jan 14 2016

A368100 Numbers of which it is possible to choose a different prime factor of each prime index.

Original entry on oeis.org

1, 3, 5, 7, 11, 13, 15, 17, 19, 23, 29, 31, 33, 35, 37, 39, 41, 43, 47, 51, 53, 55, 59, 61, 65, 67, 69, 71, 73, 77, 79, 83, 85, 87, 89, 91, 93, 95, 97, 101, 103, 107, 109, 111, 113, 119, 123, 127, 129, 131, 137, 139, 141, 143, 145, 149, 151, 155, 157, 161, 163
Offset: 1

Views

Author

Gus Wiseman, Dec 12 2023

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.

Examples

			The prime indices of 2849 are {4,5,12}, with prime factors {{2,2},{5},{2,2,3}}, and of the two choices (2,5,2) and (2,5,3) the latter has all different terms, so 2849 is in the sequence.
The terms together with their prime indices of prime indices begin:
   1: {}
   3: {{1}}
   5: {{2}}
   7: {{1,1}}
  11: {{3}}
  13: {{1,2}}
  15: {{1},{2}}
  17: {{4}}
  19: {{1,1,1}}
  23: {{2,2}}
  29: {{1,3}}
  31: {{5}}
  33: {{1},{3}}
  35: {{2},{1,1}}
  37: {{1,1,2}}
  39: {{1},{1,2}}
		

Crossrefs

The complement is A355529, odd A355535, binary A367907.
Positions of positive terms in A367771.
The version for binary indices is A367906, positive positions in A367905.
For a unique choice we have A368101, binary A367908.
The version for divisors instead of factors is A368110, complement A355740.
A058891 counts set-systems, covering A003465, connected A323818.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100], Select[Tuples[prix/@prix[#]], UnsameQ@@#&]!={}&]

A367908 Numbers n such that there is only one way to choose a different binary index of each binary index of n.

Original entry on oeis.org

1, 2, 3, 5, 6, 8, 9, 10, 11, 13, 14, 17, 19, 21, 22, 24, 26, 28, 34, 35, 37, 38, 40, 41, 44, 49, 50, 56, 67, 69, 70, 73, 74, 81, 88, 98, 104, 128, 129, 130, 131, 133, 134, 136, 137, 138, 139, 141, 142, 145, 147, 149, 150, 152, 154, 156, 162, 163, 165, 166, 168
Offset: 1

Views

Author

Gus Wiseman, Dec 11 2023

Keywords

Comments

Also BII-numbers of set-systems (sets of nonempty sets) satisfying a strict version of the axiom of choice in exactly one way.
A binary index of n (row n of A048793) is any position of a 1 in its reversed binary expansion. A set-system is a finite set of finite nonempty sets. We define the set-system with BII-number n to be obtained by taking the binary indices of each binary index of n. Every finite set of finite nonempty sets has a different BII-number. For example, 18 has reversed binary digits (0,1,0,0,1), and since the binary indices of 2 and 5 are {2} and {1,3} respectively, the BII-number of {{2},{1,3}} is 18.
The axiom of choice says that, given any set of nonempty sets Y, it is possible to choose a set containing an element from each. The strict version requires this set to have the same cardinality as Y, meaning no element is chosen more than once.

Examples

			The set-system {{1},{1,2},{1,3}} with BII-number 21 satisfies the axiom in exactly one way, namely (1,2,3), so 21 is in the sequence.
The terms together with the corresponding set-systems begin:
   1: {{1}}
   2: {{2}}
   3: {{1},{2}}
   5: {{1},{1,2}}
   6: {{2},{1,2}}
   8: {{3}}
   9: {{1},{3}}
  10: {{2},{3}}
  11: {{1},{2},{3}}
  13: {{1},{1,2},{3}}
  14: {{2},{1,2},{3}}
  17: {{1},{1,3}}
  19: {{1},{2},{1,3}}
  21: {{1},{1,2},{1,3}}
  22: {{2},{1,2},{1,3}}
		

Crossrefs

These set-systems are counted by A367904.
Positions of 1's in A367905, firsts A367910, sorted firsts A367911.
If there is at least one choice we get A367906, counted by A367902.
If there are no choices we get A367907, counted by A367903.
If there are multiple choices we get A367909, counted by A367772.
The version for MM-numbers of multiset partitions is A368101.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A058891 counts set-systems, covering A003465, connected A323818.
A059201 counts covering T_0 set-systems.
A070939 gives length of binary expansion.
A096111 gives product of binary indices.
A326031 gives weight of the set-system with BII-number n.
A368098 counts unlabeled multiset partitions for axiom, complement A368097.
BII-numbers: A309314 (hyperforests), A326701 (set partitions), A326703 (chains), A326704 (antichains), A326749 (connected), A326750 (clutters), A326751 (blobs), A326752 (hypertrees), A326754 (covers), A326783 (uniform), A326784 (regular), A326788 (simple), A330217 (achiral).

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Select[Range[100], Length[Select[Tuples[bpe/@bpe[#]], UnsameQ@@#&]]==1&]
  • Python
    from itertools import count, islice, product
    def bin_i(n): #binary indices
        return([(i+1) for i, x in enumerate(bin(n)[2:][::-1]) if x =='1'])
    def a_gen(): #generator of terms
        for n in count(1):
            p = list(product(*[bin_i(k) for k in bin_i(n)]))
            x,c = len(p),0
            for j in range(x):
                if len(set(p[j])) == len(p[j]): c += 1
                if j+1 == x and c == 1: yield(n)
    A367908_list = list(islice(a_gen(), 100)) # John Tyler Rascoe, Feb 10 2024

Formula

A367771 Number of ways to choose a different prime index of each prime index of 2n + 1.

Original entry on oeis.org

1, 1, 1, 2, 0, 1, 2, 1, 1, 3, 0, 2, 0, 0, 2, 1, 1, 2, 3, 1, 1, 2, 0, 2, 0, 1, 4, 1, 0, 1, 3, 0, 1, 1, 2, 3, 2, 0, 2, 2, 0, 1, 1, 1, 4, 2, 1, 3, 2, 0, 2, 3, 0, 3, 1, 1, 3, 0, 0, 2, 0, 1, 0, 1, 1, 5, 0, 0, 2, 2, 2, 2, 2, 0, 2, 4, 0, 1, 1, 0, 4, 2, 1, 2, 2, 0, 4
Offset: 0

Views

Author

Gus Wiseman, Dec 12 2023

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.

Examples

			The prime indices of prime indices of 427 = 2*213 + 1 are {{1,1},{1,2,2}}, with four ways to choose (1,2), so a(213) = 4.
The prime indices of prime indices of 1469 = 2*734 + 1 are {{1,2},{1,2,3}}, with four choices (1,2), (1,3), (2,1), (2,3), so a(734) = 4.
		

Crossrefs

The "extended" version below includes alternating zeros at even positions.
Extended positions of zeros are A355529, binary A367907.
The extended version for binary indices is A367905.
Extended positions of nonzeros are A368100, binary A367906.
Extended positions of ones are A368101, binary A367908.
The extended version without distinctness is A355741, for multisets A355744.
A058891 counts set-systems, covering A003465, connected A323818.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[Tuples[prix/@prix[2n+1]], UnsameQ@@#&]],{n,0,100}]

A370584 Number of subsets of {1..n} such that only one set can be obtained by choosing a different prime factor of each element.

Original entry on oeis.org

1, 1, 2, 4, 6, 12, 18, 36, 48, 68, 104, 208, 284, 568, 888, 1296, 1548, 3096, 3968, 7936, 10736, 15440, 24008, 48016, 58848, 73680, 114368, 132608, 176240, 352480, 449824, 899648, 994976, 1399968, 2160720, 2859584, 3296048, 6592096, 10156672, 14214576, 16892352
Offset: 0

Views

Author

Gus Wiseman, Feb 26 2024

Keywords

Comments

For example, the only choice of a different prime factor of each element of (4,5,6) is (2,5,3).

Examples

			The a(0) = 1 through a(6) = 18 subsets:
  {}  {}  {}   {}     {}     {}       {}
          {2}  {2}    {2}    {2}      {2}
               {3}    {3}    {3}      {3}
               {2,3}  {4}    {4}      {4}
                      {2,3}  {5}      {5}
                      {3,4}  {2,3}    {2,3}
                             {2,5}    {2,5}
                             {3,4}    {2,6}
                             {3,5}    {3,4}
                             {4,5}    {3,5}
                             {2,3,5}  {3,6}
                             {3,4,5}  {4,5}
                                      {4,6}
                                      {2,3,5}
                                      {2,5,6}
                                      {3,4,5}
                                      {3,5,6}
                                      {4,5,6}
		

Crossrefs

For divisors instead of factors we have A051026, cf. A368110, A355740.
The version for set-systems is A367904, ranks A367908.
Multisets of this type are ranked by A368101, cf. A368100, A355529.
For existence we have A370582, differences A370586.
For nonexistence we have A370583, differences A370587.
Maximal sets of this type are counted by A370585.
The version for partitions is A370594, cf. A370592, A370593.
For binary indices instead of factors we have A370638, cf. A370636, A370637.
The version for factorizations is A370645, cf. A368414, A368413.
For unlabeled multiset partitions we have A370646, cf. A368098, A368097.
A006530 gives greatest prime factor, least A020639.
A027746 lists prime factors, A112798 indices, length A001222.
A355741 counts ways to choose a prime factor of each prime index.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]], Length[Union[Sort/@Select[Tuples[If[#==1, {},First/@FactorInteger[#]]&/@#], UnsameQ@@#&]]]==1&]],{n,0,10}]

Extensions

More terms from Jinyuan Wang, Mar 28 2025

A370594 Number of integer partitions of n such that only one set can be obtained by choosing a different prime factor of each part.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 0, 3, 3, 4, 3, 4, 5, 5, 8, 10, 11, 7, 14, 13, 19, 23, 24, 20, 30, 33, 40, 47, 49, 55, 53, 72, 80, 90, 92, 110, 110, 132, 154, 169, 180, 201, 218, 246, 281, 302, 323, 348, 396, 433, 482, 530, 584, 618, 670, 754, 823, 903, 980, 1047, 1137
Offset: 0

Views

Author

Gus Wiseman, Feb 29 2024

Keywords

Examples

			The partition (10,6,4) has unique choice (5,3,2) so is counted under a(20).
The a(0) = 1 through a(12) = 5 partitions:
()  .  (2)  (3)  (4)  (5)    .  (7)    (8)    (9)    (6,4)    (11)   (6,6)
                      (3,2)     (4,3)  (5,3)  (5,4)  (7,3)    (7,4)  (7,5)
                                (5,2)  (6,2)  (6,3)  (5,3,2)  (8,3)  (10,2)
                                              (7,2)           (9,2)  (5,4,3)
                                                                     (7,3,2)
		

Crossrefs

The version for set-systems is A367904, ranks A367908.
Multisets of this type are ranked by A368101, cf. A368100, A355529.
The version for subsets is A370584, cf. A370582, A370583, A370586, A370587.
Maximal sets of this type are counted by A370585.
For existence we have A370592.
For nonexistence we have A370593.
For divisors instead of factors we have A370595.
For subsets and binary indices we have A370638, cf. A370636, A370637.
The version for factorizations is A370645, cf. A368414, A368413.
These partitions have ranks A370647.
A006530 gives greatest prime factor, least A020639.
A027746 lists prime factors, A112798 indices, length A001222.
A355741 counts ways to choose a prime factor of each prime index.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n], Length[Union[Sort/@Select[Tuples[If[#==1, {},First/@FactorInteger[#]]&/@#], UnsameQ@@#&]]]==1&]],{n,0,30}]

A370638 Number of subsets of {1..n} such that a unique set can be obtained by choosing a different binary index of each element.

Original entry on oeis.org

1, 2, 4, 6, 12, 19, 30, 45, 90, 147, 230, 343, 504, 716, 994, 1352, 2704, 4349, 6469, 9162, 12585, 16862, 22122, 28617, 36653, 46431, 58075, 72097, 88456, 107966, 130742, 157647, 315294, 494967, 704753, 950080, 1234301, 1565165, 1945681, 2387060, 2890368, 3470798
Offset: 0

Views

Author

Gus Wiseman, Mar 09 2024

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.

Examples

			The set {3,4} has binary indices {{1,2},{3}}, with two choices {1,3}, {2,3}, so is not counted under a(4).
The a(0) = 1 through a(5) = 19 subsets:
  {}  {}   {}     {}     {}       {}
      {1}  {1}    {1}    {1}      {1}
           {2}    {2}    {2}      {2}
           {1,2}  {1,2}  {4}      {4}
                  {1,3}  {1,2}    {1,2}
                  {2,3}  {1,3}    {1,3}
                         {1,4}    {1,4}
                         {2,3}    {1,5}
                         {2,4}    {2,3}
                         {1,2,4}  {2,4}
                         {1,3,4}  {4,5}
                         {2,3,4}  {1,2,4}
                                  {1,2,5}
                                  {1,3,4}
                                  {1,3,5}
                                  {2,3,4}
                                  {2,3,5}
                                  {2,4,5}
                                  {3,4,5}
		

Crossrefs

Set systems of this type are counted by A367904, ranks A367908.
A version for MM-numbers of multisets is A368101.
For prime indices we have A370584.
This is the unique version of A370636, complement A370637.
The maximal case is A370640, differences A370641.
Factorizations of this type are counted by A370645.
The case A370818 is the restriction to A000225.
A048793 lists binary indices, A000120 length, A272020 reverse, A029931 sum.
A058891 counts set-systems, A003465 covering, A323818 connected.
A070939 gives length of binary expansion.
A096111 gives product of binary indices.

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Table[Length[Select[Subsets[Range[n]],Length[Union[Sort /@ Select[Tuples[bpe/@#],UnsameQ@@#&]]]==1&]],{n,0,10}]

Formula

a(2^n - 1) = A370818(n).

Extensions

More terms from Jinyuan Wang, Mar 28 2025

A382201 MM-numbers of sets of sets with distinct sums.

Original entry on oeis.org

1, 2, 3, 5, 6, 10, 11, 13, 15, 17, 22, 26, 29, 30, 31, 33, 34, 39, 41, 43, 47, 51, 55, 58, 59, 62, 65, 66, 67, 73, 78, 79, 82, 83, 85, 86, 87, 93, 94, 101, 102, 109, 110, 113, 118, 123, 127, 129, 130, 134, 137, 139, 141, 145, 146, 149, 155, 157, 158, 163, 165
Offset: 1

Views

Author

Gus Wiseman, Mar 21 2025

Keywords

Comments

First differs from A302494 in lacking 143, corresponding to the multiset partition {{1,2},{3}}.
Also products of prime numbers of squarefree index such that the factors all have distinct sums of prime indices.
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, sum A056239. The multiset of multisets with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.

Examples

			The terms together with their prime indices of prime indices begin:
   1: {}
   2: {{}}
   3: {{1}}
   5: {{2}}
   6: {{},{1}}
  10: {{},{2}}
  11: {{3}}
  13: {{1,2}}
  15: {{1},{2}}
  17: {{4}}
  22: {{},{3}}
  26: {{},{1,2}}
  29: {{1,3}}
  30: {{},{1},{2}}
  31: {{5}}
  33: {{1},{3}}
  34: {{},{4}}
  39: {{1},{1,2}}
		

Crossrefs

Set partitions of this type are counted by A275780.
Twice-partitions of this type are counted by A279785.
For just sets of sets we have A302478.
For distinct blocks instead of block-sums we have A302494.
For equal instead of distinct sums we have A302497.
For just distinct sums we have A326535.
For normal multiset partitions see A326519, A326533, A326537, A381718.
Factorizations of this type are counted by A381633. See also A001055, A045778, A050320, A050326, A321455, A321469, A382080.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],And@@SquareFreeQ/@prix[#]&&UnsameQ@@Total/@prix/@prix[#]&]

Formula

Equals A302478 /\ A326535.

A370645 Number of integer factorizations of n into unordered factors > 1 such that only one set can be obtained by choosing a different prime factor of each factor.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 4, 1, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 3, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 5, 1, 1, 2, 2, 1, 1, 1, 4, 1, 1, 1, 3, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Mar 01 2024

Keywords

Comments

All of these factorizations are co-balanced (A340596).

Examples

			The factorization f = (3*6*10) has prime factor choices (3,2,2), (3,3,2), (3,2,5), and (3,3,5), of which only (3,2,5) has all different parts, so f is counted under a(180).
The a(n) factorizations for n = 2, 12, 24, 36, 72, 120, 144, 180, 288:
  (2)  (2*6)  (3*8)   (4*9)   (8*9)   (3*5*8)   (2*72)   (4*5*9)   (3*96)
       (3*4)  (4*6)   (6*6)   (2*36)  (4*5*6)   (3*48)   (5*6*6)   (4*72)
              (2*12)  (2*18)  (3*24)  (2*3*20)  (4*36)   (2*3*30)  (6*48)
                      (3*12)  (4*18)  (2*5*12)  (6*24)   (2*5*18)  (8*36)
                              (6*12)  (2*6*10)  (8*18)   (2*6*15)  (9*32)
                                      (3*4*10)  (9*16)   (2*9*10)  (12*24)
                                                (12*12)  (3*4*15)  (16*18)
                                                         (3*5*12)  (2*144)
                                                         (3*6*10)
		

Crossrefs

Multisets of this type are ranked by A368101, see also A368100, A355529.
For nonexistence we have A368413, complement A368414.
Subsets of this type are counted by A370584, see also A370582, A370583.
Maximal sets of this type are counted by A370585.
The version for partitions is A370594, see also A370592, A370593.
Subsets of this type are counted by A370638, see also A370636, A370637.
For unlabeled multiset partitions we have A370646, also A368098, A368097.
A001055 counts factorizations, strict A045778.
A006530 gives greatest prime factor, least A020639.
A007716 counts non-isomorphic multiset partitions, connected A007718.
A027746 lists prime factors, A112798 indices, length A001222.
A058891 counts set-systems, unlabeled A000612, connected A323818.
A283877 counts non-isomorphic set-systems, connected A300913.
A355741 counts ways to choose a prime factor of each prime index.
For set-systems see A367902-A367908.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join @@ Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]], {d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n], Length[Union[Sort/@Select[Tuples[First /@ FactorInteger[#]&/@#], UnsameQ@@#&]]]==1&]],{n,100}]

A370647 Numbers such that only one set can be obtained by choosing a different prime factor of each prime index.

Original entry on oeis.org

1, 3, 5, 7, 11, 15, 17, 19, 23, 31, 33, 35, 39, 41, 51, 53, 55, 59, 65, 67, 69, 77, 83, 85, 87, 91, 93, 95, 97, 103, 109, 111, 119, 123, 127, 129, 131, 155, 157, 161, 165, 169, 177, 179, 183, 185, 187, 191, 201, 203, 205, 209, 211, 213, 217, 227, 235, 237, 241
Offset: 1

Views

Author

Gus Wiseman, Mar 06 2024

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.

Examples

			The prime indices of 91 are {4,6}, with only choice {2,3}, so 91 is in the sequence.
The terms together with their prime indices begin:
     1: {}        53: {16}      109: {29}
     3: {2}       55: {3,5}     111: {2,12}
     5: {3}       59: {17}      119: {4,7}
     7: {4}       65: {3,6}     123: {2,13}
    11: {5}       67: {19}      127: {31}
    15: {2,3}     69: {2,9}     129: {2,14}
    17: {7}       77: {4,5}     131: {32}
    19: {8}       83: {23}      155: {3,11}
    23: {9}       85: {3,7}     157: {37}
    31: {11}      87: {2,10}    161: {4,9}
    33: {2,5}     91: {4,6}     165: {2,3,5}
    35: {3,4}     93: {2,11}    169: {6,6}
    39: {2,6}     95: {3,8}     177: {2,17}
    41: {13}      97: {25}      179: {41}
    51: {2,7}    103: {27}      183: {2,18}
		

Crossrefs

For nonexistence we have A355529, count A370593.
For binary instead of prime indices we have A367908, counted by A367904.
For existence we have A368100, count A370592.
For a sequence instead of set of factors we have A368101.
The version for subsets is A370584, see also A370582, A370583.
Maximal sets of this type are counted by A370585.
Partitions of this type are counted by A370594.
For subsets and binary indices we have A370638.
The version for factorizations is A370645, see also A368414, A368413.
For divisors instead of factors we have A370810, counted by A370595.
A006530 gives greatest prime factor, least A020639.
A027746 lists prime factors, A112798 indices, length A001222.
A355741 counts ways to choose a prime factor of each prime index.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Length[Union[Sort /@ Select[Tuples[prix/@prix[#]],UnsameQ@@#&]]]==1&]
Showing 1-10 of 12 results. Next