%I#30 2016年3月28日08:02:35

%S 1,2,4,3,2,2,4,15,2,2,6,3,10,3,26,3,4,2,2,2,15,26,7,4,2,2,6,2,2,10,10，

%T 2,20,4,2,3,4,4,6,6,4,10,2,14,16,12,3,9,10,6,24,3,4_6,2,2,2,18，

%U 6、6、3、14、5、16、9、18、3、2、2、4、3、10、6

%N最小的k，使得k*N^3+1是素数。

%e a（1）=1，因为为了使k*（1^3）+1成为最小素数，k必须是1（1*（1*3）+1=2）。

%e a（2）=2，因为为了使k*（2^3）+1成为最小素数，k必须是2（2*（2*3）+1=17）。

%e a（3）=4，因为为了使k*（3^3）+1成为最小素数，k必须是4（4*（3*3）+1=109）。

%tsk[n_]：=模[{k=1，n3=n^3}，While[！素数Q[k*n3+1]，k++]；k] ；阵列[sk，80]（*哈维·P·戴尔，2014年8月27日*）

%t表[SelectFirst[Range[10^2]，PrimeQ[#n^3+1]&]，{n，76}]（*_Michael De Vlieger_，2016年3月27日，第10版*）

%o（Python）

%o导入交响乐

%o来自sympy import isprime

%o定义f（n）：

%哦。。对于范围（1,10**3）中的k：

%o。。。。如果是素数（k*（n**3）+1）：

%o。。。。。。返回k

%o n=1

%o当n<10**3时：

%o。。打印（f（n））

%o。。n+=1

%o（PARI）a（n）={k=1；while（！isprime（k*n^3+1），k++）；k；}\\_Altug Alkan_，2016年3月26日

%Y参考A034693、A035092。

%K nonn公司

%O 1,2号机组

%2014年3月6日，阿德里克·奥尔