%I#12 2019年3月6日08:36:46
%S 1,2,6,4,10,6,42,8,18,10,66,12156,42,60,16,34,18228,20,42,66,92,24,
%电话300156108,84116,60310,32,66,34420,3622228156,40246,42172,
%电话:132180,922820,4858830024156212108660168228116590,60
%N a(N)是属于A306263的N的最小正倍数(如果有的话),或者a(N)=-1。
%C每n>0,a(n)>0吗?
%F a(A306263(n))=任意n>0时为A306262(n)。
%e对于n=7:
%e-7的除数是:1,7,
%e-相应的汉明权重为:1、3、,
%e-3不除以7,
%e-3×7的除数是:1,3,7,21,
%e-相应的汉明权重为:1、2、3、3、,
%e-2不除以3*7,
%e-2*3*7的除数是:1,2,3,6,7,14,21,42,
%e-相应的汉明权重为:1,1,2,2,3,3,
%e-他们都除以2*3*7,
%e-因此a(7)=2*3*7=42。
%t使用[{s=Select[Range[3000],使用[{k=#},AllTrue[Divisors@k,Mod[k,DigitCount[#,2,1]]==0&]]&]},Table[SelectFirst[s,Mod[#,n]==0&]/。k_/;MissingQ@k->-1,{n,60}]](*Michael De Vlieger_,2019年3月5日*)
%o(PARI)a(n)=而(1,my(m=n);对于div(m,d,m=lcm(m,hammingweight(d)););如果(n==m,返回(n),n=m))
%Y参考A000120、A306263。
%K nonn,基础
%O 1,2号机组
%2019年3月3日,A Remy Sigrist
|