整除问题(九度教程第 55 题)
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;
}