NOJ 1166 数字实力 筛素数+因子个数求和

数字实力

时间限制(普通/Java) :  1000 MS/ 3000 MS          运行内存限制 : 65536 KByte
总提交 : 406            测试通过 : 77 

比赛描述

每一个数字都有自己的实力，我们可以从不同的角度来描述数字的实力。我们做出如下定义：一个数字n的实力是它所拥有的因数个数。例如：6的实力是4 ，因为它有4个因数（1，2，3，6）。 现在我们需要知道从1开始到某一个数n中，哪个数字拥有最强的实力。

输入

一个正整数n ( 0< n < 500000)。

输出

两个数字m和s。m为从1到n中具有最强实力且最小的数字。s为m的实力。

样例输入

8

样例输出

6 4

提示

1到8中,6与8都具有最大的实力4。然而你应该输出较小的6。

题目来源

计算机学院/软件学院第二届ACM程序设计大赛

    题意：给出n，求出1到n中因子最多的数极其因子的个数。

    分析：数据很大，暴力肯定超时。首先得知道求取一个数因子个数的方法，在此给出结论，若一个数存在质因子a1、a2…ak，且a1^b1*a2*b2…ak^bk=n，则其因子数为(b1+1)*(b2+1)*…*(bk+1)，具体解释可参见http://blog.sina.com.cn/s/blog_818d3d9301017436.html ，叙述非常易懂。所以该题需要转移到求取质因子，如果直接遍历判断质数的话又会超时，在此需要筛选出质数数组，在质数数组的情况下进行遍历判断可大大节省时间。所以该题需要用到筛素数的算法，版本非常多，在此不赘述。见AC代码：

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
using namespace std;
const int maxn=500005;
bool is[maxn];
int prm[maxn];
int getprm(int n)//筛素数 
{
	int i, j, k = 0;
	int s, e = (int)(sqrt(0.0 + n) + 1);
	memset(is, 1, sizeof(is));
	prm[k++] = 2;
	is[0] = is[1] = 0;
	for (i = 4; i < n; i += 2) is[i] = 0;
	for (i = 3; i < e; i += 2)
		if (is[i])
		{
			prm[k++] = i;
			for (s = i * 2, j = i * i; j < n; j += s)
				is[j] = 0;
		}
	for ( ; i < n; i += 2)
		if (is[i])
			prm[k++] = i;
	return k;
}
int count(int n) //求一个数的因子个数 
{
	int cnt=0,res=1;
	for(int i=0; prm[i]*prm[i]<=n; i++)
		if(n%prm[i]==0)
		{
			cnt=0;//注意cnt每次需要重置为0   
			while(n%prm[i]==0)
			{
				cnt++;
				n/=prm[i];
			}
			res*=(cnt+1);
		}
	if(n!=1)//如果最后剩下非1质因子   如7之类的  
		res*=2;
	return res;
}
int main()
{
	int n;
	getprm(maxn);
	while(~scanf("%d",&n))
	{
		int res=0,num;
		for(int i=1; i<=n; i++)
			if(count(i)>res)
			{
				res=count(i);
				num=i;
			}
		printf("%d %d\n",num,res);
	}
}

    特记下，以备后日回顾。