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 31 results. Next

A381455 Number of multisets that can be obtained by taking the sum of each block of a multiset partition of the prime indices of n into a multiset of constant multisets.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 1, 5, 1, 2, 1, 2, 1, 1, 1, 3, 2, 1, 3, 2, 1, 1, 1, 7, 1, 1, 1, 4, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 5, 2, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 2, 1, 1, 2, 11, 1, 1, 1, 2, 1, 1, 1, 6, 1, 1, 2, 2, 1, 1, 1, 5, 5, 1, 1, 2, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Mar 06 2025

Keywords

Comments

First differs from A000688 at a(144) = 9, A000688(144) = 10.
First differs from A295879 at a(128) = 15, A295879(128) = 13.
Also the number of multisets that can be obtained by taking the sums of prime indices of each factor in a factorization of n into prime powers > 1.
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.
A multiset partition can be regarded as an arrow in the ranked poset of integer partitions. For example, we have {{1},{1,2},{1,3},{1,2,3}}: {1,1,1,1,2,2,3,3} -> {1,3,4,6}, or (33221111) -> (6431) (depending on notation).
Multisets of constant multisets are generally not transitive. For example, we have arrows: {{1,1},{2}}: {1,1,2} -> {2,2} and {{2,2}}: {2,2} -> {4}, but there is no multiset of constant multisets {1,1,2} -> {4}.

Examples

			The prime indices of 36 are {1,1,2,2}, with the following 4 partitions into a multiset of constant multisets:
  {{1,1},{2,2}}
  {{1},{1},{2,2}}
  {{2},{2},{1,1}}
  {{1},{1},{2},{2}}
with block-sums: {2,4}, {1,1,4}, {2,2,2}, {1,1,2,2}, which are all different, so a(36) = 4.
The prime indices of 144 are {1,1,1,1,2,2}, with the following 10 partitions into a multiset of constant multisets:
  {{2,2},{1,1,1,1}}
  {{1},{2,2},{1,1,1}}
  {{2},{2},{1,1,1,1}}
  {{1,1},{1,1},{2,2}}
  {{1},{1},{1,1},{2,2}}
  {{1},{2},{2},{1,1,1}}
  {{2},{2},{1,1},{1,1}}
  {{1},{1},{1},{1},{2,2}}
  {{1},{1},{2},{2},{1,1}}
  {{1},{1},{1},{1},{2},{2}}
with block-sums: {4,4}, {1,3,4}, {2,2,4}, {2,2,4}, {1,1,2,4}, {1,2,2,3}, {2,2,2,2}, {1,1,1,1,4}, {1,1,2,2,2}, {1,1,1,1,2,2}, of which 9 are distinct, so a(144) = 9.
The a(n) partitions for n = 4, 8, 16, 32, 36, 64, 72, 128:
  (2)   (3)    (4)     (5)      (42)    (6)       (43)     (7)
  (11)  (21)   (22)    (32)     (222)   (33)      (322)    (43)
        (111)  (31)    (41)     (411)   (42)      (421)    (52)
               (211)   (221)    (2211)  (51)      (2221)   (61)
               (1111)  (311)            (222)     (4111)   (322)
                       (2111)           (321)     (22111)  (331)
                       (11111)          (411)              (421)
                                        (2211)             (511)
                                        (3111)             (2221)
                                        (21111)            (3211)
                                        (111111)           (4111)
                                                           (22111)
                                                           (31111)
                                                           (211111)
                                                           (1111111)
		

Crossrefs

Before taking sums we had A000688.
Positions of 1 are A005117.
There is a chain from the prime indices of n to a singleton iff n belongs to A300273.
The lower version is A381453.
For distinct blocks we have A381715, before sum A050361.
For distinct block-sums we have A381716, before sums A381635 (zeros A381636).
Other multiset partitions of prime indices:
- For multiset partitions (A001055) see A317141 (upper), A300383 (lower).
- For strict multiset partitions (A045778) see A381452.
- For set multipartitions (A050320) see A381078 (upper), A381454 (lower).
- For set systems (A050326) see A381441 (upper).
- For strict multiset partitions with distinct sums (A321469) see A381637.
- For set systems with distinct sums (A381633) see A381634, A293243.
More on multiset partitions into constant blocks: A006171, A279784, A295935.
A000041 counts integer partitions, strict A000009.
A000040 lists the primes.
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.

Programs

  • Mathematica
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
    sqfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[sqfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],PrimePowerQ]}]];
    Table[Length[Union[Sort[hwt/@#]&/@sqfacs[n]]],{n,100}]

Formula

a(s) = 1 for any squarefree number s.
a(p^k) = A000041(k) for any prime p.

A381871 Numbers whose prime indices cannot be partitioned into constant blocks having a common sum.

Original entry on oeis.org

6, 10, 14, 15, 18, 20, 21, 22, 24, 26, 28, 30, 33, 34, 35, 38, 39, 42, 44, 45, 46, 50, 51, 52, 54, 55, 56, 57, 58, 60, 62, 65, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 82, 84, 85, 86, 87, 88, 90, 91, 92, 93, 94, 95, 96, 98, 99, 100, 102, 104, 105, 106, 110
Offset: 1

Views

Author

Gus Wiseman, Mar 13 2025

Keywords

Comments

First differs from A383100 in lacking 108.
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.
Also numbers that cannot be written as a product of prime powers with equal sums of prime indices.
Partitions of this type are counted by A381993.

Examples

			The terms together with their prime indices begin:
    6: {1,2}
   10: {1,3}
   14: {1,4}
   15: {2,3}
   18: {1,2,2}
   20: {1,1,3}
   21: {2,4}
   22: {1,5}
   24: {1,1,1,2}
   26: {1,6}
   28: {1,1,4}
   30: {1,2,3}
		

Crossrefs

Constant blocks: A000688, A006171, A279784, A295935, A381453 (lower), A381455 (upper).
Constant blocks with distinct sums: A381635, A381716.
For distinct instead of equal sums we have A381636, counted by A381717.
Partitions of this type are counted by A381993, complement A383093.
These are the positions of 0 in A381995.
A001055 counts multiset partitions of prime indices, strict A045778.
A050361 counts multiset partitions into distinct constant blocks.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A317141 counts coarsenings of prime indices, refinements A300383.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]&/@sps[Range[Length[mset]]]];
    Select[Range[100],Select[mps[prix[#]],SameQ@@Total/@#&&And@@SameQ@@@#&]=={}&]

A302492 Products of any power of 2 with prime numbers of prime-power index, i.e., prime numbers p of the form p = prime(q^k), for q prime, k >= 1.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 27, 28, 30, 31, 32, 33, 34, 35, 36, 38, 40, 41, 42, 44, 45, 46, 48, 49, 50, 51, 53, 54, 55, 56, 57, 59, 60, 62, 63, 64, 66, 67, 68, 69, 70, 72, 75, 76, 77, 80, 81, 82, 83
Offset: 1

Views

Author

Gus Wiseman, Apr 08 2018

Keywords

Comments

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

Examples

			Entry A302242 describes a correspondence between positive integers and multiset multisystems. In this case it gives the following sequence of multiset multisystems.
01: {}
02: {{}}
03: {{1}}
04: {{},{}}
05: {{2}}
06: {{},{1}}
07: {{1,1}}
08: {{},{},{}}
09: {{1},{1}}
10: {{},{2}}
11: {{3}}
12: {{},{},{1}}
14: {{},{1,1}}
15: {{1},{2}}
16: {{},{},{},{}}
17: {{4}}
18: {{},{1},{1}}
19: {{1,1,1}}
20: {{},{},{2}}
21: {{1},{1,1}}
22: {{},{3}}
23: {{2,2}}
24: {{},{},{},{1}}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],Or[#===1,And@@PrimePowerQ/@PrimePi/@DeleteCases[FactorInteger[#][[All,1]],2]]&]
  • PARI
    ok(n)={!#select(p->p<>2&&!isprimepower(primepi(p)), factor(n)[,1])} \\ Andrew Howroyd, Aug 26 2018

A381715 Number of multisets that can be obtained by taking the sum of each block of a multiset partition of the prime indices of n into distinct constant blocks.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 10 2025

Keywords

Comments

First differs from A050361 at a(1728) = 7, A050361(1728) = 8.
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 1728 are {1,1,1,1,1,1,2,2,2}, with multiset partitions into distinct constant blocks:
  {{2,2,2},{1,1,1,1,1,1}}
  {{1},{2,2,2},{1,1,1,1,1}}
  {{2},{2,2},{1,1,1,1,1,1}}
  {{1,1},{2,2,2},{1,1,1,1}}
  {{1},{2},{2,2},{1,1,1,1,1}}
  {{1},{1,1},{1,1,1},{2,2,2}}
  {{2},{1,1},{2,2},{1,1,1,1}}
  {{1},{2},{1,1},{2,2},{1,1,1}}
with sums:
  {6,6}
  {1,5,6}
  {2,4,6}
  {2,4,6}
  {1,2,4,5}
  {1,2,3,6}
  {2,2,4,4}
  {1,2,2,3,4}
of which 7 are distinct, so a(1728) = 7.
		

Crossrefs

Without distinct blocks (A000688) we have A381455, lower (A355731) A381453.
More on multiset partitions into constant blocks: A006171, A279784, A295935.
Positions of terms > 1 are A046099.
Before taking sums we had A050361.
For equal instead of distinct blocks we have A362421.
For strict instead of constant blocks we have A381441, before sums A050326.
For just distinct blocks we have A381452, before sums A045778.
For distinct sums we have A381716, before sums A381635, zeros A381636.
A001055 counts multiset partitions, see A317141 (upper), A300383 (lower).
A003963 gives product of prime indices.
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}]]]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]&/@sps[Range[Length[mset]]]];
    Table[Length[Union[Sort[Total/@#]&/@Select[mps[prix[n]],UnsameQ@@#&&And@@SameQ@@@#&]]],{n,100}]

A356065 Squarefree numbers whose prime indices are all prime-powers.

Original entry on oeis.org

1, 3, 5, 7, 11, 15, 17, 19, 21, 23, 31, 33, 35, 41, 51, 53, 55, 57, 59, 67, 69, 77, 83, 85, 93, 95, 97, 103, 105, 109, 115, 119, 123, 127, 131, 133, 155, 157, 159, 161, 165, 177, 179, 187, 191, 201, 205, 209, 211, 217, 227, 231, 241, 249, 253, 255, 265, 277
Offset: 1

Views

Author

Gus Wiseman, Jul 25 2022

Keywords

Examples

			105 has prime indices {2,3,4}, all three of which are prime-powers, so 105 is in the sequence.
		

Crossrefs

The multiplicative version (factorizations) is A050361, non-strict A000688.
Heinz numbers of the partitions counted by A054685, with 1's A106244, non-strict A023894, non-strict with 1's A023893.
Counting twice-partitions of this type gives A279786, non-strict A279784.
Counting twice-factorizations gives A295935, non-strict A296131.
These are the odd products of distinct elements of A302493.
Allowing prime index 1 gives A302496, non-strict A302492.
The case of primes (instead of prime-powers) is A302590, non-strict A076610.
These are the squarefree positions of 1's in A355741.
This is the squarefree case of A355743, complement A356066.
A001222 counts prime-power divisors.
A005117 lists the squarefree numbers.
A034699 gives maximal prime-power divisor.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
A355742 chooses a prime-power divisor of each prime index.

Programs

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

Formula

Intersection of A005117 and A355743.

A381995 Number of ways to partition the prime indices of n into constant blocks with a common sum.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 19 2025

Keywords

Comments

Also the number of factorizations of n into prime powers > 1 with equal sums of prime indices.

Examples

			The prime indices of 144 are {1,1,1,1,2,2}, with the following 2 multiset partitions into constant blocks with a common sum:
  {{2,2},{1,1,1,1}}
  {{2},{2},{1,1},{1,1}}
so a(144) = 2.
		

Crossrefs

For just constant blocks we have A000688.
Twice-partitions of this type are counted by A279789.
For just a common sum we have A321455.
For distinct instead of equal sums we have A381635.
Positions of 0 are A381871, counted by A381993.
MM-numbers of these multiset partitions are A382215.
A001055 counts factorizations, strict A045778.
A050361 counts factorizations into distinct prime powers.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A317141 counts coarsenings of prime indices, refinements A300383.
A353864 counts rucksack partitions, ranked by A353866.
Cf. A279784, A295935, A381453 (lower), A381455 (upper).

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]& /@ sps[Range[Length[set]]]];
    Table[Length[Select[mps[prix[n]], SameQ@@Total/@#&&And@@SameQ@@@#&]],{n,100}]

Formula

A323774(n) = Sum_{A056239(k)=n} a(k). Gus Wiseman, Apr 25 2025

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

Views

Author

Gus Wiseman, Nov 30 2017

Keywords

Comments

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

Examples

			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).
		

Crossrefs

Programs

  • Mathematica
    Table[DivisorSum[GCD@@FactorInteger[n][[All,2]],PartitionsP],{n,100}]
  • PARI
    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

Formula

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

Extensions

More terms from Antti Karttunen, Jul 29 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

Views

Author

Gus Wiseman, Nov 29 2017

Keywords

Comments

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

Examples

			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
		

Crossrefs

Programs

  • Maple
    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
  • Mathematica
    Table[Sum[DivisorSigma[0,d],{d,Divisors[GCD@@FactorInteger[n][[All,2]]]}],{n,100}]

Formula

a(A175082(k)) = 1, a(A093771(k)) = 3.
a(n) = Sum_{d|A052409(n)} A000005(d).

A381991 Numbers whose prime indices have a unique multiset partition into constant multisets with distinct sums.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79
Offset: 1

Views

Author

Gus Wiseman, Mar 22 2025

Keywords

Comments

Also numbers with a unique factorization into prime powers with 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.

Examples

			The prime indices of 270 are {1,2,2,2,3}, and there are two multiset partitions into constant multisets with distinct sums: {{1},{2},{3},{2,2}} and {{1},{3},{2,2,2}}, so 270 is not in the sequence.
The prime indices of 300 are {1,1,2,3,3}, of which there are no multiset partitions into constant multisets with distinct sums, so 300 is not in the sequence.
The prime indices of 360 are {1,1,1,2,2,3}, of which there is only one multiset partition into constant multisets with distinct sums: {{1},{1,1},{3},{2,2}}, so 360 is in the sequence.
The terms together with their prime indices begin:
    1: {}
    2: {1}
    3: {2}
    4: {1,1}
    5: {3}
    6: {1,2}
    7: {4}
    9: {2,2}
   10: {1,3}
   11: {5}
   13: {6}
   14: {1,4}
   15: {2,3}
   17: {7}
   18: {1,2,2}
   19: {8}
   20: {1,1,3}
   21: {2,4}
   22: {1,5}
   23: {9}
   24: {1,1,1,2}
   25: {3,3}
		

Crossrefs

For distinct blocks instead of block-sums we have A004709, counted by A000726.
Twice-partitions of this type are counted by A279786.
MM-numbers of these multiset partitions are A326535 /\ A355743.
These are the positions of 1 in A381635.
For no choices we have A381636 (zeros of A381635), counted by A381717.
For strict instead of constant blocks we have A381870, counted by A382079.
Partitions of this type (unique into constant with distinct) are counted by A382301.
Normal multiset partitions of this type are counted by A382203.
A001055 counts multiset partitions, see A317141 (upper), A300383 (lower), A265947.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.

Programs

  • Mathematica
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
    pfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[pfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],PrimePowerQ]}]];
    Select[Range[100],Length[Select[pfacs[#],UnsameQ@@hwt/@#&]]==1&]

A295923 Number of twice-factorizations of n where the first factorization is constant, i.e., type (P,R,P).

Original entry on oeis.org

1, 1, 1, 3, 1, 2, 1, 4, 3, 2, 1, 4, 1, 2, 2, 10, 1, 4, 1, 4, 2, 2, 1, 7, 3, 2, 4, 4, 1, 5, 1, 8, 2, 2, 2, 13, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 12, 3, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 29, 2, 5, 1, 4, 2, 5, 1, 16, 1, 2, 4, 4, 2, 5, 1, 12, 10, 2, 1, 11
Offset: 1

Views

Author

Gus Wiseman, Nov 30 2017

Keywords

Comments

a(n) is also the number of ways to choose a perfect divisor d|n and then a sequence of log_d(n) factorizations of d.

Examples

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

Crossrefs

Programs

  • Mathematica
    postfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[postfacs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    a[n_]:=Sum[Length[postfacs[n^(1/g)]]^g,{g,Divisors[GCD@@FactorInteger[n][[All,2]]]}];
    Array[a,50]
Previous Showing 11-20 of 31 results. Next