%I#16 2023年5月19日04:24:32
%S 6,12,10,18,14,24,20,22,15,36,40,26,48,28,30,21,34,54,45,38,50,42,44,
%电话:60,46,72,56,33,80,52,96,98,58,39,90,62,84,66,75,68,70108,63,74120,
%U 76,51,78100144,82126,57,86,35,88150,92,94147162112132160102,69104106192110196135
%N a(N)是最大的k<A126706(N),因此rad(k)=rad(A126708(N)),其中rad(N)=A007947(N)。
%C A024619的置换。
%C设m=A126706(n),R_m是数字k的序列,使得rad(k)=rad(m)。a(n)给出了R_m中m的前身。
%H Michael De Vlieger,n的表,n=1..10000的a(n)</a>
%H Michael De Vlieger,a(n)的对数散点图,n=1..2^16。
%H Michael De Vlieger,a(n)的对数散点图,n=1..2^12,用绿色显示无平方数,并用大的浅蓝色点高亮显示A286708中的数字k。
%e A126706(1)=12;使rad(k)=rad(12)=6的最小k<12是a(1)=6。
%e A126706(2)=18;使rad(k)=rad(18)=6的最小k<18是a(2)=12。
%e A126706(3)=20;最小的k<20,使得rad(k)=rad(20)=10是a(3)=10。
%t rad[x_]:=倍数@系数整数[x][[All,1]];表[Function[r,SelectFirst[Range[m-1,1,-1],rad[#]==r&]][rad[m]],{m,Select[Range[225],Nor[SquareFreeQ[#],PrimePowerQ[#]]&]}]
%Y参考A007947、A126706、A360529。
%K nonn公司
%O 1,1号机组
%A _米歇尔·德弗里格,2023年5月1日