#include <iostream>
#include <cstdio>
#include <cmath>

using namespace std;

unsigned long long A,B;
long long numere;
int n,k,P[80000],desc[80000];
bool X[1000005];

void functie(int in,int sf,int cate,int orig,unsigned long long div)
{
    if(cate==0)
    {
        if(orig%2==0)
        {
            numere=numere+A/div;
        }
        else numere=numere-A/div;
    }
    else
    {
        int i;
        for(i=in;i<=sf-cate+1;i++)
        {
            functie(i+1,sf,cate-1,orig,div*desc[i]);
        }
    }
}
int main()
{
    int i,j;
    FILE *f=fopen("pinex.in","r");
    FILE *g=fopen("pinex.out","w");
    fscanf(f,"%d",&n);
    int m=1000001;
    for(i=2;i<=m;i++)
    {
        if(X[i]==false)
        {
            for(j=2;j<=m/i;j++)
            {
                X[i*j]=true;
            }
        }
    }
    for(i=2;i<=m;i++)
    {
        if(!X[i])
            k++,P[k]=i;
    }
    for(i=1;i<=n;i++)
    {
        fscanf(f,"%I64d%I64d",&A,&B);
        numere=A;
        j=1;
        desc[0]=0;
        while(j<=k || B!=1)
        {
            if(B%P[j]==0)
            {
                desc[0]++;
                desc[desc[0]]=P[j];
                while(B%P[j]==0)
                    B/=P[j];
            }
            j++;
        }
        if(B!=1)
        {
            desc[0]++;
            desc[desc[0]]=B;
        }
        for(j=1;j<=desc[0];j++)
            functie(1,desc[0],j,j,1);
        fprintf(g,"%I64d\n",numere);
    }
    return 0;
}