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

最新推荐文章于 2024-06-17 01:08:30 发布

原创最新推荐文章于 2024-06-17 01:08:30 发布 · 773 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#PE #110 #分式个数 #数论 #dfs

Project Euler 同时被 3 个专栏收录

81 篇文章

订阅专栏

ACM_搜索

56 篇文章

订阅专栏

ACM_其他数论

48 篇文章

订阅专栏

本文解析了PE110问题的解决方案，该问题是PE108的加强版，主要探讨了分式1/a+1/b=1/n的个数与n的约数之间的关系，并给出了一种递归算法来寻找满足条件的最小n值。

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$ 。

我们要求的，只是倒过来求而已。

#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;
}