Даны натуральные числа p и q. Получить все делители числа q, взаимнопростые с p

x=int(input('x='))
y=int(input('y='))
i=2
k=0
while i<=y:
    if y%i==0 and x%i!=0:
        print (i,end=' ')
        k=1
    i=i+1
if k==0:
    print('таких делителей нет')

def gcd(p, delitel):
    return p if delitel == 0 else gcd(delitel, p % delitel);
 
if gcd(p, delitel) == 1:
    print(delitel);

x=int(input('x='))
y=int(input('y='))
def nod(x,y):
    while x!=0 and y!=0:
        if x>y:
            x=x%y
        else:
            y=y%x 
    nod=x+y            
i=2
k=0
while i<=y:
    if y%i==0 and nod(i,x)==1:
        print (i,end=' ')
        k=1
    i=i+1
if k==0:
    print('таких делителей нет'

p=23;
q = 9;
def gcd(a, b):
    return a if b == 0 else gcd(b, a % b)
 
for delitel in range (1, q+1):
 
    if q%delitel == 0 and  gcd(p, delitel)==1:
        print(delitel);

p=int(input('p='))
q=int(input('q='))
def gcd(a, b):
    return a if b == 0 else gcd(b, a % b)
 
for delitel in range (1, q+1):
 
    if q%delitel == 0 and  gcd(p, delitel)==1:
        print(delitel);

#Ввод данных
p=int(input('p='));
q=int(input('q='));
 
#Алгоритм Евклида для нахождения НОД
def gcd(a, b):
    return a if b == 0 else gcd(b, a % b);
 
#Ищем делители для q
for delitel in range (1, q+1):
 
    if q%delitel == 0 and  gcd(p, delitel)==1: #Если нашли делитель, то проверяем взаимопростой ли он с p
        print(delitel);

 for delitel in range (2, q+1):

  for delitel in range (2, q)