Pagini recente » Cod sursa (job #882624) | Cod sursa (job #1629920) | Cod sursa (job #2056279) | Cod sursa (job #588068) | Cod sursa (job #1640038)
#include <fstream>
#include <math.h>
#define ll long long
#define MAX_P 1000005
using namespace std;
ifstream fin("pinex.in");
ofstream fout("pinex.out");
bool ciur[MAX_P];
int prime[80000], t, curr;
long long A, B, fact[80], Sol;
void read()
{
fin>>A>>B;
}
void prec()
{
prime[++prime[0]]=2;
for(int i=3;i<MAX_P;i+=2)
{
if(ciur[i]==0)
{
for(int j=i+i+i;j<MAX_P;j+= i << 1)
ciur[j]=1;
prime[++prime[0]]=i;
}
}
}
void solve()
{
fact[0]=curr=0;
while(B>1)
{
curr++;
if(B%prime[curr]==0)
{
fact[++fact[0]]=prime[curr];
while(B%prime[curr]==0)
B/=prime[curr];
if(B>1 && prime[curr]>sqrt(B))
{
fact[++fact[0]]=B;
B=1;
}
}
}
Sol=A;
for(int i=1;i<(1 << fact[0]);i++)
{
ll nr=0, prod=1;
for(int j=0;j<fact[0];j++)
{
if((1 << j) & i)
{
prod=1LL*prod*fact[j+1];
nr++;
}
}
if(nr%2)
nr=-1;
else
nr=1;
Sol=Sol+1LL*nr*A/prod;
}
fout<<Sol<<'\n';
}
int main()
{
prec();
fin>>t;
while(t--)
{
read();
solve();
}
return 0;
}