Cod sursa(job #3218950)

Utilizator Theo20067Cismaru Theodor-Alexe Theo20067 Data 28 martie 2024 21:45:01
Problema Principiul includerii si excluderii Scor 50
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.36 kb
#include <fstream>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cmath>
using namespace std;
ifstream fin ("pinex.in");
ofstream fout("pinex.out");
long long t,F[1100002],sol,k,P,A,B;
long long v[20],D[20];
vector <long long> prime;
void precalc()
{
    F[0]=2;
    F[1]=2;
    for(int i=2;i<=1100000;i++)
    {
        if(!F[i])
        {
            F[i]=1;
            prime.push_back(i);
            for(int j=i+i;j<=1100000;j+=i)
                F[j]=2;
        }
    }
}
void backtrack(int pas,int ant)
{
    if(ant+1<=k)
    {
        for(int i=ant+1;i<=k;i++)
        {
            v[pas]=D[i];
            P*=D[i];
            backtrack(pas+1,i);
            if(pas%2)
                sol+=A/P;
            else
                sol-=A/P;
            P/=D[i];
        }
    }
}
void solve()
{
    fin>>A>>B;
    int ind=0;
    k=0;
    while(B!=1&&ind<prime.size())
    {
        if(B%prime[ind]==0)
        {
            D[++k]=prime[ind];
            while(B%prime[ind]==0)
                B/=prime[ind];
        }
        ind++;
        if(prime[ind]>sqrt(B))
        {
            D[++k]=B;
            break;
        }
    }
    P=1;
    sol=0;
    backtrack(1,0);
    fout<<A-sol<<"\n";
}
int main()
{
    precalc();
    fin>>t;
    while(t--)
        solve();
    return 0;
}