整除问题（九度教程第 55 题）

整除问题（九度教程第 55 题）

时间限制：1 秒 内存限制：32 兆 特殊判题：否

1.题目描述：

给定 n， a 求最大的 k，使 n！可以被 a^k 整除但不能被 a^(k+1)整除。
输入：
两个整数 n(2<=n<=1000)， a(2<=a<=1000)
输出：
一个整数.
样例输入：
6 10
样例输出：
1
来源：
2011 年上海交通大学计算机研究生机试真题

2.基本思路

对n!和a进行质因数分解，利用数组cntn存储n!的分解结果，cntn[i]表示素因子i的指数，利用数组cnta存储a的分解结果，cnta[i]表示素因子i的指数.值得注意的是a的素因子必定也是n!的素因子，因此只需要对cnta[i]中素数遍历一遍，然后求得min{cntn[i]/cnta[i]}便是所求解的答案。

3.代码实现

#include <iostream>
#define N 1001

using namespace std;

bool mark[N];//false表示未标记，是素数
int prime[N];
int primeSize=0;
int cntn[N],cnta[N];//cntn[i]表示分解后素数i的个数

void init(){
    for(int i=0;i<N;i++){
        mark[i]=false;
    }
    for(int i=2;i<N;i++){
        if(mark[i]==true)
            continue;//跳过被标记的元素
        else{
            prime[primeSize++]=i;
            for(int j=i*i;j<N;j++){
                if(j%i==0 && mark[j]==false)
                    mark[j]=true;
            }
        }
    }
}
int main()
{
    int n,a;

    init();
//    for(int i=0;i<primeSize;i++)
//        printf("%d ",prime[i]);

    while(scanf("%d%d",&n,&a)!=EOF){
        for(int i=0;i<N;i++)//清空计数器
            cnta[i]=cntn[i]=0;

        //对a进行素因数分解
        for(int i=0;i<primeSize;i++){//遍历素因数
                while(a%prime[i]==0){//求解num中因素prime[i]的个数
                    cnta[prime[i]]+=1;
                    a/=prime[i];
                }
                if(a==1)
                    break;
        }

        //对n!进行素因数分解
        for(int k=2;k<=n;k++){//遍历n阶乘
            int num = k;//存储k值，进行操作
            for(int i=0;i<primeSize;i++){//遍历素因数
                while(num%prime[i]==0){//求解num中因素prime[i]的个数
                    cntn[prime[i]]+=1;
                    num/=prime[i];
                }
                if(num==1)
                    break;
            }
        }
        int min=10000;
        for(int i=2;i<N;i++){
//            printf("cntn[%d]=%d cnta[%d]=%d\n",i,cntn[i],i,cnta[i]);
            if(cnta[i]!=0&&cntn[i]/cnta[i]<min){

                min=cntn[i]/cnta[i];
            }
        }
        printf("%d\n",min);
    }
    return 0;
}