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.

Showing 1-1 of 1 results.

A056676 Number of non-unitary but squarefree divisors of binomial(n,floor(n/2)). Also number of nonsquarefree but unitary divisors of binomial(n,floor(n/2)).

Original entry on oeis.org

0, 0, 0, 0, 0, 2, 0, 0, 4, 6, 0, 8, 8, 8, 8, 16, 0, 16, 0, 16, 32, 32, 0, 32, 48, 48, 56, 56, 96, 96, 64, 128, 128, 192, 256, 384, 384, 384, 512, 768, 512, 512, 512, 512, 448, 448, 768, 896, 896, 896, 896, 896, 768, 768, 2048, 2048, 4096, 4096, 2048, 2048, 2048, 2048
Offset: 1

Views

Author

Labos Elemer, Aug 10 2000

Keywords

Examples

			For n = 14, binomial(14,7) = 3432 has 32 divisors, 16 unitary, 16 squarefree. The size of overlap is 8. The complementary parts are: non-unitary/squarefree set ={2,6,22,26,66,78,286,828}, while the unitary/not squarefree set of equal size is {8,24,88,104,264,312,1144,3432}. So a(14) = 8.
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := If[e == 1, 2, 1]; a[1] = 0; a[n_] := 2^Length[fct = FactorInteger[Binomial[n, Floor[n/2]]]] - Times @@ f @@@ fct; Array[a, 100] (* Amiram Eldar, Oct 04 2024 *)
  • PARI
    a(n) = {my(f = factor(binomial(n, n\2)), e = f[, 2]); 2^omega(f) - prod(i = 1, #e, if(e[i] == 1, 2, 1)); } \\ Amiram Eldar, Oct 04 2024

Formula

a(n) = A039593(n) - A000005(A055231(x)) = A039593(n) - A000005(A007913(x)/A055229(x)), where x = A001405(n) = binomial(n, floor(n/2)).
a(n) = A039593(n) - A056673(n). - Sean A. Irvine, May 02 2022
a(n) = A056674(A001405(n)). - Amiram Eldar, Oct 04 2024
Showing 1-1 of 1 results.