A174725 -id:A174725 - cp's OEIS frontend

A174726 a(n) = (A002033(n-1) - A008683(n))/2.

0, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 1, 4, 1, 4, 1, 4, 1, 1, 1, 10, 1, 1, 2, 4, 1, 7, 1, 8, 1, 1, 1, 13, 1, 1, 1, 10, 1, 7, 1, 4, 4, 1, 1, 24, 1, 4, 1, 4, 1, 10, 1, 10, 1, 1, 1, 22, 1, 1, 4, 16, 1, 7, 1, 4, 1, 7, 1, 38, 1, 1, 4, 4, 1

Offset: 1

Mats Granvik, Mar 28 2010

nonn

a(n) is the number of permutation matrices with a negative contribution to the determinant that is the Möbius function. See A174725 for how the determinant is defined. - Mats Granvik, May 26 2017
From Gus Wiseman, Jan 04 2021: (Start)
Also the number of ordered factorizations of n into an odd number of factors > 1. The unordered case is A339890. For example, the a(n) factorizations for n = 8, 12, 24, 30, 32, 36 are:
  (8)      (12)     (24)     (30)     (32)         (36)
  (2*2*2)  (2*2*3)  (2*2*6)  (2*3*5)  (2*2*8)      (2*2*9)
           (2*3*2)  (2*3*4)  (2*5*3)  (2*4*4)      (2*3*6)
           (3*2*2)  (2*4*3)  (3*2*5)  (2*8*2)      (2*6*3)
                    (2*6*2)  (3*5*2)  (4*2*4)      (2*9*2)
                    (3*2*4)  (5*2*3)  (4*4*2)      (3*2*6)
                    (3*4*2)  (5*3*2)  (8*2*2)      (3*3*4)
                    (4*2*3)           (2*2*2*2*2)  (3*4*3)
                    (4*3*2)                        (3*6*2)
                    (6*2*2)                        (4*3*3)
                                                   (6*2*3)
                                                   (6*3*2)
                                                   (9*2*2)
(End)

Mats Granvik, Table of n, a(n) for n = 1..10000

The even version is A174725.
The unordered case is A339890, with even version A339846.
A001055 counts factorizations, with strict case A045778.
A074206 counts ordered factorizations, with strict case A254578.
A251683 counts ordered factorizations by product and length.
A340102 counts odd-length factorizations into odd factors.
Other cases of odd length:
- A024429 counts set partitions of odd length.
- A027193 counts 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.
- A332304 counts strict compositions of odd length.
Cf. A002033, A024430, A027187, A050320, A052841, A058695, A160786, A316439.

ordfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#1,d]&)/@ordfacs[n/d],{d,Rest[Divisors[n]]}]];
Table[Length[Select[ordfacs[n],OddQ@*Length]],{n,100}] (* Gus Wiseman, Jan 04 2021 *)

a(n) = (A002033(n-1) - A008683(n))/2. - Mats Granvik, May 26 2017

For n > 0, a(n) + A174725(n) = A074206(n). - Gus Wiseman, Jan 04 2021

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

Elliott Line, May 13 2016

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)

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

Alois P. Heinz, Table of n, a(n) for n = 1..10000
Matthew C. Lettington and Karl Michael Schmidt, Divisor Functions and the Number of Sum Systems, arXiv:1910.02455 [math.NT], 2019.
S. Harry White, Reversible Squares

Cf. A111588, A131514, A360098.
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.
Cf. A062312, A347437, A347439, A347440, A347442, A347456, A347459, A347464.

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

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

@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

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

A323910 Dirichlet inverse of the deficiency of n, A033879.

Original entry on oeis.org

1, -1, -2, 0, -4, 4, -6, 0, -1, 6, -10, 2, -12, 8, 10, 0, -16, 1, -18, 2, 14, 12, -22, 4, -3, 14, -2, 2, -28, -16, -30, 0, 22, 18, 26, 4, -36, 20, 26, 4, -40, -24, -42, 2, 4, 24, -46, 8, -5, -1, 34, 2, -52, 0, 42, 4, 38, 30, -58, 2, -60, 32, 6, 0, 50, -40, -66, 2, 46, -40, -70, 12, -72, 38, 2, 2, 62, -48, -78, 8, -4, 42, -82, -2, 66, 44, 58, 4, -88, 2, 74, 2

Offset: 1

Antti Karttunen, Feb 12 2019

sign

Antti Karttunen, Table of n, a(n) for n = 1..20000
Jon Maiga, Computer-generated formulas for A323910, Sequence Machine.

Cf. A033879, A323911, A323912, A359549 (parity of terms).

Sequences that appear in the convolution formulas: A002033, A008683, A023900, A055615, A046692, A067824, A074206, A174725, A191161, A327960, A328722, A330575, A345182, A349341, A346246, A349387.

b[n_] := 2 n - DivisorSigma[1, n];
a[n_] := a[n] = If[n == 1, 1, -Sum[b[n/d] a[d], {d, Most@ Divisors[n]}]];
Array[a, 100] (* Jean-François Alcover, Feb 17 2020 *)

up_to = 16384;
DirInverse(v) = { my(u=vector(#v)); u[1] = (1/v[1]); for(n=2, #v, u[n] = -sumdiv(n, d, if(dA033879(n) = (2*n-sigma(n));
v323910 = DirInverse(vector(up_to,n,A033879(n)));
A323910(n) = v323910[n];

a(1) = 1, and for n > 1, a(n) = -Sum_{d|n, dA033879(n/d) * a(d).
From Antti Karttunen, Nov 14 2024: (Start)
Following convolution formulas have been conjectured for this sequence by Sequence Machine, with each one giving the first 10000 terms correctly:
a(n) = Sum_{d|n} A046692(d)*A067824(n/d).
a(n) = Sum_{d|n} A055615(d)*A074206(n/d).
a(n) = Sum_{d|n} A023900(d)*A174725(n/d).
a(n) = Sum_{d|n} A008683(d)*A323912(n/d).
a(n) = Sum_{d|n} A191161(d)*A327960(n/d).
a(n) = Sum_{d|n} A328722(d)*A330575(n/d).
a(n) = Sum_{d|n} A345182(d)*A349341(n/d).
a(n) = Sum_{d|n} A346246(d)*A349387(n/d).
a(n) = Sum_{d|n} A002033(d-1)*A055615(n/d).
(End)

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

Gus Wiseman, Jan 08 2021

nonn

			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

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.
Cf. A000700, A027187, A078408, A174725, A236914.

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

A322353 Number of factorizations of n into distinct semiprimes; a(1) = 1 by convention.

Original entry on oeis.org

1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 1, 1, 0, 2, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 2, 1, 1, 1, 1, 0, 2, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0

Offset: 1

Antti Karttunen, Dec 06 2018

nonn

A semiprime (A001358) is a product of any two prime numbers. In the even case, these factorizations have A001222(n)/2 factors. - Gus Wiseman, Dec 31 2020

Records 1, 2, 3, 4, 5, 9, 13, 15, 17, ... occur at 1, 60, 210, 840, 1260, 4620, 27720, 30030, 69300, ...

			a(4) = 1, as there is just one way to factor 4 into distinct semiprimes, namely as {4}.
From _Gus Wiseman_, Dec 31 2020: (Start)
The a(n) factorizations for n = 60, 210, 840, 1260, 4620, 12600, 18480:
  4*15   6*35    4*6*35    4*9*35    4*15*77    4*6*15*35    4*6*10*77
  6*10   10*21   4*10*21   4*15*21   4*21*55    4*6*21*25    4*6*14*55
         14*15   4*14*15   6*10*21   4*33*35    4*9*10*35    4*6*22*35
                 6*10*14   6*14*15   6*10*77    4*9*14*25    4*10*14*33
                           9*10*14   6*14*55    4*10*15*21   4*10*21*22
                                     6*22*35    6*10*14*15   4*14*15*22
                                     10*14*33                6*10*14*22
                                     10*21*22
                                     14*15*22
(End)

Unlabeled multiset partitions of this type are counted by A007717.
The version for partitions is A112020, or A101048 without distinctness.
The non-strict version is A320655.
Positions of zeros include A320892.
Positions of nonzero terms are A320912.
The case of squarefree factors is A339661, or A320656 without distinctness.
Allowing prime factors gives A339839, or A320732 without distinctness.
A322661 counts loop-graphs, ranked by A320461.
A001055 counts factorizations, with strict case A045778.
A001358 lists semiprimes, with squarefree case A006881.
A027187 counts partitions of even length, ranked by A028260.
A037143 lists primes and semiprimes.
A338898/A338912/A338913 give the prime indices of semiprimes.
A339846 counts even-length factorizations, with ordered version A174725.
Cf. A001221, A006125, A006129, A028260, A320893, A338915, A339841.

Table[Count[Subsets[Select[Divisors[n], PrimeOmega[#] == 2 &]], ?(Times @@ # == n &)], {n, 105}] (* _Michael De Vlieger, Dec 11 2020 *)

A322353(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((2==bigomega(d)&&(d<=m)), s += A322353(n/d, d-1))); (s)); \\ Antti Karttunen, Dec 10 2020

a(n) = Sum_{d|n} (-1)^A001222(d) * A339839(n/d). - Gus Wiseman, Dec 31 2020

A340785 Number of factorizations of 2n into even factors > 1.

Original entry on oeis.org

1, 2, 1, 3, 1, 2, 1, 5, 1, 2, 1, 4, 1, 2, 1, 7, 1, 3, 1, 4, 1, 2, 1, 7, 1, 2, 1, 4, 1, 3, 1, 11, 1, 2, 1, 6, 1, 2, 1, 7, 1, 3, 1, 4, 1, 2, 1, 12, 1, 3, 1, 4, 1, 3, 1, 7, 1, 2, 1, 7, 1, 2, 1, 15, 1, 3, 1, 4, 1, 3, 1, 12, 1, 2, 1, 4, 1, 3, 1, 12, 1, 2, 1, 7, 1

Offset: 1

Gus Wiseman, Jan 30 2021

nonn

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

Antti Karttunen, Table of n, a(n) for n = 1..20000

Note: A-numbers of Heinz-number sequences are in parentheses below.
The version for partitions is A035363 (A066207).
The odd version is A340101.
The even length case is A340786.
- Factorizations -
A001055 counts factorizations, with strict case A045778.
A340653 counts balanced factorizations.
A340831/A340832 count factorizations with odd maximum/minimum.
A316439 counts factorizations by product and length
A340102 counts odd-length factorizations of odd numbers into odd factors.
- Even -
A027187 counts partitions of even length/maximum (A028260/A244990).
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A236913 counts partitions of even length and sum.
A340601 counts partitions of even rank (A340602).
Cf. A001147, A001222, A050320, A066208, A160786, A174725, A320655, A320656, A339890, A340851.
Even bisection of A349906.

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],Select[#,OddQ]=={}&]],{n,2,100,2}]

A349906(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&!(d%2), s += A349906(n/d, d))); (s));
A340785(n) = A349906(2*n); \\ Antti Karttunen, Dec 13 2021

a(n) = A349906(2*n). - 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

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

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

			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)

Vincenzo Librandi, Table of n, a(n) for n = 0..200
J.-C. Aval, V. Féray, J.-C. Novelli, J.-Y. Thibon, Quasi-symmetric functions as polynomial functions on Young diagrams, arXiv preprint arXiv:1312.2727, 2013
Gottfried Helms, Discussion of a problem concerning summing of like powers

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.
Cf. A000700, A026424, A027187, A028260, A078408, A174725, A236914, A340101.

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

Table[Sum[Binomial[n, k](-1)^(n-k)Sum[i! StirlingS2[k, i], {i, 1, k}], {k, 0, n}], {n, 0, 20}]

a(n)=if(n<0,0,n!*polcoeff(subst(y/(1-y^2),y,exp(x+x*O(x^n))-1),n))

{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 */

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

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

A191161 Hypersigma(n), definition 2: sum of the divisors of n plus the recursive sum of the divisors of the proper divisors.

Original entry on oeis.org

1, 4, 5, 12, 7, 22, 9, 32, 19, 30, 13, 72, 15, 38, 37, 80, 19, 90, 21, 96, 47, 54, 25, 208, 39, 62, 65, 120, 31, 178, 33, 192, 67, 78, 65, 316, 39, 86, 77, 272, 43, 222, 45, 168, 147, 102, 49, 560, 67, 174, 97, 192, 55

Offset: 1

Alonso del Arte, May 26 2011

In wanting to ensure the definition was not arbitrary, I initially thought that 1s had to stop the recursion. But as T. D. Noe showed me, this doesn't have to be the case: the 1s can be included in the recursion.

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Jon Maiga, Computer-generated formulas for A191161, Sequence Machine.

Cf. A000203, A191150, A202687, A255242, A378211 (Dirichlet inverse).

Sequences that appear in the convolution formulas: A000010, A000203, A007429, A038040, A060640, A067824, A074206, A174725, A253249, A323910, A323912, A330575.

hsTD[n_] := hsTD[n] = Module[{d = Divisors[n]}, Total[d] + Total[hsTD /@ Most[d]]]; Table[hsTD[n], {n, 100}] (* From T. D. Noe *)

a(n)=sumdiv(n,d,if(dCharles R Greathouse IV, Dec 20 2011

a(n) = sigma(n) + sum_{d | n, d < n} a(d). - Charles R Greathouse IV, Dec 20 2011
From Antti Karttunen, Nov 22 2024: (Start)
Following formulas were conjectured by Sequence Machine:
For n > 1, a(n) = A191150(n) + A074206(n).
a(n) = A330575(n) + A255242(n) = 2*A255242(n) + n = 2*A330575(n) - n.
a(n) = Sum_{d|n} A330575(d).
a(n) = Sum_{d|n} d*A067824(n/d).
a(n) = Sum_{d|n} A000203(d)*A074206(n/d).
a(n) = Sum_{d|n} A007429(d)*A174725(n/d).
a(n) = Sum_{d|n} A000010(d)*A253249(n/d).
a(n) = Sum_{d|n} A038040(d)*A323912(n/d).
a(n) = Sum_{d|n} A060640(d)*A323910(n/d).
(End)

A340786 Number of factorizations of 4n into an even number of even factors > 1.

Original entry on oeis.org

1, 1, 1, 3, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 6, 1, 3, 1, 4, 2, 2, 1, 6, 2, 2, 2, 4, 1, 4, 1, 7, 2, 2, 2, 7, 1, 2, 2, 6, 1, 4, 1, 4, 3, 2, 1, 10, 2, 3, 2, 4, 1, 4, 2, 6, 2, 2, 1, 8, 1, 2, 3, 12, 2, 4, 1, 4, 2, 4, 1, 10, 1, 2, 3, 4, 2, 4, 1, 10, 3, 2, 1, 8, 2, 2

Offset: 1

Gus Wiseman, Jan 31 2021

nonn

			The a(n) factorizations for n = 6, 12, 24, 36, 60, 80, 500:
  4*6   6*8      2*48      2*72      4*60      4*80          40*50
  2*12  2*24     4*24      4*36      6*40      8*40          4*500
        4*12     6*16      6*24      8*30      10*32         8*250
        2*2*2*6  8*12      8*18      10*24     16*20         10*200
                 2*2*4*6   12*12     12*20     2*160         20*100
                 2*2*2*12  2*2*6*6   2*120     2*2*2*40      2*1000
                           2*2*2*18  2*2*2*30  2*2*4*20      2*2*10*50
                                     2*2*6*10  2*2*8*10      2*2*2*250
                                               2*4*4*10      2*10*10*10
                                               2*2*2*2*2*10

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

Note: A-numbers of Heinz-number sequences are in parentheses below.
Positions of ones are 1 and A000040, or A008578.
A version for partitions is A027187 (A028260).
Allowing odd length gives A108501 (bisection of A340785).
Allowing odd factors gives A339846.
An odd version is A340102.
- Factorizations -
A001055 counts factorizations, with strict case A045778.
A316439 counts factorizations by product and length.
A340101 counts factorizations into odd factors.
A340653 counts balanced factorizations.
A340831/A340832 count factorizations with odd maximum/minimum.
- Even -
A027187 counts partitions of even maximum (A244990).
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A236913 counts partitions of even length and sum (A340784).
A340601 counts partitions of even rank (A340602).
Cf. A001147, A001222, A001749, A035363, A050320, A066207, A066208, A160786, A174725, A320655, A320656, A339890, A340851.

g:= proc(n, m, p)
 option remember;
 local F,r,x,i;
 # number of factorizations of n into even factors > m with number of factors == p (mod 2)
 if n = 1 then if p = 0 then return 1 else return 0 fi fi;
 if m > n  or n::odd then return 0 fi;
 F:= sort(convert(select(t -> t > m and t::even, numtheory:-divisors(n)),list));
 r:= 0;
 for x in F do
   for i from 1 while n mod x^i = 0 do
     r:= r + procname(n/x^i, x, (p+i) mod 2)
 od od;
 r
end proc:
f:= n -> g(4*n, 1, 0):
map(f, [$1..100]); # Robert Israel, Mar 16 2023

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[Select[facs[4n],EvenQ[Length[#]]&&Select[#,OddQ]=={}&]],{n,100}]

A340786aux(n, m=n, p=0) = if(1==n, (0==p), my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&!(d%2), s += A340786aux(n/d, d, 1-p))); (s));
A340786(n) = A340786aux(4*n); \\ Antti Karttunen, Apr 14 2022

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

Paolo P. Lava, Feb 19 2015

nonn

a(n) = 1 if n is prime.

			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.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Paolo P. Lava)
Jon Maiga, Computer-generated formulas for A255242, Sequence Machine.

Cf. A001047, A001065, A001787, A006516, A016127, A016130, A016135, A255243, A330575.

Sequences that appear in the convolution formulas: A000203, A001065, A051953, A062790, A067824, A074206, A157658, A174725, A174726, A191161, A211779, A245211, A323910, A253249.

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

f[s_] := Flatten[Most[Divisors[#]] & /@ s]; a[n_] := Total@Flatten[FixedPointList[ f, {n}]] - n; Array[a, 100] (* Amiram Eldar, Apr 06 2019 *)

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

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)

cp's OEIS Frontend

A174726 a(n) = (A002033(n-1) - A008683(n))/2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

SageMath

Formula

A323910 Dirichlet inverse of the deficiency of n, A033879.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

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

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

A322353 Number of factorizations of n into distinct semiprimes; a(1) = 1 by convention.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A340785 Number of factorizations of 2n into even factors > 1.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Sage

Formula