%I#12 2023年12月18日12:22:12

%S 2,9,28,811761562528865617921377811792024414062530720，

%电话：7971615311296430467211492992381469726562525123486784401，

%电话：38797312242137805625283115525960464477539062518454937651684605176023258369126464664570822595126814724185877482557133644888891

%N a（N）是第一个数==1（mod N），它是N个素数的乘积，用重数计数。

%C a（n）是第一个数字k==1（mod n），因此A001222（k）=n。

%C A053669（n）^n<=a（n）<=A034694（n）。

%C如果n在A007694中，则a（n）=A053669（n）^n。

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

%e a（4）=81，因为81==1（mod 4）和81=3^4是4个素数的乘积，用重数计算，并且没有更小的数起作用。

%pf:=proc（n）使用priqueue；局部p，x，四月，v；

%p初始化（Aprimes）；

%p p：=2；

%p，当n模p=0时，做p:=下一素数（p）od：

%p插入（[-p^n，p，0]，Aprimes）；

%p do公司

%p v：=提取物（Aprimes）；

%p x：=-v[1]；

%p如果x mod n=1，则返回x fi；

%p如果v[3]<n则

%p插入（[v[1]，v[2]，v[3]+1]，四月）；

%p p:=下一素数（v[2]）；

%p，而n模p=0做p:=下一素数（p）od；

%px:=x*（p/v[2]）^（n-v[3]）；

%p插入（[-x，p，v[3]]，Aprimes）；

%p fi；

%p od；

%p端程序：

%p f（1）：=2：

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

%Y请参阅A001222、A007694、A034694和A053669。

%K nonn公司

%O 1,1号机组

%以色列罗贝尔特，2023年12月17日