cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A048243 Number of non-unitary divisors of binomial(n, floor(n/2)).

This page as a plain text file.
%I A048243 #18 Jul 22 2024 15:22:34
%S A048243 0,0,0,0,0,2,0,0,4,10,0,8,8,16,8,16,0,16,0,16,32,64,0,32,80,128,224,
%T A048243 320,256,352,64,128,256,512,256,640,640,1024,512,1280,512,1024,512,
%U A048243 1024,2560,3328,1280,2432,2432,3584,3584,5120,2048,2816,2048,4096,8192
%N A048243 Number of non-unitary divisors of binomial(n, floor(n/2)).
%F A048243 a(n) = A048105(A001405(n)).
%F A048243 a(n) = A000005(A001405(n)) - A034444(A001405(n)). [corrected by _Amiram Eldar_, Jul 22 2024]
%e A048243 For n = 10, binomial(10,5) = 252 = 4*9*7 has 18 divisors, 8 are unitary and the residual 10 are non-unitary; thus a(10) = 18 - 8 = 10.
%t A048243 Table[Function[k, DivisorSum[k, 1 &, ! CoprimeQ[#, k/#] &]]@ Binomial[n, Ceiling[n/2]], {n, 57}] (* _Michael De Vlieger_, Jun 29 2017 *)
%t A048243 f[n_] := DivisorSigma[0, n] - 2^PrimeNu[n]; Table[f[Binomial[n, Floor[n/2]]], {n, 1, 60}] (* _Amiram Eldar_, Jul 22 2024 *)
%o A048243 (PARI) a(n) = apply(x -> numdiv(x) - 2^omega(x), binomial(n, n\2)); \\ _Amiram Eldar_, Jul 22 2024
%Y A048243 Cf. A000005, A001405, A034444, A048105.
%K A048243 nonn
%O A048243 1,6
%A A048243 _Labos Elemer_