%I#119 2023年11月24日09:00:10

%第4,3,4,1570,11464,54,60140251932,17194,153612,04746,1540150页，

%电话：70601388009461603348093135,036661826058,190510168369824，

%电话：41887557246，0182394640771417894341036078，131452015870194942150447002419443871308730

%N使N*k^N+/-1是孪生素数的最小数k，如果不存在这样的数，则a（N）=0。

%如果n在A014574中，则C a（n）=1。

%C如果a（n）=0，则n在A097764中。

%C如果a（n）>1，则A367566（n）除以a（n_Jon E.Schoenfield_2023年11月23日

%H Jon E.Schoenfield，n表，n=1..200的a（n）（米歇尔·马库斯的前77个术语）

%e1*1^1+1（2）和1*1^1-1（0）不是素数。1*2^1+1（3）和1*2^1-1（1）并不都是素数。1*3^1+1（4）和1*3*1-1（2）不是素数。1*4^1+1（5）和1*4*1-1（3）都是质数。因此，a（1）=4。

%t zeroQ[n_]：=模块[{f=FactorInteger[n]}，pow=GCD@@f[[；；，2]]；n>4&&AnyTrue[Divisors[pow]，#>1&&Divisible[n，#]&]]；

%t a[n_，kmax_]：=模[{k=1}，如果[zeroQ[n]，0，While[k<=kmax&&！和@@PrimeQ[n*k^n+{-1，1}]，k++]；如果[k<kmax，k，-1]]；表[a[n，10^6]，{n，1，25}]（*_Amiram Eldar_，2023年11月18日，如果搜索限制超过kmax*，则返回-1）

%o（PARI）bot（n）=对于（k=1，10^5，如果（ispseudoprime（n*k^n-1），如果（ispseudoprime（n*k^n+1），return（k）））；

%o n=1；而（n<100，打印1（bot（n），“，”）；n+=1）

%o（PARI）a（n）=如果（（n==16）||（n==27）||；我的（k=1）；while（！（ispseudoprime（n*k^n-1）&&ispseudop时间（n*k^n+1）），k++）；k、 \\_米歇尔·马库斯，2023年11月18日

%Y参见A097764、A014574、A239938、A239666、A367566。

%K非n

%O 1,1号机组

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

%E a（46）摘自Giovanni Resta，2014年3月31日