from math import sqrt def div(n): keep = set() for i in range(1,int(sqrt(n))+1): if n % i == 0: keep.add(i) keep.add(n/i) return sorted(list(keep))[::-1] top = 10**6 a = 1 final = [] while a < top: flag = True if str(a) != str(a)[::-1]: q = div(a)[1:] if len(q) > 1: r = str(q[0]) s = str(q[1]) if len(r) > 1 and len(s) > 1 and r == r[::-1] and s == s[::-1]: for x in q[2:]: if str(x) == str(x)[::-1] and len(str(x)) > 1: flag = False break else: flag = False else: flag = False else: flag = False if flag: final.append(a) a += 1 with open("A074889b.txt","w") as g: for y,z in enumerate(final): g.write(str(y+1)+" "+str(z)+"\n")