OFFSET
1,6
COMMENTS
a(n) > 1 if and only if n is in A049094.
LINKS
Antti Karttunen (terms 1..310) & Hans Havermann, Table of n, a(n) for n = 1..1206
FORMULA
EXAMPLE
2^42 - 1 = 3^2 * 7^2 * 43 * 127 * 337 * 5419, so a(42) = 3*7 = 21.
MAPLE
a:= n-> mul(i[1]^iquo(i[2], 2), i=ifactors(2^n-1)[2]):
seq(a(n), n=1..105); # Alois P. Heinz, Apr 29 2016
MATHEMATICA
a[n_] := Sqrt[(2^n-1)/Times @@ Power @@@ ({#[[1]], Mod[#[[2]], 2]}& /@ FactorInteger[2^n -1])];
Array[a, 105] (* Jean-François Alcover, Jan 03 2022 *)
PROG
(PARI) a(n)=core(2^n-1, 1)[2]
CROSSREFS
KEYWORD
nonn
AUTHOR
Charles R Greathouse IV, Apr 26 2016
EXTENSIONS
More terms from Antti Karttunen, Sep 23 2018
STATUS
approved