OFFSET
1,1
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..1200
EXAMPLE
d(10) = d(1) + d(9) = d(3) + d(7) = d(5) + d(5) = 4 and 10 is the least number with 3 partitions of two numbers with this property: therefore a(3) = 10;
d(126) = d(21) + d(105) = d(22) + d(104) = d(28) + d(98) = d(38) + d(33) = d(40) + d(86) = d(50) + d(76) = d(63) + d(63) = 12 and 126 is the least number with 7 partitions of two numbers with this property: therefore a(7) = 126.
MAPLE
with(numtheory): P:=proc(q) local a, h, k, n; for h from 1 to q do for n from 2*h to q do
a:=0; for k from 1 to trunc(n/2) do if tau(n)=tau(k)+tau(n-k) then a:=a+1; fi; od;
if a=h then print(n); break; fi; od; od; end: P(10^6);
MATHEMATICA
nn = 10^3; Table[SelectFirst[Range@ nn, Function[k, With[{e = DivisorSigma[0, k]}, Count[Transpose@ {Range[k - 1, Ceiling[k/2], -1], Range@ Floor[k/2]}, x_ /; Total@ DivisorSigma[0, x] == e] == n]]], {n, 54}] (* Michael De Vlieger, Apr 06 2016 *)
PROG
(PARI) isok(k, n) = {my(nb = 0, tau = numdiv(k)); for (j=1, k\2, if (numdiv(j)+numdiv(k-j) == tau, nb++); if (nb > n, return (0)); ); nb == n; }
a(n) = {k=2; while (!isok(k, n), k++); k; } \\ Michel Marcus, Apr 07 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paolo P. Lava, Apr 06 2016
STATUS
approved