[交流] Euler Project Q12 欧拉工程第十二题已有9人参与

# -*- coding: utf-8 -*-

# 初始化一个质数表
i = 0
primes = range(2, 100)
while i < len(primes):
    primes = [x for x in primes if x == primes[i] or x % primes[i] != 0]
    i = i + 1

def numbersOfFactors(n):
    """计算n的因子数
    使用了sudo说的计算方法
    加上一个质数分解
    """
    factorNumbers = 1
    factorTimes = 0
    i = 0
    while i < len(primes):
        if n % primes[i] == 0:
            factorTimes += 1
            n /= primes[i]
            continue
        else:
            if factorTimes != 0:
                factorNumbers *= factorTimes + 1
                factorTimes = 0
            i = i + 1
            if n == 1:
                break
    return factorNumbers

# 从第三个三角数开始算，只算第奇数个三角数
n = 3
while numbersOfFactors(n*(n+1)/2) < 500:
    n += 2

# 为防止对偶数个三角数的遗漏，有最后这个计算
if numbersOfFactors(n*(n-1)/2) > 500:
    print n*(n-1)/2
else:
    print n*(n+1)/2

#include
#include

int main(int args, char* argv[])
{
        long result=1,trinum=0,curnum=1;
        int n=0,i=0,big=0;
        int stop = 500;                //停止要求的除数个数,改成5可以测试28
        while(n         {
                n = 0;
                trinum = 0;
                // 计算当前三角数
                for(i=1;i<=curnum;i++)
                        trinum += i;

                // 从1开始除,计算除数个数
                for(i=1;i<=trinum;i++)
                {
                        if(trinum%i==0) // 如果整除,n+1
                                n += 1;
                }
               
                if(n>=stop) //如果超过stop个,保存结果跳出循环,这里是500个
                {
                        result = trinum;
                        break;
                }

                //如果没超过stop个,记录当前最大的个数并输出
                //这个if可以不要,因为我运行以为电脑死机了
                //所以加了这个判断,找到更接近stop的除数个数时输出,让我有点盼头
                if(n>big)  
                {
                        big = n;
                        printf("Current number: %i [%d]\n",trinum,big);
                }

                //自加自然数增加1,trinum = sum(1:curnum)
                curnum++;
        }
       
        // 打印结果
        printf("%i\n",result);
       
        system("PAUSE");
        return 0;
        }

% ans =
%       76576500

Originally posted by libralibra at 2011-05-17 05:38:26:
对此题我相当无语,难度在20题以内,我也不知道具体程序跑了多久,反正我吃了个泡面回来,打印出结果了,应该超过10分钟了

c的,matlab对此表示鸭梨很大

[code] #include <stdio.h>
#include <stdlib.h ...

1. 建立素数表prime[N]，N取一个比较大的值
2. 对x=(1+i)*i/2进行质因数分解：
    x=k_1^m_1 * k_2^m_2 ... * k_n^m_n
    其中^表示指数，k_1 ... k_n表示素数
3. 计算x因子总数：
    divisors = (m_1 + 1) * (m_2 + 1) * ... * (m_n + 1)
4. 判断divisors是否超过500，如果没有，跳到第2步，尝试下一个x；如果已经超过500，则输出x

x/=prime[i]