You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
28 lines
788 B
28 lines
788 B
modulus = 4002409555221667393417789825735904156556882819939007885332058136124031650490837864442687629129015664037894272559787
|
|
|
|
assert(modulus.is_prime())
|
|
|
|
Fp = GF(modulus)
|
|
|
|
generator = Fp(0);
|
|
for i in range(0, 20):
|
|
i = Fp(i);
|
|
neg_i = Fp(-i)
|
|
if not(i.is_primitive_root() or neg_i.is_primitive_root()):
|
|
continue
|
|
elif i.is_primitive_root():
|
|
assert(i.is_primitive_root());
|
|
print("Generator: %d" % i)
|
|
generator = i
|
|
break
|
|
else:
|
|
assert(neg_i.is_primitive_root());
|
|
print("Generator: %d" % neg_i)
|
|
generator = neg_i
|
|
break
|
|
|
|
|
|
two_adicity = valuation(modulus - 1, 2);
|
|
trace = (modulus - 1) / 2**two_adicity;
|
|
two_adic_root_of_unity = generator^trace
|
|
print("2-adic Root of Unity: %d " % two_adic_root_of_unity)
|