Статья Автор: Лебедев Дмитрий Алексеевич

Простые числа_файлы и не только

Создаем решето

from time import process_time as prt
N = 10000000
t0 = prt()
A = [0, 0] + [1]*(N-1)
k = 2
while k * k <= N :
    if A[k]:
        for i in range (k*k, N+1, k): A[i] = 0
    k += 1
S =[1] + sorted({i for i in range(N+1) if A[i]})

print(prt() - t0)
print(*S[:20])

from time import process_time as prt
t0 = prt()
fname = 'pr_10000000.ch'
pr = [int(s) for s in open(fname)]
print(len(pr),prt() - t0)
A, Z = [0] *10, [0]*10
AZ = []
for x in pr :
  y = int(str(x)[::-1])
  a =int(str(x)[0])
  z =int(str(x)[-1])
  A[a] += 1
  Z[z] += 1
  if x == y : AZ.append(x)
  
print(*[(i, int(A[i]*1000/len(pr))) for i in range(10)])  
print(*[(i, int(Z[i]*1000/len(pr))) for i in range(10)])  
print(len(AZ))
print(AZ[-20])

Прикрепленные файлы
pr_1000000.ch
pr_10000000.ch

Печать