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.

A351087 After initial 0, numbers k such that A327860(k) is a multiple of k.

This page as a plain text file.
%I A351087 #9 Feb 05 2022 17:24:18
%S A351087 0,1,7,8,16,25,80,91,161,175,203,343,392,440,539,588,616,968,1100,
%T A351087 1225,1331,1333,1337,1375,1400,2556,3025,3773,5112,5875,6655,7280,
%U A351087 7668,8281,8575,8600,9604,9800,10224,10820,10868,11011,11492,12595,12980,13013,15925,17303,17576,17875,20449,22308,23677,24067,24167
%N A351087 After initial 0, numbers k such that A327860(k) is a multiple of k.
%t A351087 Select[Range[25000], Block[{i, m, n = #, p}, m = i = 1; While[n > 0, p = Prime[i]; m *= p^Mod[n, p]; n = Quotient[n, p]; i++]; # == GCD[#, If[m < 2, 0, m Total[#2/#1 & @@@ FactorInteger[m]]]]] &] (* _Michael De Vlieger_, Feb 04 2022 *)
%o A351087 (PARI)
%o A351087 A327860(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= (p^e); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); };
%o A351087 isA351087(n) = (!n || !(A327860(n)%n));
%Y A351087 Fixed points of A351083.
%Y A351087 Cf. A003415, A276086, A327860.
%Y A351087 Cf. A328110 (a subsequence).
%K A351087 nonn
%O A351087 1,3
%A A351087 _Antti Karttunen_, Feb 03 2022