Morris' Blog 我在摸索我存在的意義、生命的意義、內涵及價值。 沒有真正神手般的厲害，套用模式是大部分的情況。 沒那麼厲害的我存在的意義為何，襯托出神的存在？ 有時候不能怪別人不給我機會，我已經開始省思， 為什麼我沒有給別人機會的原因，而是全要靠別人。 有時候我認為我與眾不同，那是因為我做不到很多人能辦到的事情。 我聽不到、玩不起來、說得不夠明瞭、理解得不夠多， 我缺少的語文能力，就是一切一切無能的判定。 為此，我不曉得能做些對別人有幫助的事情。 我知道我自己是個不服輸的性格，但現在我不得不認輸。 我沒有能力，因此我必須認輸。 我到底來這世界作什麼？處處充滿不確定性， 而我總是無法出類拔萃，想當個一般人似乎也回不了頭了。 越強大的挑戰需要越強力的支持，而我的支持是什麼？ 最近嘗試與人談話，但都是失敗的結局，沒人能接受， 不需要言語，不需要溝通，一切所想都表達不出來。

49愛的鼓勵 6訂閱站台

Problem D: Factovisors

The factorial function, n! is defined thus for n a non-negative integer:

   0! = 1
   n! = n * (n-1)!   (n > 0)

We say that a divides b if there exists an integer k such that

   k*a = b

The input to your program consists of several lines, each containing two non-negative integers, n and m, both less than 2^31. For each input line, output a line stating whether or not m divides n!, in the format shown below.

Sample Input

6 9
6 27
20 10000
20 100000
1000 1009

Output for Sample Input

9 divides 6!
27 does not divide 6!
10000 divides 20!
100000 does not divide 20!
1009 does not divide 1000!

#include <stdio.h>
int p[5500], pt = 0;
void sieve() {
    int i, j, m[50000] = {};
    for(i = 2; i < 50000; i++)
        if(m[i] == 0) {
            p[pt++] = i;
            for(j = i+i; j < 50000; j += i)
                m[j] = 1;
        }
}
int solve(int n, int m) {
    int i;
    for(i = 0; i < pt && p[i]*p[i] <= m; i++) {
        if(m%p[i] == 0) {
            int cnt = 0;
            while(m%p[i] == 0)
                cnt++, m /= p[i];
            long long tmp = p[i], cnt2 = 0;
            while(tmp <= n) {
                cnt2 += n/tmp;
                tmp *= p[i];
            }
            if(cnt2 < cnt)  return 0;
        }
    }
    if(m != 1) {
        if(n < m)   return 0;
    }
    return 1;
}
int main() {
    sieve();
    int n, m;
    while(scanf("%d %d", &n, &m) == 2) {
        if(solve(n, m))
            printf("%d divides %d!\n", m, n);
        else
            printf("%d does not divide %d!\n", m, n);
    }
    return 0;
}