A244918
Primes p where the digital sum is equal to 68.
Original entry on oeis.org
59999999, 69999899, 69999989, 78998999, 88989899, 88999979, 89699999, 89799989, 89989799, 89989979, 89997899, 89997989, 89999699, 89999969, 97889999, 98699999, 98879999, 98899799, 98979989, 98988899, 98989889, 98997989, 98998979, 98999969
Offset: 1
69999899 is a prime with sum of the digits = 68, hence belongs to the sequence.
Cf. Primes p where the digital sum is equal to k: 2, 11 and 101 for k=2;
A062339 (k=4),
A062341 (k=5),
A062337 (k=7),
A062343 (k=8),
A107579 (k=10),
A106754 (k=11),
A106755 (k=13),
A106756 (k=14),
A106757 (k=16),
A106758 (k=17),
A106759 (k=19),
A106760 (k=20),
A106761 (k=22),
A106762 (k=23),
A106763 (k=25),
A106764 (k=26),
A048517 (k=28),
A106766 (k=29),
A106767 (k=31),
A106768 (k=32),
A106769 (k=34),
A106770 (k=35),
A106771 (k=37),
A106772 (k=38),
A106773 (k=40),
A106774 (k=41),
A106775 (k=43),
A106776 (k=44),
A106777 (k=46),
A106778 (k=47),
A106779 (k=49),
A106780 (k=50),
A106781 (k=52),
A106782 (k=53),
A106783 (k=55),
A106784 (k=56),
A106785 (k=58),
A106786 (k=59),
A106787 (k=61),
A107617 (k=62),
A107618 (k=64),
A107619 (k=65),
A106807 (k=67), this sequence (k=68),
A181321 (k=70).
-
[p: p in PrimesUpTo(100000000) | &+Intseq(p) eq 68];
-
Select[Prime[Range[10000000]], Total[IntegerDigits[#]]==68 &]
-
# see code in A107579: the same code can be used to produce this sequence, by giving the initial term p = 6*10**7-1, for digit sum 68. - M. F. Hasler, Mar 16 2022
A062339
Primes whose sum of digits is 4.
Original entry on oeis.org
13, 31, 103, 211, 1021, 1201, 2011, 3001, 10111, 20011, 20101, 21001, 100003, 102001, 1000003, 1011001, 1020001, 1100101, 2100001, 10010101, 10100011, 20001001, 30000001, 101001001, 200001001, 1000000021, 1000001011, 1000010101, 1000020001, 1000200001, 1002000001, 1010000011
Offset: 1
3001 is a prime with sum of digits = 4, hence belongs to the sequence.
Cf. Primes p with digital sum equal to k: {2, 11 and 101} for k=2; this sequence (k=4),
A062341 (k=5),
A062337 (k=7),
A062343 (k=8),
A107579 (k=10),
A106754 (k=11),
A106755 (k=13),
A106756 (k=14),
A106757 (k=16),
A106758 (k=17),
A106759 (k=19),
A106760 (k=20),
A106761 (k=22),
A106762 (k=23),
A106763 (k=25),
A106764 (k=26),
A048517 (k=28),
A106766 (k=29),
A106767 (k=31),
A106768 (k=32),
A106769 (k=34),
A106770 (k=35),
A106771 (k=37),
A106772 (k=38),
A106773 (k=40),
A106774 (k=41),
A106775 (k=43),
A106776 (k=44),
A106777 (k=46),
A106778 (k=47),
A106779 (k=49),
A106780 (k=50),
A106781 (k=52),
A106782 (k=53),
A106783 (k=55),
A106784 (k=56),
A106785 (k=58),
A106786 (k=59),
A106787 (k=61),
A107617 (k=62),
A107618 (k=64),
A107619 (k=65),
A106807 (k=67),
A244918 (k=68),
A181321 (k=70).
Cf.
A020449 (primes with digits 0 and 1),
A036953 (primes with digits <= 2),
A106100 (primes with largest digit = 2),
A069663,
A069664 (smallest resp. largest n-digit prime with minimum digit sum).
-
[p: p in PrimesUpTo(800000000) | &+Intseq(p) eq 4]; // Vincenzo Librandi, Jul 08 2014
-
N:= 20: # to get all terms < 10^N
B[1]:= {1}:
B[2]:= {2}:
B[3]:= {3}:
A:= {}:
for d from 2 to N do
B[4]:= map(t -> 10*t+1,B[3]) union map(t -> 10*t+3,B[1]);
B[3]:= map(t -> 10*t, B[3]) union map(t -> 10*t+1,B[2]) union map(t -> 10*t+2,B[1]);
B[2]:= map(t -> 10*t, B[2]) union map(t -> 10*t+1,B[1]);
B[1]:= map(t -> 10*t, B[1]);
A:= A union select(isprime,B[4]);
od:
sort(convert(A,list)); # Robert Israel, Dec 28 2015
-
Union[FromDigits/@Select[Flatten[Table[Tuples[{0,1,2,3},k],{k,9}],1],PrimeQ[FromDigits[#]]&&Total[#]==4&]] (* Jayanta Basu, May 19 2013 *)
FromDigits/@Select[Tuples[{0,1,2,3},10],Total[#]==4&&PrimeQ[FromDigits[#]]&] (* Harvey P. Dale, Jul 23 2025 *)
-
for(a=1,20,for(b=0,a,for(c=0,b,if(isprime(k=10^a+10^b+10^c+1),print1(k", "))))) \\ Charles R Greathouse IV, Jul 26 2011
-
select( {is_A062339(p,s=4)=sumdigits(p)==s&&isprime(p)}, primes([1,10^7])) \\ 2nd optional parameter for similar sequences with other digit sums. M. F. Hasler, Mar 09 2022
-
A062339_upto_length(L,s=4,a=List(),u=[10^(L-k)|k<-[1..L]])=forvec(d=[[1,L]|i<-[1..s]], isprime(p=vecsum(vecextract(u,d))) && listput(a,p),1); Vecrev(a) \\ M. F. Hasler, Mar 09 2022
Corrected and extended by Larry Reeves (larryr(AT)acm.org), Jul 06 2001
A213354
Primes p with digit sums s(p) and s(s(p)) also prime, but s(s(s(p))) not prime.
Original entry on oeis.org
59899999, 69899899, 69899989, 69979999, 69997999, 69999799, 77899999, 78997999, 78998989, 78999889, 78999979, 79699999, 79879999, 79889899, 79979899, 79979989, 79988899, 79989979, 79996999, 79997899, 79997989, 79999789, 79999879, 79999987
Offset: 1
59899999 and s(59899999) = 5+9+8+9+9+9+9+9 = 67 and s(s(59899999)) = s(67) = 6+7 = 13 are all primes, but s(s(s(59899999))) = s(13) = 1+3 = 4 is not prime. No smaller prime has this property, so a(1) = 59899999 = A213355(3).
-
Select[Prime[Range[5000000]], PrimeQ[Apply[Plus, IntegerDigits[#]]] && PrimeQ[Apply[Plus, IntegerDigits[Apply[Plus, IntegerDigits[#]]]]] && ! PrimeQ[Apply[Plus, IntegerDigits[Apply[Plus, IntegerDigits[Apply[Plus, IntegerDigits[#]]]]]]] &]
A107349
Indices of primes with digit sum = 67.
Original entry on oeis.org
3556550, 4112560, 4112568, 4116958, 4117950, 4118050, 4554051, 4614334, 4614393, 4614437, 4614444, 4652931, 4662825, 4663397, 4668326, 4668331, 4668793, 4668843, 4669228, 4669275, 4669281
Offset: 1
a(1) = 3556550 because prime(3556550) = 59899999 is the first prime with digit sum = 67.
Showing 1-4 of 4 results.
Comments