Pagini recente » Cod sursa (job #2451195) | Cod sursa (job #1892615) | Cod sursa (job #2373220) | Cod sursa (job #2329837) | Cod sursa (job #1081410)
#include <fstream>
using namespace std;
ifstream cin("pinex.in");
ofstream cout("pinex.out");
long long q, x, y, m, f[20], bk[20], s, k, p[100000];
bool c[1000000];
void ciur()
{
long long n, r, i, j;
n=1000000;
r=1000;
c[2]=1;
for(i=3; i<=n; i+=2) c[i]=1;
for(i=3; i<=n; i+=2)
if(c[i] && i<r)
for(j=i*i; j<=n; j+=i)
c[j]=0;
for(i=2; i<n; i++)
if(c[i])
p[++k]=i;
}
void fact(long long n)
{
long long d, i;
m=0;
for(i=1; p[i]<=n && i<=k; i++)
if(n%p[i]==0)
f[++m]=p[i];
}
void aduna(long long k, long long x)
{
long long i, p;
p=1;
for(i=1; i<=k; i++)
p*=f[bk[i]];
if(k & 1)
s+=x/p;
else
s-=x/p;
}
long long pinex(long long x)
{
long long i, k;
k=1;
i=1;
s=0;
while(k)
if(k>m || i>m)
{
k--;
i=bk[k]+1;
}
else
{
bk[k]=i;
aduna(k, x);
k++;
i++;
}
return x-s;
}
int main()
{
cin>>q;
ciur();
while(q--)
{
cin>>x>>y;
fact(y);
cout<<pinex(x)<<'\n';
}
return 0;
}