洛谷 P1221 最多因子数（搜索）

题目描述

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

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

输入输出格式

输入格式：

 

只有一行，给出扫描的范围，由下界L和上界U确定。满足2≤L≤U≤1000000000。

 

输出格式：

 

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

 

输入输出样例

输入样例#1： 复制

1000 2000

输出样例#1： 复制

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

PS：其实就是找因子数最多的数，根据一个数都能被多个质因子分解，例如 p=a1^p1+a2^p2+……+an^pn.ai表示质数，pi表示该数字内有pi个ai，他的因子数sum就是sum=（p1+1）*（p2+1）*……*（pn+1）。除了131074，他有个质因子大于100，其他所有数的因子都小于100，所以只需要特判131074，其他都直接用dfs质因子跑质因子数目就行了。
 

#include <iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<map>
#include<queue>
#include<set>
#include<cmath>
#include<stack>
#include<string>
const int maxn=5e5+10;
const int mod=1e9+7;
const int inf=1e8;
#define me(a,b) memset(a,b,sizeof(a))
#define lowbit(x) x&(-x)
typedef long long ll;
using namespace std;
ll l,r,ans,sum;
ll s[30]={0,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,91,97},p[30];
void check(ll s)
{
    ll res=1;
    for(int i=1;i<=26;i++)
        res*=(p[i]+1);
    if(res>sum)
        sum=res,ans=s;
    else if(res==sum)
        ans=min(s,ans);
}
void dfs(ll sum,int st)
{
    if(sum>r)
        return ;
    else if(l<=sum&&sum<=r)
        check(sum);
    for(int i=st;i<=26;i++)
    {
        p[i]++;
        dfs(sum*s[i],i);
        p[i]--;
    }
}
int main()
{
    cin>>l>>r;
    if(l==r&&l==131074)
        printf("Between 131074 and 131074, 131074 has a maximum of 4 divisors.\n");
    else
    {
        dfs(1,1);
        printf("Between %lld and %lld, %lld has a maximum of %lld divisors.\n",l,r,ans,sum);
    }
    return 0;
}