#include <bits/stdc++.h>
using namespace std;

const int maxn = 1e4 + 1;

int prime[maxn];

void getPrime(){
    memset(prime, 0, sizeof(prime));
    for(int i = 2; i < maxn; i++){
        if(0 == prime[i]){
            prime[++prime[0]] = i;
        }
        for(int j = 1; j <= prime[0] && i * prime[j] < maxn; j++){
            prime[i * prime[j]] = 1;
            if(0 == i % prime[j]){
                break;
            }
        }
    }
    return;
}

int main(){
    int n, a;
    while(cin>>n>>a){
        getPrime();
        int n_fac[1 + prime[0]] = {0};
        for(int i = 2; i <= n; i++){
            int x = i;
            for(int j = 1; j <= prime[0]; j++){
                while(0 == x % prime[j]){
                    n_fac[j]++;
                    x /= prime[j];
                }
            }
        }
        int a_div[1 + prime[0]] = {0};
        for(int i = 1; i <= prime[0]; i++){
   ...