编写一个测试程序，用多线程计算1万以内素数之和与完全数之和的乘积。

由于1万以内的素数和完全数数量比较少，使用多线程不会带来明显的性能优势，反而会增加程序的复杂度。不过，为了演示多线程的使用方法，我们还是编写一个多线程的程序。

首先，编写计算素数和完全数的函数：

def is_prime(n):
    if n < 2:
        return False
    for i in range(2, int(n ** 0.5) + 1):
        if n % i == 0:
            return False
    return True


def sum_of_primes(start, end):
    return sum(filter(is_prime, range(start, end)))


def is_perfect(n):
    return n == sum(filter(lambda x: n % x == 0, range(1, n)))


def sum_of_perfects(start, end):
    return sum(filter(is_perfect, range(start, end)))

然后，我们使用多线程计算素数和完全数之和的乘积：

from threading import Thread

def calculate(start, end, result):
    result[0] = sum_of_primes(start, end) * sum_of_perfects(start, end)

def main():
    result = [0]
    threads = []
    for i in range(4):
        t = Thread(target=calculate, args=(i * 2500 + 1, (i + 1) * 2500, result))
        threads.append(t)
        t.start()
    for t in threads:
        t.join()
    print(result[0])

if __name__ == '__main__':
    main()

我们将1万个数分成4个区间，并创建4个线程分别计算每个区间内素数和完全数之和的乘积，最后将结果相加得到最终结果。注意，为了避免线程之间对共享变量result的并发访问，我们使用一个列表来存储结果，这样每个线程都可以通过修改列表元素来修改结果。