cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A073753 a(n) = A073752(A073752(n)), where A073752(n) = gcd(n/spf(n), n/gpf(n)), with spf(n) as the smallest and gpf(n) as the greatest prime factor of n (see A020639, A006530).

This page as a plain text file.
%I A073753 #12 Dec 03 2021 04:48:03
%S A073753 1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,4,1,1,1,1,1,1,1,2,1,1,3,1,1,1,1,8,1,1,
%T A073753 1,1,1,1,1,2,1,1,1,1,1,1,1,4,1,1,1,1,1,3,1,2,1,1,1,1,1,1,1,16,1,1,1,1,
%U A073753 1,1,1,2,1,1,1,1,1,1,1,4,9,1,1,1,1,1,1,2,1,3,1,1,1,1,1,8,1,1,1,1,1,1,1,2
%N A073753 a(n) = A073752(A073752(n)), where A073752(n) = gcd(n/spf(n), n/gpf(n)), with spf(n) as the smallest and gpf(n) as the greatest prime factor of n (see A020639, A006530).
%C A073753 a(n) = A062378(n) for n<90, but a(90)=3 <> A062378(90)=1.
%H A073753 Antti Karttunen, <a href="/A073753/b073753.txt">Table of n, a(n) for n = 1..10000</a>
%F A073753 a(n) = A073752(A073752(n)).
%t A073753 A073752[n_] := With[{f = FactorInteger[n][[All, 1]]}, GCD[n/f[[1]], n/f[[-1]]]];
%t A073753 a[n_] := A073752[A073752[n]];
%t A073753 Array[a, 100] (* _Jean-François Alcover_, Dec 03 2021 *)
%Y A073753 Cf. A006530, A020639, A073752.
%Y A073753 Differs from A062378 for the first time at n=90.
%K A073753 nonn
%O A073753 1,8
%A A073753 _Reinhard Zumkeller_, Aug 07 2002
%E A073753 Name edited by _Antti Karttunen_, Sep 23 2017