Cod sursa(job #1081295)

Utilizator Iustin_BulimarFMI Iustin Bulimar Iustin_Bulimar Data 13 ianuarie 2014 15:05:42
Problema Principiul includerii si excluderii Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.27 kb
#include <fstream>
using namespace std;
ifstream cin("pinex.in");
ofstream cout("pinex.out");

long long q, a, p, m, f[25], b[25], s, sgn, pr, prim[10000];
bool c[10000001];

void ciur()
{
    int 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])
            prim[++pr]=i;
}

void fact(long long n)
{
    long long d;
    for(d=1; prim[d]<=n; d++)
    {
        if(n%prim[d]==0)
            f[++m]=prim[d];
    }
}

void aduna(long long k, long long n)
{
    long long i, p=1;
    for(i=1; i<=k; i++)
        p*=b[i];
    s+=sgn*(n/p);
    sgn*=-1;
}

bool cond(long long k)
{
    if(b[k-1]>=b[k])
        return 0;
    return 1;
}

void pinex(long long k, long long n)
{
    long long i;
    for(i=1; i<=m; i++)
    {
        b[k]=f[i];
        if(cond(k))
        {
            aduna(k, n);
            pinex(k+1, n);
        }
    }
}

int main()
{
    cin>>q;
    ciur();
    while(q--)
    {
        cin>>a>>p;
        m=0;
        fact(p);
        s=0;
        sgn=1;
        pinex(1, a);
        cout<<a-s<<'\n';
    }
    return 0;
}