login
Numbers k such that k^2 - 1 divides binomial(2k,k).
6

%I #43 Mar 04 2024 11:55:14

%S 2,16,21,29,43,46,67,78,89,92,105,111,141,154,157,171,188,191,205,210,

%T 211,221,229,232,239,241,267,277,300,309,313,316,323,326,331,346,369,

%U 379,415,421,430,436,441,443,451,460,461,465,469,477,484,494,497,528

%N Numbers k such that k^2 - 1 divides binomial(2k,k).

%C Is a(n) asymptotic to c*n with 9 < c < 10?

%C A subset of A004782: numbers k such that 2(2k-3)!/(k!(k-1)!) is an integer.

%C Equivalently, numbers k such that k-1 divides A000108(k), the k-th Catalan number. - _M. F. Hasler_, Nov 11 2015

%C The data does not appear to support the conjectured asymptote statement (neither the constant nor being linear). - _Bill McEachen_, Feb 26 2024

%H Charles R Greathouse IV, <a href="/A081767/b081767.txt">Table of n, a(n) for n = 1..10000</a>

%H Max Alekseyev, <a href="http://home.gwu.edu/~maxal/gpscripts/">PARI/GP Scripts for Miscellaneous Math Problems</a>, sect. III: Binomial coefficients modulo integers, binomod.gp (v.1.4, 11/2015).

%t Select[Range[2,600],Divisible[Binomial[2#,#],#^2-1]&] (* _Harvey P. Dale_, May 11 2013 *)

%o (PARI) for(n=2, 999, binomial(2*n, n)%(n^2-1)||print1(n", ")) \\ _M. F. Hasler_, Nov 11 2015

%o (PARI) is_A081767(n)=!binomod(2*n, n, n^2-1) \\ Using binomod.gp by _Max Alekseyev_, cf. links. - _M. F. Hasler_, Nov 11 2015

%Y Subsequence of A094575 and of A004782.

%Y Cf. A000108.

%K nonn

%O 1,1

%A _Benoit Cloitre_, Apr 09 2003