欧拉工程第12题第一个拥有超过500个约数的三角形数是多少

最新推荐文章于 2023-07-20 10:43:15 发布

Stone Chen

最新推荐文章于 2023-07-20 10:43:15 发布

阅读量1.2k

点赞数 1

分类专栏：我的欧拉工程之路文章标签：解题算法欧拉工程我的欧拉工程之路

本文链接：https://blog.csdn.net/chen595572/article/details/52830770

版权

我的欧拉工程之路专栏收录该内容

31 篇文章 2 订阅

订阅专栏

题目

三角形数序列是由对自然数的连加构造而成的。所以第七个三角形数是 $1+2+3+4+5+6+7=28$ . 那么三角形数序列中的前十个是：
$1, 3, 6, 10, 15, 21, 28, 36, 45, 55, \dots$ $1, 3, 6, 10, 15, 21, 28, 36, 45, 55, \ldots$
下面我们列出前七个三角形数的约数：
1: 1
3: 1,3
6: 1,2,3,6
10: 1,2,5,10
15: 1,3,5,15
21: 1,3,7,21
28: 1,2,4,7,14,28
可以看出28是第一个拥有超过5个约数的三角形数。
那么第一个拥有超过500个约数的三角形数是多少？

解题方法

第n个三角形数实质上就是等差数列{ $a_{n}=n$ }的前n项和 $S_{n}$ ：

S n = n ( n + 1 ) 2

$S_{n}=\frac {n(n+1)}{2}$

任何一个正整数都可以表示成几个素数的次方的乘积

假设 $P_{n}$ 表示第n个素数，那么任意正整数可以通过下面的式子获得：

N u m = P k 1 1 P k 2 2 P k 3 3 \dots P k n n, n \in N + ， k n \in N

$Num=P_{1}^{k_{1}}P_{2}^{k_{2}}P_{3}^{k_{3}} \cdots P_{n}^{k_{n}}, n\in N^{+}，k_{n}\in N$

而一个正整数拥有的约数个数为：

d i v i s o r = (k 1 + 1) (k 2 + 1) \dots (k n + 1)

$divisor=(k_{1}+1)(k_{2}+1)\cdots (k_{n}+1)$

程序

程序中用到的Prime类是我自定义的工具类，因为在做欧拉工程的题目遇到很多素数相关的题目，所以我实现了一个Prime工具类，方便解题。具体的源码以及用法参考我的另一篇文章——《Java工具类素数类》

public static void solve() {
    int[] primeArray = Prime.getPrimeArrayByLength(2000);
    int num = 0;
    for (int i = 1;; i++) {
        // Sn的值
        num += i;
        int temp = num;
        // Sn拥有的约数个数
        int divisor = 1;
        int j = 0;
        for (; temp != 1; j++) {
            int count = 0;
            // 求kn的值
            while (temp % primeArray[j] == 0) {
                temp /= primeArray[j];
                count++;
            }
            divisor *= count + 1;
        }

        if (divisor > 500) {
            System.out.println("num " + num + " is first num whose divisors over 500(" + divisor + "）");
            break;
        }
    }
}