|
%I #29 Sep 08 2022 08:46:01
%S 11,13,17,19,31,41,61,71,101,103,107,109,113,127,131,137,139,149,151,
%T 157,163,167,173,179,181,191,193,197,199,211,241,251,271,281,311,313,
%U 317,331,401,419,421,431,461,491,521,541,571,601,613,617,619,631,641,661
%N Primes containing a digit 1.
%C Subsequence of A011531, A062634, A092911 and A092912.
%C Supersequence of A106101, A045707 and A030430.
%C Complement of A208271 with respect to A011531.
%H Vincenzo Librandi, <a href="/A208270/b208270.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/De#decimal_expansion">Index entries for sequences related to decimal expansion of n</a>
%F a(n) ~ n log n since the sequence contains almost all primes. - _Charles R Greathouse IV_, Mar 04 2012
%t Select[Prime[Range[124]], MemberQ[IntegerDigits[#], 1] &](* _Jayanta Basu_, Apr 01 2013 *)
%t Select[Prime[Range[200]],DigitCount[#,10,1]>0&] (* _Harvey P. Dale_, Dec 15 2020 *)
%o (PARI) forprime(p=2,1e3,s=vecsort(eval(Vec(Str(p))),,8);if(s[1]==1||(s[1]==0&&s[2]==1),print1(p", "))) \\ _Charles R Greathouse IV_, Mar 04 2012
%o (PARI) is_A208270(n)=isprime(n)&setsearch(Set(Vec(Str(n))),1) \\ _M. F. Hasler_, Mar 05 2012
%o (Magma) [p: p in PrimesUpTo(400) | 1 in Intseq(p)]; // _Vincenzo Librandi_, Apr 29 2019
%Y Cf. A208271 (nonprimes containing a digit 1), A011531 (numbers containing a digit 1).
%Y Complement of A038603 in A000040. - _M. F. Hasler_, Mar 05 2012
%K nonn,base,easy
%O 1,1
%A _Jaroslav Krizek_, Mar 04 2012
|
