Pagini recente » Cod sursa (job #1879385) | Cod sursa (job #3197620) | Cod sursa (job #580859) | Cod sursa (job #488598) | Cod sursa (job #2874747)
#include <bits/stdc++.h>
#define ll long long
#define mx 1000001
using namespace std;
ifstream f("pinex.in");
ofstream g("pinex.out");
ll a[25],sol[25],prim[100005];
bool ciur[1000005];
ll A,B,t,k,sum,nrp;
void desc(ll nr)
{
int j=1;
k=0;
memset(a,0,sizeof(a));
while(prim[j]*prim[j]<=nr&&nr>1)
{
if(nr%prim[j]==0)
{
k++;
a[k]=prim[j];
while(nr%prim[j]==0) nr=nr/prim[j];
}
j++;
}
if(nr>1)
{
k++;
a[k]=nr;
}
}
void verif(int nr)
{
ll prod=1;
for(int i=1; i<=nr; i++)
{
prod*=a[sol[i]];
}
if(nr%2==0) sum+=A/prod;
else sum-=A/prod;
}
void submultimi(int n)
{
int i,j;
for(i=1; i<(1<<n); i++)
{
int p=0;
for(j=0; j<n; j++)
if( i & (1<<j)) sol[++p]=j+1;
verif(p);
}
}
int main()
{
int i,j;
f>>t;
for(i=2; i<=mx; i++)
{
if(!ciur[i])
{
prim[++nrp]=i;
for(j=i+i; j<=mx; j=j+i) ciur[j]=1;
}
}
while(t)
{
f>>A>>B;
desc(B);
sum=A;
submultimi(k);
g<<sum<<'\n';
t--;
}
return 0;
}