【 在 feng321 的大作中提到: 】
: 谢谢回复，我们可以假设这个素数列表是递增不重复的。这样程序就可以不考虑重复的情况了。
: “判断规则加上p-1整除N-1则在x，y，z不相同不能成立（找不大）”————这句话什么意思？满足p|N,(p-1)|(N-1),(p+1)|(N+1) 的素数，找不到？不可能的！我已经找到一个了，并且发给了iMx
: 好奇，您是怎么用gpu运算的啊？代码看不出你是用gpu运算啊

三个质数不相同的话， 程序没有返回没有找到结果，再大的指数， 三个相乘也超出32位整数范围了，
除非程序实现长整数乘法加法器，不过，我没数学证明过
--
FROM 124.126.0.*

python可以直接计算大整数乘积，上万位都没问题。我好奇你的代码是如何利用Gpu算的？
【 在 poggy 的大作中提到: 】
:
: 三个质数不相同的话， 程序没有返回没有找到结果，再大的指数， 三个相乘也超出32位整数范围了，
: 除非程序实现长整数乘法加法器，不过，我没数学证明过
--
FROM 120.242.253.*

【 在 feng321 的大作中提到: 】
: python可以直接计算大整数乘积，上万位都没问题。我好奇你的代码是如何利用Gpu算的？

这个很简单， 首先，需要一块计算显卡，一般用Nvidia的容易些， 安装好专用驱动和CUDA库，
python端， 选择一种支持显卡计算的python库， 比如numba， 或者cupy之类的
一般安装tenorflow或者pytorch做深度学习计算也都会需要这些步骤。
这里用的numba的cuda
一个生成素数的程序，从1000开始，搜索24000个整数：

from numba import cuda
import numpy as np
import math
from time import time

@cuda.jit
def gpu_prime_gen(a, b, result, n):
    idx = cuda.threadIdx.x + cuda.blockDim.x * cuda.blockIdx.x    
    if idx < n :
        if(b[idx]%2==0):
            result[idx] = 0
        else:
            result[idx] = 3
            while(result[idx]*result[idx] <= a[idx]):
                if((b[idx] % result[idx]) == 0):
                    result[idx] = 0
                    break
                else:
                    result[idx] += 2
            else:
                result[idx] = b[idx]
def main():
    n = 8000
    x = np.arange(1000,1000+n*3).astype(np.int32)
    y = x+1

    start = time()
    x_device = cuda.to_device(x)
    y_device = cuda.to_device(y)
    out_device = cuda.device_array(n*3)

    threads_per_block = 1024
    blocks_per_grid = math.ceil(n*3 / threads_per_block)

    # 使用默认流
    gpu_prime_gen[blocks_per_grid, threads_per_block](x_device, y_device, out_device, n*3)
    gpu_result = out_device.copy_to_host()    
    cuda.synchronize()
    gpu_result = gpu_result[gpu_result>0]
    print("gpu vector gpu_prime_gen time " + str(time() - start)+" got prime "+ str(len(gpu_result)))
    print(gpu_result[:1024].astype(np.int32))

if __name__ == "__main__":
    main()
--
修改:poggy FROM 124.126.0.*
FROM 124.126.0.*

【 在 feng321 的大作中提到: 】
: 这是什么语法？ else不是跟if一块用吗？这样用也可以？
: [upload=1][/upload]

建议好好学学python基础， 别整什么并行计算， GPU运算这些，
我感觉你的基础还需要加强， 是不是个学生， 慢慢来， 不能急于求成。
--
FROM 124.126.0.*

首先，应该是a+1整除bc-1，你搞反了
其次，你的循环说明我看不懂。开头你已经说了，a,b,c都是素数
【 在 cafitren 的大作中提到: 】
: 假设三个素数是a、b、c，你的要求是abc+1整除a+1且abc-1整除a-1。
: abc+1 = (a+1)bc-bc+1 => bc-1整除a+1
: abc-1 = (a-1)bc+bc-1 => bc-1整除a-1
: ...................
--
FROM 120.242.251.*