九度OJ 1104 整除问题 （求阶乘的质因数个数）

题目描述：

给定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;
}