%I#32 2023年5月20日14:49:52

%S 0,0,0,1,0,0，0,0'，0,0，

%温度0,4,0,6,2,6,0,12,0,5,3,7,0,14,0,15,2,10,1,0,4,13,0,19,0,12,7,7,1，

%U 13,0,14,3,11,0,13,9,8,2,21,0,19,7,21,0,18,2,13,9,22,0,21,16,10,16

%N k＜N的个数，使得gcd（k，N）＞1，gcd（N^2-1，k）＝1，rad（k）不除N。

%A272619第n行中与（n-1）*（n+1）互素的项数。

%H Michael De Vlieger，n的表，a（n）表示n=1..10000</a>

%H Michael De Vlieger，a（n）的散点图，n=1..2*10^5。

%H Michael De Vlieger，a（n）的对数散点图，n=1..2*10^5。

%F a（n）<=A243822（n）。

%e设S（n）=A272619的第n行。

%e a（p）=0，因为S（p）为空。

%e a（4）=0，因为S（4）为空。

%e a（6）=0，因为S（6）为空。

%e a（8）=0，因为S（8）={6}，但gcd（6，（8+1））=3。

%e a（10）=0，因为S（10）={6}，但gcd（6，（10-1））=3。

%e a（12）=1，因为S（12）={10}，gcd（10143）=1。

%e a（16）=1，因为S（16）={6，10，12，14}，但只有14是这样的gcd（14，255）=1。

%e a（18）=3，因为S（18）={10,14,15}，并且它们都不与323共享素因子。

%e a（20）=0，因为S（20）={6，12，14，15，18}，所有这些都与21共用一个因子。

%t表格[Count[Range[k]，_？（Nor[CoprimQ[#，k]，GCD[k^2-1，#]>1，可除[k，Times@@FactorInteger[#][[All，1]]]&）]，{k，120}]

%o（PARI）rad（n）=因子回复（因子（n）[，1]）；\\A007947号

%o a（n）=总和（k=1，n-1，（gcd（k，n）>1）&（gcd）（n^2-1，k）==1）&（n%rad（k））；\\_米歇尔·马库斯，2023年5月20日

%Y参见A045763、A243823、A272619。

%K nonn公司

%O 1,18号

%A _Michael De Vlieger_，2023年5月19日