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-7 of 7 results.

A350847 Number of even parts in the conjugate of the integer partition with Heinz number n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Mar 14 2022

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so a(n) counts even prime indices of n.

Crossrefs

Positions of first appearances are A001248.
The triangular version is A116482.
Positions of zeros are A346635.
Subtracting from the number of odd conjugate parts gives A350941.
Subtracting from the number of odd parts gives A350942.
Subtracting from the number of even parts gives A350950.
There are four statistics:
- A257991 = # of odd parts, conjugate A344616.
- A257992 = # of even parts, conjugate A350847 (this sequence).
There are six possible pairings of statistics:
- A325698: # of even parts = # of odd parts, counted by A045931.
- A349157: # of even parts = # of odd conjugate parts, counted by A277579.
- A350848: # of even conj parts = # of odd conj parts, counted by A045931.
- A350943: # of even conjugate parts = # of odd parts, counted by A277579.
- A350944: # of odd parts = # of odd conjugate parts, counted by A277103.
- A350945: # of even parts = # of even conjugate parts, counted by A350948.
There are three possible double-pairings of statistics:
- A350946, counted by A351977.
- A350949, counted by A351976.
- A351980, counted by A351981.
The case of all four statistics equal is A350947, counted by A351978.
A056239 adds up prime indices, counted by A001222, row sums of A112798.
A122111 represents partition conjugation using Heinz numbers.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Table[Count[conj[primeMS[n]],_?EvenQ],{n,100}]

Formula

a(n) = A344616(n) - A350941(n).
a(n) = A257992(A122111(n)).

A347458 Number of factorizations of n^2 with integer alternating product.

Original entry on oeis.org

1, 2, 2, 4, 2, 6, 2, 8, 4, 6, 2, 17, 2, 6, 6, 15, 2, 17, 2, 16, 6, 6, 2, 41, 4, 6, 8, 16, 2, 31, 2, 27, 6, 6, 6, 56, 2, 6, 6, 39, 2, 31, 2, 17, 17, 6, 2, 90, 4, 17, 6, 17, 2, 41, 6, 39, 6, 6, 2, 105, 2, 6, 17, 48, 6, 31, 2, 17, 6, 31, 2, 148, 2, 6, 17, 17, 6, 32, 2, 86, 15, 6, 2, 107, 6, 6, 6, 40, 2, 109, 6, 17
Offset: 1

Views

Author

Gus Wiseman, Sep 21 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.
The even-length case, the case of alternating product 1, and the case of alternating sum 0 are all counted by A001055.

Examples

			The a(2) = 2 through a(8) = 8 factorizations:
  4     9     16        25    36        49    64
  2*2   3*3   4*4       5*5   6*6       7*7   8*8
              2*2*4           2*2*9           2*4*8
              2*2*2*2         2*3*6           4*4*4
                              3*3*4           2*2*16
                              2*2*3*3         2*2*4*4
                                              2*2*2*2*4
                                              2*2*2*2*2*2
		

Crossrefs

Positions of 2's are A000040, squares A001248.
The restriction to powers of 2 is A344611.
This is the restriction to perfect squares of A347437.
The nonsquared even-length version is A347438.
The reciprocal version is A347459, non-squared A347439.
The additive version (partitions) is the even bisection of A347446.
The nonsquared ordered version is A347463.
The case of alternating product 1 in the ordered version is A347464.
Allowing any alternating product gives A347466.
A000290 lists squares, complement A000037.
A001055 counts factorizations.
A046099 counts factorizations with no alternating permutations.
A071321 gives the alternating sum of prime factors of n (reverse: A071322).
A273013 counts ordered factorizations of n^2 with alternating product 1.
A347460 counts possible alternating products of factorizations.
A339846 counts even-length factorizations.
A339890 counts odd-length factorizations.
A347457 ranks partitions with integer alternating product.
Apparently, A006881 gives the positions of 6's. - Antti Karttunen, Oct 22 2023

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[facs[n^2],IntegerQ[altprod[#]]&]],{n,100}]
  • PARI
    A347437(n, m=n, ap=1, e=0) = if(1==n, if(e%2, 1==denominator(ap), 1==numerator(ap)), sumdiv(n, d, if((d>1)&&(d<=m), A347437(n/d, d, ap * d^((-1)^e), 1-e))));
    A347458(n) = A347437(n*n); \\ Antti Karttunen, Oct 22 2023

Formula

a(2^n) = A344611(n).
a(n) = A347437(n^2).

Extensions

Data section extended up to a(92) by Antti Karttunen, Oct 22 2023

A347459 Number of factorizations of n^2 with integer reciprocal alternating product.

Original entry on oeis.org

1, 1, 1, 3, 1, 4, 1, 6, 3, 4, 1, 11, 1, 4, 4, 12, 1, 11, 1, 12, 4, 4, 1, 28, 3, 4, 6, 12, 1, 19, 1, 22, 4, 4, 4, 38, 1, 4, 4, 29, 1, 21, 1, 12, 11, 4, 1, 65, 3, 11, 4, 12, 1, 29, 4, 29, 4, 4, 1, 71, 1, 4, 11, 40, 4, 22, 1, 12, 4, 18, 1, 107, 1, 4, 11, 12, 4, 22, 1, 66, 12, 4, 1, 76, 4, 4, 4, 30, 1, 71, 4, 12, 4, 4, 4, 141
Offset: 1

Views

Author

Gus Wiseman, Sep 22 2021

Keywords

Comments

We define the reciprocal alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^i).
A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.
All such factorizations have even length.
Image appears to be 1, 3, 4, 6, 11, ... , missing some numbers such as 2, 5, 7, 8, 9, ...
The case of alternating product 1, the case of alternating sum 0, and the reverse version are all counted by A001055.

Examples

			The a(2) = 1 through a(10) = 4 factorizations:
    2*2  3*3  2*8      5*5  6*6      7*7  8*8          9*9      2*50
              4*4           2*18          2*32         3*27     5*20
              2*2*2*2       3*12          4*16         3*3*3*3  10*10
                            2*2*3*3       2*2*2*8               2*2*5*5
                                          2*2*4*4
                                          2*2*2*2*2*2
		

Crossrefs

Positions of 1's are 1 and A000040, squares A001248.
The additive version (partitions) is A000041, the even bisection of A119620.
Partitions of this type are ranked by A028982 and A347451.
The restriction to powers of 2 is A236913, the even bisection of A027187.
The nonsquared nonreciprocal even-length version is A347438.
This is the restriction to perfect squares of A347439.
The nonreciprocal version is A347458, non-squared A347437.
A000290 lists squares, complement A000037.
A001055 counts factorizations.
A046099 counts factorizations with no alternating permutations.
A273013 counts ordered factorizations of n^2 with alternating product 1.
A347460 counts possible alternating products of factorizations.
A339846 counts even-length factorizations.
A339890 counts odd-length factorizations.
A347457 ranks partitions 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]]}]];
    recaltprod[q_]:=Product[q[[i]]^(-1)^i,{i,Length[q]}];
    Table[Length[Select[facs[n^2],IntegerQ[recaltprod[#]]&]],{n,100}]
  • PARI
    A347439(n, m=n, ap=1, e=0) = if(1==n, !(e%2) && 1==denominator(ap), sumdiv(n, d, if(d>1 && d<=m, A347439(n/d, d, ap * d^((-1)^e), 1-e))));
    A347459(n) = A347439(n^2); \\ Antti Karttunen, Jul 28 2024

Formula

a(2^n) = A236913(n).
a(n) = A347439(n^2).

Extensions

Data section extended up to a(96) by Antti Karttunen, Jul 28 2024

A347466 Number of factorizations of n^2.

Original entry on oeis.org

1, 2, 2, 5, 2, 9, 2, 11, 5, 9, 2, 29, 2, 9, 9, 22, 2, 29, 2, 29, 9, 9, 2, 77, 5, 9, 11, 29, 2, 66, 2, 42, 9, 9, 9, 109, 2, 9, 9, 77, 2, 66, 2, 29, 29, 9, 2, 181, 5, 29, 9, 29, 2, 77, 9, 77, 9, 9, 2, 269, 2, 9, 29, 77, 9, 66, 2, 29, 9, 66, 2, 323, 2, 9, 29, 29
Offset: 1

Views

Author

Gus Wiseman, Sep 23 2021

Keywords

Comments

A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.

Examples

			The a(1) = 1 through a(8) = 11 factorizations:
  ()  (4)    (9)    (16)       (25)   (36)       (49)   (64)
      (2*2)  (3*3)  (2*8)      (5*5)  (4*9)      (7*7)  (8*8)
                    (4*4)             (6*6)             (2*32)
                    (2*2*4)           (2*18)            (4*16)
                    (2*2*2*2)         (3*12)            (2*4*8)
                                      (2*2*9)           (4*4*4)
                                      (2*3*6)           (2*2*16)
                                      (3*3*4)           (2*2*2*8)
                                      (2*2*3*3)         (2*2*4*4)
                                                        (2*2*2*2*4)
                                                        (2*2*2*2*2*2)
		

Crossrefs

Positions of 2's are the primes (A000040), which have squares A001248.
The restriction to powers of 2 is A058696.
The additive version (partitions) is A072213.
The case of integer alternating product is A347459, nonsquared A347439.
A000290 lists squares, complement A000037.
A001055 counts factorizations.
A339846 counts even-length factorizations.
A339890 counts odd-length factorizations.
A347050 = factorizations with alternating permutation, complement A347706.

Programs

  • Maple
    b:= proc(n, k) option remember; `if`(n>k, 0, 1)+`if`(isprime(n), 0,
          add(`if`(d>k, 0, b(n/d, d)), d=numtheory[divisors](n) minus {1, n}))
        end:
    a:= proc(n) option remember; b((l-> mul(ithprime(i)^l[i], i=1..nops(l)))(
          sort(map(i-> i[2], ifactors(n^2)[2]), `>`))$2)
        end:
    seq(a(n), n=1..76);  # Alois P. Heinz, Oct 14 2021
  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[facs[n^2]],{n,25}]
  • PARI
    A001055(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A001055(n/d, d))); (s));
    A347466(n) = A001055(n^2); \\ Antti Karttunen, Oct 13 2021

Formula

a(n) = A001055(A000290(n)).

A342768 a(n) = A342767(n, n).

Original entry on oeis.org

1, 2, 3, 8, 5, 12, 7, 32, 27, 20, 11, 48, 13, 28, 45, 128, 17, 108, 19, 80, 63, 44, 23, 192, 125, 52, 243, 112, 29, 180, 31, 512, 99, 68, 175, 432, 37, 76, 117, 320, 41, 252, 43, 176, 405, 92, 47, 768, 343, 500, 153, 208, 53, 972, 275, 448, 171, 116, 59, 720
Offset: 1

Views

Author

Rémy Sigrist, Apr 02 2021

Keywords

Comments

This sequence has similarities with A087019.
These are the positions of first appearances of each positive integer in A346701, and also in A346703. - Gus Wiseman, Aug 09 2021

Examples

			For n = 42:
- 42 = 2 * 3 * 7, so:
          2 3 7
        x 2 3 7
        -------
          2 3 7
        2 3 3
    + 2 2 2
    -----------
      2 2 3 3 7
- hence a(42) = 2 * 2 * 3 * 3 * 7 = 252.
		

Crossrefs

The sum of prime indices of a(n) is 2*A056239(n) - A061395(n).
The version for even indices is A129597(n) = 2*a(n) for n > 1.
The sorted version is A346635.
These are the positions of first appearances in A346701 and in A346703.
A001221 counts distinct prime factors.
A001222 counts prime factors with multiplicity.
A027193 counts partitions of odd length, ranked by A026424.
A209281 adds up the odd bisection of standard compositions (even: A346633).
A346697 adds up the odd bisection of prime indices (reverse: A346699).

Programs

  • Mathematica
    Table[n^2/FactorInteger[n][[-1,1]],{n,100}] (* Gus Wiseman, Aug 09 2021 *)
  • PARI
    See Links section.

Formula

a(n) = n iff n = 1 or n is a prime number.
a(p^k) = p^(2*k-1) for any k > 0 and any prime number p.
A007947(a(n)) = A007947(n).
A001222(a(n)) = 2*A001222(n) - 1 for any n > 1.
From Gus Wiseman, Aug 09 2021: (Start)
A001221(a(n)) = A001221(n).
If g = A006530(n) is the greatest prime factor of n, then a(n) = n^2/g.
a(n) = A129597(n)/2.
(End)

A129597 Central diagonal of array A129595.

Original entry on oeis.org

1, 4, 6, 16, 10, 24, 14, 64, 54, 40, 22, 96, 26, 56, 90, 256, 34, 216, 38, 160, 126, 88, 46, 384, 250, 104, 486, 224, 58, 360, 62, 1024, 198, 136, 350, 864, 74, 152, 234, 640, 82, 504, 86, 352, 810, 184, 94, 1536, 686, 1000, 306, 416, 106, 1944, 550, 896, 342
Offset: 1

Views

Author

Antti Karttunen, May 01 2007, based on Marc LeBrun's Jan 11 2006 message on SeqFan mailing list

Keywords

Comments

These are the positions of first appearances of each positive integer in A346704. - Gus Wiseman, Oct 16 2021

Crossrefs

a(n) = A129595(n,n).
The sum of prime indices of a(n) is 2*A056239(n) - A061395(n) + 1 for n > 1.
The version for odd indices is A342768(n) = a(n)/2 for n > 1.
Except the first term, the sorted version is 2*A346635.
These are the positions of first appearances in A346704.
A001221 counts distinct prime factors.
A001222 counts prime factors with multiplicity.
A027187 counts partitions of even length, ranked by A028260.
A346633 adds up the even bisection of standard compositions (odd: A209281).
A346698 adds up the even bisection of prime indices (reverse: A346699).

Programs

  • Mathematica
    Table[If[n==1,1,2*n^2/FactorInteger[n][[-1,1]]],{n,100}] (* Gus Wiseman, Aug 10 2021 *)
  • PARI
    A129597(n) = if(1==n, n, my(f=factor(n)); (2*n*n)/f[#f~, 1]); \\ Antti Karttunen, Oct 16 2021

Formula

From Gus Wiseman, Aug 10 2021: (Start)
For n > 1, A001221(a(n)) = A099812(n).
If g = A006530(n) is the greatest prime factor of n > 1, then a(n) = 2n^2/g.
a(n) = A100484(A000720(n)) = 2n iff n is prime.
a(n > 1) = 2*A342768(n).
(End)

A352143 Numbers whose prime indices and conjugate prime indices are all odd.

Original entry on oeis.org

1, 2, 5, 8, 11, 17, 20, 23, 31, 32, 41, 44, 47, 59, 67, 68, 73, 80, 83, 92, 97, 103, 109, 124, 125, 127, 128, 137, 149, 157, 164, 167, 176, 179, 188, 191, 197, 211, 227, 233, 236, 241, 257, 268, 269, 272, 275, 277, 283, 292, 307, 313, 320, 331, 332, 347, 353
Offset: 1

Views

Author

Gus Wiseman, Mar 18 2022

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798, sum A056239, length A001222.
A number's prime signature is the sequence of positive exponents in its prime factorization, which is row n of A124010, length A001221, sum A001222.
These are the Heinz numbers of integer partitions whose parts and conjugate parts are all odd. They are counted by A053253.

Examples

			The terms together with their prime indices begin:
   1: {}
   2: {1}
   5: {3}
   8: {1,1,1}
  11: {5}
  17: {7}
  20: {1,1,3}
  23: {9}
  31: {11}
  32: {1,1,1,1,1}
  41: {13}
  44: {1,1,5}
  47: {15}
  59: {17}
  67: {19}
  68: {1,1,7}
  73: {21}
  80: {1,1,1,1,3}
		

Crossrefs

The restriction to primes is A031368.
These partitions appear to be counted by A053253.
The even version is A066207^2.
For even instead of odd conjugate parts we get A066208^2.
The first condition alone (all odd indices) is A066208, counted by A000009.
The second condition alone is A346635, counted by A000009.
A055922 counts partitions with odd multiplicities, ranked by A268335.
A066207 = indices all even, counted by A035363 (complement A086543).
A109297 = same indices as exponents, counted by A114640.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A162642 counts odd prime exponents, even A162641.
A238745 gives the Heinz number of the conjugate prime signature.
A257991 counts odd indices, even A257992.
A258116 ranks strict partitions with all odd parts, even A258117.
A351979 = odd indices and even multiplicities, counted by A035457.
A352140 = even indices and odd multiplicities, counted by A055922 aerated.
A352141 = even indices and even multiplicities, counted by A035444.
A352142 = odd indices and odd multiplicities, counted by A117958.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Select[Range[100],And@@OddQ/@primeMS[#]&&And@@OddQ/@conj[primeMS[#]]&]

Formula

Intersection of A066208 and A346635.
Showing 1-7 of 7 results.