（eular）最大公约（1552）

#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;

int n;
int eular(int n)
{
int ret=1,i;
for (i=2;i*i<=n;i++)
if (n%i==0)
{
n/=i,ret*=i-1;
while (n%i==0)
n/=i,ret*=i;
}
if (n>1)
ret*=n-1;
return ret;
}

int main()
{
int i,j,k;
int t,m;
cin>>t;
while (t--)
{
int ans=0;
scanf("%d%d",&n,&m);
if (m==1)
{
printf("%d\n",n);
continue;
}
if (n<m)
{
printf("0\n");
continue;
}
k=sqrt((double)n);
for (i=2;i<=k;i++) if (n%i==0)
{
if (i>=m)
ans+=eular(n/i);
if (n/i>m && n/i!=i)
ans+=eular(i);
}
printf("%d\n",ans+1);
}
return 0;
}