Factor a^(3y)+1