Как ускорить обработку обработку данных в Пайтон?!?!

f=open("input.txt","r")
a=f.read().split()
c1=int(a[0])
c2=int(a[1])
mas=[]
for i in range(c1):
    mas.append([])
    for j in range(c2):
        mas[i].append((i+1)*(j+1))
black=0
for i in range(c1):
    for j in range(c2):
        black=black+1
red=[]
for i in range(c1):
    for j in range(c2):
        if mas[i][j]%2==0:
            red.append(mas[i][j])
green=[]
for i in range(c1):
    for j in range(c2):
        if mas[i][j]%3==0:
            green.append(mas[i][j])
for i in range(len(green)):
    if green[i] in red:
        b=red.index(green[i])
        red.pop(b)
blue=[]
for i in range(c1):
    for j in range(c2):
        if mas[i][j]%5==0:
            blue.append(mas[i][j])
for i in range(len(blue)):
    if blue[i] in red:
        b1=red.index(blue[i])
        red.pop(b1)
for i in range(len(blue)):
    if blue[i] in green:
         b2=green.index(blue[i])
         green.pop(b2)
red_k=len(red)
blue_k=len(blue)
green_k=len(green)
black=black-red_k-green_k-blue_k
with open("output.txt","w") as f:
    f.write("RED:"+str(red_k))
    f.write("\n")
    f.write("GREEN:"+str(green_k))
    f.write("\n")
    f.write("BLUE:"+str(blue_k))
    f.write("\n")
    f.write("BLACK:"+str(black))

def coloring(rank):
    black, red, green, blue = rank * rank, 0, 0, 0
    for row in range(1,rank+1):
        for column in range(1,rank+1):
            a = row * column
            if a % 5 == 0: blue += 1
            elif (a % 5 != 0)  and (a % 3 == 0):green += 1
            elif (a % 5 != 0)  and (a % 3 != 0) and (a % 2 == 0):red += 1
    black = black - red - green - blue
    return red, green, blue, black

for rank in [10,100,1000]:
    start = time.time()    
    print(
        'RED:{}\nGREEN:{}\nBLUE:{}\nBLACK:{}'.format(   
        *coloring(rank))
    )
    stop = time.time() 
    print("ELAPSED %s sec FOR RANK %s " % (stop-start, rank)) 
    print('*' * 25)