|
%I #14 Jul 08 2024 09:00:07
%S 1,2,3,11,7,9,10,11,14,19,12,21,28,19,24,66,54,21,20,58,24,25,32,70,
%T 44,28,39,55,91,97,44,66,35,94,39,80,86,47,129,66,45,75,100,58,129,75,
%U 56,70,68,74,178,62,66,340,58,75,209,97,93,124,115,101,138,66,84,75,111,94
%N a(n) is the smallest k>=n such that the number of partitions of k is a multiple of n, or -1 if no such k exists.
%C I do not know if a(n) exists for all n. First term which is currently unknown is a(2219) (a(2219) > 11600, while a(2218) = 2602).
%e a(4) is 11 because 11 is the smallest number for which P(11) is divisible by 4, where P() is the partition function.
%p for i from 2 while i < 30000 do for j from i while j < 1000000000 do c := numbpart(j); if (c mod i = 0) then print(i,j); break; end if; end do; end do;
%p # alternative
%p A086146 := proc(n)
%p local k ;
%p for k from n do
%p if combinat[numbpart](k) mod n =0 then
%p return k;
%p end if;
%p end do:
%p end proc: # _R. J. Mathar_, Apr 22 2013
%t kmax = 10^9;
%t a[n_] := Module[{k}, For[k = n, k <= kmax, k++, If[Divisible[ PartitionsP[k], n], Return[k]]]] /. Null -> -1;
%t Table[a[n], {n, 1, 68}] (* _Jean-François Alcover_, Jul 08 2024 *)
%Y Cf. A046641.
%K nonn
%O 1,2
%A _Jim Nastos_, Aug 25 2003
%E More terms from _Ray Chandler_, Sep 16 2003
|
