Cod sursa(job #3143438)

Utilizator ezluciPirtac Eduard ezluci Data 30 iulie 2023 11:03:24
Problema Principiul includerii si excluderii Scor 60
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.35 kb
using namespace std;
#ifdef EZ
   #include "./ez/ez.h"
   const string FILE_NAME = "test";
#else
   #include <bits/stdc++.h>
   const string FILE_NAME = "pinex";
#endif
#define mp make_pair
#define mt make_tuple
#define ll long long
#define pb push_back
#define eb emplace_back
#define fi first
#define se second
#define cin fin
#define cout fout
ifstream fin (FILE_NAME + ".in");
ofstream fout (FILE_NAME + ".out");

bool prime[1'000'001];
vector<int> primes;

void genPrimes()
{
   memset(prime, 1, 1'000'001);
   prime[0] = prime[1] = 0;

   for (int i = 2; i*i <= 1e6; ++i)
      if (prime[i])
         for (int j = i*i; j <= 1e6; j += i)
            prime[j] = 0;
   
   for (int x = 2; x <= 1e6; ++x)
      if (prime[x])
         primes.pb(x);
}

int main()
{
   genPrimes();

   int m;
   cin >> m;
   while (m--)
   {
      ll a, b;
      cin >> a >> b;

      vector<int> divi;
      
      for (int d : primes)
         if (d > b)
            break;
         else if (b % d == 0)
            divi.pb(d);
      
      ll ans = a;
      for (int x = 1; x < (1 << divi.size()); ++x)
      {
         ll nr = 1;
         for (int j = 0; j < divi.size(); ++j)
            if (x & 1<<j)
               nr *= divi[j];
         
         if (__builtin_popcount(x) % 2 == 1)
            ans -= a / nr;
         else
            ans += a / nr;
      }

      cout << ans << '\n';
   }
}