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 51-60 of 98 results. Next

A317885 Number of series-reduced free pure achiral multifunctions with one atom and n positions.

Original entry on oeis.org

1, 0, 0, 1, 1, 1, 2, 3, 4, 7, 9, 14, 21, 32, 45, 69, 103, 153, 224, 338, 500, 746, 1107, 1645, 2447, 3652, 5413, 8052, 11993, 17834, 26500, 39447, 58655, 87240, 129772, 193001, 287034, 427014, 635048, 944501, 1404910, 2089633, 3107864, 4622670, 6875533
Offset: 1

Views

Author

Gus Wiseman, Aug 09 2018

Keywords

Comments

A series-reduced free pure achiral multifunction (SRAM) is either (case 1) the leaf symbol "o", or (case 2) a nonempty and non-unitary expression of the form h[g, ..., g] where h and g are SRAMs. The number of positions in a SRAM is the number of brackets [...] plus the number of o's.

Examples

			The a(10) = 7 SRAMs:
  o[o[o,o],o[o,o]]
  o[o,o][o,o][o,o]
  o[o,o][o,o,o,o,o]
  o[o,o,o][o,o,o,o]
  o[o,o,o,o][o,o,o]
  o[o,o,o,o,o][o,o]
  o[o,o,o,o,o,o,o,o]
		

Crossrefs

Programs

  • Mathematica
    a[n_]:=If[n==1,1,Sum[a[k]*Sum[a[d],{d,Most[Divisors[n-k-1]]}],{k,n-2}]];
    Array[a,12]
  • PARI
    seq(n)={my(p=O(x)); for(n=1, n, p = x + p*x*sum(k=2, n-2, subst(p + O(x^(n\k+1)), x, x^k)) + O(x*x^n)); Vec(p)} \\ Andrew Howroyd, Aug 19 2018
    
  • PARI
    seq(n)={my(v=vector(n)); v[1]=1; for(n=2, #v, v[n]=sum(i=1, n-2, v[i]*sumdiv(n-i-1, d, if(dAndrew Howroyd, Aug 19 2018

Formula

a(1) = 1; a(n > 1) = Sum_{0 < k < n - 1} a(k) * Sum_{d|(n - k - 1), d < n - k - 1} a(d).

Extensions

Terms a(17) and beyond from Andrew Howroyd, Aug 19 2018

A342492 Number of compositions of n with weakly increasing first quotients.

Original entry on oeis.org

1, 1, 2, 4, 7, 11, 17, 26, 37, 52, 73, 95, 125, 163, 208, 261, 330, 407, 498, 607, 734, 881, 1056, 1250, 1480, 1738, 2029, 2359, 2742, 3160, 3635, 4169, 4760, 5414, 6151, 6957, 7861, 8858, 9952, 11148, 12483, 13934, 15526, 17267, 19173, 21252, 23535, 25991
Offset: 0

Views

Author

Gus Wiseman, Mar 16 2021

Keywords

Comments

Also called log-concave-up compositions.
The first quotients of a sequence are defined as if the sequence were an increasing divisor chain, so for example the first quotients of (6,3,1) are (1/2,1/3).

Examples

			The composition (4,2,1,2,3) has first quotients (1/2,1/2,2,3/2) so is not counted under a(12), even though the first differences (-2,-1,1,1) are weakly increasing.
The a(1) = 1 through a(6) = 17 compositions:
  (1)  (2)    (3)      (4)        (5)          (6)
       (1,1)  (1,2)    (1,3)      (1,4)        (1,5)
              (2,1)    (2,2)      (2,3)        (2,4)
              (1,1,1)  (3,1)      (3,2)        (3,3)
                       (1,1,2)    (4,1)        (4,2)
                       (2,1,1)    (1,1,3)      (5,1)
                       (1,1,1,1)  (2,1,2)      (1,1,4)
                                  (3,1,1)      (2,1,3)
                                  (1,1,1,2)    (2,2,2)
                                  (2,1,1,1)    (3,1,2)
                                  (1,1,1,1,1)  (4,1,1)
                                               (1,1,1,3)
                                               (2,1,1,2)
                                               (3,1,1,1)
                                               (1,1,1,1,2)
                                               (2,1,1,1,1)
                                               (1,1,1,1,1,1)
		

Crossrefs

The weakly decreasing version is A069916.
The version for differences instead of quotients is A325546.
The strictly increasing version is A342493.
The unordered version is A342497, ranked by A342523.
The strict unordered version is A342516.
A000005 counts constant compositions.
A000009 counts strictly increasing (or strictly decreasing) compositions.
A000041 counts weakly increasing (or weakly decreasing) compositions.
A000929 counts partitions with all adjacent parts x >= 2y.
A001055 counts factorizations.
A002843 counts compositions with all adjacent parts x <= 2y.
A003238 counts chains of divisors summing to n - 1 (strict: A122651).
A074206 counts ordered factorizations.
A167865 counts strict chains of divisors > 1 summing to n.

Programs

  • Maple
    b:= proc(n, q, l) option remember; `if`(n=0, 1, add(
         `if`(q=0 or q>=l/j, b(n-j, l/j, j), 0), j=1..n))
        end:
    a:= n-> b(n, 0$2):
    seq(a(n), n=0..50);  # Alois P. Heinz, Mar 25 2021
  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],LessEqual@@Divide@@@Reverse/@Partition[#,2,1]&]],{n,0,15}]
    (* Second program: *)
    b[n_, q_, l_] := b[n, q, l] = If[n == 0, 1, Sum[
         If[q == 0 || q >= l/j, b[n - j, l/j, j], 0], {j, 1, n}]];
    a[n_] := b[n, 0, 0];
    a /@ Range[0, 50] (* Jean-François Alcover, May 19 2021, after Alois P. Heinz *)

Extensions

a(21)-a(47) from Alois P. Heinz, Mar 25 2021

A342493 Number of compositions of n with strictly increasing first quotients.

Original entry on oeis.org

1, 1, 2, 3, 6, 8, 11, 16, 22, 28, 39, 49, 61, 77, 93, 114, 140, 169, 198, 233, 276, 321, 381, 439, 509, 591, 678, 774, 883, 1007, 1147, 1300, 1465, 1641, 1845, 2068, 2317, 2590, 2881, 3193, 3549, 3928, 4341, 4793, 5282, 5813, 6401, 7027, 7699, 8432, 9221, 10076
Offset: 0

Views

Author

Gus Wiseman, Mar 16 2021

Keywords

Comments

The first quotients of a sequence are defined as if the sequence were an increasing divisor chain, so for example the first quotients of (6,3,1) are (1/2,1/3).

Examples

			The composition (3,1,1,2) has first quotients (1/3,1,2) so is counted under a(7).
The a(1) = 1 through a(7) = 16 compositions:
  (1)  (2)    (3)    (4)      (5)      (6)        (7)
       (1,1)  (1,2)  (1,3)    (1,4)    (1,5)      (1,6)
              (2,1)  (2,2)    (2,3)    (2,4)      (2,5)
                     (3,1)    (3,2)    (3,3)      (3,4)
                     (1,1,2)  (4,1)    (4,2)      (4,3)
                     (2,1,1)  (1,1,3)  (5,1)      (5,2)
                              (2,1,2)  (1,1,4)    (6,1)
                              (3,1,1)  (2,1,3)    (1,1,5)
                                       (3,1,2)    (2,1,4)
                                       (4,1,1)    (2,2,3)
                                       (2,1,1,2)  (3,1,3)
                                                  (3,2,2)
                                                  (4,1,2)
                                                  (5,1,1)
                                                  (2,1,1,3)
                                                  (3,1,1,2)
		

Crossrefs

The version for differences instead of quotients is A325547.
The weakly increasing version is A342492.
The strictly decreasing version is A342494.
The unordered version is A342498, ranked by A342524.
The strict unordered version is A342517.
A000005 counts constant compositions.
A000009 counts strictly increasing (or strictly decreasing) compositions.
A000041 counts weakly increasing (or weakly decreasing) compositions.
A001055 counts factorizations.
A003238 counts chains of divisors summing to n - 1 (strict: A122651).
A074206 counts ordered factorizations.
A167865 counts strict chains of divisors > 1 summing to n.
A274199 counts compositions with all adjacent parts x < 2y.

Programs

  • Maple
    b:= proc(n, q, l) option remember; `if`(n=0, 1, add(
         `if`(q=0 or q>l/j, b(n-j, l/j, j), 0), j=1..n))
        end:
    a:= n-> b(n, 0$2):
    seq(a(n), n=0..55);  # Alois P. Heinz, Mar 25 2021
  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],Less@@Divide@@@Reverse/@Partition[#,2,1]&]],{n,0,15}]
    (* Second program: *)
    b[n_, q_, l_] := b[n, q, l] = If[n == 0, 1, Sum[
         If[q == 0 || q > l/j, b[n - j, l/j, j], 0], {j, 1, n}]];
    a[n_] := b[n, 0, 0];
    a /@ Range[0, 55] (* Jean-François Alcover, May 19 2021, after Alois P. Heinz *)

Extensions

a(21)-a(51) from Alois P. Heinz, Mar 18 2021

A342520 Number of strict integer partitions of n with distinct first quotients.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 4, 4, 6, 8, 10, 12, 13, 16, 20, 25, 30, 37, 42, 50, 57, 65, 80, 93, 108, 127, 147, 170, 198, 225, 258, 297, 340, 385, 448, 499, 566, 647, 737, 832, 937, 1064, 1186, 1348, 1522, 1701, 1916, 2157, 2402, 2697, 3013, 3355, 3742, 4190, 4656, 5191
Offset: 0

Views

Author

Gus Wiseman, Mar 20 2021

Keywords

Comments

Also the number of reversed strict integer partitions of n with distinct first quotients.
The first quotients of a sequence are defined as if the sequence were an increasing divisor chain, so for example the first quotients of (6,3,1) are (1/2,1/3).

Examples

			The strict partition (12,10,5,2,1) has first quotients (5/6,1/2,2/5,1/2) so is not counted under a(30), even though the first differences (-2,-5,-3,-1) are distinct.
The a(1) = 1 through a(13) = 16 partitions (A..D = 10..13):
  1   2   3    4    5    6     7    8     9     A      B      C     D
          21   31   32   42    43   53    54    64     65     75    76
                    41   51    52   62    63    73     74     84    85
                         321   61   71    72    82     83     93    94
                                    431   81    91     92     A2    A3
                                    521   432   532    A1     B1    B2
                                          531   541    542    543   C1
                                          621   631    632    642   643
                                                721    641    651   652
                                                4321   731    732   742
                                                       821    741   751
                                                       5321   831   832
                                                              921   841
                                                                    A21
                                                                    5431
                                                                    7321
		

Crossrefs

The version for differences instead of quotients is A320347.
The non-strict version is A342514 (ranking: A342521).
The equal instead of distinct version is A342515.
The non-strict ordered version is A342529.
The version for strict divisor chains is A342530.
A000041 counts partitions (strict: A000009).
A001055 counts factorizations (strict: A045778, ordered: A074206).
A003238 counts chains of divisors summing to n - 1 (strict: A122651).
A167865 counts strict chains of divisors > 1 summing to n.
A342086 counts strict chains of divisors with strictly increasing quotients.
A342098 counts (strict) partitions with all adjacent parts x > 2y.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&UnsameQ@@Divide@@@Partition[#,2,1]&]],{n,0,30}]

A342523 Heinz numbers of integer partitions with weakly increasing first quotients.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 31, 32, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 51, 52, 53, 55, 56, 57, 58, 59, 61, 62, 63, 64, 65, 66, 67, 68, 69, 71, 73, 74, 76
Offset: 1

Views

Author

Gus Wiseman, Mar 23 2021

Keywords

Comments

Also called log-concave-up partitions.
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
The first quotients of a sequence are defined as if the sequence were an increasing divisor chain, so for example the first quotients of (6,3,1) are (1/2,1/3).

Examples

			The prime indices of 60 are {1,1,2,3}, with first quotients (1,2,3/2), so 60 is not in the sequence.
Most small numbers are in the sequence, but the sequence of non-terms together with their prime indices begins:
   18: {1,2,2}
   30: {1,2,3}
   36: {1,1,2,2}
   50: {1,3,3}
   54: {1,2,2,2}
   60: {1,1,2,3}
   70: {1,3,4}
   72: {1,1,1,2,2}
   75: {2,3,3}
   90: {1,2,2,3}
   98: {1,4,4}
  100: {1,1,3,3}
		

Crossrefs

The version counting strict divisor chains is A057567.
For multiplicities (prime signature) instead of quotients we have A304678.
For differences instead of quotients we have A325360 (count: A240026).
These partitions are counted by A342523 (strict: A342516, ordered: A342492).
The strictly increasing version is A342524.
The weakly decreasing version is A342526.
A000041 counts partitions (strict: A000009).
A000929 counts partitions with adjacent parts x >= 2y.
A001055 counts factorizations (strict: A045778, ordered: A074206).
A003238 counts chains of divisors summing to n - 1 (strict: A122651).
A167865 counts strict chains of divisors > 1 summing to n.
A318991/A318992 rank reversed partitions with/without integer quotients.
A342086 counts strict chains of divisors with strictly increasing quotients.

Programs

  • Mathematica
    primeptn[n_]:=If[n==1,{},Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Select[Range[100],LessEqual@@Divide@@@Reverse/@Partition[primeptn[#],2,1]&]

A342530 Number of strict chains of divisors ending with n and having distinct first quotients.

Original entry on oeis.org

1, 2, 2, 3, 2, 6, 2, 6, 3, 6, 2, 12, 2, 6, 6, 9, 2, 12, 2, 12, 6, 6, 2, 28, 3, 6, 6, 12, 2, 26, 2, 14, 6, 6, 6, 31, 2, 6, 6, 28, 2, 26, 2, 12, 12, 6, 2, 52, 3, 12, 6, 12, 2, 28, 6, 28, 6, 6, 2, 66, 2, 6, 12, 25, 6, 26, 2, 12, 6, 26, 2, 76, 2, 6, 12, 12, 6, 26
Offset: 1

Views

Author

Gus Wiseman, Mar 25 2021

Keywords

Comments

The first quotients of a sequence are defined as if the sequence were an increasing divisor chain, so for example the quotients of (6,3,1) are (1/2,1/3).

Examples

			The a(1) = 1 through a(12) = 12 chains (reversed):
  1  2    3    4    5    6      7    8      9    10      11    12
     2/1  3/1  4/1  5/1  6/1    7/1  8/1    9/1  10/1    11/1  12/1
               4/2       6/2         8/2    9/3  10/2          12/2
                         6/3         8/4         10/5          12/3
                         6/2/1       8/2/1       10/2/1        12/4
                         6/3/1       8/4/1       10/5/1        12/6
                                                               12/2/1
                                                               12/3/1
                                                               12/4/1
                                                               12/4/2
                                                               12/6/1
                                                               12/6/2
Not counted under a(12) are: 12/4/2/1, 12/6/2/1, 12/6/3, 12/6/3/1.
		

Crossrefs

The version for weakly increasing first quotients is A057567.
The version for equal first quotients is A169594.
The case of chains starting with 1 is A254578.
The version for strictly increasing first quotients is A342086.
A001055 counts factorizations (strict: A045778, ordered: A074206).
A067824 counts strict chains of divisors ending with n.
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts strict chains of divisors.
A334997 counts chains of divisors of n by length.
A342495/A342529 count compositions with equal/distinct quotients.
A342496/A342514 count partitions with equal/distinct quotients.
A342515/A342520 count strict partitions with equal/distinct quotients.
A342522/A342521 rank partitions with equal/distinct quotients.

Programs

  • Mathematica
    cmi[n_]:=Prepend[Prepend[#,n]&/@Join@@cmi/@Most[Divisors[n]],{n}];
    Table[Length[Select[cmi[n],UnsameQ@@Divide@@@Partition[#,2,1]&]],{n,100}]

Formula

a(n) = Sum_{d|n} A254578(d). - Ridouane Oudra, Jun 17 2025

A343339 Numbers with no prime index dividing all the other prime indices, but with a prime index divisible by all the other prime indices.

Original entry on oeis.org

195, 555, 585, 915, 957, 975, 1295, 1335, 1665, 1695, 1755, 2193, 2265, 2343, 2535, 2585, 2715, 2745, 2775, 2871, 2925, 3115, 3345, 3367, 3729, 3765, 3885, 4005, 4209, 4215, 4575, 4755, 4875, 4995, 5085, 5265, 5285, 5385, 5457, 5467, 5709, 5955, 6205, 6215
Offset: 1

Views

Author

Gus Wiseman, Apr 15 2021

Keywords

Comments

Alternative name: Numbers > 1 whose smallest prime index does not divide all the other prime indices, but whose greatest prime index is divisible by all the other 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.
Also Heinz numbers of partitions with greatest part divisible by all the others, but smallest part not dividing all the others (counted by A343344). The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

Examples

			The sequence of terms together with their prime indices begins:
     195: {2,3,6}        2585: {3,5,15}       4575: {2,3,3,18}
     555: {2,3,12}       2715: {2,3,42}       4755: {2,3,66}
     585: {2,2,3,6}      2745: {2,2,3,18}     4875: {2,3,3,3,6}
     915: {2,3,18}       2775: {2,3,3,12}     4995: {2,2,2,3,12}
     957: {2,5,10}       2871: {2,2,5,10}     5085: {2,2,3,30}
     975: {2,3,3,6}      2925: {2,2,3,3,6}    5265: {2,2,2,2,3,6}
    1295: {3,4,12}       3115: {3,4,24}       5285: {3,4,36}
    1335: {2,3,24}       3345: {2,3,48}       5385: {2,3,72}
    1665: {2,2,3,12}     3367: {4,6,12}       5457: {2,7,28}
    1695: {2,3,30}       3729: {2,5,30}       5467: {4,5,20}
    1755: {2,2,2,3,6}    3765: {2,3,54}       5709: {2,5,40}
    2193: {2,7,14}       3885: {2,3,4,12}     5955: {2,3,78}
    2265: {2,3,36}       4005: {2,2,3,24}     6205: {3,7,21}
    2343: {2,5,20}       4209: {2,9,18}       6215: {3,5,30}
    2535: {2,3,6,6}      4215: {2,3,60}       6475: {3,3,4,12}
		

Crossrefs

The first condition alone gives A342193.
The second condition alone gives the complement of A343337.
The partitions with these Heinz numbers are counted by A343344.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A056239 adds up prime indices, row sums of A112798.
A067824 counts strict chains of divisors starting with n.
A253249 counts strict chains of divisors.
A339564 counts factorizations with a selected factor.

Programs

  • Mathematica
    Select[Range[2,1000],With[{p=PrimePi/@First/@FactorInteger[#]},And@@IntegerQ/@(Max@@p/p)&&!And@@IntegerQ/@(p/Min@@p)]&]

Formula

Complement of A343337 in A342193.

A343653 Number of non-singleton pairwise coprime nonempty sets of divisors > 1 of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 2, 0, 1, 1, 0, 0, 2, 0, 2, 1, 1, 0, 3, 0, 1, 0, 2, 0, 7, 0, 0, 1, 1, 1, 4, 0, 1, 1, 3, 0, 7, 0, 2, 2, 1, 0, 4, 0, 2, 1, 2, 0, 3, 1, 3, 1, 1, 0, 13, 0, 1, 2, 0, 1, 7, 0, 2, 1, 7, 0, 6, 0, 1, 2, 2, 1, 7, 0, 4, 0, 1, 0, 13, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Apr 25 2021

Keywords

Comments

First differs from A066620 at a(210) = 36, A066620(210) = 35.

Examples

			The a(n) sets for n = 6, 12, 24, 30, 36, 60, 72, 96:
  {2,3}  {2,3}  {2,3}  {2,3}    {2,3}  {2,3}    {2,3}  {2,3}
         {3,4}  {3,4}  {2,5}    {2,9}  {2,5}    {2,9}  {3,4}
                {3,8}  {3,5}    {3,4}  {3,4}    {3,4}  {3,8}
                       {5,6}    {4,9}  {3,5}    {3,8}  {3,16}
                       {2,15}          {4,5}    {4,9}  {3,32}
                       {3,10}          {5,6}    {8,9}
                       {2,3,5}         {2,15}
                                       {3,10}
                                       {3,20}
                                       {4,15}
                                       {5,12}
                                       {2,3,5}
                                       {3,4,5}
		

Crossrefs

The case of pairs is A089233.
The version with 1's, empty sets, and singletons is A225520.
The version for subsets of {1..n} is A320426.
The version for strict partitions is A337485.
The version for compositions is A337697.
The version for prime indices is A337984.
The maximal case with 1's is A343652.
The version with empty sets is a(n) + 1.
The version with singletons is A343654(n) - 1.
The version with empty sets and singletons is A343654.
The version with 1's is A343655.
The maximal case is A343660.
A018892 counts pairwise coprime unordered pairs of divisors.
A048691 counts pairwise coprime ordered pairs of divisors.
A048785 counts pairwise coprime ordered triples of divisors.
A051026 counts pairwise indivisible subsets of {1..n}.
A100565 counts pairwise coprime unordered triples of divisors.
A305713 counts pairwise coprime non-singleton strict partitions.
A343659 counts maximal pairwise coprime subsets of {1..n}.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Rest[Divisors[n]]],CoprimeQ@@#&]],{n,100}]

A255242 Calculate the aliquot parts of a number n and take their sum. Then repeat the process calculating the aliquot parts of all the previous aliquot parts and add their sum to the previous one. Repeat the process until the sum to be added is zero. Sequence lists these sums.

Original entry on oeis.org

0, 1, 1, 4, 1, 8, 1, 12, 5, 10, 1, 30, 1, 12, 11, 32, 1, 36, 1, 38, 13, 16, 1, 92, 7, 18, 19, 46, 1, 74, 1, 80, 17, 22, 15, 140, 1, 24, 19, 116, 1, 90, 1, 62, 51, 28, 1, 256, 9, 62, 23, 70, 1, 136, 19, 140, 25, 34, 1, 286, 1, 36, 61, 192, 21, 122, 1, 86, 29, 114
Offset: 1

Views

Author

Paolo P. Lava, Feb 19 2015

Keywords

Comments

a(n) = 1 if n is prime.

Examples

			The aliquot parts of 8 are 1, 2, 4 and their sum is 7.
Now, let us calculate the aliquot parts of 1, 2 and 4:
1 => 0;  2 => 1;  4 => 1, 2.  Their sum is 0 + 1 + 1 + 2 = 4.
Let us calculate the aliquot parts of 1, 1, 2:
1 => 0;  1 = > 0; 2 => 1. Their sum is 1.
We have left 1: 1 => 0.
Finally, 7 + 4 + 1 = 12. Therefore a(8) = 12.
		

Crossrefs

Programs

  • Maple
    with(numtheory): P:=proc(q) local a,b,c,k,n,t,v;
    for n from 1 to q do b:=0; a:=sort([op(divisors(n))]); t:=nops(a)-1;
    while add(a[k],k=1..t)>0 do b:=b+add(a[k],k=1..t); v:=[];
    for k from 2 to t do c:=sort([op(divisors(a[k]))]); v:=[op(v),op(c[1..nops(c)-1])]; od;
    a:=v; t:=nops(a); od; print(b); od; end: P(10^3);
  • Mathematica
    f[s_] := Flatten[Most[Divisors[#]] & /@ s]; a[n_] := Total@Flatten[FixedPointList[ f, {n}]] - n; Array[a, 100] (* Amiram Eldar, Apr 06 2019 *)
  • PARI
    ali(n) = setminus(divisors(n), Set(n));
    a(n) = my(list = List(), v = [n]); while (#v, my(w = []); for (i=1, #v, my(s=ali(v[i])); for (j=1, #s, w = concat(w, s[j]); listput(list, s[j]));); v = w;); vecsum(Vec(list)); \\ Michel Marcus, Jul 15 2023

Formula

a(1) = 0.
a(2^k) = k*2^(k-1) = A001787(k), for k>=1.
a(n^k) = (n^k-2^k)/(n-2), for n odd prime and k>=1.
In particular:
a(3^k) = A001047(k-1);
a(5^k) = A016127(k-1);
a(7^k) = A016130(k-1);
a(11^k) = A016135(k-1).
From Antti Karttunen, Nov 22 2024: (Start)
a(n) = A330575(n) - n.
Also, following formulas were conjectured by Sequence Machine:
a(n) = (A191161(n)-n)/2.
a(n) = Sum_{d|n} A001065(d)*A074206(n/d). [Compare to David A. Corneth's Apr 13 2020 formula for A330575]
a(n) = Sum_{d|n} A051953(d)*A067824(n/d).
a(n) = Sum_{d|n} A000203(d)*A174726(n/d).
a(n) = Sum_{d|n} A062790(d)*A253249(n/d).
a(n) = Sum_{d|n} A157658(d)*A191161(n/d).
a(n) = Sum_{d|n} A174725(d)*A211779(n/d).
a(n) = Sum_{d|n} A245211(d)*A323910(n/d).
(End)

A317099 Number of series-reduced planted achiral trees whose leaves span an initial interval of positive integers appearing with multiplicities an integer partition of n.

Original entry on oeis.org

1, 3, 4, 9, 8, 19, 16, 35, 35, 54, 57, 113, 102, 155, 189, 279, 298, 447, 491, 702, 813, 1063, 1256, 1759, 1967, 2542, 3050, 3902, 4566, 5882, 6843, 8676, 10205, 12612, 14908, 18608, 21638, 26510, 31292, 38150, 44584, 54185, 63262, 76308, 89371, 106818, 124755
Offset: 1

Views

Author

Gus Wiseman, Aug 01 2018

Keywords

Comments

In these trees, achiral means that all branches directly under any given node that is not a leaf or a cover of leaves are equal, and series-reduced means that every node that is not a leaf or a cover of leaves has at least two branches.

Examples

			The a(4) = 9 trees:
  (1111), ((11)(11)), (((1)(1))((1)(1))), ((1)(1)(1)(1)),
  (1112),
  (1122), ((12)(12)),
  (1123),
  (1234).
The a(6) = 19 trees:
  (111111), ((111)(111)), (((1)(1)(1))((1)(1)(1))), ((11)(11)(11)), (((1)(1))((1)(1))((1)(1))), ((1)(1)(1)(1)(1)(1)),
  (111112),
  (111122), ((112)(112)),
  (111123),
  (111222), ((12)(12)(12)),
  (111223),
  (111234),
  (112233), ((123)(123)),
  (112234),
  (112345),
  (123456).
		

Crossrefs

Programs

  • Mathematica
    b[n_]:=1+Sum[b[n/d],{d,Rest[Divisors[n]]}];
    a[n_]:=Sum[b[GCD@@Length/@Split[ptn]],{ptn,IntegerPartitions[n]}];
    Array[a,30]
Previous Showing 51-60 of 98 results. Next