-
题目描述:
-
给定n,a求最大的k,使n!可以被a^k整除但不能被a^(k+1)整除。
-
输入:
-
两个整数n(2<=n<=1000),a(2<=a<=1000)
-
输出:
-
一个整数.
-
样例输入:
-
6 10
-
样例输出:
-
1
身体欠佳....分析不想写了,最重要的一点就是n!的质因数的求法,要转为n的质因数 的求法,详见机试指南P82,
直接上代码:
#include <stdio.h>
int Prime[1001];
int PrimeSize;
bool mark[1001];
int cnt1[1001];//存储所有素数中,n的质因数的指数
int cnt2[1001];//存储所有素数中,a的质因数的指数
void Init(){
int i;
PrimeSize=0;
for (i=1;i<=1000;i++) mark[i]=false;
for(i=2;i<=1000;i++){
if(mark[i]==true) continue;
Prime[PrimeSize++]=i;
int j;
for(j=i*i;j<=1000;j+=i)
{
mark[j]=true;
}
}
}
int main(int argc, char** argv) {
int n,a;
Init();
while(scanf("%d%d",&n,&a)==2){
int i;
//初始化
for(i=0;i<=PrimeSize;i++)
{
cnt1[i]=cnt2[i]=0;
}
//统计n!的该质因数的指数
for(i=0;i<PrimeSize;i++)
{
int t=n;
while(t){
cnt1[i]+=t/Prime[i];
t=t/Prime[i];//将n/(p^k)转换为n/p/p/p.....
}
}
int ans =12312312;//为遍历最小做准备
//统计a的质因数的指数
for(i=0;i<PrimeSize;i++)
{
while(a%Prime[i]==0)
{
cnt2[i]++;
a/=Prime[i];
}
if(cnt2[i]==0) continue;//如果该素数不能整除a,则跳过,不影响a对n的整除性
//计算同一个质因数n对应的指数是a对应指数的多少倍,取最小的即为所求
if(cnt1[i]/cnt2[i]<ans)
{
ans=cnt1[i]/cnt2[i];
}
}
printf("%d\n",ans);
}
return 0;
}