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.

A063770 Numbers k such that Sum_{j=1..k} sigma(j) divides Product_{j=1..k} phi(j).

This page as a plain text file.
%I A063770 #16 Feb 28 2024 19:08:01
%S A063770 1,28,52,53,55,63,76,159,166,176,219,230,289,302,303,318,321,327,348,
%T A063770 360,364,365,381,383,402,417,430,434,438,444,451,452,454,465,469,478,
%U A063770 504,512,522,530,531,557,559,570,572,584,595,613,629,631,641,645,684
%N A063770 Numbers k such that Sum_{j=1..k} sigma(j) divides Product_{j=1..k} phi(j).
%H A063770 Harry J. Smith, <a href="/A063770/b063770.txt">Table of n, a(n) for n = 1..1000</a>
%t A063770 Module[{nn=700,s,p},s=Accumulate[DivisorSigma[1,Range[nn]]];p=FoldList[ Times,EulerPhi[ Range[nn]]]; Position[ Thread[{s,p}],_?(Mod[#[[2]],#[[1]]]==0&),1,Heads->False]]//Flatten (* _Harvey P. Dale_, Feb 28 2024 *)
%o A063770 (PARI) j=[]; for(n=1,1300,a=sum(k=1,n, sigma(k)); b=prod(k=1,n,eulerphi(k)); if(Mod(b,a)==0,j=concat(j,n))); j
%o A063770 (PARI) { n=0; s=0; p=1; for (m=1, 10^9, s+=sigma(m); p*=eulerphi(m); if(p%s == 0, write("b063770.txt", n++, " ", m); if (n==1000, break)) ) } \\ _Harry J. Smith_, Aug 30 2009
%Y A063770 Cf. A001088, A024916.
%K A063770 nonn
%O A063770 1,2
%A A063770 _Jason Earls_, Aug 15 2001