fzu 1564 Combination 组合数是否包含因数

对于本题的因数m，只需要将其化素因子的乘积形式，之后只需要判定每一个素因子p在组合数中的个数是否大于或等于在因数m中的个数，若是则输出Yes，若存在某一个不是则输出No。至于如何判断一个素因子在组合数中的个数，如下所示：

//获得素因子p在1*2*··n中的个数
int getNum(int n, int p) {
    int res;
    res = 0;
    while (n) {
        res += n / p;
        n /= p;
    }
    return res;
}

代码如下：

#include<stdio.h>
#include<string.h>
#include<math.h>
#define nmax 46345
int flag[nmax], prime[nmax], pfactor[nmax], cpfactor[nmax];
int plen, cplen;
//打素数表
void mkprime() {
    int i, j;
    memset(flag, -1, sizeof(flag));
    for (i = 2, plen = 0; i < nmax; i++) {
        if (flag[i]) {
            prime[plen++] = i;
        }
        for (j = 0; (j < plen) && (i * prime[j] < nmax); j++) {
            flag[i * prime[j]] = 0;
            if (i % prime[j] == 0) {
                break;
            }
        }
    }
}
//将K化成素因子的乘积的形式
void findpFactor(int k) {
    int i, te, cnt;
    te = (int) sqrt(k * 1.0);
    for (i = 0, cplen = 0; (i < plen) && (prime[i] <= te); i++) {
        if (k % prime[i] == 0) {
            cnt = 0;
            while (k % prime[i] == 0) {
                cnt++;
                k /= prime[i];
            }
            pfactor[cplen] = prime[i];
            cpfactor[cplen++] = cnt;
        }
    }
    if (k > 1) {
        pfactor[cplen] = k;
        cpfactor[cplen++] = 1;
    }
}
//获得素因子p在1*2*··n中的个数
int getNum(int n, int p) {
    int res;
    res = 0;
    while (n) {
        res += n / p;
        n /= p;
    }
    return res;
}
void solve(int n, int m) {
    int i, temp;
    for (i = 0; i < cplen; i++) {
        temp = getNum(n, pfactor[i]) - getNum(m, pfactor[i])
                - getNum(n - m, pfactor[i]);
        if (temp < cpfactor[i]) {
            puts("No");
            return;
        }
    }
    puts("Yes");
}
int main() {
#ifndef ONLINE_JUDGE
    freopen("data.in", "r", stdin);
#endif
    int t, n, m, k;
    mkprime();
    scanf("%d", &t);
    while (t--) {
        scanf("%d %d %d", &n, &m, &k);
        findpFactor(k);
        solve(n, m);
    }
    return 0;
}