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.

A203182 Primes p such that A008472(p-1) = A008472(p+1), where A008472 = sum of distinct primes dividing n.

This page as a plain text file.
%I A203182 #11 Sep 22 2013 14:51:13
%S A203182 3,18913,24733,29633,32429,42719,45751,46103,61409,117991,149351,
%T A203182 171529,174019,176017,223099,294893,326369,363691,421727,423503,
%U A203182 434237,472631,658579,678077,686423,706841,735901,770059,771629,906949,936827,937571,1073447,1256029
%N A203182 Primes p such that A008472(p-1) = A008472(p+1), where A008472 = sum of distinct primes dividing n.
%C A203182 Conjecture: the sequence is infinite.
%H A203182 Donovan Johnson, <a href="/A203182/b203182.txt">Table of n, a(n) for n = 1..2000</a>
%e A203182 18913 is in the sequence because:
%e A203182 sum of the distinct prime divisors of 18912 = 2+3+197 = 202;
%e A203182 sum of the distinct prime divisors of 18914 = 2+7+193 = 202.
%p A203182 with(numtheory):for n from 1 to 100000 do:p:=ithprime(n):p1:=p-1: p2:=p+1:t1:=ifactors(p1)[2]; t11 := sum(t1[i][1], i=1..nops(t1)):t2:=ifactors(p2)[2]; t22 := sum(t2[i][1], i=1..nops(t2)):if t11=t22 then printf(`%d, `,p):else fi:od:
%t A203182 Select[Prime[Range[100000]],Total[Transpose[FactorInteger[#-1]][[1]]] == Total[Transpose[FactorInteger[#+1]][[1]]]&] (* _Harvey P. Dale_, Sep 22 2013 *)
%Y A203182 Cf. A008472, A190680, A086711.
%K A203182 nonn
%O A203182 1,1
%A A203182 _Michel Lagneau_, Dec 30 2011