#include <fstream>
using namespace std;


#define ll long long
#define max_p 1000005

ll i, j, n, B, M, k, A;
int fprim[7330000], fact[6330000];
bool prim[1000005];

ifstream in("pinex.in");
ofstream out("pinex.out");

void sieve()
{
     int i, j;
     for(i = 2; i < max_p; i++)
             if(!prim[i])
                        for(j = i + i; j < max_p; j += i) prim[j] = 1;
     for(j = 2; j < max_p; j++) if(v[j] == 0) fprim[++k] = j;
}

void solve()
{
     ll ret, nr = 0, prod, i = 0, ok = 1, sol = 0;
     while(B > 1)
     {
             i++;
             if(B % fprim[i] == 0)
             {
                  fact[++nr] = fprim[i];
                  while(B % fprim[i] == 0) B /= fprim[i];
             }
             if(fprim[i] * fprim[i] > B && B > 1)
             {
                         fact[++nr] = B;
                         B = 1;
             }
     }
     for(int i = 1; i < (1 << nr); i++)
     {
             prod = 1, ok = 1;
             for(int j = 0; j < nr; j++)
             {
                     if(i & (1 << j))
                     {
                          prod *= fact[j + 1];
                          ok++;
                     }
             }
             ret = A / prod;
             if(ok & 1) sol -= ret;
             else sol += ret;
     }
     out << A - sol << "\n";
}

int main()
{
    int i;
    sieve();
    in >> M;
    for(i = 0; i < M; i++)
    {
          in >> A >> B;
          solve();
    }
    return 0;
}