P1221 最多因子数题解

被回文题号吸引过来的

整体数据范围：$ 2\le L\le U\le10^9 $。

完全可以爆搜。

PS：其实就是找因子数最多的数。根据一个数都能被多个质因子分解，例如 $p=a{1}^{p1}+a{2}^{p2}+\dots \dots +a{n}^{pn.ai}$ 表示质数，$pi$ 表示该数字内有 $pi$ 个 $ai$，他的因子数 $sum$ 就是 $sum=(p1+1)\ast (p2+1)\ast \dots \dots \ast (pn+1)$。除了131074，他有个质因子大于100，其他所有数的因子都小于100，所以只需要特判131074，其他都直接用dfs质因子跑质因子数目就行了。

Code

#include<bits/stdc++.h>
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;
}