fermats(n) := makelist(2^(2^i)+1, i, 0, n)$
fermats(9);
map(ifactors, fermats(7));