%I#11 2017年3月6日17:20:09
%S 125125525101100911099112551137912326129551437915033,
%电话1530316325174822110925105310073150333011352133712737921,
%电话:413034401149393195136753491630136941370319710577201372517743417701181767
%复合数N,使N中的不同数字和N的适当除数中的不同数相同。
%C由于1总是n的除数,因此这是A011531的子序列_米歇尔·马库斯,2014年2月12日
%H Giovanni Resta,n的表格,n=1..10000的a(n)</a>
%e125在序列中是因为125的适当除数是{1,5,25},其数字与125相同。
%t Pn[n]:=排序[DeleteDuplicates[Flatten[InterDigits[Take[Divisions[n],DivisiorSigma[0,n]-1]]]];Fn[n_]:=排序[DeleteDuplicates[IntegerDigits[n]]];lst={};Do[If[!PrimeQ[n]&&Pn[n]===Fn[n],AppendTo[lst,n]],{n,2,10^5}];第一次
%t选择[范围[10^5]!PrimeQ@#&&Union@IntegerDigits@#=Union@@IntegerDigits/@Most@Divisiors@#&](*_Giovanni Resta_,2014年2月12日*)
%o(PARI)isok(n)={my(digs=[]);fordiv(n,d,if(d!=n,digs=concat(diggs,digits(d)));(n!=1)&&!isprime(n)&&vecsort(digs,,8)==vecsort
%Y参考A032741、A035141。
%Y有关变体,请参见A282755。
%K nonn,基础
%O 1,1号机组
%2014年2月12日,拉格瑙市
|