login
Number of strict partitions of n having standard deviation sigma > 1.
4

%I #7 Mar 14 2014 13:26:30

%S 0,0,0,0,1,1,3,4,5,8,10,12,16,20,24,30,36,43,52,62,73,87,102,119,140,

%T 163,189,220,254,293,338,388,445,510,583,665,758,862,979,1111,1258,

%U 1423,1608,1814,2045,2302,2588,2907,3262,3656,4094,4580,5118,5715,6376

%N Number of strict partitions of n having standard deviation sigma > 1.

%C Regarding standard deviation, see Comments at A238616.

%F a(n) + A239141(n) = A000009(n) for n >=1.

%F G.f.: Product_{m>=1} (1+x^m) -1 +(x^5+x^4+x^3+2*x^2+x+1)*x / ((x-1)*(x^2+x+1)). - _Alois P. Heinz_, Mar 14 2014

%e The standard deviations of the strict partitions of 9 are 0., 3.5, 2.5, 1.5, 2.16025, 0.5, 1.63299, 0.816497, so that a(9) = 5.

%t z = 30; g[n_] := Select[IntegerPartitions[n], Max[Length /@ Split@#] == 1 &]; s[t_] := s[t] = Sqrt[Sum[(t[[k]] - Mean[t])^2, {k, 1, Length[t]}]/Length[t]]

%t Table[Count[g[n], p_ /; s[p] < 1], {n, z}] (* A239140 *)

%t Table[Count[g[n], p_ /; s[p] <= 1], {n, z}] (* A239141 *)

%t Table[Count[g[n], p_ /; s[p] == 1], {n, z}] (* periodic 01 *)

%t Table[Count[g[n], p_ /; s[p] > 1], {n, z}] (* A239142 *)

%t Table[Count[g[n], p_ /; s[p] >= 1], {n, z}] (* A239143 *)

%t t[n_] := t[n] = N[Table[s[g[n][[k]]], {k, 1, PartitionsQ[n]}]]

%t ListPlot[Sort[t[30]]] (*plot of st.dev's of strict partitions of 30*)

%t (* _Peter J. C. Moses_, Mar 03 2014 *)

%Y Cf. A239140, A239141, A239142, A000009, A238616.

%K nonn,easy

%O 1,7

%A _Clark Kimberling_, Mar 11 2014