|
%I #20 Oct 13 2023 10:02:26
%S 16,104,5368,40880,2532160,20390552,1364393896,788854912240,
%T 6641649422408,4056611764783760,296528425830656800,
%U 2544627654221217656,188573151481968108424,121907205457107043376080
%N Fermat quotient of the n-th prime with base 3.
%H Takashi Agoh and Ladislav Skula, <a href="http://dx.doi.org/10.1016/j.jnt.2008.06.001">The fourth power of the Fermat quotient</a>, J. Numb. Theory 128 (2008) 2865-2873.
%F a(n) = (3^(p-1)-1)/p, where p=A000040(n).
%F a(n) = A046211(A000040(n)), for n >= 3. - _Amiram Eldar_, Oct 13 2023
%p A146211:= n-> map (p-> (3^(p-1)-1)/p, ithprime(n)):
%p seq (A146211(n), n=3..16); # _Jani Melik_, Jan 24 2010
%t Table[(3^(p - 1) - 1)/p, {p, Prime[Range[3, 16]]}] (* _Amiram Eldar_, Oct 13 2023 *)
%o (PARI) a(n) = my(p=prime(n)); (3^(p-1)-1)/p; \\ _Michel Marcus_, Oct 13 2023
%Y Subsequence of A046211. Cf. A007663.
%K easy,nonn
%O 3,1
%A _R. J. Mathar_, Oct 28 2008
| 