BZOJ1053

最新推荐文章于 2020-02-16 14:11:41 发布

Le_ballon_rouge

最新推荐文章于 2020-02-16 14:11:41 发布

阅读量1.4k

点赞数

分类专栏：搜索深搜数学文章标签：数学约数个数定理质因数分解搜索

本文链接：https://blog.csdn.net/Le_ballon_rouge/article/details/47909363

版权

数学同时被 3 个专栏收录

7 篇文章 0 订阅

订阅专栏

搜索

4 篇文章 0 订阅

订阅专栏

深搜

2 篇文章 0 订阅

订阅专栏

传送门：BZOJ1053

我喜欢数学题，我讨厌结论题……

约数个数定理：
$记\Phi (p)为p不平凡的约数的个数，对\forall q \in \mathbb N^{*} ,q>1,作Lagrange分解q=\prod p^{a},则\Phi (a)=\prod (a+1)$

对于这道题，我们会注意到a应是随p递减的，且2*3*5*7*11*-…*23*29>2000000000，于是我们只要直接搜索这11个质数在分解式中的指数即可。

代码上的小细节见下。

最近眼睛稍微好受一点点了……

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;

int prime[11]={1,2,3,5,7,11,13,17,19,23,29};
int n;
int ans,ans_cnt;
long long m;

void Dfs(int where,long long num,int size,int cnt)
{
    if(where==11){
        if(cnt>ans_cnt&&num>ans)
            ans=num,ans_cnt=cnt;
        if(ans>=num&&cnt>=ans_cnt)
            ans=num,ans_cnt=cnt;
        return;
    }
    int t=1;
    for(int i=0;i<=size;i++){
        if(num*t>=n)
            break;
        Dfs(where+1,num*t,i,cnt*(i+1));
        t*=prime[where];
    }
}

void Readdata()
{
    freopen("loli.in","r",stdin);
    scanf("%d",&n);
}

void Solve()
{
    long long a=1;
    Dfs(0,a,30,1);
    printf("%d\n",ans);
}

void Close()
{
    fclose(stdin);
    fclose(stdout);
}

int main()
{
    Readdata();
    Solve();
    Close();
    return 0;
}