欧拉计划第12题

高度可约的三角形数

三角形数数列是通过逐个加上自然数来生成的。例如，第7个三角形数是 1 + 2 + 3 + 4 + 5 + 6 + 7 = 28。三角形数数列的前十项分别是：

1, 3, 6, 10, 15, 21, 28, 36, 45, 55, …

让我们列举出前七个三角形数的所有约数：

 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个约数的三角形数是多少？

 

C代码演示

#include<iostream>
#include<cmath>
using namespace std;
#define MAX_N 10000
int f[MAX_N +5],cnt[MAX_N+5];
int prime[MAX_N +5];

int main(){
    f[1]=1;
    for(int i=2;2*i<=MAX_N;i++){
        if(!prime[i]){
            prime[++prime [0]]=i;
            f[i]=2;
            cnt[i]=1;
        }
        for(int j=1;prime[j]*i<=MAX_N;j++){
            prime[prime[j]*i]=1;
            if(i%prime[j]==0) {
                f[prime[j]*i]=f[i]/(cnt[i]+1)*(cnt[i]+2);
                cnt[prime[j]*i]=cnt[i]+1;
                break;
            }else{
                f[prime[j]*i]=f[prime[j]]*f[i];
                cnt[prime[j]*i]=1;
            }
        }
    }
    for(int i=2;i<=100;i++){
        printf("F[%d]=%d\n",i,f[i]);
    }

    return 0;
}