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

A080469 Composite n such that binomial(3*n,n)==3^n (mod n).

Original entry on oeis.org

36, 57, 121, 132, 552, 8397, 7000713, 9692541, 36294723, 564033861
Offset: 1

Views

Author

Benoit Cloitre, Oct 15 2003

Keywords

Comments

If p is prime, binomial(3*p,p)==3^p (mod p)
No other terms below 10^9.
A subsequence of A109641. The terms a(n) with n=2, 6, 7, 8, 9, 10 are of the form 3^k*p where p is prime and k=1, 3, 2, 5, 6, 7, respectively. It is tempting to conjecture that there are (infinitely many?) more terms of that form. - M. F. Hasler, Nov 11 2015

Examples

			57 is a term because binomial(3*57, 57) = 12039059761216294940321619222324879408784636200 mod 57 = 27 == 3^57 mod 57.
		

Crossrefs

Programs

  • Mathematica
    Do[If[ !PrimeQ[n], k = Binomial[3*n, n]; m = 3^n; If[Mod[k, n] == Mod[m, n], Print[n]]], {n, 1, 70000}] (* Ryan Propper, Aug 12 2005 *)
  • PARI
    forcomposite(n=1,1e9, binomod(3*n,n,n)==Mod(3,n)^n && print1(n",")) \\ Cf. Alekseyev link. - M. F. Hasler, Nov 14 2015

Extensions

One more term a(6) from Ryan Propper, Aug 12 2005
Four new terms a(7)-a(10) added by Max Alekseyev, Nov 05 2009

A084699 Composite integers j such that binomial(2*j,j) == 2^j (mod j).

Original entry on oeis.org

12, 30, 56, 424, 992, 16256, 58288, 119984, 356992, 1194649, 9973504, 12327121, 13141696, 22891184, 67100672, 233850649
Offset: 1

Views

Author

Benoit Cloitre, Oct 15 2003

Keywords

Comments

If p is prime, binomial(2*p,p) == 2^p (mod p).
a(17) > 10^9.
From Gabriel Guedes and Ricardo Machado, Nov 16 2023: (Start)
Theorem. Let j = (2^k)*p, where p is an odd prime and k is in N; then binomial(2*j,j) == 2^j (mod j) if and only if p satisfies the following conditions:
a) p divides binomial(2^(k+1),2^k) - 2^(2^k);
b) p has at least k 1's in its binary expansion.
Theorem. If m is an even perfect number then j = 2m satisfies the congruence binomial(2*j,j) == 2^j (mod j). See A000396.
Theorem. Let j = p^2 with p a prime number. Then p is a Wieferich prime if and only if binomial(2*j,j) == 2^j (mod j). See A001220. (End)
Contains 17179738112 and 274877382656 (from Guedes-Machado paper). - Michael De Vlieger, Nov 22 2023
Contains 3386741824, 750984028672, 33029195197184, 1145067923695616, 422612863956511744. - Ricardo Machado, Nov 23 2023
Contains 84385517065596416, 62648180117928433664, 273984397779878971648, 36506097537257040703232. - Max Alekseyev, Dec 07 2023

Crossrefs

Contains A139256 as a subsequence.

Programs

  • PARI
    lista(nn) = {forcomposite(n=1, nn, if (binomod(2*n, n, n) == Mod(2, n)^n, print1(n, ", ")));} \\ Michel Marcus, Dec 06 2013 and Dec 03 2023

Extensions

More terms from David Wasserman, Jan 03 2005
a(11)-a(16) from Max Alekseyev, Aug 05 2011
Showing 1-2 of 2 results.