Cod sursa(job #755681)

Utilizator visanrVisan Radu visanr Data 6 iunie 2012 18:29:49
Problema Principiul includerii si excluderii Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.58 kb
#include <fstream>
#include <iostream>
using namespace std;


#define ll long long
#define max_p 1000005

ll i, j, n, B, M, k, A;
int fprim[7330000], fact[6330000];
bool prim[1000005];

ifstream in("pinex.in");
ofstream out("pinex.out");

void sieve()
{
     int i, j;
     for(i = 2; i < max_p; i++)
     {
             if(!prim[i])
             {
                        fprim[++k] = i;
                        for(j = i + i; j < max_p; j += i) prim[j] = true;
             }
     }
}

void solve()
{
     ll prod, sol = 0;
     int i = 0, nr = 0, ok = 1, j;
     while(B > 1)
     {
             i++;
             if(B % fprim[i] == 0)
             {
                  fact[++nr] = fprim[i];
                  while(B % fprim[i] == 0) B /= fprim[i];
             }
             if(fprim[i] * fprim[i] > B && B > 1)
             {
                         fact[++nr] = B;
                         B = 1;
             }
     }
     for(i = 1; i < (1 << nr); i++)
     {
             prod = 1, ok = 1;
             for(j = 0; j < nr; j++)
             {
                     if(i & (1 << j))
                     {
                          prod *= fact[j + 1];
                          ok++;
                     }
             }
             if(ok & 1) sol -= A / prod;
             else sol += A / prod;
     }
     out << A - sol << "\n";
}

int main()
{
    int i;
    sieve();
    in >> M;
    for(i = 0; i < M; i++)
    {
          in >> A >> B;
          solve();
    }
   // cin >> i;
    return 0;
}