#include <bits/stdc++.h>

using namespace std;

long long a,b;
char c[100005];
int diviz[100005],k,fprim[100005];

void ciur()
{
    for (int i=2;i<=100005;++i)
    {
        if (!c[i])
            for (int j=i+i;j<=100005;++j)
                c[j]=1;
        diviz[++k]=i;
    }
}
void solve()
{
    memset(fprim,0,sizeof(fprim));
    int total=a;
    for (int i=1;diviz[i]*diviz[i]<=b;++i)
    {
        if (b%diviz[i]==0)
            {
                fprim[++fprim[0]]=diviz[i];
                while (b%diviz[i]==0)
                    b/=diviz[i];
            }
    }
    if (b>1)
        fprim[++fprim[0]]=b;
    for (int i=1;i<(1<<fprim[0]);++i)///asa iau toate posibilitatile de divizori
    {
        int nr=0,t=1;
        for (int j=0;j<fprim[0];++j)
            if (i & (1<<j)) ///daca e in componenta binara a lui i
        {
            ++nr;
            t*=fprim[j+1];
        }
        if (nr%2==1)
            total-=a/t;
        else
            total+=a/t;
    }
    printf("%d\n",total);
}
int main()
{
    freopen("pinex.in","r",stdin);
    freopen("pinex.out","w",stdout);
    int n;
    scanf("%d",&n);
    ciur();
    for (int i=1;i<=n;++i)
    {
        scanf("%lld%lld",&a,&b);
        solve();
    }
    return 0;
}