%I#28 2023年3月9日06:24:46

%编号2047158417038006058911341293609485152551437178110761055201，

%电话：413413211299212211217547859878787235851228370670192243914382191，

%电话：29341325713808066913211064404682551259205381381347461340596158825,6590614118359482705

%反对偶向下读取的N平方数组A（N，k）：对于k，N>=2，具有k个不同素因子的最小基-N强费马伪素。

%C数组A（n，k）的开头如下：

%C k=2 3 4 5 6

%编号：2047 15841 800605 293609485 10761055201

%账号=3:703 8911 152551 41341321 12283706701

%Cn=4:341 4371 129921 9224391 2592053871

%编号：781 24211 4382191 381347461 9075517561

%电话：217 29341 3405961 557795161 333515107081

%H Daniel Suteu，n表，n（n）表示n=2..137</a>

%o（PARI）

%o strong_check（p，base，e，r）=我的（tv=估价（p-1，2））；tv>e&&Mod（基本，p）^（（p-1）>>（tv-e））==r；

%o strong_fermat_psp（A，B，k，base）=A=最大值（A，vecprod（质数（k）））；（f（m，l，lo，k，e，r）=my（list=list（））；my（hi=sqrtnint（B\m，k））；如果（lo>hi，返回（列表））；如果（k==1，对于步长（p=升力（1/Mod（m，l））），hi，l，如果（isprimepower（p）&&gcd（m*base，p）==1&strong_check（p，base，e，r），my（n=m*p）；如果（n>=A&&（n-1）%znorder（Mod（base，p））==0，listput（list，n））），则表示prime（p=lo，hi，base%p==0&&next；strong_check（p，base，e，r）||next；my（z=znorder（Mod（base，p）））；gcd（m，z）==1||下一个；my（q=p，v=m*p）；而（v<=B，list=concat（list，f（v，lcm（l，z），p+1，k-1，e，r））；q*=p；Mod（基数，q）^z==1||中断；v*=p））；列表）；我的（res=f（1，1，2，k，0，1））；对于（v=0，logint（B，2），res=concat（res，f（1，1，2，k，v，-1））；向量排序（Set（res））；

%o T（n，k）=如果（n<2，返回（））；my（x=vecprod（素数（k）），y=2*x）；而（1，my（v=strong_fermat_psp（x，y，k，n））；如果（#v>=1，返回（v[1]））；x=y+1；y=2*x）；

%o打印表（n，k）=对于（x=2，n，对于（y=2，k，打印1（T（x，y），“，”））；打印（“”）；

%o表示（k=2,9，表示（n=2,k，打印1（T（n，k-n+2）“，”））；

%Y参见A001262、A180065（第n=2行）、A271873。

%K nonn，表

%O 2.1个

%2023年3月4日，安妮尔·苏图