login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A271382 Least k with precisely n partitions k = x + y satisfying d(k) = d(x) + d(y), where d(k) is the number of divisors of k. 2

%I #13 Apr 09 2016 16:26:43

%S 2,14,10,26,44,45,126,68,99,104,162,117,98,124,232,164,148,200,260,

%T 333,231,244,248,297,273,284,315,406,332,345,385,430,344,399,388,436,

%U 429,488,465,495,472,525,561,555,621,556,632,604,652,712,536,693,735,675

%N Least k with precisely n partitions k = x + y satisfying d(k) = d(x) + d(y), where d(k) is the number of divisors of k.

%H Paolo P. Lava, <a href="/A271382/b271382.txt">Table of n, a(n) for n = 1..1200</a>

%e 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;

%e 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.

%p 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

%p 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;

%p if a=h then print(n); break; fi; od; od; end: P(10^6);

%t 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 *)

%o (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;}

%o a(n) = {k=2; while (!isok(k, n), k++); k;} \\ _Michel Marcus_, Apr 07 2016

%Y Cf. A000005, A211224, A271384.

%K nonn,easy

%O 1,1

%A _Paolo P. Lava_, Apr 06 2016

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 10:22 EDT 2024. Contains 371268 sequences. (Running on oeis4.)