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.

A193990 Number of distinct prime factors <= n of binomial(2*n,n).

This page as a plain text file.
%I A193990 #11 Oct 31 2015 05:30:19
%S A193990 0,1,1,1,2,2,2,3,2,1,3,2,3,3,3,3,4,5,5,6,5,4,4,4,4,4,3,5,7,5,5,6,7,6,
%T A193990 6,5,5,6,5,5,6,6,7,7,7,7,8,9,8,8,7,6,8,8,6,6,5,6,7,8,9,9,9,9,9,9,11,
%U A193990 10,9,8,9,10,11,11,11,10,11,11,11,12,12,11
%N A193990 Number of distinct prime factors <= n of binomial(2*n,n).
%H A193990 T. D. Noe, <a href="/A193990/b193990.txt">Table of n, a(n) for n = 1..10000</a>
%p A193990 a := n -> nops(numtheory:-factorset(binomial(2*n,n)) minus select(isprime, {$n+1..2*n})): seq(a(n), n=1..82); # _Peter Luschny_, Oct 31 2015
%t A193990 Table[Length[Select[Transpose[FactorInteger[Binomial[2 n, n]]][[1]], # <= n &]], {n, 100}]
%o A193990 (PARI) valp(n, p)=my(s); while(n\=p, s+=n); s
%o A193990 a(n)=my(s); forprime(p=2, n, if(valp(2*n, p)>2*valp(n, p), s++)); s \\ _Charles R Greathouse IV_, May 25 2013
%Y A193990 Cf. A067434 (number of prime factors in binomial(2*n,n)), A108954.
%K A193990 nonn
%O A193990 1,5
%A A193990 _T. D. Noe_, Aug 17 2011