Pagini recente » Cod sursa (job #1275208) | Cod sursa (job #2174264) | Cod sursa (job #1284684) | Cod sursa (job #1324830) | Cod sursa (job #2174591)
#include <cstdio>
#include <bitset>
using namespace std;
long long n, a, b, d[10005000], k, prim[10005000];
bitset <10005000>c;
void ciur()
{
for(int i=2;i<=10005000;i++)
if(c[i]==0)
{
for(int j=i+i;j<=10005000;j+=i)
c[j]=1;
prim[++prim[0]]=i;
}
}
void diviz()
{
int i=3;
k=0;
/**if(b%2==0)
{
d[k++]=2;
while(b%2==0)
b/=2;
}**/
for(i=1;prim[i]*prim[i]<=b;i++)
{
if(b%prim[i]==0)
{
while(b%prim[i]==0)
b/=prim[i];
d[k++]=prim[i];
}
}
if(b!=1)
d[k++]=b;
}
void zero()
{
for(int j=0;j<k;j++)
{
d[j]=0;
}
k=0;
}
long long subsir()
{
long long m=1<<k;
long long s=0;
for(int i=1;i<m;i++)
{
long long p=1;
long long nr=0;
for(int j=0;j<k;j++)
{
if((i&(1<<j))!=0)
{
p*=d[j];
nr++;
}
}
if(nr%2==0)
s-=(a/p);
else s+=(a/p);
}
return s;
}
int main()
{
freopen("pinex.in", "r", stdin);
freopen("pinex.out", "w", stdout);
scanf("%d", &n);
ciur();
for(int i=1;i<=n;i++)
{
scanf("%d %d", &a, &b);
diviz();
printf("%d", a-subsir());
zero();
printf("\n");
}
return 0;
}