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.

A278961 Triangle read by rows: row n consists of k, 1<=k<=n, such that binomial(n,k) is divisible by gcd(n,k).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 3, 1, 2, 3, 4, 1, 5, 1, 2, 3, 4, 5, 6, 1, 2, 3, 5, 6, 7, 1, 2, 3, 4, 5, 6, 7, 8, 1, 3, 4, 6, 7, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 5, 6, 7, 10, 11, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 3, 5, 9, 11, 13, 1, 2, 4, 7, 8, 11, 13, 14, 1, 2
Offset: 1

Views

Author

Robert Israel, Dec 02 2016

Keywords

Comments

All k coprime to n are always included, in particular 1 and n-1.
Row n contains k if and only if it contains n-k.
A014847 consists of k such that row 2k contains k.
If n is prime or the square of a prime, row n contains all numbers from 1 to n-1. This is not true for higher powers: row p^r does not contain any multiples of p^(r-1) if r > 2.
Prime p is in row n>p if and only if the p-adic order of n is not 1.

Examples

			Row 8 contains 2 because gcd(8,2)=2 divides binomial(8,2) = 28, but not 4 because gcd(8,4)=4 does not divide binomial(8,4)= 70.
		

Crossrefs

Cf. A014847.

Programs

  • Maple
    f:= proc(n,m) if binomial(n,m) mod igcd(n,m) = 0 then m else NULL fi end proc:
    seq(seq(f(n,m),m=1..n),n=1..40);
  • Mathematica
    Table[If[Divisible[Binomial[n,k],GCD[n,k]],k,Nothing],{n,20},{k,n}]//Flatten (* Harvey P. Dale, Dec 04 2022 *)