A099908 C(2n-1,n-1) mod n^4.

0, 3, 10, 35, 126, 462, 1716, 2339, 4627, 2378, 1332, 4238, 2198, 5148, 1260, 57635, 14740, 85026, 61732, 64410, 100509, 163716, 158172, 171918, 93876, 309780, 148969, 444220, 268280, 370712, 29792, 532771, 652200, 938386, 816466, 907874

Offset: 1

Author

Henry Bottomley, Oct 29 2004

Keywords

• nonn

Comments

a(16843)=a(2124679)=1 meaning that 16843 and 2124679 are Wolstenholme primes A088164.

Examples

			a(11) =352716 mod 1461 =1332.
		

Links

• Chai Wah Wu, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A088164, A088218, A099905, A099906, A099907.

Programs

Mathematica

Table[Mod[Binomial[2n-1,n-1],n^4],{n,40}] (* Harvey P. Dale, Dec 12 2021 *)

Python

from _future_ import division
A099908_list, b = [], 1
for n in range(1,10001):
    A099908_list.append(b % n**4)
    b = b*2*(2*n+1)//(n+1) # Chai Wah Wu, Jan 26 2016