Количество пар элементов последовательности, разность которых четна и хотя бы одно из чисел делится на 31

with open('C:/Users/333/Downloads/17 (2).txt') as f:
    s = [int(x) for x in f]
    n = []
    for i in range(len(s) -1):
        if (s[i] % 31 == 0 or s[i + 1] % 31 == 0) and (abs(s[i]) - (abs(s[i + 1])) % 2 == 0):
            n.append(s[i] + s[i + 1])
print(len(n), max(n))

num_even = num_31_even_1 = num_31_even_2 = 0
num_odd = num_31_odd_1 = num_31_odd_2 = 0
k_odd = k_even = k_31_odd = k_31_even = 0
 
with open('17.txt') as f:
    for num in f:
        num = int(num)
        if num%31:
            if num%2:
                num_odd = max(num_odd, num)
                k_odd += 1
            else:
                num_even = max(num_even, num)
                k_even += 1
        elif num%2:
            if num > num_31_odd_1:
                num_31_odd_2 = num_31_odd_1
                num_31_odd_1 = num
            elif num > num_31_odd_2:
                num_31_odd_2 = num 
            k_31_odd += 1
        else:
            if num > num_31_even_1:
                num_31_even_2 = num_31_even_1
                num_31_even_1 = num
            elif num > num_31_even_2:
                num_31_even_2 = num 
            k_31_even += 1
 
sum_odd = sum(sorted((num_odd, num_31_odd_1, num_31_odd_2))[1:])
sum_even = sum(sorted((num_even, num_31_even_1, num_31_even_2))[1:])
cnt = (k_odd + k_31_odd - 1)*k_31_odd + (k_even + k_31_even - 1)*k_31_even
print(cnt, max(sum_odd, sum_even))

from itertools import combinations as cmb
res = [x + y for x, y in cmb(map(int, open('17.txt').readlines()), 2) if not (abs(x - y) % 2 or x % 31 and y % 31)]
print(len(res))
print(max(res))

cnt = k_odd*k_31_odd + (k_31_odd - 1)*k_31_odd//2 + \
      k_even*k_31_even + (k_31_even - 1)*k_31_even//2