Cod sursa(job #1081410)

Utilizator Iustin_BulimarFMI Iustin Bulimar Iustin_Bulimar Data 13 ianuarie 2014 16:54:31
Problema Principiul includerii si excluderii Scor 60
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, x, y, m, f[20], bk[20], s, k, p[100000];
bool c[1000000];

void ciur()
{
    long long 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])
            p[++k]=i;
}

void fact(long long n)
{
    long long d, i;
    m=0;
    for(i=1; p[i]<=n && i<=k; i++)
        if(n%p[i]==0)
            f[++m]=p[i];
}

void aduna(long long k, long long x)
{
    long long i, p;
    p=1;
    for(i=1; i<=k; i++)
        p*=f[bk[i]];
    if(k & 1)
        s+=x/p;
    else
        s-=x/p;
}

long long pinex(long long x)
{
    long long i, k;
    k=1;
    i=1;
    s=0;
    while(k)
        if(k>m || i>m)
        {
            k--;
            i=bk[k]+1;
        }
        else
        {
            bk[k]=i;
            aduna(k, x);
            k++;
            i++;
        }
        return x-s;
}

int main()
{
    cin>>q;
    ciur();
    while(q--)
    {
        cin>>x>>y;
        fact(y);
        cout<<pinex(x)<<'\n';
    }
    return 0;
}