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

A340102 Number of factorizations of 2n + 1 into an odd number of odd factors > 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 30 2020

Keywords

Examples

			The factorizations for 2n + 1 = 135, 225, 315, 405, 675, 1155, 1215:
  135      225      315      405         675         1155      1215
  3*5*9    5*5*9    5*7*9    5*9*9       3*3*75      3*5*77    3*5*81
  3*3*15   3*3*25   3*3*35   3*3*45      3*5*45      3*7*55    3*9*45
           3*5*15   3*5*21   3*5*27      3*9*25      5*7*33    5*9*27
                    3*7*15   3*9*15      5*5*27      3*11*35   9*9*15
                             3*3*3*3*5   5*9*15      5*11*21   3*15*27
                                         3*15*15     7*11*15   3*3*135
                                         3*3*3*5*5             3*3*3*5*9
                                                               3*3*3*3*15
		

Crossrefs

The version for partitions is A160786, ranked by A300272.
The not necessarily odd-length version is A340101.
A000009 counts partitions into odd parts, ranked by A066208.
A001055 counts factorizations, with strict case A045778.
A027193 counts partitions of odd length, ranked by A026424.
A058695 counts partitions of odd numbers, ranked by A300063.
A316439 counts factorizations by product and length.

Programs

  • Maple
    g:= proc(n, k, t) option remember; `if`(n>k, 0, t)+
          `if`(isprime(n), 0, add(`if`(d>k, 0, g(n/d, d, 1-t)),
              d=numtheory[divisors](n) minus {1, n}))
        end:
    a:= n-> `if`(n=0, 0, g(2*n+1$2, 1)):
    seq(a(n), n=0..100);  # Alois P. Heinz, Dec 30 2020
  • 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],OddQ[Length[#]]&&OddQ[Times@@#]&]],{n,1,100,2}];

A273013 Number of different arrangements of nonnegative integers on a pair of n-sided dice such that the dice can add to any integer from 0 to n^2-1.

Original entry on oeis.org

1, 1, 1, 3, 1, 7, 1, 10, 3, 7, 1, 42, 1, 7, 7, 35, 1, 42, 1, 42, 7, 7, 1, 230, 3, 7, 10, 42, 1, 115, 1, 126, 7, 7, 7, 393, 1, 7, 7, 230, 1, 115, 1, 42, 42, 7, 1, 1190, 3, 42, 7, 42, 1, 230, 7, 230, 7, 7, 1, 1158, 1, 7, 42, 462, 7, 115, 1, 42, 7, 115, 1, 3030
Offset: 1

Views

Author

Elliott Line, May 13 2016

Keywords

Comments

The set of b values (see formula), and therefore also a(n), depends only on the prime signature of n. So, for example, a(24) will be identical to a(n) of any other n which is also of the form p_1^3*p_2, (e.g., 40, 54, 56).
The value of b_1 will always be 1. When n is prime, the only nonzero b will be b_1, so therefore a(n) will be 1.
In any arrangement, both dice will have a 0, and one will have a 1 (here called the lead die). To determine any one of the actual arrangements to numbers on the dice, select one of the permutations of divisors (for the lead die), then select another permutation that is either the same length as that of the lead die, or one less. For example, if n = 12, we might select 2*3*2 for the lead die, and 3*4 for the second die. These numbers effectively tell you when to "switch track" when numbering the dice, and will uniquely result in the numbering: (0,1,6,7,12,13,72,73,78,79,84,85; 0,2,4,18,20,22,36,38,40,54,56,58).
a(n) is the number of (unordered) pairs of polynomials c(x) = x^c_1 + x^c_2 + ... + x^c_n, d(x) = x^d_1 + x^d_2 + ... + x^d_n with nonnegative integer exponents such that c(x)*d(x) = (x^(n^2)-1)/(x-1). - Alois P. Heinz, May 13 2016
a(n) is also the number of principal reversible squares of order n. - S. Harry White, May 19 2018
From Gus Wiseman, Oct 29 2021: (Start)
Also the number of ordered factorizations of n^2 with alternating product 1. This follows from the author's formula. Taking n instead of n^2 gives a(sqrt(n)) if n is a perfect square, otherwise 0. Here, an ordered factorization of n is a sequence of positive integers > 1 with product n, and the alternating product of a sequence (y_1,...,y_k) is Product_i y_i^((-1)^(i-1)). For example, the a(1) = 1 through a(9) = 3 factorizations are:
() (22) (33) (44) (55) (66) (77) (88) (99)
(242) (263) (284) (393)
(2222) (362) (482) (3333)
(2233) (2244)
(2332) (2442)
(3223) (4224)
(3322) (4422)
(22242)
(24222)
(222222)
The even-length case is A347464.
(End)

Examples

			When n = 4, a(n) = 3; the three arrangements are (0,1,2,3; 0,4,8,12), (0,1,4,5; 0,2,8,10), (0,1,8,9; 0,2,4,6).
When n = 5, a(n) = 1; the sole arrangement is (0,1,2,3,4; 0,5,10,15,20).
		

Crossrefs

Positions of 1's are 1 and A000040.
A000290 lists squares, complement A000037.
A001055 counts factorizations, ordered A074206.
A119620 counts partitions with alternating product 1, ranked by A028982.
A339846 counts even-length factorizations, ordered A174725.
A339890 counts odd-length factorizations, ordered A174726.
A347438 counts factorizations with alternating product 1.
A347460 counts possible alternating products of factorizations.
A347463 counts ordered factorizations with integer alternating product.
A347466 counts factorizations of n^2.

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[Select[Join@@Permutations/@facs[n^2], altprod[#] == 1&]],{n, 30}]
    (* Gus Wiseman, Oct 29 2021 *)
    (* or *)
    ofc[n_,k_] := If[k > PrimeOmega[n], 0, If[k == 0 && n == 1, 1, Sum[ofc[n/d, k-1],{d, Rest[Divisors[n]]}]]];
    Table[If[n == 1, 1, Sum[ofc[n, x]^2 + ofc[n, x]*ofc[n, x+1], {x, n}]],{n, 30}]
    (* Gus Wiseman, Oct 29 2021, based on author's formula *)
  • PARI
    A273013aux(n, k=0, t=1) = if(1==n, (1==t), my(s=0); fordiv(n, d, if((d>1), s += A273013aux(n/d, 1-k, t*(d^((-1)^k))))); (s));
    A273013(n) = A273013aux(n^2); \\ Antti Karttunen, Oct 30 2021
    
  • SageMath
    @cached_function
    def r(m,n):
        if n==1:
            return(1)
        divList = divisors(m)[:-1]
        return(sum(r(n,d) for d in divList))
    def A273013(n):
        return(r(n,n)) # William P. Orrick, Feb 19 2023

Formula

a(n) = b_1^2 + b_2^2 + b_3^2 + ... + b_1*b_2 + b_2*b_3 + b_3*b_4 + ..., where b_k is the number of different permutations of k divisors of n to achieve a product of n. For example, for n=24, b_3 = 9 (6 permutations of 2*3*4 and 3 permutations of 2*2*6).
a(n) = r(n,n) where r(m,1) = 1 and r(m,n) = Sum_{d|m,dWilliam P. Orrick, Feb 19 2023

A340854 Numbers that cannot be factored into factors > 1, the least of which is odd.

Original entry on oeis.org

1, 2, 4, 6, 8, 10, 14, 16, 20, 22, 26, 28, 32, 34, 38, 44, 46, 52, 58, 62, 64, 68, 74, 76, 82, 86, 88, 92, 94, 104, 106, 116, 118, 122, 124, 128, 134, 136, 142, 146, 148, 152, 158, 164, 166, 172, 178, 184, 188, 194, 202, 206, 212, 214, 218, 226, 232, 236, 244
Offset: 1

Views

Author

Gus Wiseman, Feb 04 2021

Keywords

Comments

Consists of 1 and all numbers that are even and have no odd divisor 1 < d <= n/d.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}              44: {1,1,5}          106: {1,16}
      2: {1}             46: {1,9}            116: {1,1,10}
      4: {1,1}           52: {1,1,6}          118: {1,17}
      6: {1,2}           58: {1,10}           122: {1,18}
      8: {1,1,1}         62: {1,11}           124: {1,1,11}
     10: {1,3}           64: {1,1,1,1,1,1}    128: {1,1,1,1,1,1,1}
     14: {1,4}           68: {1,1,7}          134: {1,19}
     16: {1,1,1,1}       74: {1,12}           136: {1,1,1,7}
     20: {1,1,3}         76: {1,1,8}          142: {1,20}
     22: {1,5}           82: {1,13}           146: {1,21}
     26: {1,6}           86: {1,14}           148: {1,1,12}
     28: {1,1,4}         88: {1,1,1,5}        152: {1,1,1,8}
     32: {1,1,1,1,1}     92: {1,1,9}          158: {1,22}
     34: {1,7}           94: {1,15}           164: {1,1,13}
     38: {1,8}          104: {1,1,1,6}        166: {1,23}
For example, the factorizations of 88 are (2*2*2*11), (2*2*22), (2*4*11), (2*44), (4*22), (8*11), (88), none of which has odd minimum, so 88 is in the sequence.
		

Crossrefs

The version looking at greatest factor is A000079.
The version for twice-balanced is A340656, with complement A340657.
These factorization are counted by A340832.
The complement is A340855.
A033676 selects the maximum inferior divisor.
A038548 counts inferior divisors.
A055396 selects the least prime index.
- Factorizations -
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
A339890 counts factorizations of odd length.
A340653 counts balanced factorizations.
- Odd -
A000009 counts partitions into odd parts.
A024429 counts set partitions of odd length.
A026424 lists numbers with odd Omega.
A066208 lists Heinz numbers of partitions into odd parts.
A067659 counts strict partitions of odd length (A030059).
A174726 counts ordered factorizations of odd length.

Programs

  • Mathematica
    Select[Range[100],Function[n,n==1||EvenQ[n]&&Select[Rest[Divisors[n]],OddQ[#]&&#<=n/#&]=={}]]

A340855 Numbers that can be factored into factors > 1, the least of which is odd.

Original entry on oeis.org

3, 5, 7, 9, 11, 12, 13, 15, 17, 18, 19, 21, 23, 24, 25, 27, 29, 30, 31, 33, 35, 36, 37, 39, 40, 41, 42, 43, 45, 47, 48, 49, 50, 51, 53, 54, 55, 56, 57, 59, 60, 61, 63, 65, 66, 67, 69, 70, 71, 72, 73, 75, 77, 78, 79, 80, 81, 83, 84, 85, 87, 89, 90, 91, 93, 95
Offset: 1

Views

Author

Gus Wiseman, Feb 04 2021

Keywords

Comments

These are numbers that are odd or have an odd divisor 1 < d <= n/d.

Examples

			The sequence of terms together with their prime indices begins:
     3: {2}          27: {2,2,2}      48: {1,1,1,1,2}
     5: {3}          29: {10}         49: {4,4}
     7: {4}          30: {1,2,3}      50: {1,3,3}
     9: {2,2}        31: {11}         51: {2,7}
    11: {5}          33: {2,5}        53: {16}
    12: {1,1,2}      35: {3,4}        54: {1,2,2,2}
    13: {6}          36: {1,1,2,2}    55: {3,5}
    15: {2,3}        37: {12}         56: {1,1,1,4}
    17: {7}          39: {2,6}        57: {2,8}
    18: {1,2,2}      40: {1,1,1,3}    59: {17}
    19: {8}          41: {13}         60: {1,1,2,3}
    21: {2,4}        42: {1,2,4}      61: {18}
    23: {9}          43: {14}         63: {2,2,4}
    24: {1,1,1,2}    45: {2,2,3}      65: {3,6}
    25: {3,3}        47: {15}         66: {1,2,5}
For example, 72 is in the sequence because it has three suitable factorizations: (3*3*8), (3*4*6), (3*24).
		

Crossrefs

The version looking at greatest factor is A057716.
The version for twice-balanced is A340657, with complement A340656.
These factorization are counted by A340832.
The complement is A340854.
A033676 selects the maximum inferior divisor.
A038548 counts inferior divisors, listed by A161906.
A055396 selects the least prime index.
- Factorizations -
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
A339890 counts factorizations of odd length.
A340653 counts balanced factorizations.
- Odd -
A000009 counts partitions into odd parts.
A024429 counts set partitions of odd length.
A026424 lists numbers with odd Omega.
A066208 lists Heinz numbers of partitions into odd parts.
A067659 counts strict partitions of odd length (A030059).
A174726 counts ordered factorizations of odd length.
A332304 counts strict compositions of odd length.
A340692 counts partitions of odd rank.

Programs

  • Mathematica
    Select[Range[100],Function[n,n>1&&(OddQ[n]||Select[Rest[Divisors[n]],OddQ[#]&&#<=n/#&]!={})]]

A340385 Number of integer partitions of n into an odd number of parts, the greatest of which is odd.

Original entry on oeis.org

1, 0, 2, 0, 3, 1, 6, 3, 10, 7, 18, 15, 30, 28, 51, 50, 82, 87, 134, 145, 211, 235, 331, 375, 510, 586, 779, 901, 1172, 1366, 1750, 2045, 2581, 3026, 3778, 4433, 5476, 6430, 7878, 9246, 11240, 13189, 15931, 18670, 22417, 26242, 31349, 36646, 43567, 50854
Offset: 1

Views

Author

Gus Wiseman, Jan 08 2021

Keywords

Examples

			The a(3) = 2 through a(10) = 7 partitions:
  3     5       321   7         332     9           532
  111   311           322       521     333         541
        11111         331       32111   522         721
                      511               531         32221
                      31111             711         33211
                      1111111           32211       52111
                                        33111       3211111
                                        51111
                                        3111111
                                        111111111
		

Crossrefs

Partitions of odd length are counted by A027193, ranked by A026424.
Partitions with odd maximum are counted by A027193, ranked by A244991.
The Heinz numbers of these partitions are given by A340386.
Other cases of odd length:
- A024429 counts set partitions of odd length.
- A067659 counts strict partitions of odd length.
- A089677 counts ordered set partitions of odd length.
- A166444 counts compositions of odd length.
- A174726 counts ordered factorizations of odd length.
- A332304 counts strict compositions of odd length.
- A339890 counts factorizations of odd length.
A000009 counts partitions into odd parts, ranked by A066208.
A026804 counts partitions whose least part is odd.
A058695 counts partitions of odd numbers, ranked by A300063.
A072233 counts partitions by sum and length.
A101707 counts partitions with odd rank.
A160786 counts odd-length partitions of odd numbers, ranked by A300272.
A340101 counts factorizations into odd factors.
A340102 counts odd-length factorizations into odd factors.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],OddQ[Length[#]*Max[#]]&]],{n,30}]

A340832 Number of factorizations of n into factors > 1 with odd least factor.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Feb 04 2021

Keywords

Examples

			The a(n) factorizations for n = 45, 108, 135, 180, 252:
  (45)     (3*36)     (135)      (3*60)     (3*84)
  (5*9)    (9*12)     (3*45)     (5*36)     (7*36)
  (3*15)   (3*4*9)    (5*27)     (9*20)     (9*28)
  (3*3*5)  (3*6*6)    (9*15)     (5*6*6)    (3*3*28)
           (3*3*12)   (3*5*9)    (3*3*20)   (3*4*21)
           (3*3*3*4)  (3*3*15)   (3*4*15)   (3*6*14)
                      (3*3*3*5)  (3*5*12)   (3*7*12)
                                 (3*6*10)   (3*3*4*7)
                                 (3*3*4*5)
		

Crossrefs

Positions of 0's are A340854.
Positions of nonzero terms are A340855.
The version for partitions is A026804.
Odd-length factorizations are counted by A339890.
The version looking at greatest factor is A340831.
- Factorizations -
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
A340101 counts factorizations into odd factors, odd-length case A340102.
A340607 counts factorizations with odd length and greatest factor.
A340653 counts balanced factorizations.
- Odd -
A000009 counts partitions into odd parts.
A026424 lists numbers with odd Omega.
A027193 counts partitions of odd length.
A058695 counts partitions of odd numbers (A300063).
A066208 lists numbers with odd-indexed prime factors.
A067659 counts strict partitions of odd length (A030059).
A174726 counts ordered factorizations of odd length.
A244991 lists numbers whose greatest prime index is odd.
A340692 counts partitions of odd rank.

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],OddQ@*Min]],{n,100}]
  • PARI
    A340832(n, m=n, fc=1) = if(1==n, (m%2)&&!fc, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A340832(n/d, d, 0*fc))); (s)); \\ Antti Karttunen, Dec 13 2021

Extensions

Data section extended up to 108 terms by Antti Karttunen, Dec 13 2021

A089677 Exponential convolution of A000670(n), with A000670(0)=0, with the sequence of all ones alternating in sign.

Original entry on oeis.org

0, 1, 1, 7, 37, 271, 2341, 23647, 272917, 3543631, 51123781, 811316287, 14045783797, 263429174191, 5320671485221, 115141595488927, 2657827340990677, 65185383514567951, 1692767331628422661, 46400793659664205567, 1338843898122192101557
Offset: 0

Views

Author

Mario Catalani (mario.catalani(AT)unito.it), Jan 03 2004

Keywords

Comments

Stirling transform of A005212(n)=[1,0,6,0,120,0,5040,...] is a(n)=[1,1,7,37,271,...]. - Michael Somos, Mar 04 2004
Occurs also as first column of a matrix-inversion occurring in a sum-of-like-powers problem. Consider the problem for any fixed natural number m>2 of finding solutions to sum(k=1,n,k^m) = (k+1)^m. Erdos conjectured that there are no solutions for n,m>2. Let D be the matrix of differences of D[m,n] := sum(k=1,n,k^m) - (k+1)^m. Then the generating functions for the rows of this matrix D constitute a set of polynomials in n (for varying n along columns) and the m-th polynomial defining the m-th row. Let GF_D be the matrix of the coefficients of this set of polynomials. Then the present sequence is the (unsigned) second column of GF_D^-1. - Gottfried Helms, Apr 01 2007

Examples

			From _Gus Wiseman_, Jan 06 2021: (Start)
a(n) is the number of ordered set partitions of {1..n} into an odd number of blocks. The a(1) = 1 through a(3) = 7 ordered set partitions are:
  {{1}}  {{1,2}}  {{1,2,3}}
                  {{1},{2},{3}}
                  {{1},{3},{2}}
                  {{2},{1},{3}}
                  {{2},{3},{1}}
                  {{3},{1},{2}}
                  {{3},{2},{1}}
(End)
		

Crossrefs

Ordered set partitions are counted by A000670.
The case of (unordered) set partitions is A024429.
The complement (even-length ordered set partitions) is counted by A052841.
A058695 counts partitions of odd numbers, ranked by A300063.
A101707 counts partitions of odd positive rank.
A160786 counts odd-length partitions of odd numbers, ranked by A300272.
A340102 counts odd-length factorizations into odd factors.
A340692 counts partitions of odd rank.
Other cases of odd length:
- A027193 counts partitions of odd length.
- A067659 counts strict partitions of odd length.
- A166444 counts compositions of odd length.
- A174726 counts ordered factorizations of odd length.
- A332304 counts strict compositions of odd length.
- A339890 counts factorizations of odd length.

Programs

  • Maple
    h := n -> add(combinat:-eulerian1(n,k)*2^k,k=0..n):
    a := n -> (h(n)-(-1)^n)/2: seq(a(n),n=0..20); # Peter Luschny, Jul 09 2015
  • Mathematica
    Table[Sum[Binomial[n, k](-1)^(n-k)Sum[i! StirlingS2[k, i], {i, 1, k}], {k, 0, n}], {n, 0, 20}]
  • PARI
    a(n)=if(n<0,0,n!*polcoeff(subst(y/(1-y^2),y,exp(x+x*O(x^n))-1),n))
    
  • PARI
    {a(n)=polcoeff(sum(m=0,n,(2*m+1)!*x^(2*m+1)/prod(k=1,2*m+1,1-k*x+x*O(x^n))),n)} /* Paul D. Hanna, Jul 20 2011 */
    
  • Sage
    def A089677_list(len):  # with a(0)=1
        e, r = [1], [1]
        for i in (1..len-1):
            for k in range(i-1, -1, -1): e[k] = (e[k]*i)//(i-k)
            r.append(-sum(e[j]*(-1)^(i-j) for j in (0..i-1)))
            e.append(sum(e))
        return r
    A089677_list(21) # Peter Luschny, Jul 09 2015

Formula

E.g.f.: (exp(x)-1)/(exp(x)*(2-exp(x))).
O.g.f.: Sum_{n>=0} (2*n+1)! * x^(2*n+1) / Product_{k=1..2*n+1} (1-k*x). - Paul D. Hanna, Jul 20 2011
a(n)=Sum(Binomial(n, k)(-1)^(n-k)Sum(i! Stirling2(k, i), i=1, ..k), k=0, .., n).
a(n) = (A000670(n)-(-1)^n)/2. - Vladeta Jovovic, Jan 17 2005
a(n) ~ n! / (4*(log(2))^(n+1)). - Vaclav Kotesovec, Feb 25 2014
a(n) = Sum_{k=0..floor(n/2)} (2*k+1)!*Stirling2(n, 2*k+1). - Peter Luschny, Sep 20 2015

A340831 Number of factorizations of n into factors > 1 with odd greatest factor.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Feb 04 2021

Keywords

Examples

			The a(n) factorizations for n = 45, 108, 135, 180, 252:
  (45)      (4*27)        (135)       (4*45)        (4*63)
  (5*9)     (2*6*9)       (3*45)      (12*15)       (12*21)
  (3*15)    (3*4*9)       (5*27)      (4*5*9)       (4*7*9)
  (3*3*5)   (2*2*27)      (9*15)      (2*2*45)      (6*6*7)
            (2*2*3*9)     (3*5*9)     (2*6*15)      (2*2*63)
            (2*2*3*3*3)   (3*3*15)    (3*4*15)      (2*6*21)
                          (3*3*3*5)   (2*2*5*9)     (3*4*21)
                                      (3*3*4*5)     (2*2*7*9)
                                      (2*2*3*15)    (2*3*6*7)
                                      (2*2*3*3*5)   (3*3*4*7)
                                                    (2*2*3*21)
                                                    (2*2*3*3*7)
		

Crossrefs

Positions of 0's are A000079.
The version for partitions is A027193.
The version for prime indices is A244991.
The version looking at length instead of greatest factor is A339890.
The version that also has odd length is A340607.
The version looking at least factor is A340832.
- Factorizations -
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
A340101 counts factorizations into odd factors, odd-length case A340102.
A340653 counts balanced factorizations.
- Odd -
A000009 counts partitions into odd parts.
A024429 counts set partitions of odd length.
A026424 lists numbers with odd Omega.
A058695 counts partitions of odd numbers.
A066208 lists numbers with odd-indexed prime factors.
A067659 counts strict partitions of odd length (A030059).
A174726 counts ordered factorizations of odd length.
A340692 counts partitions of odd rank.

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],OddQ@*Max]],{n,100}]
  • PARI
    A340831(n, m=n, fc=1) = if(1==n, !fc, my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&(!fc||(d%2)), s += A340831(n/d, d, 0*fc))); (s)); \\ Antti Karttunen, Dec 13 2021

Extensions

Data section extended up to 108 terms by Antti Karttunen, Dec 13 2021

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)

A347048 Number of even-length ordered factorizations of n with integer alternating product.

Original entry on oeis.org

1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 3, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 4, 0, 0, 0, 7, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 6, 1, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 11, 0, 0, 0, 1, 0, 0, 0, 11, 0, 0, 1, 1, 0, 0, 0, 6, 3, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 8, 0, 1, 1, 7, 0, 0, 0, 1, 0
Offset: 1

Views

Author

Gus Wiseman, Oct 10 2021

Keywords

Comments

An ordered factorization of n is a 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)).

Examples

			The a(n) ordered factorizations for n = 16, 32, 36, 48, 64, 96:
  4*4       8*4       6*6       12*4      8*8           24*4
  8*2       16*2      12*3      24*2      16*4          48*2
  2*2*2*2   2*2*4*2   18*2      2*2*6*2   32*2          3*2*8*2
            4*2*2*2   2*2*3*3   3*2*4*2   2*2*4*4       4*2*6*2
                      2*3*3*2   4*2*3*2   2*2*8*2       6*2*4*2
                      3*2*2*3   6*2*2*2   2*4*4*2       8*2*3*2
                      3*3*2*2             4*2*2*4       12*2*2*2
                                          4*2*4*2       2*2*12*2
                                          4*4*2*2
                                          8*2*2*2
                                          2*2*2*2*2*2
		

Crossrefs

Positions of 0's are A005117 \ {2}.
The restriction to powers of 2 is A027306.
Heinz numbers of partitions of this type are A028260 /\ A347457.
Positions of 3's appear to be A030514.
Positions of 1's are 1 and A082293.
Allowing non-integer alternating product gives A174725, unordered A339846.
The odd-length version is A347049.
The unordered version is A347438, reverse A347439.
Allowing any length gives A347463.
Partitions of this type are counted by A347704, reverse A035363.
A001055 counts factorizations (strict A045778, ordered A074206).
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A119620 counts partitions with alternating product 1, ranked by A028982.
A273013 counts ordered factorizations of n^2 with alternating product 1.
A339890 counts odd-length factorizations, ordered A174726.
A347050 = factorizations with alternating permutation, complement A347706.
A347437 = factorizations with integer alternating product, reverse A347442.
A347446 = partitions with integer alternating product, reverse A347445.
A347460 counts possible alternating products of factorizations.

Programs

  • Mathematica
    ordfacs[n_]:=If[n<=1,{{}},Join@@Table[Prepend[#,d]&/@ordfacs[n/d],{d,Rest[Divisors[n]]}]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Table[Length[Select[ordfacs[n],EvenQ[Length[#]]&&IntegerQ[altprod[#]]&]],{n,100}]
  • PARI
    A347048(n, m=n, ap=1, e=0) = if(1==n,!(e%2) && 1==numerator(ap), sumdiv(n, d, if(d>1, A347048(n/d, d, ap * d^((-1)^e), 1-e)))); \\ Antti Karttunen, Jul 28 2024

Formula

a(n) = A347463(n) - A347049(n).

Extensions

Data section extended up to a(105) by Antti Karttunen, Jul 28 2024
Previous Showing 11-20 of 22 results. Next