%我#32 2024年2月17日04:04:44

%S 2,4,6,6,10,12,14,12,12,20,22,28,26,28,30,20,34,24,38,30,42,44,46,36，

%电话：30,52,36,42,58,60,62,40,66,68,70,42,74,76,78,60,82,84,86,66,60,92,94，

%U 60,56,60102,78106,72110,814116118,90122124,84,72

%N a（N）是N除以k^2的最小k>N。

%C A073353（n）<=a（n）<=2*n。任何除n的素数也必须除a（n），因为n除（2*n）^2。

%H Peter Kagey，n的表格，n=1..5000的a（n）</a>

%F a（n）=平方（n*A072905（n））。

%F a（n）=A019554（n）*（A000188（n）+1）。

%F和{k=1..n}a（k）~c*n^2/2，其中c=1+zeta（3）/zeta（2）=1+A253905=1.73076296940143849872….-_Amiram Eldar，2024年2月17日

%e a（12）=18，因为12除以18^2，但12不除以13^2、14^2、15^2、16^2或17^2。

%tlk[n_]：=模[{k=n+1}，While[！可除[k^2，n]，k++]；k] ；阵列[lk，70]（*哈维·P·达尔，2017年11月5日*）

%o（红宝石）

%o定义a（n）

%o（n+1..2*n）.查找{|k|k**2%n==0}

%o端

%o（PARI）

%o a（n）=对于（k=n+1,2*n，如果（k^2%n==0，返回（k））

%o矢量（100，n，a（n））\\德雷克·奥尔，2015年2月6日

%o（PARI）a（n）=我的（t=系数（n）[，1]）；对于步骤（k=n+t，2*n，t，if（k^2%n==0，return（k）））\\_Charles R Greathouse IV_，2015年2月7日

%o（哈斯克尔）

%o a254732 n=头部[k|k<-[n+1..]，mod（k^2）n==0]

%o——Reinhard Zumkeller，2015年2月7日

%o（Python）

%o定义A254732（n）：

%o。。。。k=n+1

%o。。。。而pow（k，2，n）：

%o。。。。。。。。k+=1

%o。。。。return k#_柴华武_，2015年2月15日

%Y参见A254733（类似，带有k^3）、A254734（类似，带k^4）、A073353（类似的，极限m->无穷大k^m）。

%Y参考A253905。

%K nonn，简单

%O 1,1号机组

%A _彼得·凯吉，2015年2月6日