login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A202425 Number of partitions of n into parts having pairwise common factors but no overall common factor. 7

%I

%S 1,0,0,0,0,0,2,0,0,0,3,0,3,0,0,1,6,0,5,0,2,2,9,0,8,2,4,3,16,0,22,5,6,

%T 5,19,2,35,8,14,6,44,4,55,13,16,19,64,6,82,17,39,31,108,10,105,40,66,

%U 46,161,14,182,61,97,72,207,37,287,85,144,93,357,59

%N Number of partitions of n into parts having pairwise common factors but no overall common factor.

%H Alois P. Heinz, <a href="/A202425/b202425.txt">Table of n, a(n) for n = 31..251</a>

%F a(n > 0) = A328672(n) - 1. - _Gus Wiseman_, Nov 04 2019

%e a(31) = 1: [6,10,15] = [2*3,2*5,3*5].

%e a(37) = 2: [6,6,10,15], [10,12,15].

%e a(41) = 3: [6,10,10,15], [6,15,20], [6,14,21].

%e a(47) = 6: [6,6,10,10,15], [10,10,12,15], [6,6,15,20], [12,15,20], [6,6,14,21], [12,14,21].

%e a(49) = 5: [6,6,6,6,10,15], [6,6,10,12,15], [10,12,12,15], [6,10,15,18], [10,15,24].

%p with(numtheory):

%p w:= (m, h)-> mul(`if`(j>=h, 1, j), j=factorset(m)):

%p b:= proc(n, i, g, s) option remember; local j, ok, si;

%p if n<0 then 0

%p elif n=0 then `if`(g>1, 0, 1)

%p elif i<2 or member(1, s) then 0

%p else ok:= evalb(i<=n);

%p si:= map(x->w(x, i), s);

%p for j in s while ok do ok:= igcd(i, j)>1 od;

%p b(n, i-1, g, si) +`if`(ok, add(b(n-t*i, i-1, igcd(i, g),

%p si union {w(i,i)} ), t=1..iquo(n, i)), 0)

%p fi

%p end:

%p a:= n-> b(n, n, 0, {}):

%p seq(a(n), n=31..100);

%t w[m_, h_] := Product[If[j >= h, 1, j], {j, FactorInteger[m][[All, 1]]}]; b[n_, i_, g_, s_] := b[n, i, g, s] = Module[{j, ok, si}, Which[n<0, 0, n == 0, If[g>1, 0, 1], i<2 || MemberQ[s, 1], 0, True, ok = (i <= n); si = w[#, i]& /@ s; Do[If[ok, ok = (GCD[i, j]>1)], {j, s}]; b[n, i-1, g, si] + If[ok, Sum[b[n-t*i, i-1, GCD[i, g], si ~Union~ {w[i, i]}], {t, 1, Quotient[n, i]}], 0]]]; a[n_] := b[n, n, 0, {}]; Table[a[n], {n, 31, 100}] (* _Jean-Fran├žois Alcover_, Feb 16 2017, translated from Maple *)

%t Table[Length[Select[IntegerPartitions[n],GCD@@#==1&&And@@(GCD[##]>1&)@@@Tuples[#,2]&]],{n,0,40}] (* _Gus Wiseman_, Nov 04 2019 *)

%Y Cf. A202385, A018783.

%Y The version with only distinct parts compared is A328672.

%Y The Heinz numbers of these partitions are A328868.

%Y The strict case is A202385, which is essentially the same as A318715.

%Y The version for non-isomorphic multiset partitions is A319759.

%Y The version for set-systems is A326364.

%Y Intersecting partitions are A200976.

%Y Cf. A000837, A305148, A305843, A306006, A318717, A328673, A328867.

%K nonn

%O 31,7

%A _Alois P. Heinz_, Dec 19 2011

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified February 21 20:39 EST 2020. Contains 332111 sequences. (Running on oeis4.)