v=pareval(vector(1000,i,()->factor(2^i+1)[1,1]));
vecmin(v)