A303362 -id:A303362 - cp's OEIS frontend

A316467 Matula-Goebel numbers of locally stable rooted identity trees, meaning no branch is a subset of any other branch of the same root.

1, 2, 3, 5, 11, 15, 31, 33, 47, 55, 93, 127, 137, 141, 155, 165, 211, 257, 341, 381, 411, 465, 487, 633, 635, 709, 771, 773, 811, 907, 977, 1023, 1055, 1285, 1297, 1397, 1457, 1461, 1507, 1621, 1705, 1905, 2127, 2293, 2319, 2321, 2433, 2621, 2721, 2833, 2931

Offset: 1

Gus Wiseman, Jul 04 2018

nonn

A prime index of n is a number m such that prime(m) divides n. A number belongs to this sequence iff it is squarefree, its distinct prime indices are pairwise indivisible, and its prime indices also belong to this sequence.

			165 = prime(2)*prime(3)*prime(5) belongs to the sequence because it is squarefree, the indices {2,3,5} are pairwise indivisible, and each of them already belongs to the sequence.
Sequence of locally stable rooted identity trees preceded by their Matula-Goebel numbers begins:
    1: o
    2: (o)
    3: ((o))
    5: (((o)))
   11: ((((o))))
   15: ((o)((o)))
   31: (((((o)))))
   33: ((o)(((o))))
   47: (((o)((o))))
   55: (((o))(((o))))
   93: ((o)((((o)))))
  127: ((((((o))))))
  137: (((o)(((o)))))
  141: ((o)((o)((o))))
  155: (((o))((((o)))))
  165: ((o)((o))(((o))))

Cf. A000081, A004111, A007097, A276625, A277098, A285572, A285573, A302796, A303362, A304713, A316468, A316469, A316471, A316474, A316476, A316494.

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
ain[n_]:=And[Select[Tuples[primeMS[n],2],UnsameQ@@#&&Divisible@@#&]=={},SquareFreeQ[n],And@@ain/@primeMS[n]];
Select[Range[100],ain]

A316474 Number of locally stable rooted identity trees with n nodes, meaning no branch is a subset of any other branch of the same root.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 3, 5, 8, 14, 23, 42, 73, 133, 241, 442, 812, 1508, 2802, 5247, 9842, 18554, 35045, 66453, 126249

Offset: 1

Gus Wiseman, Jul 04 2018

			The a(9) = 8 locally stable rooted identity trees:
((((((((o))))))))
(((((o)((o))))))
((((o)(((o))))))
(((o)((((o))))))
((((o))(((o)))))
((o)(((((o))))))
((o)((o)((o))))
(((o))((((o)))))

Gus Wiseman, The a(10) = 14 locally stable rooted identity trees with 10 nodes.

Cf. A000081, A004111, A276625, A277098, A285572, A285573, A303362, A304713, A316467, A316471, A316471, A316475.

strut[n_]:=strut[n]=If[n===1,{{}},Select[Join@@Function[c,Union[Sort/@Tuples[strut/@c]]]/@IntegerPartitions[n-1],UnsameQ@@#&&Select[Tuples[#,2],UnsameQ@@#&&Complement@@#=={}&]=={}&]];
Table[Length[strut[n]],{n,20}]

A316494 Matula-Goebel numbers of locally disjoint rooted identity trees, meaning no branch overlaps any other branch of the same root.

Original entry on oeis.org

1, 2, 3, 5, 6, 10, 11, 13, 15, 22, 26, 29, 30, 31, 33, 41, 47, 55, 58, 62, 66, 79, 82, 93, 94, 101, 109, 110, 113, 123, 127, 137, 141, 143, 145, 155, 158, 165, 179, 186, 202, 205, 211, 218, 226, 246, 254, 257, 271, 274, 282, 286, 290, 293, 310, 317, 327, 330

Offset: 1

Gus Wiseman, Jul 04 2018

nonn

A prime index of n is a number m such that prime(m) divides n. A number is in the sequence iff either it is equal to 1, it is a prime number whose prime index already belongs to the sequence, or its prime indices are pairwise coprime, distinct, and already belong to the sequence.

			The sequence of all locally disjoint rooted identity trees preceded by their Matula-Goebel numbers begins:
   1: o
   2: (o)
   3: ((o))
   5: (((o)))
   6: (o(o))
  10: (o((o)))
  11: ((((o))))
  13: ((o(o)))
  15: ((o)((o)))
  22: (o(((o))))
  26: (o(o(o)))
  29: ((o((o))))
  30: (o(o)((o)))
  31: (((((o)))))

Cf. A000081, A004111, A007097, A276625, A277098, A302696, A303362, A304713, A316467, A316471, A316474, A316495.

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[1000],Or[#==1,And[SquareFreeQ[#],Or[PrimeQ[#],CoprimeQ@@primeMS[#]],And@@#0/@primeMS[#]]]&]

A322437 Number of unordered pairs of factorizations of n into factors > 1 where no factor of one divides any factor of the other.

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, 0, 0, 0, 0, 0, 0, 1, 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, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2

Offset: 1

Gus Wiseman, Dec 08 2018

nonn

First differs from A322438 at a(144) = 3, A322438(144) = 4.
From Antti Karttunen, Dec 11 2020: (Start)
Zeros occur on numbers that are either of the form p^k, or q * p^k, or p*q*r, for some primes p, q, r, and exponent k >= 0. [Note also that in all these cases, when x > 1, A307408(x) = 2+A307409(x) = 2 + (A001222(x) - 1)*A001221(x) = A000005(x)].
Proof:
  It is easy to see that for such numbers it is not possible to obtain two such distinct factorizations, that no factor of the other would not divide some factor of the other.
  Conversely, the complement set of above is formed of such composites n that have at least one unitary divisor that is either of the form
  (1) p^x * q^y, with x, y >= 2,
  or
  (2) p^x * q^y * r^z, with x >= 2, and y, z >= 1,
  or
  (3) p^x * q^y * r^z * s^w, with x, y, z, w >= 1,
  where p, q, r, s are distinct primes. Let's indicate with C the remaining portion of k coprime to p, q, r and s (which could be also 1). Then in case (1) we can construct two factorizations, the first having factors (p*q*C) and (p^(x-1) * q^(y-1)), and the second having factors (p^x * C) and (q^y) that are guaranteed to satisfy the condition that no factor in the other factorization divides any of the factors of the other factorization. For case (2) pairs like {(p * q^y * C), (p^(x-1) * r^z)} and {(p^x * C), (q^y * r^z)}, and for case (3) pairs like {(p^x * q^y * C), (r^z * s^w)} and {(p^x * r^z * C), (q^y * s^w)} offer similar examples, therefore a(n) > 0 for all such cases.
(End)

			The a(120) = 2 pairs of such factorizations:
   (6*20)|(8*15)
   (8*15)|(10*12)
The a(144) = 3 pairs of factorizations:
   (6*24)|(9,16)
   (8*18)|(12*12)
   (9*16)|(12*12)
The a(210) = 3 pairs of factorizations:
   (6*35)|(10*21)
   (6*35)|(14*15)
  (10*21)|(14*15)
[Note that 210 is the first squarefree number obtaining nonzero value]
The a(240) = 4 pairs of factorizations:
   (6*40)|(15*16)
   (8*30)|(12*20)
  (10*24)|(15*16)
  (12*20)|(15*16)
The a(1728) = 14 pairs of factorizations:
    (6*6*48)|(27*64)
   (6*12*24)|(27*64)
     (6*288)|(27*64)
    (8*8*27)|(12*12*12)
  (12*12*12)|(27*64)
  (12*12*12)|(32*54)
    (12*144)|(27*64)
    (12*144)|(32*54)
    (16*108)|(24*72)
     (18*96)|(27*64)
     (24*72)|(27*64)
     (24*72)|(32*54)
     (27*64)|(36*48)
     (32*54)|(36*48)

Antti Karttunen, Table of n, a(n) for n = 1..65537
Index entries for sequences computed from exponents in factorization of n

Cf. A000688, A001055, A285572, A285573, A294068, A303362, A304713, A305149, A305150, A305193, A316476, A317144, A322435, A322436, A322441, A322442, A339626.

Cf. also comments A307408 and A307409.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Select[Subsets[facs[n],{2}],And[!Or@@Divisible@@@Tuples[#],!Or@@Divisible@@@Reverse/@Tuples[#]]&]],{n,100}]

factorizations(n, m=n, f=List([]), z=List([])) = if(1==n, listput(z,Vec(f)); z, my(newf); fordiv(n, d, if((d>1)&&(d<=m), newf = List(f); listput(newf,d); z = factorizations(n/d, d, newf, z))); (z));
is_ndf_pair(fac1,fac2) = { for(i=1,#fac1,for(j=1,#fac2,if(!(fac1[i]%fac2[j])||!(fac2[j]%fac1[i]),return(0)))); (1); };
number_of_ndfpairs(z) = sum(i=1,#z,sum(j=i+1,#z,is_ndf_pair(z[i],z[j])));
A322437(n) = number_of_ndfpairs(Vec(factorizations(n))); \\ Antti Karttunen, Dec 10 2020

For n > 0, a(A002110(n)) = A322441(n)/2 = A339626(n). - Antti Karttunen, Dec 10 2020

Data section extended up to a(120) and more examples added by Antti Karttunen, Dec 10 2020

A328677 Numbers whose distinct prime indices are relatively prime and pairwise indivisible.

Original entry on oeis.org

2, 4, 8, 15, 16, 32, 33, 35, 45, 51, 55, 64, 69, 75, 77, 85, 93, 95, 99, 119, 123, 128, 135, 141, 143, 145, 153, 155, 161, 165, 175, 177, 187, 201, 205, 207, 209, 215, 217, 219, 221, 225, 245, 249, 253, 255, 256, 265, 275, 279, 287, 291, 295, 297, 309, 323

Offset: 1

Gus Wiseman, Oct 30 2019

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. Stable numbers are listed in A316476.

			The sequence of terms together with their prime indices begins:
    2: {1}
    4: {1,1}
    8: {1,1,1}
   15: {2,3}
   16: {1,1,1,1}
   32: {1,1,1,1,1}
   33: {2,5}
   35: {3,4}
   45: {2,2,3}
   51: {2,7}
   55: {3,5}
   64: {1,1,1,1,1,1}
   69: {2,9}
   75: {2,3,3}
   77: {4,5}
   85: {3,7}
   93: {2,11}
   95: {3,8}
   99: {2,2,5}
  119: {4,7}

These are the Heinz numbers of the partitions counted by A328676.
Numbers whose prime indices are relatively prime are A289509.
Partitions whose distinct parts are pairwise indivisible are A305148.
The version for binary indices (instead of prime indices) is A328671.
Cf. A000837, A056239, A112798, A285573, A289508, A303362, A304713, A327393, A328460.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
Select[Range[100],GCD@@primeMS[#]==1&&stableQ[primeMS[#],Divisible]&]

Intersection of A316476 and A289509.

A350840 Number of strict integer partitions of n with no adjacent parts of quotient 2.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 2, 4, 5, 6, 7, 8, 10, 13, 17, 19, 22, 25, 30, 35, 43, 52, 60, 70, 81, 93, 106, 122, 142, 166, 190, 216, 249, 287, 325, 371, 420, 479, 543, 617, 695, 784, 888, 1000, 1126, 1266, 1420, 1594, 1792, 2008, 2247, 2514, 2809, 3135, 3496, 3891, 4332

Offset: 0

Gus Wiseman, Jan 20 2022

nonn

			The a(1) = 1 through a(13) = 13 partitions (A..D = 10..13):
  1   2   3   4    5    6    7    8     9     A     B     C     D
              31   32   51   43   53    54    64    65    75    76
                   41        52   62    72    73    74    93    85
                             61   71    81    82    83    A2    94
                                  431   432   91    92    B1    A3
                                        531   532   A1    543   B2
                                              541   641   651   C1
                                                    731   732   643
                                                          741   652
                                                          831   751
                                                                832
                                                                931
                                                                5431

The version for subsets of prescribed maximum is A045691.
The double-free case is A120641.
The non-strict case is A350837, ranked by A350838.
An additive version (differences) is A350844, non-strict A350842.
The non-strict complement is counted by A350846, ranked by A350845.
Versions for prescribed quotients:
  = 2:  A154402, sets A001511.
  != 2: A350840 (this sequence), sets A045691.
  >= 2: A000929, sets A018819.
  <= 2: A342095, non-strict A342094.
  < 2:  A342097, non-strict A342096, sets A045690.
  > 2:  A342098, sets A040039.
A000041 = integer partitions.
A000045 = sets containing n with all differences > 2.
A003114 = strict partitions with no successions, ranked by A325160.
A116931 = partitions with no successions, ranked by A319630.
A116932 = partitions with differences != 1 or 2, strict A025157.
A323092 = double-free integer partitions, ranked by A320340.
A350839 = partitions with gaps and conjugate gaps, ranked by A350841.
Cf. A003000, A018819, A303362, A323093, A323094, A337135, A342191.

Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&And@@Table[#[[i-1]]/#[[i]]!=2,{i,2,Length[#]}]&]],{n,0,30}]

A317102 Powerful numbers whose distinct prime multiplicities are pairwise indivisible.

Original entry on oeis.org

1, 4, 8, 9, 16, 25, 27, 32, 36, 49, 64, 72, 81, 100, 108, 121, 125, 128, 169, 196, 200, 216, 225, 243, 256, 288, 289, 343, 361, 392, 432, 441, 484, 500, 512, 529, 625, 648, 675, 676, 729, 800, 841, 864, 900, 961, 968, 972, 1000, 1024, 1089, 1125, 1152, 1156, 1225

Offset: 1

Gus Wiseman, Aug 01 2018

nonn

A number is powerful if its prime multiplicities are all greater than 1.

			144 = 2^4 * 3^2 is not in the sequence because 4 and 2 are not pairwise indivisible.

Robert Israel, Table of n, a(n) for n = 1..3000

Cf. A001694, A056239, A112798, A118914, A124010, A285572, A285573, A303362, A304713, A316475, A317101, A317616.

filter:= proc(n) local L,i,j,q;
  L:= convert(map(t -> t[2],ifactors(n)[2]),set);
  if min(L) = 1 then return false fi;
  for j from 2 to nops(L) do
    for i from 1 to j-1 do
      q:= L[i]/L[j];
      if q::integer or (1/q)::integer then return false fi;
  od od;
  true
end proc:
select(filter, [$4..10000]); # Robert Israel, Jun 23 2019

Select[Range[1000],And[Max@@Last/@FactorInteger[#]>=2,Select[Tuples[Last/@FactorInteger[#],2],And[UnsameQ@@#,Divisible@@#]&]=={}]&]

Definition corrected and a(1)=1 inserted by Robert Israel, Jun 23 2019

A318727 Number of integer compositions of n where adjacent parts are indivisible (either way) and the last and first part are also indivisible (either way).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 5, 3, 5, 13, 9, 23, 15, 37, 45, 63, 115, 131, 207, 265, 415, 603, 823, 1251, 1673, 2521, 3519, 5147, 7409, 10449, 15225, 21497, 31285, 44719, 64171, 92315, 131619, 190085, 271871, 391189, 560979, 804265, 1155977, 1656429, 2381307, 3414847

Offset: 1

Gus Wiseman, Sep 02 2018

nonn

			The a(10) = 13 compositions:
  (10)
  (7,3) (3,7) (6,4) (4,6)
  (5,3,2) (5,2,3) (3,5,2) (3,2,5) (2,5,3) (2,3,5)
  (3,2,3,2) (2,3,2,3)

Andrew Howroyd, Table of n, a(n) for n = 1..200

Cf. A000740, A008965, A167606, A285573, A296302, A303362, A304713, A316476, A318726.

Table[Select[Join@@Permutations/@IntegerPartitions[n],!MatchQ[#,({_,x_,y_,_}/;Divisible[x,y]||Divisible[y,x])|({y_,_,x_}/;Divisible[x,y]||Divisible[y,x])]&]//Length,{n,20}]

b(n,k,pred)={my(M=matrix(n,n)); for(n=1, n, M[n,n]=pred(k,n); for(j=1, n-1, M[n,j]=sum(i=1, n-j, if(pred(i,j), M[n-j,i], 0)))); sum(i=1, n, if(pred(i,k), M[n,i], 0))}
a(n)={1 + sum(k=1, n-1, b(n-k, k, (i,j)->i%j<>0&&j%i<>0))} \\ Andrew Howroyd, Sep 08 2018

a(21)-a(28) from Robert Price, Sep 07 2018

Terms a(29) and beyond from Andrew Howroyd, Sep 08 2018

A323054 Number of strict integer partitions of n with no 1's such that no part is a power of any other part.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 1, 3, 3, 4, 4, 6, 6, 8, 9, 12, 13, 16, 19, 21, 25, 30, 36, 40, 47, 53, 63, 71, 83, 94, 107, 121, 140, 159, 180, 204, 233, 260, 296, 334, 377, 421, 474, 532, 598, 668, 750, 835, 933, 1038, 1163, 1292, 1435, 1597, 1771, 1966, 2180, 2421, 2673

Offset: 0

Gus Wiseman, Jan 04 2019

nonn

			The a(2) = 1 through a(13) = 8 strict integer partitions (A = 10, B = 11, C = 12, D = 13):
  (2)  (3)  (4)  (5)   (6)  (7)   (8)   (9)   (A)    (B)    (C)    (D)
                 (32)       (43)  (53)  (54)  (64)   (65)   (75)   (76)
                            (52)  (62)  (63)  (73)   (74)   (84)   (85)
                                        (72)  (532)  (83)   (A2)   (94)
                                                     (92)   (543)  (A3)
                                                     (632)  (732)  (B2)
                                                                   (643)
                                                                   (652)

Fausto A. C. Cariboni, Table of n, a(n) for n = 0..350

Cf. A001597, A007916, A025147, A052410, A078374, A087897, A120641, A275972, A303362, A323053, A323087, A323088.

stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
Table[Length[Select[IntegerPartitions[n],And[FreeQ[#,1],UnsameQ@@#,stableQ[#,IntegerQ[Log[#1,#2]]&]]&]],{n,30}]

A323094 Number of strict integer partitions of n where no part is 2^k times any other part, for any k > 0.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 2, 4, 4, 4, 5, 7, 8, 10, 12, 12, 15, 17, 20, 24, 27, 33, 35, 41, 48, 54, 61, 69, 79, 87, 101, 113, 128, 144, 159, 181, 201, 225, 251, 281, 311, 347, 388, 428, 477, 525, 579, 643, 712, 788, 868, 954, 1051, 1155, 1272, 1398, 1534, 1682, 1840, 2016

Offset: 0

Gus Wiseman, Jan 04 2019

nonn

			The a(1) = 1 through a(12) = 8 strict integer partitions (A = 10, B = 11, C = 12):
  (1)  (2)  (3)  (4)   (5)   (6)   (7)   (8)   (9)    (A)    (B)    (C)
                 (31)  (32)  (51)  (43)  (53)  (54)   (64)   (65)   (75)
                                   (52)  (62)  (72)   (73)   (74)   (93)
                                   (61)  (71)  (531)  (91)   (83)   (A2)
                                                      (532)  (92)   (B1)
                                                             (A1)   (543)
                                                             (731)  (651)
                                                                    (732)

Fausto A. C. Cariboni, Table of n, a(n) for n = 0..400

Cf. A018819, A087897, A120641, A276431, A303362, A323054, A323092, A323093.

stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&stableQ[#,IntegerQ[Log[2,#1/#2]]&]&]],{n,30}]

cp's OEIS Frontend

A316467 Matula-Goebel numbers of locally stable rooted identity trees, meaning no branch is a subset of any other branch of the same root.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A316474 Number of locally stable rooted identity trees with n nodes, meaning no branch is a subset of any other branch of the same root.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

A316494 Matula-Goebel numbers of locally disjoint rooted identity trees, meaning no branch overlaps any other branch of the same root.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A322437 Number of unordered pairs of factorizations of n into factors > 1 where no factor of one divides any factor of the other.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A328677 Numbers whose distinct prime indices are relatively prime and pairwise indivisible.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A350840 Number of strict integer partitions of n with no adjacent parts of quotient 2.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

A317102 Powerful numbers whose distinct prime multiplicities are pairwise indivisible.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Extensions

A318727 Number of integer compositions of n where adjacent parts are indivisible (either way) and the last and first part are also indivisible (either way).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica