NewRingFp(P: INT): RING
/**/ P := NextPrime(1000); /**/ Fp := NewRingFp(P); /**/ Use Fp[x]; /**/ product([x-i | i in 1..P]); x^1009 - x