P1221 最多因子数

# 最多因子数

## 题目描述

数学家们喜欢各种类型的有奇怪特性的数。例如，他们认为 $945$ 是一个有趣的数，因为它是第一个所有约数之和大于本身的奇数。

为了帮助他们寻找有趣的数，你将写一个程序扫描一定范围内的数，并确定在此范围内约数个数最多的那个数。不幸的是，这个数和给定的范围的都比较大，用简单的方法寻找可能需要较多的运行时间。所以请确定你的算法能在几秒内完成最大范围内的扫描。

## 输入格式

只有一行，给出扫描的范围，由下界 $L$ 和上界 $U$ 确定。满足 $2\le L\le U\le10^9$。

## 输出格式

对于给定的范围，输出该范围内约数个数 $D$ 最多的数 $P$。若有多个，则输出最小的那个。请输出 $\texttt{Between }L\texttt{ and }U\texttt{, }P\texttt{ has a maximum of }D\texttt{ divisors.}$，其中 $L,U,P,D$ 的含义同前面所述。

## 样例 #1

### 样例输入 #1

```
1000 2000
```

### 样例输出 #1

```
Between 1000 and 2000, 1680 has a maximum of 40 divisors.
```

//god with me
//#pragma GCC optimize(1)
//#pragma GCC optimize(2)
//#pragma GCC optimize(3,"Ofast","inline")
//#include <bits/stdc++.h>
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#define inf 0x7fffffff
//#define ll long long
#define int long long
//#define double long double
//#define double long long
#define re int
//#define void inline void
#define eps 1e-8
//#define mod 1e9+7
#define ls(p) p<<1
#define rs(p) p<<1|1
#define pi acos(-1.0)
#define pb push_back
#define mk make_pair
#define P pair < int , int >
#define pd pair <double,double>
#define x first
#define y second
#define db double
using namespace std;
const int mod=4933;
//const int inf=1e18;
const int M=1e8;
const int N=1e5+5;//??????.???? 4e8

//反素数
int l,r;
int ans,maxn;
int prime[]={2,3,5,7,11,13,17,19,23,29};
int divide(int n)
{
    int num=1;
    for(int i=2;i<=sqrt(n);i++)
    {
        if(n%i==0)
        {
            num++;
            int temp=n/i;
            if(n%temp==0&&temp!=i)num++;         //求反就是跟他相积为n的另一个数，判断是不是平方
        }
    }
    if(n!=1)num++;
    return num;
}

void dfs(int deep,int arr,int cur,int num)  //质数层数  当前质数最大质数   因子个数  当前的质数之积
{
    if(cur>maxn||(cur==maxn&&num<ans))
    {
        maxn=cur,ans=num;
    }
    if(deep>8)return;    //原因是质数单个乘积超范围
    for(int i=1;i<=arr;i++)           //从1开始因为后面可能是0
    {
        num*=prime[deep];
        if(num>r) return;
        dfs(deep+1,i,(i+1)*cur,num);
    }
}
void solve()
{
    cin>>l>>r;
    if(r-l<100000)
    {
        for(int i=l;i<=r;i++)
        {
            int cnt=0;
            cnt=divide(i);
            if(cnt>maxn)
            {
                maxn=cnt;
                ans=i;
            }
        }
        cout<<"Between "<<l<<" and "<<r<<", "<<ans<<" has a maximum of "<<maxn<<" divisors.";

    }
    else
    {
        dfs(0,31,1,1);
        cout<<"Between "<<l<<" and "<<r<<", "<<ans<<" has a maximum of "<<maxn<<" divisors.";
    }
}

signed main()
{
    int T=1;
//    cin>>T;
    for(re index=1;index<=T;index++)
    {
        solve();
    }
    return 0;
}