Cod sursa(job #2279341)

Utilizator adimiclaus15Miclaus Adrian Stefan adimiclaus15 Data 9 noiembrie 2018 14:20:14
Problema Principiul includerii si excluderii Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.38 kb
#include <iostream>
#include <fstream>
#define NMAX 1000001
using namespace std;
ifstream f("pinex.in");
ofstream g("pinex.out");
long long t,a,b,nr,v[101],nrc,rez,p,prim[100000],nrp;
bool ciur[NMAX];
int i,j;
int main()
{
    for(i=2;i<=NMAX;i++)
    {
        if(ciur[i]==0)
        {
            nrp++;
            prim[nrp]=i;
            for(j=2;j<=NMAX/i;j++)
            {
                ciur[i*j]=1;
            }
        }
    }
    f>>t;
    while(t--)
    {
        f>>a>>b;
        nr=0;
        for(i=1;prim[i]*prim[i]<=b;i++)
        {
            if(b%prim[i]==0)
            {
                nr++;
                v[nr]=prim[i];
                while(b%prim[i]==0)
                {
                    b=b/prim[i];
                }
            }
        }
        if(b>1)
        {
            nr++;
            v[nr]=b;
        }
        rez=a;
        for(i=1;i<(1LL<<nr);i++)
        {
            nrc=0;
            p=1;
            for(j=0;j<nr;j++)
            {
                if(((1LL<<j)&i)!=0)
                {
                    nrc++;
                    p=p*v[j+1];
                }
            }
            if(nrc%2==1)
            {
                rez=rez-a/p;
            }
            else
            {
                rez=rez+a/p;
            }
        }
        g<<rez<<"\n";
    }
    return 0;
}