login
a(n) = A356867(n) / gcd(n, A356867(n)).
5

%I #8 Sep 15 2023 23:09:21

%S 1,1,1,5,4,1,10,1,1,7,14,5,25,10,4,25,16,1,35,7,10,125,40,1,4,16,1,11,

%T 22,7,55,11,14,35,8,5,49,49,25,35,140,10,250,20,4,245,196,25,625,4,16,

%U 125,64,1,7,55,35,275,88,7,350,56,10,343,98,125,875,70,40,125,160,1,1225,196,4,3125,400,16,1000,8,1,13,26

%N a(n) = A356867(n) / gcd(n, A356867(n)).

%C Denominator of n / A356867(n).

%H Antti Karttunen, <a href="/A365465/b365465.txt">Table of n, a(n) for n = 1..19683</a>

%F a(n) = A356867(n) / A365463(n).

%o (PARI)

%o up_to = 19683;

%o A356867list(up_to) = { my(v=vector(up_to),met=Map(),h=0,ak); for(i=1,#v,if(1==vecsum(digits(i,3)), v[i] = i; h = i, ak = v[i-h]; forprime(p=2,,if(3!=p && !mapisdefined(met,p*ak), v[i] = p*ak; break))); mapput(met,v[i],i)); (v); };

%o v356867 = A356867list(up_to);

%o A356867(n) = v356867[n];

%o A365465(n) = (A356867(n)/gcd(n, A356867(n)));

%Y Cf. A356867, A365462, A365463, A365464 (numerators).

%Y Cf. also A364502.

%K nonn

%O 1,4

%A _Antti Karttunen_, Sep 15 2023