%I #14 Aug 21 2013 17:59:10
%S 61,8205,3440347021,7080447489597,171336855102372210685,
%T 1747517658865390518778893,610345691966794096778276272763149,
%U 49983985045539556672075839852554462798428935229
%N a(n)=(3/2)*(1/p)*(2*p+1)*(3^p+1)*B(2*p) where p=prime(n) and where B(k) denotes the k-th Bernoulli number.
%C a(n)==1 mod (prime(n))
%t Table[p=Prime[n];3/(2p) (2p+1)(3^p+1)BernoulliB[2p],{n,3,10}] (* _Harvey P. Dale_, Aug 21 2013 *)
%o (PARI) a(n)=3/2/prime(n)*(2*prime(n)+1)*(3^prime(n)+1)*bernfrac(2*prime(n))
%K nonn
%O 3,1
%A _Benoit Cloitre_, Feb 11 2004