cp's OEIS Frontend

A323891 a(n) is the number of partitions of 72*n + 42 into 10 odd squares.

%I A323891 #21 Sep 19 2022 06:18:42
%S A323891 2,9,22,41,68,106,154,212,285,368,477,598,741,898,1076,1286,1524,1785,
%T A323891 2068,2379,2741,3131,3554,4002,4497,5044,5644,6274,6939,7653,8445,
%U A323891 9295,10186,11117,12113,13192,14355,15556,16807,18147,19570,21089,22673,24300,26029,27865,29821,31822,33894,36088
%N A323891 a(n) is the number of partitions of 72*n + 42 into 10 odd squares.
%D A323891 Laurențiu Panaitopol, Alexandru Gica, Arithmetic problems and number theory, Ed. Gil, Zalău, (2006), ch. 14, p. 85, pr. 32. (in Romanian).
%H A323891 Robert Israel, <a href="/A323891/b323891.txt">Table of n, a(n) for n = 0..2000</a>
%e A323891 For n=0, 72*0+42 = 42 = 25+9+1+1+1+1+1+1+1+1 = 9+9+9+9+1+1+1+1+1+1, so a(0)=2.
%e A323891 For n=1, 72*1+42 = 114 = 81+25+1+1+1+1+1+1+1+1 = 81+9+9+9+1+1+1+1+1+1 = 49+49+9+1+1+1+1+1+1+1 = 49+25+25+9+1+1+1+1+1+1 = 49+25+9+9+9+9+1+1+1+1 = 49+9+9+9+9+9+9+9+1+1 = 25+25+25+25+9+1+1+1+1+1 = 25+25+25+9+9+9+9+1+1+1 = 25+25+9+9+9+9+9+9+9+1, so a(1)=9.
%p A323891 S:= proc(n, k, m)
%p A323891    option remember;
%p A323891    local p,j;
%p A323891    if k = 0 then if n = 0 then return 1 else return 0 fi
%p A323891    elif m < 1 then return 0
%p A323891    elif n < k then return 0
%p A323891    elif n > k*m^2 then return 0
%p A323891    fi;
%p A323891    if m^2 > n then
%p A323891      p:= floor(sqrt(n));
%p A323891      if p::even then p:= p-1 fi;
%p A323891      return procname(n, k, p)
%p A323891    fi;
%p A323891    add(procname(n-j*m^2,k-j,m-2), j=0..n/m^2)
%p A323891 end proc:
%p A323891 seq(S(72*n+42, 10, 72*n+42), n=0..100); # _Robert Israel_, Feb 24 2019
%t A323891 a[n_] := IntegerPartitions[72n+42, {10}, Select[ Range[1, 72n+42, 2], IntegerQ@Sqrt@#&]] // Length;
%t A323891 Table[a[n], {n, 0, 100}] (* _Jean-François Alcover_, Sep 19 2022 *)
%o A323891 (Magma) [#RestrictedPartitions(72*n+42, 10, {(2*d+1)^2:d in [0..100]}): n in [0..100]];
%Y A323891 Cf. A000041, A001156, A016754, A025425, A025434, A033461, A035294, A078406, A090677, A167661, A167700.
%K A323891 nonn
%O A323891 0,1
%A A323891 _Marius A. Burtea_, Feb 12 2019