输出1000~9999之间的素数对，如（1031,1301）

num=1000
a=list()
b=list()
while num<=9999:
    j=2
    counter=0
    while j<=num/2:
        if num%j==0:
           counter=counter+1
           break
        j=j+1
    if counter==0:
        #a是素数的列表
        a.append(num)
    num=num+1
Len=len(a)
i=0
while i<Len:
    str1=str(int(a[i]/1000))
    str2=str(int((a[i]%1000)/100))
    str3=str(int((a[i]%100)/10))
    str4=str(int(a[i]%10))
    CgNum=int(str4+str3+str2+str1) #求反过来的数,并判断是否为素数
    k=2
    count=0
    while k<=CgNum/2:
        if CgNum%k==0:
            count=count+1
            break
        if count==0:
            if CgNum in a:
                b.extend([[a[i],CgNum]])
                a.remove(CgNum)
                Len=Len-1 #每移走一个数，长度减一，防止超出范围
                break
        k=k+1
    i=i+1
print(b) 
print(f'共有{len(b)}个素数对')