九度OJ 1104 整除问题

题目地址：http://ac.jobdu.com/problem.php?pid=1104

题目描述：

给定n，a求最大的k，使n！可以被a^k整除但不能被a^(k+1)整除。

输入：

两个整数n(2<=n<=1000)，a(2<=a<=1000)

输出：

一个整数.

样例输入：

6 10

样例输出：

1

来源：

2011年上海交通大学计算机研究生机试真题

#include <stdio.h>
#include <string.h>
#include <math.h>

int IsPrime (int n){
    if (n <= 1)
        return 0;
    int sq = (int)sqrt((double)n);
    while (sq >= 2){
        if (n % sq == 0)
            break;
        --sq;
    }
    return (sq >= 2) ? 0 : 1;
} 

void Initialize(int Prime[], int n, int * primeSize){
    int index = 1;
    int num = 3; 

    Prime[0] = 2;
    while (num < n){
        if (IsPrime (num)){
            Prime[index] = num;
            ++index;
        }
        num += 2;
    }
    *primeSize = index;
}

int main(void){
    int n, a;
    int Prime[200];
    int primeSize;
    int cnt1[200], cnt2[200];
    int tmp;
    int i;

    Initialize(Prime, 1000, &primeSize);
    while (scanf ("%d%d", &n, &a) != EOF){
        memset (cnt1, 0, sizeof(cnt1));
        memset (cnt2, 0, sizeof(cnt2));
        for (i=0; i<primeSize; ++i){
            tmp = n;
            while (tmp){
                cnt1[i] += tmp / Prime[i];
                tmp /= Prime[i];
            }
        }
        int ans = 10000000;
        for (i=0; i<primeSize; ++i){
            while (a % Prime[i] == 0){
                ++cnt2[i];
                a /= Prime[i];
            }
            if (cnt2[i] == 0)
                continue;
            if (cnt1[i] / cnt2[i] < ans)
                ans = cnt1[i] / cnt2[i];
        }
        printf ("%d\n", ans);
    }

    return 0;
}

参考资料：2013年王道论坛计算机考研机试指南