Cod sursa(job #2794617)

Utilizator db_123Balaban David db_123 Data 5 noiembrie 2021 10:51:35
Problema Divizori Primi Scor 45
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 2.5 kb
#include <fstream>
#include <bitset>

using namespace std;

ifstream cin("divprim.in");
ofstream cout("divprim.out");

const int DIM=1000000;
const int DIM_MATRIX=380000;

int    v[DIM+1];
bitset<DIM+1> ciur;
int    ma[8][DIM_MATRIX];
int    t,n,k;

int main()
{
    int left,right,mid,sol,j;
    for(int i=0 ; i < 8;i++)
        ma[i][0]=1; /// counter of numbers

    v[2]=1;
    for(int i = 4;i<=DIM;i+=2)
    {
        v[i]++;
        ciur[i]=1;
    }
    for(int i=3;i<=DIM;i+=2)
    {
        if(ciur[i]==0)
        {
            v[i]++;
            ciur[i]=1;

            for(int j=2*i; j<=DIM; j+=i)
            {
                v[j]++;
                ciur[j]=1;
            }
        }
    }
    v[1]=1;
    for(int i = 1;i<=DIM;i++)
    {
        if(v[i]<=7)
            ma[v[i]][ma[v[i]][0]++]=i;
    }

//    for(int i=2; i<=DIM; i++)
//    {
//        if(v[i]==0)
//        {
//            for(int j=i; j<=DIM; j+=i)
//            {
//                v[j]++;
//            }
//        }
//        if(v[i]<=7)
//            ma[v[i]][ma[v[i]][0]++]=i;
//    }

//    for(int i=0; i<=1000; i++)
//    {
//        cout <<  "( "<< v[i] << ") ";
//    }
//    return 0;

/*

    for(int i=0; i<=100; i++)
    {
        cout <<  ma[1][i] << " ";
    }
    cout << endl;
    for(int i=0; i<=100; i++)
    {
        cout <<  ma[2][i] << " ";
    }
    cout << endl;
    for(int i=0; i<=100; i++)
    {
        cout <<  ma[3][i] << " ";
    }
    cout << endl;
    for(int i=0; i<=100; i++)
    {
        cout <<  ma[4][i] << " ";
    }
    cout << endl;
    for(int i=0; i<=100; i++)
    {
        cout <<  ma[5][i] << " ";
    }
    cout << endl;
    for(int i=0; i<=100; i++)
    {
        cout <<  ma[6][i] << " ";
    }
    cout << endl;
    for(int i=0; i<=100; i++)
    {
        cout <<  ma[7][i] << " ";
    }
    return 0;
*/
    cin >>t;
    for(int i=1; i<=t; i++)
    {
        cin >>n>>k;
        left = 1;
        right =  ma[k][0]-1;
        sol=0;
        if(k==0)
            right=1;
        while(left<=right)
        {
            mid = (left+ right)/2;
            if(ma[k][mid]==n)
            {
                sol =  ma[k][mid];
                break;
            }
            else if(ma[k][mid]<n)
            {
                sol = ma[k][mid];
                left = mid+1;
            }
            else
                right = mid -1;
        }
        cout <<sol<<endl;
    }
    return 0;
}