%I#14 2017年5月11日13:43:12

%S 1,5,8,14,17,34,30,44,19,51,68103,93,72196152155103192132,72，

%电话：126278349，32159,53165437976498560709237786739705282，

%电话：15940286366094824333738411301273,491321546288143334912645928256717722761189361652144211921228039027027072046255827732703

%N复合过滤器：a（N）=T（A046523（N），A109395（N）），其中T（N，k）是用作配对函数的序列A000027。

%H Antti Karttunen，n的表，n=1..10000的a（n）</a>

%H数学世界，<a href=“http://mathworld.wolfram.com/PairingFunction.html“>配对功能</a>

%F a（n）=（1/2）*（2+（（A046523（n）+A109395（n））^2）-A046523（n）-3*A109395。

%t表[（2+（#1+#2）^2-#1-3#2）/2&@@{Times@@MapIndexed[Prime[First@#2]^#1&，Sort[FactorInteger[n][[All，-1]]，Greater]]-Boole[n==1]，分母[EulerPhi[n]/n]}，{n，73}]（*Michael De Vlieger_，2017年5月4日*）

%o（PARI）

%o A109395（n）=n/gcd（n，eulerphi（n））；

%o A046523（n）={my（f=vecsort（factor（n）[，2]，4），p）；prod（i=1，#f，（p=nextprime（p+1））^f[i]）；}；\\此功能来自_Charles R Greathouse IV_，2011年8月17日

%o A286149（n）=（1/2）*（2+（（A046523（n）+A109395（n））^2）-A046523（n）-3*A109395*（n）；

%o表示（n=110000，写入（“b286149.txt”，n，“”，A286149（n））；

%o（方案）（定义（A286149 n）（*（/1 2）（+（导出（+（A046523 n）（A109395 n））2）（-（A046533 n））（-

%o（Python）

%o来自症状导入因子int，totient，gcd

%o定义T（n，m）：返回（（n+m）**2-n-3*m+2）/2

%o定义P（n）：

%o f=因子（n）

%o返回排序（[f[i]代表f中的i）

%o定义a046523（n）：

%o x=1

%o为True时：

%o如果P（n）==P（x）：返回x

%o其他：x+=1

%o def a（n）：返回T（a046523（n），n/gcd（n，totient（n）））#_Indranil Ghosh，2017年5月5日

%Y参见A000027、A046523、A109395、A285729、A286142、A286143、A286 144、A28.6152、A2186154、A28160、A286161、A287162、A284163、A283164。

%K nonn公司

%O 1,2号机组

%A _Antti Karttunen，2017年5月4日