A066808 - OEIS

%I #18 Feb 16 2025 08:32:45

%S 0,1,1,4,1,4,16,4,1,256,16,4,4081,4,16,256,1,4,261121,4,65536,256,16,

%T 4,65536,33554305,16,67108864,65536,4,16,4,1,256,16,262144,

%U 68451041281,4,16,256,65536,4,4398042316801,4,65536,35184371957761,16,4,281474976645121

%N a(n) = F(n)-1 mod 2^n+1 with F(n) = n-th Fermat number = 1+2^2^n.

%C All terms except n=12,18,25,36,42,45,48,55 result in a(n) that are powers of 2, whereas these exceptions (4081, 261121, 33554305, 68451041281, 4398042316801, 35184371957761, 281474976645121, 36020000925941761) are all odd.

%H Alois P. Heinz, <a href="/A066808/b066808.txt">Table of n, a(n) for n = 0..3323</a>

%H Chris Caldwell, <a href="https://primes.utm.edu/glossary/page.php?sort=FermatNumber">Fermat Number</a>, The Prime Glossary.

%H Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/FermatNumber.html">Fermat Number</a>.

%F F(n)-1=1 mod (2^n+1) for all n=2^k because F(n)=2+ F(1)F(2)..F(n-1)

%p a:= n-> 2&^(2^n) mod (2^n+1):

%p seq(a(n), n=0..50);  # _Alois P. Heinz_, Jul 04 2022

%t Table[ PowerMod[ 2, 2^n, 2^n+1 ], {n, 64} ]

%Y Cf. A004249, A007516, A000215, A019434.

%K nonn,look

%O 0,4

%A _Wouter Meeussen_, Jan 19 2002

%E a(0)=0 prepended by _Alois P. Heinz_, Jul 04 2022