#Python program to generate the b-file from sympy import isprime def omega(n): s=0 for i in range(1,n+1): if n%i==0 and isprime(i)==True: s+=1 return s def spf(n): #spf is the smallest prime factor of n if n==1: return 1 for i in range(1,n+1): if n%i==0 and isprime(i)==True: return i def bigomega(n): if n==1: return 0 return bigomega(n/spf(n))+1 i=1 j=1 while j<=1000: if bigomega(i)-omega(i)==1 and bigomega(i+1)-omega(i+1)==1: print str(j)+" "+str(i) j+=1 i+=1