cinta第四次作业

7-2,3,6,8

8-1,3,5

8-9：编程完成以下工作：对任意给定的一个素数 p，求出 Zp* 的最小生成元。任取一个整数 n，大于 1 小于 n 的所有素数 p，求 Zp* 的最小生成元，并求以上最小生成元集合中最大者所对应的素数 p。

# 寻找最小生成元
def is_prime(n):
    if n <= 1:
        return False
    if n <= 3:
        return True
    if n % 2 == 0 or n % 3 == 0:
        return False
    i = 5
    while i * i <= n:
        if n % i == 0 or n % (i + 2) == 0:
            return False
        i += 6
    return True

def find_root(p):
    if not is_prime(p):
        return None  # 输入不是素数，无法找到最小生成元
    if p == 2:
        return 1  # 对于素数2，1是它的最小生成元

    factors = [p - 1]  # 模p的最小生成元的指数必须是p-1
    for i in range(2, int(p**0.5) + 1):
        if (p - 1) % i == 0:
            factors.append(i)
            factors.append((p - 1) // i)

    for g in range(2, p):
        is_primitive = True
        for factor in factors:
            if pow(g, factor, p) == 1:
                is_primitive = False
                break
        if is_primitive:
            return g

    return None  # 找不到最小生成元

# 输入素数p
p = int(input("请输入一个素数 p："))
primitive_root = find_root(p)
if primitive_root is not None:
    print(f"Z_{p}* 的最小生成元是 {primitive_root}")
else:
    print(f"无法找到 Z_{p}* 的最小生成元")

# 任取一个整数 n，大于 1 小于 n 的所有素数 p，求 Zp* 的最小生成元，并求以上最小生成元集合中最大者所对应的素数 p。
def is_prime(num):
    if num <= 1:
        return False
    if num == 2:
        return True
    if num % 2 == 0:
        return False

    for i in range(3, int(num**0.5) + 1, 2):
        if num % i == 0:
            return False

    return True

def prime_factors(n):
    factors = []
    for i in range(2, n + 1):
        while n % i == 0:
            factors.append(i)
            n //= i
    return factors

def is_primitive_root(g, p):
    if pow(g, p - 1, p) != 1:
        return False

    factors = prime_factors(p - 1)
    for factor in factors:
        if pow(g, (p - 1) // factor, p) == 1:
            return False

    return True

def find_min_primitive_root(n):
    max_primitive_root = 1
    max_p = 0

    for p in range(2, n):
        if is_prime(p):
            g = 2
            while g < p:
                if is_primitive_root(g, p):
                    if g > max_primitive_root:
                        max_primitive_root = g
                        max_p = p
                    break
                g += 1

    return max_p

# 输入整数 n
n = int(input("请输入一个整数 n："))

# 找到大于1小于n的所有素数p，并计算Zp*的最小生成元
max_p = find_min_primitive_root(n)

if max_p != 0:
    print(f"Zp* 的最小生成元集合中最大的素数 p 是: {max_p}")
else:
    print(f"在范围内未找到素数 p")