CalmoSoft primes: Difference between revisions

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 w in range(len(pri)-1, 1, -1):

~~for~~ i ~~in range(len~~(pri)~~-window_size):~~

psum = sum(pri[:w])

~~if isprime(sum~~(pri~~[i:i+window_size]~~)):

for d in range(-1, len(pri)-w):

~~print(~~

if d >= 0:

psum -= pri[d]

⚫

~~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[w + 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,

⚫

') of primes less than', maxp, 'totals', sum(pri[d+1:d+w+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(pri[i:i+~~window_size~~], '\n')

print("The sequence is: ", pri[d+1:d+w+1], "\n")

return

⚫

calmo_prime_sequence()

~~calmo_prime_sequence(~~50_000_000)

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:

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, 1 , ... 49999699, 49999711, 49999739, 49999751, 49999753]

</pre>