PE 110 Diophantine reciprocals II （数论：分式个数）（加强版）（dfs）

PE 110题目
PE 108 题解

从PE 108 我们可以知道，分式 $\dfrac{1}{a}+\dfrac{1}{b}=\dfrac{1}{n}$ 的个数 就是$\dfrac{1+d(n^2)}{2}$，其中$d(n)$是n的约数的个数。、

而 PE 110 只是PE 108的加强版。用PE 108的方法 换个方式就可以了。
比如：

$n=2^{p_1}3^{p_2}5^{p_3}......$

$d(n^2)=N=\prod_{i=1}(2p_i+1)$

所以我们可以知道样例是怎么来的。

样例：

$1260=2^2 3^2 5^1 7^1$

$d(n^2)=N=5*5*3*3=225$

综合 $\dfrac{1+d(n^2)}{2}$，所以答案就是$\dfrac{225+1}{2}=113$。

我们要求的，只是倒过来求 n <script type="math/tex" id="MathJax-Element-10">n</script> 而已。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll primes[16]={2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47,53};
ll ans;
void dfs(ll down, ll upper, ll build, ll cost)
{
    if(cost > ans)return;
    if(build > 7999999) //4000000 * 2 -1
    {
        ans=min(ans, cost);
        return;
    }
    if(down>=16)return;

    for(ll i = upper;i>=1;i--)
    {
        ll tmp=1;
        for(int j = 1 ; j <= i ; j++)
        {
            tmp*=primes[down]; 
        }

        if((double) tmp*cost > (double)ans) continue;

        dfs(down + 1, i , build * ( 2 * i + 1 ), tmp * cost);
    }
}

int main()
{
    ans=1e18;
    dfs(0, 20, 1, 1);
    printf("%lld\n", ans);
    return 0;
}