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-5 of 5 results.

A030979 Numbers k such that binomial(2k,k) is not divisible by 3, 5 or 7.

Original entry on oeis.org

0, 1, 10, 756, 757, 3160, 3186, 3187, 3250, 7560, 7561, 7651, 20007, 59548377, 59548401, 45773612811, 45775397187, 237617431723407, 24991943420078301, 24991943420078302, 24991943420078307, 24991943715007536, 24991943715007537
Offset: 1

Views

Author

Shawn Godin (sgodin(AT)onlink.net)

Keywords

Comments

By Lucas's theorem, binomial(2k,k) is not divisible by a prime p iff all base-p digits of k are smaller than p/2.
Ronald L. Graham offered $1000 to the first person who could settle the question of whether this sequence is finite or infinite. He remarked that heuristic arguments show that it should be infinite, but finite if it is required that binomial(2k,k) is prime to 3, 5, 7 and 11, with k = 3160 probably the last k which has this property.
The Erdős et al. paper shows that for any two odd primes p and q there are an infinite number of k for which gcd(p*q,binomial(2k,k))=1; i.e., p and q do not divide binomial(2k,k). The paper does not deal with the case of three primes. - T. D. Noe, Apr 18 2007
Pomerance gives a heuristic suggesting that there are on the order of x^0.02595... terms up to x. - Charles R Greathouse IV, Oct 09 2015

References

  • R. K. Guy, Unsolved Problems in Number Theory, Springer, 1st edition, 1981. See section B33.

Crossrefs

Programs

  • Mathematica
    lim=10000; Intersection[Table[FromDigits[IntegerDigits[k,2],3], {k,0,lim}], Table[FromDigits[IntegerDigits[k,3],5], {k,0,lim}], Table[FromDigits[IntegerDigits[k,4],7], {k,0,lim}]] (* T. D. Noe, Apr 18 2007 *)
  • PARI
    fval(n,p)=my(s);while(n\=p,s+=n);s
    is(n)=fval(2*n,3)==2*fval(n,3) && fval(2*n,5)==2*fval(n,5) && fval(2*n,7)==2*fval(n,7) \\ Charles R Greathouse IV, Oct 09 2015

Formula

Intersection of A005836, A037453 and A037461. - T. D. Noe, Apr 18 2007

Extensions

More terms from Naohiro Nomoto, May 06 2002
Additional comments from R. L. Graham, Apr 25 2007
Additional comments and terms up 3^41 in b-file from Max Alekseyev, Nov 23 2008
Additional terms up to 10^70 in b-file from Christopher E. Thompson, Nov 06 2015

A196747 Numbers n such that 3 does not divide swing(n) = A056040(n).

Original entry on oeis.org

0, 1, 2, 6, 7, 8, 18, 19, 20, 24, 25, 26, 54, 55, 56, 60, 61, 62, 72, 73, 74, 78, 79, 80, 162, 163, 164, 168, 169, 170, 180, 181, 182, 186, 187, 188, 216, 217, 218, 222, 223, 224, 234, 235, 236, 240, 241, 242, 486, 487, 488, 492, 493, 494, 504, 505, 506, 510
Offset: 1

Views

Author

Peter Luschny, Oct 06 2011

Keywords

Crossrefs

Programs

  • Maple
    SwingExp := proc(m,n) local p, q; p := m;
    do q := iquo(n,p);
       if (q mod 2) = 1 then RETURN(1) fi;
       if q = 0 then RETURN(0) fi;
       p := p * m;
    od end:
    Search := proc(n,L) local m, i, r; m := n;
    for i in L do r := SwingExp(i,m);
       if r <> 0 then RETURN(NULL) fi
    od; n end:
    A196747_list := n -> Search(n,[3]):  # n is a search limit
  • Mathematica
    (* A naive solution *) sf[n_] := n!/Quotient[n, 2]!^2; Select[Range[0, 600], ! Divisible[sf[#], 3] &] (* Jean-François Alcover, Jun 28 2013 *)
  • PARI
    valp(n,p)=my(s); while(n\=p, s+=n); s
    is(n)=my(t=valp(n,3)); t%2==0 && 2*valp(n\2,3)==t \\ Charles R Greathouse IV, Feb 02 2016

A196748 Numbers n such that 3 and 5 do not divide swing(n) = A056040(n).

Original entry on oeis.org

0, 1, 2, 20, 24, 54, 60, 61, 62, 72, 73, 74, 504, 510, 511, 512, 560, 564, 1512, 1513, 1514, 1520, 1620, 1621, 1622, 6320, 6324, 6372, 6373, 6374, 6500, 6504, 6552, 6553, 6554, 6560, 13122, 13123, 13124, 13770, 13771, 13772, 13824, 13850, 15072, 15073, 15074
Offset: 1

Views

Author

Peter Luschny, Oct 06 2011

Keywords

Crossrefs

Programs

  • Maple
    # The function Search is defined in A196747.
    A196748_list := n -> Search(n,[3,5]):  # n is a search limit
  • Mathematica
    (* A naive solution *) sf[n_] := n!/Quotient[n, 2]!^2; Select[Range[0, 16000], !Divisible[sf[#], 3] && !Divisible[sf[#], 5] &] (* Jean-François Alcover, Jun 28 2013 *)
  • PARI
    valp(n,p)=my(s); while(n\=p, s+=n); s
    is(n)=valp(n,3)==2*valp(n\2,3) && valp(n,5)==2*valp(n\2,5) \\ Charles R Greathouse IV, Feb 02 2016

A196749 Numbers n such that 3, 5 and 7 do not divide swing(n) = A056040(n).

Original entry on oeis.org

0, 1, 2, 20, 1512, 1513, 1514, 6320, 6372, 6373, 6374, 6500, 15120, 15121, 15122, 15302, 40014, 119096754, 119096802, 91547225622, 91550794374
Offset: 1

Views

Author

Peter Luschny, Oct 06 2011

Keywords

Crossrefs

Programs

  • Maple
    # The function Search is defined in A196747.
    A196749_list := n -> Search(n,[3,5,7]):  # n is a search limit
  • PARI
    valp(n,p)=my(s); while(n\=p, s+=n); s
    is(n)=valp(n,3)==2*valp(n\2,3) && valp(n,5)==2*valp(n\2,5) && valp(n,7)==2*valp(n\2,7) \\ Charles R Greathouse IV, Feb 02 2016

A196750 Numbers n such that 3, 5, 7 and 11 do not divide swing(n) = A056040(n).

Original entry on oeis.org

0, 1, 2, 6320
Offset: 1

Views

Author

Peter Luschny, Oct 06 2011

Keywords

Comments

It is conjectured that there are no other terms.

Crossrefs

Programs

  • Maple
    # The function Search is defined in A196747.
    A196750_list := n -> Search(n,[3,5,7,11]):  # n is a search limit
  • PARI
    valp(n,p)=my(s); while(n\=p, s+=n); s
    is(n)=valp(n,3)==2*valp(n\2,3) && valp(n,5)==2*valp(n\2,5) && valp(n,7)==2*valp(n\2,7) && valp(n,11)==2*valp(n\2,11) \\ Charles R Greathouse IV, Feb 02 2016
Showing 1-5 of 5 results.