Cod sursa(job #3322176)

Utilizator Costy2345Costi Dimian Costy2345 Data 13 noiembrie 2025 01:09:17
Problema Principiul includerii si excluderii Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.68 kb
#include <bits/stdc++.h>
#define int long long
using namespace std;

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

const int VMAX = 1e6 + 2;

int a, b, m;
int ciur[VMAX];

signed main()
{
    vector<int> primes;
    ciur[1] = 1;
    for (int i = 2; i < VMAX; i++)
    {
        if (ciur[i] == 0)
        {
            primes.push_back(i);
            for (int j = i * i; j < VMAX; j += i)
            {
                ciur[j] = 1;
            }
        }
    }
    fin >> m;
    while (m--)
    {
        fin >> a >> b;
        vector<int> fact;
        int idx = 0;
        for (auto x : primes)
        {
            if(x * x > b)break;
            if (b % x == 0)
            {
                // cout << x << " ";
                fact.push_back(x);
                while(b % x == 0){
                    b /= x;
                } 
                idx++;
            }
        }
        if(b > 1){fact.push_back(b);idx++;}
        
        fact.resize(idx);
        // for(auto x: fact)
        // {
        //     cout << x << " ";
        // }
        cout << "\n";
        int ans = 0;
        for (int i = 1; i <= (1 << idx) - 1; i++)
        {
            int d = 1;
            for (int j = 0; j < idx; j++)
            {
                if (((i >> j) & 1) == 1)
                {
                    d *= fact[j];
                }
            }
            int bits = __builtin_popcount(i);
            if (bits % 2 != 0)
            {
                ans += a / d;
            }
            else
            {
                ans -= a / d;
            }
        }
        fout << a - ans << "\n";
    }
    return 0;
}