Jump to content

CalmoSoft primes: Difference between revisions

m
Line 992:
 
 
def calmo_prime_sequence(N=100maxp):
""" find the largest prime seq in primes < Nmaxp that sums to a prime """
pri = list(primerange(Nmaxp))
for window_sizew in range(len(pri)+-1, 1, -1):
forpsum i= in range(lensum(pri[:w])-window_size):
for d in range(-1, if isprime(sumlen(pri[i:i+window_size])-w):
if d >= print(0:
psum -= pri[d]
f'Longest Calmo prime seq (length {window_size}) of primes less than {N} totals {sum(pri[i:i+window_size])}:')
ifpsum window_size+= >pri[w 24:+ d]
if isprime(psum):
print('[', ', '.join(map(str, pri[i:i+6])), ', ... ', ', '.join(map(str, pri[i+window_size-6:i+window_size])), ']\n', sep='')
print('Longest Calmo prime seq (length', w,
f'Longest Calmo prime seq (length {window_size}') of primes less than', {N}maxp, 'totals', {sum(pri[id+1:id+window_sizew+1])}:')
if w > 24:
print(''.join(list(str(pri[d+1:d+6]))[:-2]), ", ... ",
''.join(list(str(pri[d-5+w:d+w]))[1:]), "\n")
else:
print("The sequence is: ", pri[id+1:id+w+window_size1], '"\n'")
return
 
 
 
calmo_prime_sequence()
calmo_prime_sequence(for pmax in [100, 500_000, 50_000_000)]:
calmo_prime_sequence(pmax)
</syntaxhighlight>{{out}}
<pre>
Longest Calmo prime seq (length 21 ) of primes less than 100 totals 953:
The sequence is: [7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89]
 
Longest Calmo prime seq (length 3001117) of primes less than 50000000 totals 72618848632313:
Longest Calmo prime seq (length 41530 ) of primes less than 500000 totals 9910236647
[7, 11, 13, 17, 19, 23, ... 49999699, 49999711, 49999739, 49999751, 49999753, 49999757]
[2, 3, 5, 7, 1 , ... 499787, 499801, 499819, 499853, 499879]
 
Longest Calmo prime seq (length 3001117 ) of primes less than 50000000 totals 72618848632313:
[7, 11, 13, 17, 19,1 23, ... 49999699, 49999711, 49999739, 49999751, 49999753,] 49999757]
</pre>
 
4,105

edits

Cookies help us deliver our services. By using our services, you agree to our use of cookies.