OFFSET
1,1
MAPLE
select(isprime, [4&^n mod 3^n$n=1..200])[]; # Alois P. Heinz, May 18 2019
MATHEMATICA
Select[Table[PowerMod[4, n, 3^n], {n, 100}], PrimeQ] (* Alonso del Arte, Jan 03 2011 *)
PROG
(PARI) terms(n) = my(i=0); for(k=0, oo, if(i>=n, break); my(x=lift(Mod(4, 3^k)^k)); if(ispseudoprime(x), print1(x, ", "); i++))
/* Print initial 7 terms as follows: */
terms(7) \\ Felix Fröhlich, May 18 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Jan 03 2011
STATUS
approved