cp's OEIS Frontend

A365465 a(n) = A356867(n) / gcd(n, A356867(n)).

%I A365465 #8 Sep 15 2023 23:09:21
%S A365465 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 A365465 22,7,55,11,14,35,8,5,49,49,25,35,140,10,250,20,4,245,196,25,625,4,16,
%U A365465 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 A365465 a(n) = A356867(n) / gcd(n, A356867(n)).
%C A365465 Denominator of n / A356867(n).
%H A365465 Antti Karttunen, <a href="/A365465/b365465.txt">Table of n, a(n) for n = 1..19683</a>
%F A365465 a(n) = A356867(n) / A365463(n).
%o A365465 (PARI)
%o A365465 up_to = 19683;
%o A365465 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 A365465 v356867 = A356867list(up_to);
%o A365465 A356867(n) = v356867[n];
%o A365465 A365465(n) = (A356867(n)/gcd(n, A356867(n)));
%Y A365465 Cf. A356867, A365462, A365463, A365464 (numerators).
%Y A365465 Cf. also A364502.
%K A365465 nonn
%O A365465 1,4
%A A365465 _Antti Karttunen_, Sep 15 2023