login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A228562 Composite numbers k that are not prime powers such that binomial(2k-1, k-1) is congruent to 1 (mod k). 11
27173, 2001341, 16024189487 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
LINKS
Richard J. McIntosh, On the converse of Wolstenholme's Theorem, Acta Arithmetica, 71 (1995), 381-389.
FORMULA
A099905(a(n)) = 1. - Jonathan Sondow, Jan 24 2016
MATHEMATICA
Select[Range[30000], PrimeNu[#] > 1 && Mod[Binomial[2# - 1, # - 1], #] == 1 &] (* Alonso del Arte, May 11 2014 *)
PROG
(PARI) N=10^9; for(n=2, N, if(Mod(binomial(2*n-1, n-1), n)==1 && !ispower(n) && !isprime(n), print1(n, ", "))); \\ Felix Fröhlich, May 11 2014
(PARI) vp(n, p)=my(s); while(n\=p, s+=n); s
is(n)=my(f=factor(n)[, 1], G); if(#f==1, return(0)); for(i=1, #f, if(vp(2*n-1, f[i]) > vp(n, f[i])+vp(n-1, f[i]), return(0))); G=prod(i=1, #f, f[i]^(log(n)\log(f[i]))); prod(i=n+1, 2*n-1, i/gcd(i, G), Mod(1, n))/prod(i=2, n-1, i/gcd(i, G), Mod(1, n))==1
forcomposite(n=4, 1e9, if(is(n), print1(n", "))) \\ Charles R Greathouse IV, May 12 2014
CROSSREFS
The odd terms of A328497.
Sequence in context: A186481 A210073 A236275 * A224466 A127411 A268358
KEYWORD
hard,nonn,more,bref
AUTHOR
Felix Fröhlich, Aug 25 2013
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 04:13 EDT 2024. Contains 371235 sequences. (Running on oeis4.)