%I#21 2022年3月23日12:49:14

%序号3,7,13,31,47,57,53,731371031311091372392572349257269，

%电话：331347389409257389251229499487509491541487353739571，

%电话：7437273078833743929827971911887569106375110138831451977125985398479671049

%N a（N）是最小素数（k），k>N，因此A036689（k）或A036690（k）是s（N）+s（N+1）+…+s（j），j<k，其中每个s（i）是A036689（i）或A036690（i）

%C a（n）是最小素数p，使得p*（p-1）或p*（p+1）是序列的和，其中对于i从n到一些j，每个项都是素数（i）*（素数（i）-1）或素数（i）*（素数（i）+1）。

%H Robert Israel，n表，n=1..2400的a（n）</a>

%e a（3）=13，因为素数（3）=5，接下来的两个素数是7和11，以及5*6+7*6+11*10=182=13*14。

%p p：=选择（isprime，[2，seq（i，i=3..10^6,2）]）：

%p R：=转换（映射（p->（p*（p-1），p*（p+1）），p），集合）：

%p f：=proc（n）局部S，T，SR，i，S；

%pS:={p[n]*（p[n]-1），p[n]*（p[n]+1）}；

%n+1 do中i的p

%p T：=[p[i]*（p[i]-1），p[i]*（p[i]+1）]；

%p S：=映射（S->（S+T[1]，S+T[2]），S）；

%p SR:=S与R相交；

%p如果SR<>{}，则

%p s：=（sqrt（1+4*min（SR））-1）/2；

%p如果isprime（s）则返回s，否则返回s+1 fi

%功率因数

%日期

%p端程序：

%p映射（f，[1..100]美元）；

%Y参考A036889、A036890。

%K nonn公司

%O 1,1号机组

%A _J.M.Bergot和_Robert Israel，2022年3月18日