divmod 3 == 0

from random import randint

print('Введите размер массива\n')
N = int(input())
mas = [int(randint(1, 1000)) for i in range(N)]
newMas = [int(i) for i in mas if i % 3 == 0]

print('Начальный массив:\n', mas)
print('Массив из чисел которые делятся на 3 без остатка:\n', newMas)