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

A196724 Number of subsets of {1..n} (including empty set) such that the pairwise products of distinct elements are all distinct.

Original entry on oeis.org

1, 2, 4, 8, 16, 32, 58, 116, 212, 416, 720, 1440, 2340, 4680, 7920, 13024, 23328, 46656, 74168, 148336, 229856, 371424, 615304, 1230608, 1780224, 3401568, 5589360, 9468504, 14397744, 28795488, 40312128, 80624256, 131388480, 206363168, 335814288, 521401536
Offset: 0

Views

Author

Alois P. Heinz, Oct 06 2011

Keywords

Comments

The number of subsets of {1..n} such that every orderless pair of (not necessarily distinct) elements has a different product is A325860(n). - Gus Wiseman, Jun 03 2019

Examples

			a(6) = 58: from the 2^6=64 subsets of {1,2,3,4,5,6} only 6 do not have all the pairwise products of elements distinct: {1,2,3,6}, {2,3,4,6}, {1,2,3,4,6}, {1,2,3,5,6}, {2,3,4,5,6}, {1,2,3,4,5,6}.
		

Crossrefs

The subset case is A196724 (this sequence).
The maximal case is A325859.
The integer partition case is A325856.
The strict integer partition case is A325855.
Heinz numbers of the counterexamples are given by A325993.

Programs

  • Maple
    b:= proc(n, s) local sn, m;
          m:= nops(s);
          sn:= [s[], n];
          `if`(n<1, 1, b(n-1, s) +`if`(m*(m+1)/2 = nops(({seq(seq(
           sn[i]*sn[j], j=i+1..m+1), i=1..m)})), b(n-1, sn), 0))
        end:
    a:= proc(n) option remember;
          b(n-1, [n]) +`if`(n=0, 0, a(n-1))
        end:
    seq(a(n), n=0..20);
  • Mathematica
    b[n_, s_] := b[n, s] = Module[{sn, m}, m = Length[s]; sn = Append[s, n]; If[n < 1, 1, b[n - 1, s] + If[m*(m + 1)/2 == Length[Union[Flatten[Table[ sn[[i]] * sn[[j]], {i, 1, m}, {j, i + 1, m + 1}]]]], b[n - 1, sn], 0]]]; a[n_] := a[n] = b[n - 1, {n}] + If[n == 0, 0, a[n - 1]]; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Jan 31 2017, translated from Maple *)
    Table[Length[Select[Subsets[Range[n]],UnsameQ@@Times@@@Subsets[#,{2}]&]],{n,0,10}] (* Gus Wiseman, Jun 03 2019 *)

Extensions

Name edited by Gus Wiseman, Jun 03 2019
a(33)-a(35) from Fausto A. C. Cariboni, Oct 05 2020

A339564 Number of ways to choose a distinct factor in a factorization of n (pointed factorizations).

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 1, 4, 2, 3, 1, 7, 1, 3, 3, 7, 1, 7, 1, 7, 3, 3, 1, 14, 2, 3, 4, 7, 1, 10, 1, 12, 3, 3, 3, 17, 1, 3, 3, 14, 1, 10, 1, 7, 7, 3, 1, 26, 2, 7, 3, 7, 1, 14, 3, 14, 3, 3, 1, 25, 1, 3, 7, 19, 3, 10, 1, 7, 3, 10, 1, 36, 1, 3, 7, 7, 3, 10, 1, 26, 7, 3
Offset: 1

Views

Author

Gus Wiseman, Apr 10 2021

Keywords

Examples

			The pointed factorizations of n for n = 2, 4, 6, 8, 12, 24, 30:
  ((2))  ((4))    ((6))    ((8))      ((12))     ((24))       ((30))
         ((2)*2)  ((2)*3)  ((2)*4)    ((2)*6)    ((3)*8)      ((5)*6)
                  (2*(3))  (2*(4))    (2*(6))    (3*(8))      (5*(6))
                           ((2)*2*2)  ((3)*4)    ((4)*6)      ((2)*15)
                                      (3*(4))    (4*(6))      (2*(15))
                                      ((2)*2*3)  ((2)*12)     ((3)*10)
                                      (2*2*(3))  (2*(12))     (3*(10))
                                                 ((2)*2*6)    ((2)*3*5)
                                                 (2*2*(6))    (2*(3)*5)
                                                 ((2)*3*4)    (2*3*(5))
                                                 (2*(3)*4)
                                                 (2*3*(4))
                                                 ((2)*2*2*3)
                                                 (2*2*2*(3))
		

Crossrefs

The additive version is A000070 (strict: A015723).
The unpointed version is A001055 (strict: A045778, ordered: A074206, listed: A162247).
Allowing point (1) gives A057567.
Choosing a position instead of value gives A066637.
The ordered additive version is A336875.
A000005 counts divisors.
A001787 count normal multisets with a selected position.
A001792 counts compositions with a selected position.
A006128 counts partitions with a selected position.
A066186 count strongly normal multisets with a selected position.
A254577 counts ordered factorizations with a selected position.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Sum[Length[Union[fac]],{fac,facs[n]}],{n,50}]

Formula

a(n) = A057567(n) - A001055(n).
a(n) = Sum_{d|n, d>1} A001055(n/d).

A347460 Number of distinct possible alternating products of factorizations of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 4, 1, 4, 1, 4, 2, 2, 1, 6, 2, 2, 3, 4, 1, 5, 1, 5, 2, 2, 2, 7, 1, 2, 2, 6, 1, 5, 1, 4, 4, 2, 1, 8, 2, 4, 2, 4, 1, 5, 2, 6, 2, 2, 1, 10, 1, 2, 4, 6, 2, 5, 1, 4, 2, 5, 1, 10, 1, 2, 4, 4, 2, 5, 1, 8, 4, 2, 1, 10, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Oct 06 2021

Keywords

Comments

We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).
A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.

Examples

			The a(n) alternating products for n = 1, 4, 8, 12, 24, 30, 36, 48, 60, 120:
  1  4  8    12   24   30    36   48    60    120
     1  2    3    6    10/3  9    12    15    30
        1/2  3/4  8/3  5/6   4    16/3  20/3  40/3
             1/3  2/3  3/10  1    3     15/4  15/2
                  3/8  2/15  4/9  3/4   12/5  24/5
                  1/6        1/4  1/3   3/5   10/3
                             1/9  3/16  5/12  5/6
                                  1/12  4/15  8/15
                                        3/20  3/10
                                        1/15  5/24
                                              2/15
                                              3/40
                                              1/30
		

Crossrefs

Positions of 1's are 1 and A000040.
Positions of 2's appear to be A001358.
Positions of 3's appear to be A030078.
Dominates A038548, the version for reverse-alternating product.
Counting only integers gives A046951.
The even-length case is A072670.
The version for partitions (not factorizations) is A347461, reverse A347462.
The odd-length case is A347708.
The length-3 case is A347709.
A001055 counts factorizations (strict A045778, ordered A074206).
A056239 adds up prime indices, row sums of A112798.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A108917 counts knapsack partitions, ranked by A299702.
A276024 counts distinct positive subset-sums of partitions, strict A284640.
A292886 counts knapsack factorizations, by sum A293627.
A299701 counts distinct subset-sums of prime indices, positive A304793.
A301957 counts distinct subset-products of prime indices.
A304792 counts distinct subset-sums of partitions.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Table[Length[Union[altprod/@facs[n]]],{n,100}]

A347461 Number of distinct possible alternating products of integer partitions of n.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 7, 10, 12, 16, 19, 23, 27, 34, 41, 49, 57, 67, 78, 91, 106, 125, 147, 166, 187, 215, 245, 277, 317, 357, 405, 460, 524, 592, 666, 740, 829, 928, 1032, 1147, 1273, 1399, 1555, 1713, 1892, 2087, 2298, 2523, 2783, 3070, 3383, 3724, 4104, 4504
Offset: 0

Views

Author

Gus Wiseman, Oct 06 2021

Keywords

Comments

We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).

Examples

			Partitions representing each of the a(7) = 10 alternating products are:
     (7) -> 7
    (61) -> 6
    (52) -> 5/2
   (511) -> 5
    (43) -> 4/3
   (421) -> 2
  (4111) -> 4
   (331) -> 1
   (322) -> 3
  (3211) -> 3/2
		

Crossrefs

The version for alternating sum is A004526.
Counting only integers gives A028310, reverse A347707.
The version for factorizations is A347460, reverse A038548.
The reverse version is A347462.
A000041 counts partitions.
A027187 counts partitions of even length.
A027193 counts partitions of odd length.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A108917 counts knapsack partitions, ranked by A299702.
A122768 counts distinct submultisets of partitions.
A126796 counts complete partitions.
A293627 counts knapsack factorizations by sum.
A301957 counts distinct subset-products of prime indices.
A304792 counts subset-sums of partitions, positive A276024, strict A284640.
A304793 counts distinct positive subset-sums of prime indices.
A325534 counts separable partitions, ranked by A335433.
A325535 counts inseparable partitions, ranked by A335448.

Programs

  • Mathematica
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Table[Length[Union[altprod/@IntegerPartitions[n]]],{n,0,30}]

A347462 Number of distinct possible reverse-alternating products of integer partitions of n.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 8, 11, 13, 17, 22, 28, 33, 42, 51, 59, 69, 84, 100, 117, 137, 163, 191, 222, 256, 290, 332, 378, 429, 489, 564, 643, 729, 819, 929, 1040, 1167, 1313, 1473, 1647, 1845, 2045, 2272, 2521, 2785, 3076, 3398, 3744, 4115, 4548, 5010, 5524, 6086
Offset: 0

Views

Author

Gus Wiseman, Oct 06 2021

Keywords

Comments

We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)). The reverse-alternating product is the alternating product of the reversed sequence.

Examples

			Partitions representing each of the a(7) = 11 reverse-alternating products:
     (7) -> 7
    (61) -> 1/6
    (52) -> 2/5
   (511) -> 5
    (43) -> 3/4
   (421) -> 2
  (4111) -> 1/4
   (331) -> 1
   (322) -> 3
  (3211) -> 2/3
  (2221) -> 1/2
		

Crossrefs

The version for non-reverse alternating sum instead of product is A004526.
Counting only integers gives A028310, non-reverse A347707.
The version for factorizations is A038548, non-reverse A347460.
The non-reverse version is A347461.
A000041 counts partitions.
A027187 counts partitions of even length.
A027193 counts partitions of odd length.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A108917 counts knapsack partitions, ranked by A299702.
A122768 counts distinct submultisets of partitions.
A126796 counts complete partitions.
A293627 counts knapsack factorizations by sum.
A301957 counts distinct subset-products of prime indices.
A304792 counts subset-sums of partitions, positive A276024, strict A284640.
A304793 counts distinct positive subset-sums of prime indices.
A325534 counts separable partitions, ranked by A335433.
A325535 counts inseparable partitions, ranked by A335448.

Programs

  • Mathematica
    revaltprod[q_]:=Product[Reverse[q][[i]]^(-1)^(i-1),{i,Length[q]}];
    Table[Length[Union[revaltprod/@IntegerPartitions[n]]],{n,0,30}]

A347705 Number of factorizations of n with reverse-alternating product > 1.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 3, 1, 2, 1, 4, 1, 2, 2, 3, 1, 4, 1, 4, 2, 2, 1, 7, 1, 2, 3, 4, 1, 5, 1, 7, 2, 2, 2, 7, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 12, 1, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 8, 2, 5, 1, 4, 2, 5, 1, 16, 1, 2, 4, 4, 2, 5, 1, 12, 3, 2, 1, 11, 2
Offset: 1

Views

Author

Gus Wiseman, Oct 12 2021

Keywords

Comments

A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.
We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)). The reverse-alternating product is the alternating product of the reversed sequence.

Examples

			The a(n) factorizations for n = 2, 6, 8, 12, 24, 30, 48, 60:
  2   6     8       12      24        30      48          60
      2*3   2*4     2*6     3*8       5*6     6*8         2*30
            2*2*2   3*4     4*6       2*15    2*24        3*20
                    2*2*3   2*12      3*10    3*16        4*15
                            2*2*6     2*3*5   4*12        5*12
                            2*3*4             2*3*8       6*10
                            2*2*2*3           2*4*6       2*5*6
                                              3*4*4       3*4*5
                                              2*2*12      2*2*15
                                              2*2*2*6     2*3*10
                                              2*2*3*4     2*2*3*5
                                              2*2*2*2*3
		

Crossrefs

Positions of 1's are A000430.
The weak version (>= instead of >) is A001055, non-reverse A347456.
The non-reverse version is A339890, strict A347447.
The version for reverse-alternating product 1 is A347438.
Allowing any integer reciprocal alternating product gives A347439.
The even-length case is A347440, also the opposite reverse version.
Allowing any integer rev-alt product gives A347442, non-reverse A347437.
The version for partitions is A347449, non-reverse A347448.
A001055 counts factorizations (strict A045778, ordered A074206).
A038548 counts possible rev-alt products of factorizations, integer A046951.
A103919 counts partitions by sum and alternating sum, reverse A344612.
A292886 counts knapsack factorizations, by sum A293627.
A347707 counts possible integer reverse-alternating products of partitions.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    revaltprod[q_]:=Product[q[[-i]]^(-1)^(i-1),{i,Length[q]}];
    Table[Length[Select[facs[n],revaltprod[#]>1&]],{n,100}]

Formula

a(n) = A001055(n) - A347438(n).

A347708 Number of distinct possible alternating products of odd-length factorizations of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 11 2021

Keywords

Comments

We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).
A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.
Note that it is sufficient to look at only length-1 and length-3 factorizations; cf. A347709.

Examples

			Representative factorizations for each of the a(180) = 7 alternating products:
  (2*2*3*3*5) -> 5
     (2*2*45) -> 45
     (2*3*30) -> 20
     (2*5*18) -> 36/5
     (2*9*10) -> 20/9
     (3*4*15) -> 45/4
        (180) -> 180
		

Crossrefs

The version for partitions is A028310, reverse A347707.
Positions of 1's appear to be A037143 \ {1}.
The even-length version for n > 1 is A072670, strict A211159.
Counting only integers appears to give A293234, with evens A046951.
This is the odd-length case of A347460, reverse A038548.
The any-length version for partitions is A347461, reverse A347462.
The length-3 case is A347709.
A001055 counts factorizations (strict A045778, ordered A074206).
A056239 adds up prime indices, row sums of A112798.
A276024 counts distinct positive subset-sums of partitions.
A292886 counts knapsack factorizations, by sum A293627.
A301957 counts distinct subset-products of prime indices.
A304792 counts distinct subset-sums of partitions.
A347050 = factorizations w/ an alternating permutation, complement A347706.
A347441 counts odd-length factorizations with integer alternating product.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Table[Length[Union[altprod/@Select[facs[n],OddQ[Length[#]]&]]],{n,100}]
  • PARI
    altprod(facs) = prod(i=1,#facs,facs[i]^((-1)^(i-1)));
    A347708aux(n, m=n, facs=List([])) = if(1==n, if((#facs)%2, altprod(facs), 0), my(newfacs, r, rats=List([])); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); r = A347708aux(n/d, d, newfacs); if(r, rats = concat(rats,r)))); (rats));
    A347708(n) = if(1==n,0,#Set(A347708aux(n))); \\ Antti Karttunen, Jan 29 2025

Formula

Conjecture: For n > 1, a(n) = 1 + A347460(n) - A038548(n) + A072670(n).

Extensions

Data section extended to a(105) by Antti Karttunen, Jan 29 2025

A294150 Number of knapsack partitions of n that are also knapsack factorizations.

Original entry on oeis.org

1, 1, 1, 2, 2, 4, 4, 6, 8, 10, 12, 13, 20, 20, 29, 30, 41, 41, 56, 53, 81, 75
Offset: 1

Views

Author

Gus Wiseman, Oct 23 2017

Keywords

Comments

a(n) is the number of finite multisets of positive integers summing to n such that every distinct submultiset has a different sum, and also every distinct submultiset has a different product.

Examples

			The a(12) = 13 partitions are:
(12),
(10 2), (9 3), (8 4), (7 5), (6 6),
(8 2 2), (7 3 2), (5 5 2), (5 4 3), (4 4 4),
(3 3 3 3),
(2 2 2 2 2 2).
		

Crossrefs

Programs

  • Mathematica
    nn=22;
    dubQ[y_]:=And[UnsameQ@@Times@@@Union[Rest@Subsets[y]],UnsameQ@@Plus@@@Union[Rest@Subsets[y]]];
    Table[Length@Select[IntegerPartitions[n],dubQ],{n,nn}]

A347709 Number of distinct rational numbers of the form x * z / y for some factorization x * y * z = n, 1 < x <= y <= z.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 14 2021

Keywords

Comments

This is also the number of distinct possible alternating products of length-3 factorizations of n, where we define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)), and where a factorization of n is a weakly increasing sequence of positive integers > 1 with product n.

Examples

			Representative factorizations for each of the a(360) = 9 alternating products:
   (2,2,90) -> 90
   (2,3,60) -> 40
   (2,4,45) -> 45/2
   (2,5,36) -> 72/5
   (2,6,30) -> 10
   (2,9,20) -> 40/9
  (2,10,18) -> 18/5
  (2,12,15) -> 5/2
   (3,8,15) -> 45/8
		

Crossrefs

Allowing factorizations of any length <= 3 gives A033273.
Positions of positive terms are A033942.
Positions of 0's are A037143.
The length-2 version is A072670.
Allowing any length gives A347460, reverse A038548.
Allowing any odd length gives A347708.
A001055 counts factorizations (strict A045778, ordered A074206).
A122179 counts length-3 factorizations.
A292886 counts knapsack factorizations, by sum A293627.
A301957 counts distinct subset-products of prime indices.
A304792 counts distinct subset-sums of partitions, positive A276024.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Table[Length[Union[altprod/@Select[facs[n],Length[#]==3&]]],{n,100}]
  • PARI
    A347709(n) = { my(rats=List([])); fordiv(n,z,my(yx=n/z); fordiv(yx, y, my(x = yx/y); if((y <= z) && (x <= y) && (x > 1), listput(rats,x*z/y)))); #Set(rats); }; \\ Antti Karttunen, Jan 29 2025

Extensions

More terms from Antti Karttunen, Jan 29 2025

A316364 Number of factorizations of n into factors > 1 such that every distinct submultiset of the factors has a different average.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 3, 1, 2, 2, 2, 1, 3, 1, 3, 2, 2, 1, 4, 1, 2, 2, 3, 1, 5, 1, 3, 2, 2, 2, 5, 1, 2, 2, 5, 1, 5, 1, 3, 3, 2, 1, 6, 1, 3, 2, 3, 1, 5, 2, 5, 2, 2, 1, 8, 1, 2, 3, 4, 2, 5, 1, 3, 2, 5, 1, 9, 1, 2, 3, 3, 2, 5, 1, 6, 2, 2, 1, 9, 2, 2, 2, 5, 1, 9, 2, 3, 2, 2, 2, 10, 1, 3, 3, 5, 1, 5, 1, 5, 4
Offset: 1

Views

Author

Gus Wiseman, Jun 30 2018

Keywords

Comments

Note that such a factorization is necessarily strict.

Examples

			The a(80) = 6 factorizations are (80), (10*8), (16*5), (20*4), (40*2), (10*4*2).
		

Crossrefs

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],UnsameQ@@Mean/@Union[Subsets[#]]&]],{n,50}]
  • PARI
    choosebybits(v,m) = { my(s=vector(hammingweight(m)),i=j=1); while(m>0,if(m%2,s[j] = v[i];j++); i++; m >>= 1); s; };
    hasdupavgs(v) = { my(avgs=Map(), k); for(i=1,(2^(#v))-1,k = (vecsum(choosebybits(v,i))/hammingweight(i)); if(mapisdefined(avgs,k),return(i),mapput(avgs,k,i))); (0); };
    A316364(n, m=n, facs=List([])) = if(1==n, (0==hasdupavgs(Vec(facs))), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A316364(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Sep 21 2018

Extensions

More terms from Antti Karttunen, Sep 21 2018
Showing 1-10 of 11 results. Next