A007304 -id:A007304 - cp's OEIS frontend

A211540 Number of ordered triples (w,x,y) with all terms in {1..n} and 2w = 3x + 4y.

0, 0, 0, 0, 0, 1, 1, 2, 3, 4, 5, 7, 8, 10, 12, 14, 16, 19, 21, 24, 27, 30, 33, 37, 40, 44, 48, 52, 56, 61, 65, 70, 75, 80, 85, 91, 96, 102, 108, 114, 120, 127, 133, 140, 147, 154, 161, 169, 176, 184, 192, 200, 208, 217, 225, 234, 243, 252, 261, 271, 280, 290

Offset: 0

Clark Kimberling, Apr 15 2012

For a guide to related sequences, see A211422.
Also the number of partitions of n+1 into three parts, where each part > 1. - Peter Woodward, May 25 2015
a(n) is also equal to the number of partitions of n+4 into three distinct parts, where each part > 1. - Giovanni Resta, May 26 2015
Number of different distributions of n+1 identical balls in 3 boxes as x,y,z where 0 < x < y < z. - Ece Uslu and Esin Becenen, Dec 31 2015
After the first three terms, partial sums of A008615. - Robert Israel, Dec 31 2015
For n >= 2, also the number of partitions of n - 2 into 3 parts. The Heinz numbers of these partitions are given by A014612. - Gus Wiseman, Oct 11 2020

			a(5) = a(6) = 1 with only one ordered triple (5,2,1). - _Michael Somos_, Feb 02 2015
a(11) = 5 Number of different distributions of 11 identical balls in 3 boxes as x,y and z where 0 < x < y < z. - _Ece Uslu_, Esin Becenen, Dec 31 2015
a(1) = a(2) = a(3) = a(4) = a(5) = 0, since with fewer than 6 identical balls there is no such distribution with 3 boxes that holds for 0 < x < y < z. - _Ece Uslu_, Esin Becenen, Dec 31 2015
G.f.: x^5 + x^6 + 2*x^7 + 3*x^8 + 4*x^9 + 5*x^10 + 7*x^11 + 8*x^12 + ...
From _Gus Wiseman_, Oct 11 2020: (Start)
The a(5) = 1 through a(15) = 14 partitions of n + 1 into three parts > 1 [Woodward] are the following (A = 10, B = 11, C = 12). The ordered version is A000217(n - 4) and the Heinz numbers are A046316.
  222  322  332  333  433  443  444  544  554  555  655
            422  432  442  533  543  553  644  654  664
                 522  532  542  552  643  653  663  754
                      622  632  633  652  662  744  763
                           722  642  733  743  753  772
                                732  742  752  762  844
                                822  832  833  843  853
                                     922  842  852  862
                                          932  933  943
                                          A22  942  952
                                               A32  A33
                                               B22  A42
                                                    B32
                                                    C22
The a(5) = 1 through a(15) = 14 partitions of n + 4 into three distinct parts > 1 [Resta] are the following (A = 10, B = 11, C = 12, D = 13, E = 14). The ordered version is A211540*6 and the Heinz numbers are A046389.
  432  532  542  543  643  653  654  754  764  765  865
            632  642  652  743  753  763  854  864  874
                 732  742  752  762  853  863  873  964
                      832  842  843  862  872  954  973
                           932  852  943  953  963  982
                                942  952  962  972  A54
                                A32  A42  A43  A53  A63
                                     B32  A52  A62  A72
                                          B42  B43  B53
                                          C32  B52  B62
                                               C42  C43
                                               D32  C52
                                                    D42
                                                    E32
The a(5) = 1 through a(15) = 14 partitions of n + 1 into three distinct parts [Uslu and Becenen] are the following (A = 10, B = 11, C = 12, D = 13). The ordered version is A211540(n)*6 and the Heinz numbers are A007304.
  321  421  431  432  532  542  543  643  653  654  754
            521  531  541  632  642  652  743  753  763
                 621  631  641  651  742  752  762  853
                      721  731  732  751  761  843  862
                           821  741  832  842  852  871
                                831  841  851  861  943
                                921  931  932  942  952
                                     A21  941  951  961
                                          A31  A32  A42
                                          B21  A41  A51
                                               B31  B32
                                               C21  B41
                                                    C31
                                                    D21
(End)

Robert Israel, Table of n, a(n) for n = 0..10000
Clark Kimberling, A Combinatorial Classification of Triangle Centers on the Line at Infinity, J. Int. Seq., Vol. 22 (2019), Article 19.5.4.
Ece Uslu and Esin Becenen, Identical Object Distributions.
Index entries for linear recurrences with constant coefficients, signature (1,1,0,-1,-1,1).

Cf. A001399, A069905, A211422.
All of the following pertain to 3-part strict partitions.
- A000009 counts these partitions of any length, with non-strict version A000041.
- A007304 gives the Heinz numbers, with non-strict version A014612.
- A101271 counts the relatively prime case, with non-strict version A023023.
- A220377 counts the pairwise coprime case, with non-strict version A307719.
- A337605 counts the pairwise non-coprime case, with non-strict version A337599.
Cf. A000217, A001840, A156040, A284825, A337453, A337483, A337484, A337563.

I:=[0,0,0,0,0,1]; [n le 6 select I[n] else Self(n-1)+Self(n-2)-Self(n-4)-Self(n-5)+Self(n-6): n in [1..70]]; // Vincenzo Librandi, Dec 31 2015

f:= gfun:-rectoproc({a(n) = a(n-1)+a(n-2)-a(n-4)-a(n-5)+a(n-6),seq(a(i)=0,i=0..4),a(5)=1},a(n),remember):
seq(f(i),i=0..100); # Robert Israel, Dec 31 2015

t[n_] := t[n] = Flatten[Table[-2 w + 3 x + 4 y, {w, n}, {x, n}, {y, n}]]
c[n_] := Count[t[n], 0]
t = Table[c[n], {n, 0, 80}]  (* A211540 *)
FindLinearRecurrence[t]
LinearRecurrence[{1, 1, 0, -1, -1, 1}, {0, 0, 0, 0, 0, 1}, 70] (* Vincenzo Librandi, Dec 31 2015 *)
Table[Length[Select[IntegerPartitions[n+1,{3}],UnsameQ@@#&]],{n,0,30}] (* Gus Wiseman, Oct 05 2020 *)

{a(n) = round( (n-2)^2 / 12 )}; / * Michael Somos, Feb 02 2015 */

concat(vector(5), Vec(x^5/(1-x-x^2+x^4+x^5-x^6) + O(x^100))) \\ Altug Alkan, Jan 10 2016

a(n) = a(n-1) + a(n-2) - a(n-4) - a(n-5) + a(n-6).
a(n) = A069905(n-2) = A001399(n-5) for n >= 5. - Alois P. Heinz, Nov 03 2012
a(n) = 3*k^2-6*k+3 (for n = 6*k-3), 3*k^2-5*k+2 (for n = 6*k-2), 3*k^2-4*k+1 (for n = 6*k-1), 3*k^2-3*k+1 (for n = 6*k), 3*k^2-2*k (for n = 6*k+1), 3*k^2-k (for n = 6*k+2). - Ece Uslu, Esin Becenen, Dec 31 2015
a(n) = A004526(n-2) + a(n-2) for n > 2. - Ece Uslu, Esin Becenen, Dec 31 2015
G.f.: x^5/(1 - x - x^2 + x^4 + x^5 - x^6). - Robert Israel, Dec 31 2015
a(n) = Sum_{k=1..floor(n/3)} floor((n-k)/2)-k. - Wesley Ivan Hurt, Apr 27 2019
From Gus Wiseman, Oct 11 2020: (Start)
a(n+2) = A069905(n) = A001399(n-3) counts 3-part partitions.
a(n-1) = A069905(n-3) = A001399(n-6) counts 3-part strict partitions.
a(n-1) = A069905(n-3) = A001399(n-6) counts 3-part partitions with no 1's.
a(n-4) = A069905(n-6) = A001399(n-9) counts 3-part strict partitions with no 1's.
A000217(n-2) counts 3-part compositions.
a(n-1)*6 = A069905(n-3)*6 = A001399(n-6)*6 counts 3-part strict compositions.
A000217(n-5) counts 3-part compositions with no 1's.
a(n-4)*6 = A069905(n-6)*6 = A001399(n-9)*6 counts 3-part strict compositions with no 1's.
(End)

A284825 Number of partitions of n into 3 parts without common divisors such that every pair of them has common divisors.

Original entry on oeis.org

1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 3, 0, 5, 0, 0, 0, 1, 0, 5, 0, 1, 0, 6, 0, 6, 0, 0, 0, 4, 0, 6, 0, 0, 0, 9, 0, 2, 1, 2, 0, 9, 0, 8, 1, 1, 0, 5, 0, 14, 0, 1, 0, 15, 0, 14, 0, 0, 1, 14, 0, 14, 0, 2, 0, 15, 0, 6, 1, 2, 1, 11, 0, 18, 1, 1, 0, 10, 0, 23

Offset: 31

Alois P. Heinz, Apr 03 2017

The Heinz numbers of these partitions are the intersection of A014612 (triples), A289509 (relatively prime), and A337694 (pairwise non-coprime). - Gus Wiseman, Oct 16 2020

			a(31) = 1: [6,10,15] = [2*3,2*5,3*5].
a(41) = 2: [6,14,21], [6,15,20].
From _Gus Wiseman_, Oct 14 2020: (Start)
Selected terms and the corresponding triples:
  a(31)=1: a(41)=2: a(59)=3:  a(77)=4:  a(61)=5:  a(71)=6:
-------------------------------------------------------------
  15,10,6  20,15,6  24,20,15  39,26,12  33,22,6   39,26,6
           21,14,6  24,21,14  42,20,15  40,15,6   45,20,6
                    35,14,10  45,20,12  45,10,6   50,15,6
                              50,15,12  28,21,12  35,21,15
                                        36,15,10  36,20,15
                                                  36,21,14
(End)

Alois P. Heinz, Table of n, a(n) for n = 31..10000

Cf. A082024, A230034, A230035.
A023023 does not require pairwise non-coprimality, with strict case A101271.
A202425 and A328672 count these partitions of any length, ranked by A328868.
A284825*6 is the ordered version.
A307719 is the pairwise coprime instead of non-coprime version.
A337599 does not require relatively primality, with strict case A337605.
A200976 and A328673 count pairwise non-coprime partitions.
A289509 gives Heinz numbers of relatively prime partitions.
A327516 counts pairwise coprime partitions, ranked by A333227.
A337694 gives Heinz numbers of pairwise non-coprime partitions.
Cf. A000217, A000741, A001399, A007304, A014612, A220377, A318716, A328679, A337604, A337666.

a:= proc(n) option remember; add(add(`if`(igcd(i, j)>1
      and igcd(i, j, n-i-j)=1 and igcd(i, n-i-j)>1 and
      igcd(j, n-i-j)>1, 1, 0), j=i..(n-i)/2), i=2..n/3)
    end:
seq(a(n), n=31..137);

a[n_] := a[n] = Sum[Sum[If[GCD[i, j] > 1 && GCD[i, j, n - i - j] == 1 && GCD[i, n - i - j] > 1 && GCD[j, n - i - j] > 1, 1, 0], {j, i, (n - i)/2} ], {i, 2, n/3}];
Table[a[n], {n, 31, 137}] (* Jean-François Alcover, Jun 13 2018, from Maple *)
stabQ[u_,Q_]:=And@@Not/@Q@@@Tuples[u,2];
Table[Length[Select[IntegerPartitions[n,{3}],GCD@@#==1&&stabQ[#,CoprimeQ]&]],{n,31,100}] (* Gus Wiseman, Oct 14 2020 *)

a(n) > 0 iff n in { A230035 }.

a(n) = 0 iff n in { A230034 }.

A140106 Number of noncongruent diagonals in a regular n-gon.

Original entry on oeis.org

0, 0, 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8, 8, 9, 9, 10, 10, 11, 11, 12, 12, 13, 13, 14, 14, 15, 15, 16, 16, 17, 17, 18, 18, 19, 19, 20, 20, 21, 21, 22, 22, 23, 23, 24, 24, 25, 25, 26, 26, 27, 27, 28, 28, 29, 29, 30, 30, 31, 31, 32, 32, 33, 33, 34, 34, 35, 35, 36, 36, 37

Offset: 1

Andrew McFarland, Jun 03 2008

Number of double-stars (diameter 3 trees) with n nodes. For n >= 3, number of partitions of n-2 into two parts. - Washington Bomfim, Feb 12 2011
Number of roots of the n-th Bernoulli polynomial in the left half-plane. - Michel Lagneau, Nov 08 2012
From Gus Wiseman, Oct 17 2020: (Start)
Also the number of 3-part non-strict integer partitions of n - 1. The Heinz numbers of these partitions are given by A285508. The version for partitions of any length is A047967, with Heinz numbers A013929. The a(4) = 1 through a(15) = 6 partitions are (A = 10, B = 11, C = 12):
  111  211  221  222  322  332  333  433  443  444  544  554
            311  411  331  422  441  442  533  552  553  644
                      511  611  522  622  551  633  661  662
                                711  811  722  822  733  833
                                          911  A11  922  A22
                                                    B11  C11
(End)

			The square (n=4) has two congruent diagonals; so a(4)=1. The regular pentagon also has congruent diagonals; so a(5)=1. Among all the diagonals in a regular hexagon, there are two noncongruent ones; hence a(6)=2, etc.

G. C. Greubel, Table of n, a(n) for n = 1..5000
Washington Bomfim, Double-star corresponding to the partition [3,7]
Index entries for linear recurrences with constant coefficients, signature (1,1,-1).
Index entries for sequences related to trees

Cf. A000554, A007304, A007997, A013929, A022998.
Cf. A047967, A129194, A235451, A285508, A321773.
A001399(n-3) = A069905(n) = A211540(n+2) counts 3-part partitions.
Essentially the same as A004526.

A140106:= func< n | n eq 1 select 0 else Floor((n-2)/2) >;
[A140106(n): n in [1..80]]; // G. C. Greubel, Feb 10 2023

with(numtheory): for n from 1 to 80 do:it:=0:
y:=[fsolve(bernoulli(n,x) , x, complex)] : for m from 1 to nops(y) do : if Re(y[m])<0 then it:=it+1:else fi:od: printf(`%d, `,it):od:

a[1]=0; a[n_?OddQ] := (n-3)/2; a[n_] := n/2-1; Array[a, 100] (* Jean-François Alcover, Nov 17 2015 *)

a(n)=if(n>1,n\2-1,0) \\ Charles R Greathouse IV, Oct 16 2015

def A140106(n): return n-2>>1 if n>1 else 0 # Chai Wah Wu, Sep 18 2023

def A140106(n): return 0 if (n==1) else (n-2)//2
[A140106(n) for n in range(1,81)] # G. C. Greubel, Feb 10 2023

a(n) = floor((n-2)/2), for n > 1, otherwise 0. - Washington Bomfim, Feb 12 2011
G.f.: x^4/(1-x-x^2+x^3). - Colin Barker, Jan 31 2012
a(n) = floor(A129194(n-1)/A022998(n)), for n > 1. - Paul Curtz, Jul 23 2017
a(n) = A001399(n-3) - A001399(n-6). Compare to A007997(n) = A001399(n-3) + A001399(n-6). - Gus Wiseman, Oct 17 2020

More terms from Joseph Myers, Sep 05 2009

A336568 Numbers that are not a product of two numbers each having distinct prime multiplicities.

Original entry on oeis.org

30, 42, 66, 70, 78, 102, 105, 110, 114, 130, 138, 154, 165, 170, 174, 182, 186, 190, 195, 210, 222, 230, 231, 238, 246, 255, 258, 266, 273, 282, 285, 286, 290, 310, 318, 322, 330, 345, 354, 357, 366, 370, 374, 385, 390, 399, 402, 406, 410, 418, 420, 426, 429

Offset: 1

Gus Wiseman, Aug 06 2020

nonn

First differs from A007304 and A093599 in having 210.
First differs from A287483 in having 222.
First differs from A350352 in having 420.
A number's prime signature (row n of A124010) is the sequence of positive exponents in its prime factorization, so a number has distinct prime multiplicities iff all the exponents in its prime signature are distinct.

			Selected terms together with their prime indices:
   660: {1,1,2,3,5}
   798: {1,2,4,8}
   840: {1,1,1,2,3,4}
  3120: {1,1,1,1,2,3,6}
  9900: {1,1,2,2,3,3,5}

A336500 has zeros at these positions.
A007425 counts divisors of divisors.
A056924 counts divisors greater than their quotient.
A074206 counts strict chains of divisors from n to 1.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A336424 counts factorizations using A130091.
A336422 counts divisible pairs of divisors, both in A130091.
A327498 is the maximum divisor with distinct prime multiplicities.
A336423 counts chains in A130091, with maximal version A336569.
A336571 counts divisor sets using A130091, with maximal version A336570.
Cf. A001055, A002033, A098859, A124010, A167865, A253249, A336420.

strsig[n_]:=UnsameQ@@Last/@FactorInteger[n]
Select[Range[100],Function[n,Select[Divisors[n],strsig[#]&&strsig[n/#]&]=={}]]

A101271 Number of partitions of n into 3 distinct and relatively prime parts.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 6, 8, 9, 12, 12, 16, 15, 21, 20, 26, 25, 33, 28, 40, 36, 45, 42, 56, 44, 65, 56, 70, 64, 84, 66, 96, 81, 100, 88, 120, 90, 133, 110, 132, 121, 161, 120, 175, 140, 176, 156, 208, 153, 220, 180, 222, 196, 261, 184, 280, 225, 270, 240, 312, 230, 341, 272

Offset: 6

Vladeta Jovovic, Dec 19 2004

The Heinz numbers of these partitions are the intersection of A289509 (relatively prime), A005117 (strict), and A014612 (triple). - Gus Wiseman, Oct 15 2020

			For n=10 we have 4 such partitions: 1+2+7, 1+3+6, 1+4+5 and 2+3+5.
From _Gus Wiseman_, Oct 13 2020: (Start)
The a(6) = 1 through a(18) = 15 triples (A..F = 10..15):
  321  421  431  432  532  542  543  643  653  654  754  764  765
            521  531  541  632  651  652  743  753  763  854  873
                 621  631  641  732  742  752  762  853  863  954
                      721  731  741  751  761  843  871  872  972
                           821  831  832  851  852  943  953  981
                                921  841  932  861  952  962  A53
                                     931  941  942  961  971  A71
                                     A21  A31  951  A51  A43  B43
                                          B21  A32  B32  A52  B52
                                               A41  B41  A61  B61
                                               B31  C31  B42  C51
                                               C21  D21  B51  D32
                                                         C32  D41
                                                         C41  E31
                                                         D31  F21
                                                         E21
(End)

Fausto A. C. Cariboni, Table of n, a(n) for n = 6..10000

Cf. A023024-A023030, A000742-A000743, A023031-A023035.
A000741 is the ordered non-strict version.
A001399(n-6) does not require relative primality.
A023022 counts pairs instead of triples.
A023023 is the not necessarily strict version.
A078374 counts these partitions of any length, with Heinz numbers A302796.
A101271*6 is the ordered version.
A220377 is the pairwise coprime instead of relatively prime version.
A284825 counts the case that is pairwise non-coprime also.
A337605 is the pairwise non-coprime instead of relatively prime version.
A008289 counts strict partitions by sum and length.
A007304 gives the Heinz numbers of 3-part strict partitions.
A307719 counts 3-part pairwise coprime partitions.
A337601 counts 3-part partitions whose distinct parts are pairwise coprime.
Cf. A000010, A000217, A000837, A007360, A014612, A055684, A289509, A332004, A337452, A337563.

m:=3: with(numtheory): g:=sum(mobius(k)*x^(m*(m+1)/2*k)/Product(1-x^(i*k),i=1..m),k=1..20): gser:=series(g,x=0,80): seq(coeff(gser,x^n),n=6..77); # Emeric Deutsch, May 31 2005

Table[Length[Select[IntegerPartitions[n,{3}],UnsameQ@@#&&GCD@@#==1&]],{n,6,50}] (* Gus Wiseman, Oct 13 2020 *)

G.f. for the number of partitions of n into m distinct and relatively prime parts is Sum(moebius(k)*x^(m*(m+1)/2*k)/Product(1-x^(i*k), i=1..m), k=1..infinity).

More terms from Emeric Deutsch, May 31 2005

A068443 Triangular numbers which are the product of two primes.

Original entry on oeis.org

6, 10, 15, 21, 55, 91, 253, 703, 1081, 1711, 1891, 2701, 3403, 5671, 12403, 13861, 15931, 18721, 25651, 34453, 38503, 49141, 60031, 64261, 73153, 79003, 88831, 104653, 108811, 114481, 126253, 146611, 158203, 171991, 188191, 218791, 226801, 258121, 269011

Offset: 1

Stephan Wagler (stephanwagler(AT)aol.com), Mar 09 2002

These triangular numbers are equal to p * (2p +- 1).
All terms belong to A006987. For n>2 all terms are odd and belong to A095147. - Alexander Adamchuk, Oct 31 2006
A156592 is a subsequence. - Reinhard Zumkeller, Feb 10 2009
Triangular numbers with exactly 4 divisors. - Jon E. Schoenfield, Sep 05 2018

			Triangular numbers begin 0, 1, 3, 6, 10, ...; 6=2*3, and 2 and 3 are two distinct primes; 10=2*5, and 2 and 5 are two distinct primes, etc. - _Vladimir Joseph Stephan Orlovsky_, Feb 27 2009
a(11) = 1891 and 1891 = 31 * 61.

Jon E. Schoenfield, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)

Cf. A000217, A005382, A005384, A006987, A095147, A001358, A005385, A006881, A007304, A066179, A111206, A157342, A157344-A157347, A157352-A157357, A164977.

Cf. A075875.

q:= n-> is(numtheory[bigomega](n)=2):
select(q, [i*(i+1)/2$i=0..1000])[];  # Alois P. Heinz, Mar 27 2024

Select[ Table[ n(n + 1)/2, {n, 1000}], Apply[Plus, Transpose[ FactorInteger[ # ]] [[2]]] == 2 &]
Select[Accumulate[Range[1000]],PrimeOmega[#]==2&] (* Harvey P. Dale, Apr 03 2016 *)

list(lim)=my(v=List());forprime(p=2,(sqrtint(lim\1*8+1)+1)\4, if(isprime(2*p-1),listput(v,2*p^2-p)); if(isprime(2*p+1), listput(v,2*p^2+p))); Vec(v) \\ Charles R Greathouse IV, Jun 13 2013

A010054(a(n))*A064911(a(n)) = 1. - Reinhard Zumkeller, Dec 03 2009

a(n) = A000217(A164977(n)). - Zak Seidov, Feb 16 2015

Edited by Robert G. Wilson v, Jul 08 2002

Definition corrected by Zak Seidov, Mar 09 2008

A179643 Products of exactly 2 distinct squares of primes and a different prime (p^2 * q^2 * r).

Original entry on oeis.org

180, 252, 300, 396, 450, 468, 588, 612, 684, 700, 828, 882, 980, 1044, 1100, 1116, 1300, 1332, 1452, 1476, 1548, 1575, 1692, 1700, 1900, 1908, 2028, 2124, 2156, 2178, 2196, 2205, 2300, 2412, 2420, 2450, 2475, 2548, 2556, 2628, 2844, 2900, 2925, 2988

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 21 2010

nonn

A050326(a(n)) = 5, subsequence of A225228. - Reinhard Zumkeller, May 03 2013

			180 = 2^2 * 3^2 * 5, 252 = 2^2 * 3^2 * 7, 300 = 2^2 * 3 * 5^2, ...

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, List of Prime Signatures
Index to sequences related to prime signature

Cf. A006881, A007304, A065036, A085986, A085987, A178739, A179642.

Cf. A179695.

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,2,2}; Select[Range[3000], f]

list(lim)=my(v=List(),t);forprime(p=2,sqrt(lim\12),forprime(q=p+1,sqrt(lim\p^2\2),t=(p*q)^2;forprime(r=2,lim\t,if(p==r||q==r,next);listput(v,t*r))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 19 2011

from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A179643(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x+sum((t:=primepi(s:=isqrt(y:=isqrt(x//r))))+(t*(t-1)>>1)-sum(primepi(y//k) for k in primerange(1, s+1)) for r in primerange(x+1))+sum(primepi(isqrt(x//p**3)) for p in primerange(integer_nthroot(x,3)[0]+1))-primepi(integer_nthroot(x,5)[0])
    return bisection(f,n,n) # Chai Wah Wu, Mar 27 2025

A337563 Number of pairwise coprime unordered triples of positive integers > 1 summing to n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 2, 0, 2, 1, 4, 0, 7, 1, 7, 3, 9, 2, 15, 3, 13, 5, 17, 4, 29, 5, 20, 8, 28, 8, 42, 8, 31, 14, 42, 10, 59, 12, 45, 21, 52, 14, 77, 17, 68, 26, 69, 19, 101, 26, 84, 34, 86, 25, 138, 28, 95, 43, 111, 36, 161, 35, 118, 52, 151

Offset: 0

Gus Wiseman, Sep 21 2020

nonn

Such partitions are necessarily strict.

The Heinz numbers of these partitions are the intersection of A005408 (no 1's), A014612 (triples), and A302696 (coprime).

			The a(10) = 1 through a(24) = 15 triples (empty columns indicated by dots, A..J = 10..19):
  532  .  543  .  743  753  754  .  765  B53  875  975  985  B75  987
          732     752       853     873       974  B73  B65  D73  B76
                            952     954       A73  D53  B74       B85
                            B32     972       B54       B83       B94
                                    B43       B72       B92       BA3
                                    B52       D43       D54       C75
                                    D32       D52       D72       D65
                                                        E53       D74
                                                        H32       D83
                                                                  D92
                                                                  F72
                                                                  G53
                                                                  H43
                                                                  H52
                                                                  J32

Fausto A. C. Cariboni, Table of n, a(n) for n = 0..10000

A055684 is the version for pairs.
A220377 allows 1's, with non-strict version A307719.
A337485 counts these partitions of any length.
A337563*6 is the ordered version.
A001399(n - 3) = A069905(n) = A211540(n + 2) counts 3-part partitions.
A002865 counts partitions with no 1's, with strict case A025147.
A007359 counts pairwise coprime partitions with no 1's.
A078374 counts relatively prime strict partitions.
A200976 and A328673 count pairwise non-coprime partitions.
A302696 ranks pairwise coprime partitions.
A302698 counts relatively prime partitions with no 1's.
A305713 counts pairwise coprime strict partitions.
A327516 counts pairwise coprime partitions.
A337452 counts relatively prime strict partitions with no 1's.
Cf. A007304, A082024, A101268, A284825, A332004, A337451, A337461, A337462, A337561, A337599, A337601, A337605.

Table[Length[Select[IntegerPartitions[n,{3}],!MemberQ[#,1]&&CoprimeQ@@#&]],{n,0,30}]

A156040 Number of compositions (ordered partitions) of n into 3 parts (some of which may be zero), where the first is at least as great as each of the others.

Original entry on oeis.org

1, 1, 3, 4, 6, 8, 11, 13, 17, 20, 24, 28, 33, 37, 43, 48, 54, 60, 67, 73, 81, 88, 96, 104, 113, 121, 131, 140, 150, 160, 171, 181, 193, 204, 216, 228, 241, 253, 267, 280, 294, 308, 323, 337, 353, 368, 384, 400, 417, 433, 451, 468, 486, 504, 523, 541, 561, 580, 600

Offset: 0

Jack W Grahl, Feb 02 2009, Feb 11 2009

For n = 1, 2 these are just the triangular numbers. a(n) is always at least 1/3 of the corresponding triangular number, since each partition of this type gives up to three ordered partitions with the same cyclical order.
An alternative definition, which avoids using parts of size 0: a(n) is the third diagonal of A184957. - N. J. A. Sloane, Feb 27 2011
Diagonal sums of the triangle formed by rows T(2, k) k = 0, 1, ..., 2m of ascending m-nomial triangles (see A004737):
1
1 2 1
1 2 3 2 1
1 2 3 4 3 2 1
1 2 3 4 5 4 3 2 1
1 2 3 4 5 6 5 4 3 2 1
- Bob Selcoe, Feb 07 2014
Arrange A004396 in rows successively shifted to the right two spaces and sum the columns:
1  1  2  3  3  4  5  5  6 ...
      1  1  2  3  3  4  5 ...
            1  1  2  3  3 ...
                  1  1  2 ...
                        1 ...
------------------------------
1  1  3  4  6  8 11 13 17 ... - L. Edson Jeffery, Jul 30 2014
a(n) is the dimension of three-dimensional (2n + 2)-homogeneous polynomial vector fields with full tetrahedral symmetry (for a given orthogonal representation), and which are solenoidal. - Giedrius Alkauskas, Sep 30 2017
Also the number of compositions of n + 3 into three parts, the first at least as great as each of the other two. Also the number of compositions of n + 4 into three parts, the first strictly greater than each of the other two. - Gus Wiseman, Oct 09 2020

			G.f. = 1 + x + 3*x^2 + 4*x^3 + 6*x^4 + 8*x^5 + 11*x^6 + 13*x^7 + 17*x^8 + 20*x^9 + ...
The a(4) = 6 compositions of 4 are: (4 0 0), (3 1 0), (3 0 1), (2 2 0), (2 1 1), (2 0 2).
From _Gus Wiseman_, Oct 05 2020: (Start)
The a(0) = 1 through a(7) = 13 triples of nonnegative integers summing to n where the first is at least as great as each of the other two are:
  (000)  (100)  (101)  (111)  (202)  (212)  (222)  (313)
                (110)  (201)  (211)  (221)  (303)  (322)
                (200)  (210)  (220)  (302)  (312)  (331)
                       (300)  (301)  (311)  (321)  (403)
                              (310)  (320)  (330)  (412)
                              (400)  (401)  (402)  (421)
                                     (410)  (411)  (430)
                                     (500)  (420)  (502)
                                            (501)  (511)
                                            (510)  (520)
                                            (600)  (601)
                                                   (610)
                                                   (700)
(End)

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Giedrius Alkauskas, Projective and polynomial superflows. I, arxiv.org/1601.06570 [math.AG] (2017), Section 4.3.
Index entries for linear recurrences with constant coefficients, signature (1,1,0,-1,-1,1).

For compositions into 4 summands see A156039; also see A156041 and A156042.
Cf. A184957, A071619 (bisection).
A001399(n-2)*2 is the strict case.
A001840(n-2) is the version with opposite relations.
A001840(n-1) is the version with strict opposite relations.
A069905 is the case with strict relations.
A014311 ranks 3-part compositions, with strict case A337453.
A014612 ranks 3-part partitions, with strict case A007304.
Cf. A000217, A001523, A211540, A218004, A220377, A337483, A337484.

a:= proc(n) local m, r; m := iquo(n, 6, 'r'); (4 +6*m +2*r) *m + [1, 1, 3, 4, 6, 8][r+1] end: seq(a(n), n=0..60); # Alois P. Heinz, Jun 14 2009

nn = 58; CoefficientList[Series[x^3/(1 - x^2)^2/(1 - x^3) + 1/(1 - x^2)^2/(1 - x), {x, 0, nn}], x] (* Geoffrey Critzer, Jul 14 2013 *)
CoefficientList[Series[(1 + x^2)/((1 + x) * (1 + x + x^2) * (1 - x)^3), {x, 0, 58}], x] (* L. Edson Jeffery, Jul 29 2014 *)
LinearRecurrence[{1, 1, 0, -1, -1, 1}, {1, 1, 3, 4, 6, 8}, 60] (* Harvey P. Dale, May 28 2015 *)
Table[Length[Select[Join@@Permutations/@IntegerPartitions[n+3,{3}],#[[1]]>=#[[2]]&&#[[1]]>=#[[3]]&]],{n,0,15}] (* Gus Wiseman, Oct 05 2020*)

{a(n) = n*(n+4)\6 + 1}; /* Michael Somos, Mar 26 2017 */

G.f.: (x^2+1) / (1-x-x^2+x^4+x^5-x^6). - Alois P. Heinz, Jun 14 2009
Slightly nicer g.f.: (1+x^2)/((1-x)*(1-x^2)*(1-x^3)). - N. J. A. Sloane, Apr 29 2011
a(n) = A007590(n+2) - A000212(n+2). - Richard R. Forberg, Dec 08 2013
a(2*n) = A071619(n+1). - L. Edson Jeffery, Jul 29 2014
a(n) = a(n-1) + a(n-2) - a(n-4) - a(n-5) + a(n-6), with a(0) = 1, a(1) = 1, a(2) = 3, a(3) = 4, a(4) = 6, a(5) = 8. - Harvey P. Dale, May 28 2015
a(n) = (n^2 + 4*n + 3)/6 + IF(MOD(n, 2) = 0, 1/2) + IF(MOD(n, 3) = 1, -1/3). - Heinrich Ludwig, Mar 21 2017
a(n) = 1 + floor((n^2 + 4*n)/6). - Giovanni Resta, Mar 21 2017
Euler transform of length 4 sequence [1, 2, 1, -1]. - Michael Somos, Mar 26 2017
a(n) = a(-4 - n) for all n in Z. - Michael Somos, Mar 26 2017
0 = a(n)*(-1 + a(n) - 2*a(n+1) - 2*a(n+2) + 2*a(n+3)) + a(n+1)*(+1 + a(n+1) + 2*a(n+2) - 2*a(n+3)) + a(n+2)*(+1 + a(n+2) - 2*a(n+3)) + a(n+3)*(-1 + a(n+3)) for all n in Z. - Michael Somos, Mar 26 2017
a(n) = round((n+1)*(n+3)/6). - Bill McEachen, Feb 16 2021
Sum_{n>=0} 1/a(n) = 3/2 + Pi^2/36 + (tan(c1)-1)*c1 + 3*c2*sinh(c2)/(1+2*cosh(c2)), where c1 = Pi/(2*sqrt(3)) and c2 = Pi*sqrt(2)/3. - Amiram Eldar, Dec 10 2022
E.g.f.: ((16 + 15*x + 3*x^2)*cosh(x) + 2*exp(-x/2)*(cos(sqrt(3)*x/2) - sqrt(3)*sin(sqrt(3)*x/2)) + (7 + 15*x + 3*x^2)*sinh(x))/18. - Stefano Spezia, Apr 05 2023

More terms from Alois P. Heinz, Jun 14 2009

A179644 Product of the 4th power of a prime and 2 different distinct primes (p^4qr).

Original entry on oeis.org

240, 336, 528, 560, 624, 810, 816, 880, 912, 1040, 1104, 1134, 1232, 1360, 1392, 1456, 1488, 1520, 1776, 1782, 1840, 1904, 1968, 2064, 2106, 2128, 2256, 2288, 2320, 2480, 2544, 2576, 2754, 2832, 2835, 2928, 2960, 2992, 3078, 3216, 3248, 3280, 3344, 3408

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 21 2010

nonn

			240=2^4*3*5,336=2^4*3*7,..810=2^3^4*5,..

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, List of Prime Signatures
Index to sequences related to prime signature

Cf. A006881, A007304, A065036, A085986, A085987, A178739, A179642, A179643.

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,1,4}; Select[Range[4000], f]
Take[Union[#[[1]]^4 #[[2]]#[[3]]&/@(Flatten[Permutations/@ Subsets[ Prime[ Range[ 20]],{3}],1])],50] (* Harvey P. Dale, Feb 07 2013 *)

list(lim)=my(v=List(),t);forprime(p=2,(lim\6)^(1/4),forprime(q=2,sqrt(lim\p^4),if(p==q,next);t=p^4*q;forprime(r=q+1,lim\t,if(p==r,next);listput(v,t*r))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 19 2011

from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A179644(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x+sum((t:=primepi(s:=isqrt(y:=x//r**4)))+(t*(t-1)>>1)-sum(primepi(y//k) for k in primerange(1, s+1)) for r in primerange(integer_nthroot(x,4)[0]+1))+sum(primepi(x//p**5) for p in primerange(integer_nthroot(x,5)[0]+1))-primepi(integer_nthroot(x,6)[0])
    return bisection(f,n,n) # Chai Wah Wu, Mar 27 2025

cp's OEIS Frontend

A211540 Number of ordered triples (w,x,y) with all terms in {1..n} and 2w = 3x + 4y.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Formula

A284825 Number of partitions of n into 3 parts without common divisors such that every pair of them has common divisors.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A140106 Number of noncongruent diagonals in a regular n-gon.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Python

SageMath

Formula

Extensions

A336568 Numbers that are not a product of two numbers each having distinct prime multiplicities.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A101271 Number of partitions of n into 3 distinct and relatively prime parts.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A068443 Triangular numbers which are the product of two primes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A179644 Product of the 4th power of a prime and 2 different distinct primes (p^4qr).