nyoj 56阶乘因式分解（一）

阶乘因式分解（一）

时间限制： 3000 ms  |  内存限制： 65535 KB

难度： 2

描述

给定两个数m,n,其中m是一个素数。

将n（0<=n<=10000）的阶乘分解质因数，求其中有多少个m。

输入

第一行是一个整数s（0<s<=100)，表示测试数据的组数
随后的s行, 每行有两个整数n，m。

输出

输出m的个数。

样例输入

2
100 5
16 2

样例输出

24
15

#include<stdio.h>
int a[10000];
void prim()  //计算10000以内的质数 
{
	int i,j,k=1;
	for(i=2;i<10000;i++) 
	{
		for(j=2;j<i;j++)
			if(i%j==0)
				break;
			if(j==i)
		        a[k++]=i;
	}
}
int f(int n, int m)
{
	int i,k=0;
	for(i=1;n>=a[i];)
		if(n%a[i]==0)
		{
			if(a[i]==m||n==m)
				k++;
			n=n/a[i];
			continue;
		}
		else
			i++;
	return k;
}
int main()
{
	int m,n,s,i,su,count;
	prim();
	scanf("%d",&s);
	while(s--)
	{
		su=count=0;
		scanf("%d%d",&n,&m);
		for(i=2;i<=n;i++)
		{
			su=f(i,m);
			count+=su;
		}
		printf("%d\n",count);
	}
	return 0;
}