A281113 -id:A281113 - cp's OEIS frontend

A294788 Number of twice-factorizations of type (Q,P,Q) and product n.

1, 1, 1, 1, 1, 3, 1, 3, 1, 3, 1, 5, 1, 3, 3, 3, 1, 5, 1, 5, 3, 3, 1, 12, 1, 3, 3, 5, 1, 12, 1, 5, 3, 3, 3, 13, 1, 3, 3, 12, 1, 12, 1, 5, 5, 3, 1, 19, 1, 5, 3, 5, 1, 12, 3, 12, 3, 3, 1, 26, 1, 3, 5, 11, 3, 12, 1, 5, 3, 12, 1, 26, 1, 3, 5, 5, 3, 12, 1, 19, 3, 3

Offset: 1

Gus Wiseman, Nov 08 2017

nonn

a(n) is the number of ways to choose a product-preserving permutation of a set partition of a factorization of n into distinct factors greater than one.

			The a(36) = 13 twice-factorizations are: (2)*(3)*(6), (2)*(3*6), (6)*(2*3), (2*3)*(6), (2*6)*(3), (2*3*6), (2)*(18), (2*18), (3)*(12), (3*12), (4)*(9), (4*9), (36).

Cf. A001055, A063834, A279790, A281113, A294786, A294787.

nn=100;
sfs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sfs[n/d],Min@@#>d&]],{d,Rest[Divisors[n]]}]];
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
Table[Total[Sum[Times@@Factorial/@Length/@Split[Sort[Times@@@f]],{f,sps[Sort[#]]}]&/@sfs[n]],{n,nn}]

A295279 Number of strict tree-factorizations of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 2, 3, 1, 4, 1, 4, 2, 2, 1, 10, 1, 2, 2, 4, 1, 8, 1, 6, 2, 2, 2, 12, 1, 2, 2, 10, 1, 8, 1, 4, 4, 2, 1, 26, 1, 4, 2, 4, 1, 10, 2, 10, 2, 2, 1, 28, 1, 2, 4, 12, 2, 8, 1, 4, 2, 8, 1, 44, 1, 2, 4, 4, 2, 8, 1, 26, 3, 2, 1

Offset: 1

Gus Wiseman, Nov 19 2017

nonn

A strict tree-factorization of n is either (case 1) the number n itself or (case 2) a set of two or more strict tree-factorizations, one of each factor in a factorization of n into distinct factors greater than one.

a(n) depends only on the prime signature of n. - Andrew Howroyd, Nov 18 2018

			The a(30) = 8 strict tree-factorizations are: 30, (2*3*5), (2*15), (2*(3*5)), (3*10), (3*(2*5)), (5*6), (5*(2*3)).
The a(36) = 12 strict tree-factorizations are: 36, (2*3*6), (2*3*(2*3)), (2*18), (2*(2*9)), (2*(3*6)), (2*(3*(2*3))), (3*12), (3*(2*6)), (3*(2*(2*3))), (3*(3*4)), (4*9).

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

Cf. A005804, A045778, A273873, A281113 A281118, A292504, A294786, A295281.

sfs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sfs[n/d],Min@@#>d&]],{d,Rest[Divisors[n]]}]];
sft[n_]:=1+Total[Function[fac,Times@@sft/@fac]/@Select[sfs[n],Length[#]>1&]];
Array[sft,100]

seq(n)={my(v=vector(n), w=vector(n)); w[1]=v[1]=1; for(k=2, n, w[k]=v[k]+1; forstep(j=n\k*k, k, -k, v[j]+=w[k]*v[j/k])); w} \\ Andrew Howroyd, Nov 18 2018

a(product of n distinct primes) = A005804(n).
a(prime^n) = A273873(n).
Dirichlet g.f.: (Zeta(s) + Product_{n >= 2}(1 + a(n)/n^s))/2.

A294786 Number of ways to choose a set partition of a factorization of n into distinct factors greater than one such that different blocks have different products.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 3, 1, 3, 1, 5, 1, 3, 3, 3, 1, 5, 1, 5, 3, 3, 1, 12, 1, 3, 3, 5, 1, 12, 1, 5, 3, 3, 3, 11, 1, 3, 3, 12, 1, 12, 1, 5, 5, 3, 1, 19, 1, 5, 3, 5, 1, 12, 3, 12, 3, 3, 1, 26, 1, 3, 5, 9, 3, 12, 1, 5, 3, 12, 1, 26, 1, 3, 5, 5, 3, 12, 1, 19, 3, 3

Offset: 1

Gus Wiseman, Nov 08 2017

nonn

			The a(36)=11 ways are:
(2)*(3)*(6),
(2)*(3*6), (2*6)*(3), (2)*(18), (3)*(12), (4)*(9),
(2*3*6), (2*18), (3*12), (4*9), (36).

Michael De Vlieger, Table of n, a(n) for n = 1..16384

Cf. A000258, A001055, A045778, A279375, A281113, A294787, A294788.

sfs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sfs[n/d],Min@@#>d&]],{d,Rest[Divisors[n]]}]];
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
Table[Length[Join@@Function[fac,Select[sps[fac],UnsameQ@@Times@@@#&]]/@sfs[n]],{n,100}]

a(product of n distinct primes) = A000258(n).

a(prime^n) = A279375(n).

A295924 Number of twice-factorizations of n of type (R,P,R).

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 1, 4, 3, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 4, 1, 1, 1, 1, 8, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 17, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1

Offset: 1

Gus Wiseman, Nov 30 2017

nonn

a(n) is the number of ways to choose an integer partition of a divisor of A052409(n).

			The a(16) = 8 twice-factorizations are (2)*(2)*(2)*(2), (2)*(2)*(2*2), (2)*(2*2*2), (2*2)*(2*2), (2*2*2*2), (4)*(4), (4*4), (16).

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

Cf. A000005, A000041, A001055, A047968, A052409, A052410, A089723, A281113, A284639, A295923, A295931, A295935, A296134.

Table[DivisorSum[GCD@@FactorInteger[n][[All,2]],PartitionsP],{n,100}]

A052409(n) = { my(k=ispower(n)); if(k, k, n>1); }; \\ From A052409
A295924(n) = if(1==n,n,sumdiv(A052409(n),d,numbpart(d))); \\ Antti Karttunen, Jul 29 2018

a(1) = 1; for n > 1, a(n) = Sum_{d|A052409(n)} A000041(d). - Antti Karttunen, Jul 29 2018

More terms from Antti Karttunen, Jul 29 2018

A296118 Number of ways to choose a factorization of each factor in a strict factorization of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 5, 2, 3, 1, 8, 1, 3, 3, 8, 1, 8, 1, 8, 3, 3, 1, 20, 2, 3, 5, 8, 1, 12, 1, 18, 3, 3, 3, 23, 1, 3, 3, 20, 1, 12, 1, 8, 8, 3, 1, 45, 2, 8, 3, 8, 1, 20, 3, 20, 3, 3, 1, 38, 1, 3, 8, 34, 3, 12, 1, 8, 3, 12, 1, 66, 1, 3, 8, 8, 3, 12, 1, 45, 8, 3

Offset: 1

Gus Wiseman, Dec 05 2017

nonn

			The a(12) = 8 twice-factorizations are (2)*(2*3), (2)*(6), (3)*(2*2), (3)*(4), (2*2*3), (2*6), (3*4), (12).

Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000

Cf. A000009, A005117, A045778, A261049, A271619, A281113, A294788, A295923, A296119, A296121.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Sum[Times@@(Length[facs[#]]&/@f),{f,Select[facs[n],UnsameQ@@#&]}],{n,100}]

A001055(n, m=n) = if(1==n, 1, sumdiv(n, d, if((d>1)&&(d<=m), A001055(n/d, d))));
A296118(n, m=n) = ((n<=m)*A001055(n) + sumdiv(n, d, if((d>1)&&(d<=m)&&(dA001055(d)*A296118(n/d, d-1)))); \\ Antti Karttunen, Oct 08 2018

Dirichlet g.f.: Product_{n > 1}(1 + A001055(n)/n^s).

A305150 Number of factorizations of n into distinct, pairwise indivisible factors greater than 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 1, 2, 1, 5, 1, 1, 2, 2, 2, 2, 1, 2, 2, 3, 1, 5, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 1, 3, 2, 3, 2, 2, 1, 6, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 3, 1, 2, 2, 2, 2, 5, 1, 3, 1, 2, 1, 6, 2, 2, 2, 3, 1, 6, 2, 2, 2, 2, 2, 4, 1, 2, 2, 2, 1, 5, 1, 3, 5

Offset: 1

Gus Wiseman, May 26 2018

nonn

			The a(60) = 6 factorizations are (3 * 4 * 5), (3 * 20), (4 * 15), (5 * 12), (6 * 10), (60). Missing from this list are (2 * 3 * 10), (2 * 5 * 6), (2 * 30).

Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..100000

Cf. A001055, A001970, A007716, A045778, A162247, A259936, A275024, A285572, A281113, A281116, A303386, A303431, A305001, A305148, A305149, A305253.

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], UnsameQ@@ # && Select[Tuples[Union[#], 2], UnsameQ@@ # && Divisible@@ # &] == {} &]], {n, 100}]

A305150(n, m=n, facs=List([])) = if(1==n, 1, my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m)&&factorback(apply(x -> (x%d),Vec(facs))), newfacs = List(facs); listput(newfacs,d); s += A305150(n/d, d-1, newfacs))); (s)); \\ Antti Karttunen, Dec 06 2018

a(n) <= A045778(n) <= A001055(n). - Antti Karttunen, Dec 06 2018

More terms from Antti Karttunen, Dec 06 2018

A295931 Number of ways to write n in the form n = (x^y)^z where x, y, and z are positive integers.

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 1, 3, 3, 1, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 3, 1, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 9, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 1

Offset: 1

Gus Wiseman, Nov 29 2017

nonn

By convention a(1) = 1.

Values can be 1, 3, 6, 9, 10, 15, 18, 21, 27, 28, 30, 36, 45, 54, 60, 63, 84, 90, etc. - Robert G. Wilson v, Dec 10 2017

			The a(256) = 10 ways are:
(2^1)^8    (2^2)^4   (2^4)^2  (2^8)^1
(4^1)^4    (4^2)^2   (4^4)^1
(16^1)^2   (16^2)^1
(256^1)^1

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

Cf. A000005, A007425, A052409, A052410, A089723, A093771, A175082, A277562, A281113, A284639, A294786, A294336, A294338, A295920, A295923, A295924, A295935.

f:= proc(n) local m,d,t;
  m:= igcd(seq(t[2],t=ifactors(n)[2]));
  add(numtheory:-tau(d),d=numtheory:-divisors(m))
end proc:
f(1):= 1:
map(f, [$1..100]); # Robert Israel, Dec 19 2017

Table[Sum[DivisorSigma[0,d],{d,Divisors[GCD@@FactorInteger[n][[All,2]]]}],{n,100}]

a(A175082(k)) = 1, a(A093771(k)) = 3.

a(n) = Sum_{d|A052409(n)} A000005(d).

A323718 Array read by antidiagonals upwards where A(n,k) is the number of k-times partitions of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 3, 1, 1, 1, 5, 6, 4, 1, 1, 1, 7, 15, 10, 5, 1, 1, 1, 11, 28, 34, 15, 6, 1, 1, 1, 15, 66, 80, 65, 21, 7, 1, 1, 1, 22, 122, 254, 185, 111, 28, 8, 1, 1, 1, 30, 266, 604, 739, 371, 175, 36, 9, 1, 1, 1, 42, 503, 1785, 2163, 1785, 672, 260, 45, 10, 1, 1

Offset: 0

Gus Wiseman, Jan 25 2019

A k-times partition of n for k > 1 is a sequence of (k-1)-times partitions, one of each part in an integer partition of n. A 1-times partition of n is just an integer partition of n, and the only 0-times partition of n is the number n itself.

			Array begins:
       k=0:   k=1:   k=2:   k=3:   k=4:   k=5:
  n=0:  1      1      1      1      1      1
  n=1:  1      1      1      1      1      1
  n=2:  1      2      3      4      5      6
  n=3:  1      3      6     10     15     21
  n=4:  1      5     15     34     65    111
  n=5:  1      7     28     80    185    371
  n=6:  1     11     66    254    739   1785
  n=7:  1     15    122    604   2163   6223
  n=8:  1     22    266   1785   8120  28413
  n=9:  1     30    503   4370  24446 101534
The A(4,2) = 15 twice-partitions:
  (4)  (31)    (22)    (211)      (1111)
       (3)(1)  (2)(2)  (11)(2)    (11)(11)
                       (2)(11)    (111)(1)
                       (21)(1)    (11)(1)(1)
                       (2)(1)(1)  (1)(1)(1)(1)

Alois P. Heinz, Rows n = 0..140, flattened

Columns: A000012 (k=0), A000041 (k=1), A063834 (k=2), A301595 (k=3).
Rows: A000027 (n=2), A000217 (n=3), A006003 (n=4).
Main diagonal gives A306187.
Cf. A001970, A055884, A096751, A144150, A196545, A281113, A289501, A290353, A300383, A323719, A327618, A327639.

b:= proc(n, i, k) option remember; `if`(n=0 or k=0 or i=1,
      1, b(n, i-1, k)+b(i$2, k-1)*b(n-i, min(n-i, i), k))
    end:
A:= (n, k)-> b(n$2, k):
seq(seq(A(d-k, k), k=0..d), d=0..14);  # Alois P. Heinz, Jan 25 2019

ptnlev[n_,k_]:=Switch[k,0,{n},1,IntegerPartitions[n],_,Join@@Table[Tuples[ptnlev[#,k-1]&/@ptn],{ptn,IntegerPartitions[n]}]];
Table[Length[ptnlev[sum-k,k]],{sum,0,12},{k,0,sum}]
(* Second program: *)
b[n_, i_, k_] := b[n, i, k] = If[n == 0 || k == 0 || i == 1, 1,
     b[n, i - 1, k] + b[i, i, k - 1]*b[n - i, Min[n - i, i], k]];
A[n_, k_] := b[n, n, k];
Table[Table[A[d - k, k], {k, 0, d}], {d, 0, 14}] // Flatten (* Jean-François Alcover, May 13 2021, after Alois P. Heinz *)

Column k is the formal power product transform of column k-1, where the formal power product transform of a sequence q with offset 1 is the sequence whose ordinary generating function is Product_{n >= 1} 1/(1 - q(n) * x^n).

A(n,k) = Sum_{i=0..k} binomial(k,i) * A327639(n,i). - Alois P. Heinz, Sep 20 2019

A302493 Prime numbers of prime-power index.

Original entry on oeis.org

2, 3, 5, 7, 11, 17, 19, 23, 31, 41, 53, 59, 67, 83, 97, 103, 109, 127, 131, 157, 179, 191, 211, 227, 241, 277, 283, 311, 331, 353, 367, 401, 419, 431, 461, 509, 547, 563, 587, 599, 617, 661, 691, 709, 719, 739, 773, 797, 859, 877, 919, 967, 991, 1009, 1031

Offset: 1

Gus Wiseman, Apr 08 2018

nonn

A prime index of n is a number m such that prime(m) divides n.

Cf. A000040, A000961, A001222, A005117, A007097, A056239, A063834, A275024, A281113, A302242, A302491, A302492.

Prime/@Select[Range[100],Or[#===1,PrimePowerQ[#]]&]

forprime(p=1, 500, if(p==2 || isprimepower(primepi(p)), print1(p, ", "))) \\ Felix Fröhlich, Apr 10 2018

a(n) = A000040(A000961(n)).

A318948 Number of ways to choose an integer partition of each factor in a factorization of n.

Original entry on oeis.org

1, 2, 3, 9, 7, 17, 15, 40, 39, 56, 56, 126, 101, 165, 197, 336, 297, 496, 490, 774, 837, 1114, 1255, 1948, 2007, 2638, 3127, 4123, 4565, 6201, 6842, 9131, 10311, 12904, 14988, 19516, 21637, 26995, 31488, 39250, 44583, 55418, 63261, 77683, 89935, 108068, 124754

Offset: 1

Gus Wiseman, Sep 05 2018

nonn

			The a(4) = 9 ways: (1+1)*(1+1), (1+1+1+1), (1+1)*(2), (2)*(1+1), (2+1+1), (2)*(2), (2+2), (3+1), (4).

Cf. A000041, A001055, A001970, A063834, A065026, A066739, A066815, A121229, A281113, A284639, A318949.

facs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@Select[facs[n/d],Min@@#1>=d&],{d,Rest[Divisors[n]]}]];
Table[Sum[Times@@PartitionsP/@fac,{fac,facs[n]}],{n,10}]

Dirichlet g.f.: Product_{n > 1} 1 / (1 - P(n) / n^s) where P = A000041. [clarified by Ilya Gutkovskiy, Oct 26 2019]

cp's OEIS Frontend

A294788 Number of twice-factorizations of type (Q,P,Q) and product n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A295279 Number of strict tree-factorizations of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A294786 Number of ways to choose a set partition of a factorization of n into distinct factors greater than one such that different blocks have different products.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A295924 Number of twice-factorizations of n of type (R,P,R).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A296118 Number of ways to choose a factorization of each factor in a strict factorization of n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A305150 Number of factorizations of n into distinct, pairwise indivisible factors greater than 1.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A295931 Number of ways to write n in the form n = (x^y)^z where x, y, and z are positive integers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula