【九度】题目1104：整除问题

时间限制：1 秒

内存限制：32 兆

特殊判题：否

提交：5282

解决：1768

题目描述：

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

输入：

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

输出：

一个整数.

样例输入：

6 10

样例输出：

1

来源： 2011年上海交通大学计算机研究生机试真
n!可以分解为多个质因数相乘的形式。设P1 P2..PN为不同的质数。 则n！=p1^e1p2^e2...pn^en 同理将b分解为b= p1^E1p2^E2...pn^En 如若n!能整除b那么，n！中的p一定包含了b中的p，并且，所有的指数都不小于b中的。 当n！能整除b的k次方并要求k为最小值，只需要求得指数的差中最小的那个。 所以重点就是将两个数字分解质因数。 b的质因数好分解，前面已经做过类似的。重点来看n！的分解： 对于质因子p来说，从1到n中，有质因子p的数字有：n/p个，每个数字有一个质因子p，一共有n/p个。（从p到p*（n/p）） 对于质因子p*p来说，从1到n中，有质因子p*p的数字有：n/p*p个，每个数字有两个质因子p，但因为前面的n/p个里面包含了其中的一个，所以一共有p因子n/p*p个。 以此类推。 n/p*p*p。。 一共有p因子n/p*p*p个。 所以对于n！求解质数因子来说，只需要对每个质数，进行类似的循环，即可得到质因数以及相应的指数。

#include "stdafx.h"
#include <string.h>
bool mark[1010];
int prime[1010];
int primesize;
void init(int x){
	for (int i = 0; i < x; i++)
	{
		mark[i]=false; //初始化标记数组
	}
	primesize=0;
	for(int i=2;i<=x;i++){
		if (mark[i]==true) //如果被标记了，则不是素数，执行下一次循环
		continue;
		prime[primesize++]=i;//将素数存入数组里面
		for (int j = i*i; j <=x; j+=i)//如果这个数字是素数，那么将其倍数的所有数字都标记为非负数
		{
			mark[j]=true;
		}
	}
}
int cnt1[1001];
int cnt2[1001];
int _tmain(int argc, _TCHAR* argv[])
{
	int n,a;
	while (scanf("%d%d",&n,&a))
	{
		init(n);
		for (int i = 0; i < primesize; i++)
		{
			printf("%d\n",prime[i]);
		}
	for(int i=0;i<primesize;i++){
	int t=n;
	while (t)  //将n!分解质因数，并使以该质因数为下标的数组单位计数
	{ cnt1[i]+=t/prime[i]; //每整除以一个质数，代表N!可以分解t/prime[i]个
	t/=prime[i];
	}
	}
      int ans=123123123;
	  for (int i = 0; i < primesize; i++)
	 { while (a%prime[i]==0)
	  {cnt2[i]++;
	   a/=prime[i];
	  }
	  if(cnt2[i]==0||cnt1[i]==0) continue; //当他为0；
	   if(cnt1[i]/cnt2[i]<ans)
		    ans=cnt1[i]/cnt2[i];           //质因数的幂的差最小的一个数
	  }
	  printf("%d",ans);
	}
	return 0;
}

未完成 。。明天继续补充。。