%I
%S 1,2,3,4,5,6,9,8,7,10,15,12,25,18,15,16,11,14,21,20,27,30,45,24,49,50,
%T 75,36,125,30,81,32,45,22,45,28,55,42,75,40,77,54,105,60,35,90,135,48,
%U 121,98,33,100,245,150,75,72,63,250,375,60,625,162,63,64,125,90,39,44,135,90,99,56,91,110,147,84,135,150,189,80,143,154,231,108,55
%N Multiplicative with a(p^e) = A005940(p^e).
%C Question: are there any other numbers n besides 1 and those in A070776, for which a(n) = A005940(n)? At least not below 2^25. This is probably easy to prove.
%H Antti Karttunen, <a href="/A324106/b324106.txt">Table of n, a(n) for n = 1..16384</a>
%H <a href="/index/Bi#binary">Index entries for sequences related to binary expansion of n</a>
%e For n = 85 = 5*17, a(85) = A005940(5) * A005940(17) = 5*11 = 55. Note that A005940(5) is obtained from the binary expansion of 51 = 4, which is "100", and A005940(17) is obtained from the binary expansion of 171 = 16, which is "1000".
%o (PARI)
%o A005940(n) = { my(p=2, t=1); n; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t }; \\ From A005940
%o A324106(n) = { my(f=factor(n)); prod(i=1, #f~, A005940(f[i,1]^f[i,2])); };
%Y Cf. A005940, A070776, A324107 (fixed points), A324108, A324109.
%K nonn,mult
%O 1,2
%A _Antti Karttunen_, Feb 15 2019
