login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A350087
a(n) = A000032(n)^A000032(n+1) mod A000032(n+2).
1
2, 1, 4, 5, 13, 4, 16, 33, 64, 137, 123, 5, 733, 1241, 786, 235, 4331, 721, 13156, 18253, 32935, 37347, 24281, 34270, 127055, 370341, 630347, 52210, 76027, 255327, 2279978, 1527009, 10107310, 13246385, 17943817, 1128698, 75596554, 135885089, 46740145, 143055674, 76508345, 850178467, 978803396
OFFSET
0,1
LINKS
EXAMPLE
a(3) = 4^7 mod 11 = 5.
MAPLE
luc:= n -> combinat:-fibonacci(n-1)+combinat:-fibonacci(n+1):
f:= n -> luc(n) &^ luc(n+1) mod luc(n+2):
map(f, [$0..50]);
MATHEMATICA
Table[PowerMod @@ LucasL[n + {0, 1, 2}], {n, 0, 50}] (* Amiram Eldar, Dec 23 2021 *)
PROG
(Python)
from gmpy2 import lucas2
def A350087(n):
a, b = lucas2(n+1)
return pow(b, a, a+b) # Chai Wah Wu, Dec 28 2021
CROSSREFS
Sequence in context: A137224 A191830 A155944 * A091232 A336398 A209337
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Dec 22 2021
STATUS
approved