%I#6 2014年9月17日20:48:42
%S 1290733393726429337146145049531956075162911790118009,
%电话:18441194132034920655209432121921267215732162213722191,
%电话:23355248312564725731260192614526163274052753728035282152828227283052960130429
%N个数x,使其所有循环置换的和等于sigma(x)的所有循环置换和Euler totiten函数phi(x)所有循环置换之和。
%C A247315和A247316的交叉口。
%C所有数字似乎都是3的倍数。
%C在a(135)与a(136)、a(1387)和a(1388)之间迈出了一大步,。。。
%H Paolo P.Lava,n表,n=1..1500的a(n)</a>
%H Paolo P.Lava,<a href=“/A247317/A247317.jpg”>序列前1500项的绘图</a>
%e 4371的循环置换之和为4371+1437+7143+3714=16667;sigma(4371)=6144,其循环置换之和为6144+4614+4461+1446=16667;φ(4371)=2760,其循环置换之和为2760+276+6027+7602=16667。
%e 24831的循环置换之和为24831+12483+31248+83124+48312=19998;sigma(24831)=37440,其循环置换之和为37440+3744+4044+44037+74403=19998;φ(24831)=15840,其循环置换之和为15840+1584+40158+84015+58401=199998。
%p with(numtheory):p:=进程(q)局部a,b,c,d,f,k,n;
%p代表n从1到q做a:=n;b: =a;c: =ilog10(a);
%p表示k从1到c做a:=(a mod 10)*10^c+trunc(a/10);b: =b+a;od;
%p a:=σ(n);d: =a;c: =ilog10(a);
%k从1到c的p do a:=(a mod 10)*10^c+trunk(a/10);d: =d+a;od;
%pa:=φ(n);f: =a;c: =ilog10(a);
%p表示k从1到c做a:=(a mod 10)*10^c+trunc(a/10);f: =f+a;od;
%p如果b=d和d=f,则打印(n);fi;od;结束:P(10^9);
%Y A000010、A000203、A247315、A247316。
%K nonn,基础,简单
%O 1,2号机组
%A _Paolo P.Lava,2014年9月12日
|