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

%I #13 Dec 04 2022 12:39:59

%S 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,

%T 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,

%U 11,12,1,3,5,9,11,13,1,2,4,7,8,11,13,14,1,2

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

%C All k coprime to n are always included, in particular 1 and n-1.

%C Row n contains k if and only if it contains n-k.

%C A014847 consists of k such that row 2k contains k.

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

%C Prime p is in row n>p if and only if the p-adic order of n is not 1.

%H Robert Israel, <a href="/A278961/b278961.txt">Table of n, a(n) for n = 1..10029</a>(rows 1 to 155 flattened)

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

%p f:= proc(n,m) if binomial(n,m) mod igcd(n,m) = 0 then m else NULL fi end proc:

%p seq(seq(f(n,m),m=1..n),n=1..40);

%t Table[If[Divisible[Binomial[n,k],GCD[n,k]],k,Nothing],{n,20},{k,n}]//Flatten (* _Harvey P. Dale_, Dec 04 2022 *)

%Y Cf. A014847.

%K nonn,tabf

%O 1,4

%A _Robert Israel_, Dec 02 2016