Efficient polynomial time algorithms computing industrial-strength primitive roots