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.

A228181 Numbers k such that sum of square of prime divisors of k equals sum of prime divisors of k+1.

This page as a plain text file.
%I A228181 #12 Sep 08 2022 08:46:05
%S A228181 12,27,385,1120,4840,9936,14500,29440,95795,105875,178904,223155,
%T A228181 341248,343343,754985,830908,1059630,1841049,2408832,2949375,3564704,
%U A228181 4934358,5368792,5500312,6695000,6805372,8332831,8846656,10126336,12956040,13157235,17254600
%N A228181 Numbers k such that sum of square of prime divisors of k equals sum of prime divisors of k+1.
%C A228181 Numbers k such that A005063(k) = A008472(k+1).
%H A228181 Amiram Eldar, <a href="/A228181/b228181.txt">Table of n, a(n) for n = 1..333</a> (terms below 10^10)
%e A228181 The prime divisors of 9936 are {2, 3, 23} and the prime divisors of 9937 are {19, 523} => 2^2 + 3^2 + 23^2 = 19 + 523 = 542, hence 9936 is in the sequence.
%t A228181 fQ[n_] := Plus @@ (First@# & /@ FactorInteger[n]^2) == Plus @@ (First@# & /@ FactorInteger[n + 1]); Select[ Range@ 100000, fQ]
%o A228181 (Magma) [k:k in [2..2500000]| &+PrimeDivisors(k+1) eq &+[PrimeDivisors(k)[i]^2: i in [1..#PrimeDivisors(k)]]]; // _Marius A. Burtea_, Feb 18 2020
%Y A228181 Cf. A000961, A005063, A008472, A164722.
%K A228181 nonn
%O A228181 1,1
%A A228181 _Michel Lagneau_, Aug 15 2013