Pagini recente » Cod sursa (job #1570457) | Cod sursa (job #1182085) | Cod sursa (job #362491) | Cod sursa (job #702570) | Cod sursa (job #2279341)
#include <iostream>
#include <fstream>
#define NMAX 1000001
using namespace std;
ifstream f("pinex.in");
ofstream g("pinex.out");
long long t,a,b,nr,v[101],nrc,rez,p,prim[100000],nrp;
bool ciur[NMAX];
int i,j;
int main()
{
for(i=2;i<=NMAX;i++)
{
if(ciur[i]==0)
{
nrp++;
prim[nrp]=i;
for(j=2;j<=NMAX/i;j++)
{
ciur[i*j]=1;
}
}
}
f>>t;
while(t--)
{
f>>a>>b;
nr=0;
for(i=1;prim[i]*prim[i]<=b;i++)
{
if(b%prim[i]==0)
{
nr++;
v[nr]=prim[i];
while(b%prim[i]==0)
{
b=b/prim[i];
}
}
}
if(b>1)
{
nr++;
v[nr]=b;
}
rez=a;
for(i=1;i<(1LL<<nr);i++)
{
nrc=0;
p=1;
for(j=0;j<nr;j++)
{
if(((1LL<<j)&i)!=0)
{
nrc++;
p=p*v[j+1];
}
}
if(nrc%2==1)
{
rez=rez-a/p;
}
else
{
rez=rez+a/p;
}
}
g<<rez<<"\n";
}
return 0;
}