记录编号 141613 评测结果 AAAAAAAAAA
题目名称 最多因子数 最终得分 100
用户昵称 Gravatar席一鸣 是否通过 通过
代码语言 C++ 运行时间 0.077 s
提交时间 2014-12-03 17:39:59 内存使用 0.33 MiB
显示代码纯文本
#include<cmath>
#include<cstdio>
#include<iostream>
using namespace std;
int a,b,c[3500]={2},m,s,u=1,w;
bool h(int n)
{
	int i,t;
	if(n==1)
		return 0;
	t=sqrt(n);
	for(i=0;i<u&&c[i]<=t;i++)
		if(!(n%c[i]))
			return 0;
	return 1;
}
void d(int p,int n,int y)
{
	int e,f,i,l=(a-1)/n,r=b/n,t;
	if(n>=a)
		if(m<y)
		{
			m=y;
			w=n;
		}
		else
			if(m==y&&w>n)
				w=n;
	if(l==r)
		return;
	if(n<s)
		for(i=l+1;i<=r;i++)
			if(h(i))
			{
				t=(y<<1);
				e=n*i;
				if (m<t)
				{
					m=t;
					w=e;
				}
				else
					if(m==t&&w>e)
						w=e;
				break;
			}
	if(m*1.0>y*pow(2.0,log(b*1.0/n)/log(c[p]*1.0)))
		return;
	if(n*c[p]>b||p>=u)
		return;
	for(i=2,f=c[p];;i++,f*=c[p])
	{
		if(f>r)
			break;
		d(p+1,n*f,y*i);
		if(f*1.0>r*1.0/c[p])
			break;
	}
	d(p+1,n,y);
}
main()
{
	freopen("divisors.in","r",stdin);
	freopen("divisors.out","w",stdout);
	int i;
	cin>>a>>b;
	if(b==1)
	{
		cout<<"Between 1 and 1, 1 has a maximum of 1 divisors.";
		return 0;
	}
	s=sqrt(b);
	for(i=3;i<=s;i++)
		if(h(i))
			c[u++]=i;
	d(0,1,1);
	cout<<"Between "<<a<<" and "<<b<<", "<<w<<" has a maximum of "<<m<<" divisors.";
}