Pagini recente » Cod sursa (job #554126) | Cod sursa (job #2817813) | Cod sursa (job #897220) | Cod sursa (job #2352030) | Cod sursa (job #755673)
Cod sursa(job #755673)
#include <fstream>
using namespace std;
#define ll long long
#define max_p 1000005
ll i, j, n, B, M, k, A;
int fprim[7330000], fact[6330000];
bool prim[1000005];
ifstream in("pinex.in");
ofstream out("pinex.out");
void sieve()
{
int i, j;
for(i = 2; i < max_p; i++)
{
if(!prim[i])
{
fprim[++k] = i;
for(j = i + i; j < max_p; j += i) prim[j] = true;
}
}
}
void solve()
{
ll ret, nr = 0, prod, i = 0, ok = 1, sol = 0;
while(B > 1)
{
i++;
if(B % fprim[i] == 0)
{
fact[++nr] = fprim[i];
while(B % fprim[i] == 0) B /= fprim[i];
}
if(fprim[i] * fprim[i] > B && B > 1)
{
fact[++nr] = B;
B = 1;
}
}
for(int i = 1; i < (1 << nr); i++)
{
prod = 1, ok = 1;
for(int j = 0; j < nr; j++)
{
if(i & (1 << j))
{
prod *= fact[j + 1];
ok++;
}
}
ret = A / prod;
if(nr & 1) sol += ret;
else sol -= ret;
}
out << A - sol << "\n";
}
int main()
{
int i;
sieve();
in >> M;
for(i = 0; i < M; i++)
{
in >> A >> B;
solve();
}
return 0;
}