OFFSET
1,2
COMMENTS
Greatest common divisor of n and 2^n - 2.
a(n)=n iff n=1 or n is prime or n is Fermat pseudoprime to base 2 or even pseudoprime to base 2. - Corrected by Thomas Ordowski, Jan 25 2016
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 3 because 2^3 - 2 = 6 and gcd(3, 6) = 3.
a(4) = 2 because 2^4 - 2 = 14 and gcd(4, 14) = 2.
MAPLE
seq(igcd(n, (2&^n - 2) mod n), n=1 .. 1000); # Robert Israel, Jan 26 2016
MATHEMATICA
Table[GCD[n, 2^n - 2], {n, 1, 59}] (* Alonso del Arte, Jul 22 2012 *)
PROG
(Java)
import java.math.BigInteger;
public class A214606 {
public static void main (String[] args) {
BigInteger c1 = BigInteger.valueOf(1);
BigInteger c2 = BigInteger.valueOf(2);
for (int n=0; n<222; n++) {
BigInteger bn=BigInteger.valueOf(n), pm2=c1.shiftLeft(n).subtract(c2);
System.out.printf("%s, ", bn.gcd(pm2).toString());
}
}
}
(PARI) a(n)=gcd(n, lift(Mod(2, n)^n-2)) \\ Charles R Greathouse IV, May 29 2014
(Magma) [GCD(n, 2^n-2): n in [1..80]]; // Vincenzo Librandi, Jan 26 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Alex Ratushnyak, Jul 22 2012
STATUS
approved