Cod sursa(job #1640038)

Utilizator GinguIonutGinguIonut GinguIonut Data 8 martie 2016 15:26:24
Problema Principiul includerii si excluderii Scor 60
Compilator cpp Status done
Runda Arhiva educationala Marime 1.4 kb
#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;
}