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.

A360979 Primes that share no digits with their digit sum.

This page as a plain text file.
%I A360979 #15 Mar 04 2023 02:06:41
%S A360979 11,13,17,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97,101,103,
%T A360979 107,113,131,151,211,223,227,229,233,239,241,251,257,263,269,277,293,
%U A360979 311,313,331,349,353,359,383,389,401,409,421,431,433,439,443,449,457,479,487,499,503,509,521,523
%N A360979 Primes that share no digits with their digit sum.
%H A360979 Robert Israel, <a href="/A360979/b360979.txt">Table of n, a(n) for n = 1..10000</a>
%e A360979 a(5) = 29 is a term because 29 is prime and 2+9 = 11 shares no digits with 29.
%p A360979 filter:= proc(n) local L,s;
%p A360979   L:= convert(n,base,10);
%p A360979   s:= convert(L,`+`);
%p A360979   convert(convert(s,base,10),set) intersect convert(L,set) = {}
%p A360979 end proc:
%p A360979 select(filter, [seq(ithprime(i),i=1..1000)]);
%t A360979 Select[Prime[Range[99]],!IntersectingQ[IntegerDigits[#],List[Total[IntegerDigits[#]]]]&] (* _Stefano Spezia_, Feb 28 2023 *)
%o A360979 (PARI) isok(p) = isprime(p) && !#setintersect(Set(digits(sumdigits(p))), Set(digits(p))); \\ _Michel Marcus_, Feb 27 2023
%Y A360979 Primes not in A158473.
%Y A360979 Cf. A007605, A007953.
%K A360979 nonn,base
%O A360979 1,1
%A A360979 _Robert Israel_, Feb 27 2023