OFFSET
1,2
COMMENTS
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
EXAMPLE
MAPLE
a:= proc(n) option remember; local k;
for k from `if`(n=1, 1, a(n-1)+1)
while (2&^k mod k <> 4&^k mod k) or
(irem(k, 2, 'r')=0 and (2&^r mod r = 4&^r mod r))
do od; k
end:
seq(a(n), n=1..50); # Alois P. Heinz, May 16 2013
MATHEMATICA
Prepend[Select[Range[2, 2100], (PowerMod[2, #, #] == PowerMod[4, #, #]) && ! (PowerMod[2, #/2, #/2] == PowerMod[4, #/2, #/2]) &], 1]
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Oct 23 2006
STATUS
approved
