CalmoSoft primes: Difference between revisions

Content added Content deleted
(Undo revision 340330 by Wherrera (talk))
Tag: Undo
Line 992: Line 992:




def calmo_prime_sequence(N=100):
def calmo_prime_sequence(maxp):
""" find the largest prime seq in primes < N that sums to a prime """
""" find the largest prime seq in primes < maxp that sums to a prime """
pri = list(primerange(N))
pri = list(primerange(maxp))
for window_size in range(len(pri)+1, 1, -1):
for win in range(len(pri)-1, 1, -1): # window size
for i in range(len(pri)-window_size):
psum = sum(pri[:win])
for bot in range(-1, len(pri)-win): # the last bottom of window
if isprime(sum(pri[i:i+window_size])):
print(
if bot >= 0:
psum -= pri[bot]
f'Longest Calmo prime seq (length {window_size}) of primes less than {N} totals {sum(pri[i:i+window_size])}:')
if window_size > 24:
psum += pri[win + bot]
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', win,
') of primes less than', maxp, 'totals', sum(pri[bot+1:bot+win+1]))
if win > 24:
print('[', ', '.join(map(str, pri[bot+1:bot+7])), ', ... ',
', '.join(map(str, pri[bot-5+win:bot+win+1])), ']\n', sep='')
else:
else:
print(pri[i:i+window_size], '\n')
print('The sequence is:', pri[bot+1:bot+win+1], '\n')
return
return




for pmax in [100, 500_000, 50_000_000]:
calmo_prime_sequence()
calmo_prime_sequence(50_000_000)
calmo_prime_sequence(pmax)
</syntaxhighlight>{{out}}
</syntaxhighlight>{{out}}
<pre>
<pre>
Longest Calmo prime seq (length 21) of primes less than 100 totals 953:
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]
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
[2, 3, 5, 7, 11, 13, ... 499787, 499801, 499819, 499853, 499879, 499883]

Longest Calmo prime seq (length 3001117 ) of primes less than 50000000 totals 72618848632313
[7, 11, 13, 17, 19, 23, ... 49999699, 49999711, 49999739, 49999751, 49999753, 49999757]
[7, 11, 13, 17, 19, 23, ... 49999699, 49999711, 49999739, 49999751, 49999753, 49999757]
</pre>
</pre>