【问题标题】:Find the first Triangular number which has 50 factors?找到第一个有 50 个因数的三角形数?
【发布时间】:2013-07-09 06:24:28
【问题描述】:

-----请求的代码修改--------

问题:计算具有 50 个因数的快速三角形序列数?

详细说明:假设有一个系列

   1 : 1
   3 : 1+2
   6 : 1+2+3 
   10 : 1+2+3+4
   15 : 1+2+3+4+5
   21 : 1+2+3+4+5+6
   28 : 1+2+3+4+5+6+7

这里的 1,3,6,10,15,21,28 是三角形系列下的数字。

让我们看看数字的因素

    Number factors         Count
    1     : 1               1              
    3     : 1,3             2
    6     : 1,2,3,6         4
    10    : 1,2,5,10        4
    15    : 1,3,5,15        4
    21    : 1,3,7,21        4
    28    : 1,2,4,7,14,28   6

这里 6 是第一个有 4 个因数的三角形数。 即使 10、15、21 也有 4 个因子,但它们不是第一个。 像这样让我们取一个数字为 2,它有 2 个因子为 1 和 2 与数字 3 相同,也有 2 个因子作为 1 和 3

但根据问题 3 将不是 2 的答案,因为 2 不在三角形序列号列表中,即使它比 3 快。

【问题讨论】:

  • 这个问题似乎是题外话,因为它是关于数学,而不是编程。
  • 是的 JNF,但我需要的是针对上述问题的适当解决方案(步骤/逻辑)。我得到了输出,但我希望可以编写更好的代码。这就是为什么要求一个比我给定逻辑更好的解决方案的原因。
  • @razlebe,Stack Overflow 上的算法问题都可以。
  • AFAIR,当您回答有关 Project Euler 的问题时,您将自动获得访问论坛主题的权限,以便向全世界分享您的解决方案并阅读其他人。这些帖子是更好逻辑的重要来源。

标签: c# algorithm number-theory


【解决方案1】:

三角形数 #2591 = 3357936 是第一个具有正好 50 个因子的三角形: 1、2、3、4、6、8、9、12、16、18、24、27、36、48、54、72、81、108、144、162、216、324、432、648、1296、 2591,5182,7773,10364,15546,20728,23319,31092,41456,46638,62184,69957,93276,124368,139914,186552,209871,279828,373104,419742,559656,839484,1119312,1678968,3357936

三角形数 #12375 = 76576500 是第一个具有至少 500 个因子(实际上是 576 个因子)的三角形:1、2、3、4、5、6、7、9、10, 11, ..., 19144125, 25525500, 38288250, 76576500

三角形数 #1569375 = 1231469730000 是第一个具有正好 500 个因子的数

解决方案代码本身很简单,只要你可以得到除数:

   public static long Solution(int factorsCount) {
      for (long i = 1; ; ++i) {
        long n = i * (i + 1) / 2;

        IList<long> factors = GetDivisors(n);

        // This code tests if a triangle number has exactly factorsCount factors
        // if you want to find out a triangle number which has at least factorsCount factors
        // change "==" comparison to ">=" one:
        // if (factors.Count >= factorsCount)  
        if (factors.Count == factorsCount) 
          return n;
      }
    }

  ...

  long solution = Solution(50);

如果你没有获取数字因子的例程,你可以使用这个:

// Get prime divisors 
private static IList<long> CoreGetPrimeDivisors(long value, IList<int> primes) {
  List<long> results = new List<long>();

  int v = 0;
  long threshould = (long) (Math.Sqrt(value) + 1);

  for (int i = 0; i < primes.Count; ++i) {
    v = primes[i];

    if (v > threshould)
      break;

    if ((value % v) != 0)
      continue;

    while ((value % v) == 0) {
      value = value / v;

      results.Add(v);
    }

    threshould = (long) (Math.Sqrt(value) + 1);
  }

  if (value > 1)
    results.Add(value);

  return results;
}

/// <summary>
/// Get prime divisors 
/// </summary>
public static IList<long> GetPrimeDivisors(long value, IList<int> primes) {
  if (!Object.ReferenceEquals(null, primes))
    return CoreGetPrimeDivisors(value, primes);

  List<long> results = new List<long>();

  while ((value % 2) == 0) {
    results.Add(2);

    value = value / 2;
  }

  while ((value % 3) == 0) {
    results.Add(3);

    value = value / 3;
  }

  while ((value % 5) == 0) {
    results.Add(5);

    value = value / 5;
  }

  while ((value % 7) == 0) {
    results.Add(7);

    value = value / 7;
  }

  int v = 0;
  long n = (long) (Math.Sqrt(value) / 6.0 + 1);
  long threshould = (long) (Math.Sqrt(value) + 1);

  for (int i = 2; i <= n; ++i) {
    v = 6 * i - 1;

    if ((value % v) == 0) {
      while ((value % v) == 0) {
        results.Add(v);

        value = value / v;
      }

      threshould = (long) (Math.Sqrt(value) + 1);
    }

    v = 6 * i + 1;

    if ((value % v) == 0) {
      while ((value % v) == 0) {
        results.Add(v);

        value = value / v;
      }

      threshould = (long) (Math.Sqrt(value) + 1);
    }

    if (v > threshould)
      break;
  }

  if (value > 1) {
    if (results.Count <= 0)
      results.Add(value);
    else if (value != results[results.Count - 1])
      results.Add(value);
  }

  return results;
}

/// <summary>
/// Get all divisors
/// </summary>
public static IList<long> GetDivisors(long value, IList<int> primes) {
  HashSet<long> hs = new HashSet<long>();

  IList<long> divisors = GetPrimeDivisors(value, primes);

  ulong n = (ulong) 1;
  n = n << divisors.Count;

  for (ulong i = 1; i < n; ++i) {
    ulong v = i;
    long p = 1;

    for (int j = 0; j < divisors.Count; ++j) {
      if ((v % 2) != 0)
        p *= divisors[j];

      v = v / 2;
    }

    hs.Add(p);
  }

  List<long> result = new List<long>();

  result.Add(1);

  var en = hs.GetEnumerator();

  while (en.MoveNext())
    result.Add(en.Current);

  result.Sort();

  return result;
}

/// <summary>
/// Get all divisors
/// </summary>
public static IList<long> GetDivisors(long value) {
  return GetDivisors(value, null);
}

【讨论】:

  • 咦,我是不是做错了什么?使用我的 50 算法,我得到#224 is 25200: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 24, 25, 28, 30, 35, 36, 40, 42, 45, 48, 50, 56, 60, 63, 70, 72, 75, 80, 84, 90, 100, 105, 112, 120, 126, 140, 144, 150, 168, 175, 180, 200, 210, 225, 240, 252, 280, 300, 315, 336, 350, 360, 400, 420, 450, 504, 525, 560, 600, 630, 700, 720, 840, 900, 1008, 1050, 1200, 1260, 1400, 1575, 1680, 1800, 2100, 2520, 2800, 3150, 3600, 4200, 5040, 6300, 8400, 12600, 25200。同样的算法产生#12375 is 76576500 500...
  • 等等...我们说的是正好 50 还是超过 50?因为 Project Euler 问题指出“第一个三角形数有超过五百个除数”
  • 只需将比较从“==”更改为“>=”(请参阅​​我的 cmets),您会发现至少有 50 或 500 个因子的三角形数。所以,三角形数#2591 = 3357936 有 50 个因子,#1569375 = 1231469730000 有 500 个因子,#224 = 25200 至少有 50 个因子(实际上是 90),#12375 = 76576500 有至少 500 个因子(实际上是 576)
【解决方案2】:

解决方案: 让我将问题分解为多个模块。

1) 找到三角级数直到一个数。

2) 将所有已识别的数字存储在整数列表中

3) 找出特定数的因子个数

4) 循环遍历三角形数列的每一项,找出每个数的因数。

5) 检查第一个计数为 50 的值,然后显示该值

6) 编写 break 语句以仅显示前 50 个数字。

程序

using System;
using System.Collections.Generic;
using System.Collections;
using System.Linq;
using System.Text;

namespace IsNumberTringularSeriesConsoleApp
{ 
    class Program
    {
        /// <summary>
        /// Listing all numbers comes under Triangular series.
        /// </summary>
        /// <param name="number"></param>
        /// <returns></returns>
        static List<int> GetTriangularNumbers(int number)
        {
            List<int> lstTriangularNumbers = new List<int>();
            int i;
            int sum = 0;
            int triangularNumber = 0;
            for (i = 1; i < number; i++)
            {
                sum = sum + i;
                triangularNumber = sum;
                lstTriangularNumbers.Add(triangularNumber);
            }
            return lstTriangularNumbers;
        }

        /// <summary>
        /// returns(count) the number of factors for each number
        /// </summary>
        /// <param name="number"></param>
        /// <returns></returns>
        public static int FactorCount(int number)
        {
            List<int> factors = new List<int>();
            int max = (int)Math.Sqrt(number);  //round down
            for (int factor = 1; factor <= max; ++factor)
            { 
                //test from 1 to the square root, or the int below it, inclusive.
                if (number % factor == 0)
                {
                    factors.Add(factor);
                    if (factor != number / factor)
                   {
                     // Don't add the square root twice!  
                        factors.Add(number / factor);
                   }
                }
            }
            return factors.Count;
        }

        static void Main(string[] args)
        {
            List<int> lstTriangularNumbers = new List<int>();
            List<int> factors = new List<int>();
            int count = 0;
            //Getting the list of numbers comes under triangular series till 5000
            lstTriangularNumbers = GetTriangularNumbers(5000);

            foreach (int number in lstTriangularNumbers)
            {
                /*
                 * Calling the FactorCount(number) function to check no of factors 
                 * available for the specific triangular number - number.
                 */
                 count = FactorCount(number);
                 //Console.WriteLine("No of factors for : " + number + " is : " + count);
                if (count == 50)
                {
                    Console.WriteLine("No of factors for first Triangular Number : " + number + " is : " + count);
                    break;
                }
            }
            Console.ReadLine();
        }
    }
}

【讨论】:

  • @HighPerformanceMark:OP 大概可以替代BigInteger,尽管可能需要一些工作
  • 老大,它给出的输出是“第一个三角形数的因子数:3357936 是 50”,这是正确的。我需要更好的代码。
  • 首先要“变得更好”:您的FactorCountGetFactors 基本上具有相同的主体。所以只需将FactorCount 设为return GetFactors(number).Count;。再一次,看看mathblog,看看其他显着减少执行时间的方法。或者你所说的“更好”的代码到底是什么意思?
  • 我会把GetTriangularNumbers改成这样的:public static IEnumerable&lt;int&gt; GetTriangularInt32() { var sum = 0; var i = 1; while (Int32.MaxValue - sum &gt;= i) { sum += i++; yield return sum; } }因为你foreachbreak不管怎样,你不必事先创建不必要的数字。
【解决方案3】:

这是我的答案

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;

namespace TriangularSeries
{
    class MyClass
    {
        static void Main(string[] args)
        {
            int result;
            TriangularSeries aSeries = new TriangularSeries();
            result = aSeries.TSeries();
            Console.WriteLine("The first Triangular Series number that has 50Factors is : " + result);
            Console.Read();
        }
    }

    //Find the Triangular Series numbers
    class TriangularSeries
    {
        public int TSeries()
        {
            int fCount = 0, T1 = 1, i = 1, T2 = 0, fval = 0;
            while (fCount != 50)
            {
                i += 1;
                T2 = T1 + i;

                fCount = CalcFactors(T1);
                fval = T1;                   
                T1 = T2;

            }
            return fval;
        }

        public int CalcFactors(int num1)
        {

            List<int> factors = new List<int>();
            int max = (int)Math.Sqrt(num1);  //round down
            for (int factor = 1; factor <= max; ++factor)
            {
                //test from 1 to the square root, or the int below it, inclusive.
                if (num1 % factor == 0)
                {
                    factors.Add(factor);
                    if (factor != num1 / factor)
                    {
                        // Don't add the square root twice!  
                        factors.Add(num1 / factor);
                    }
                }
            }
            return factors.Count;

        }
    }   
}

【讨论】:

  • 感谢您的解决方案闪亮。
【解决方案4】:

这是我的 C 语言程序

#include<stdio.h>
int i;
int num1=0,num2=1;
int a[3000];
int tri_series()         //This function finds the Triangular series numbers
{
    for(i=0;num2<=3000;i++)
    {
    num1=num1+num2;
    a[i]=num1;
    num2++;
    }
}
int main()
{
tri_series();            //Calling the function tri_series
int num,count;
    for(i=0;i<=3000;i++)
    {
      count=0;
      for(num=1;num<=a[i];num++)
      {
        if(a[i]%num==0)  //Finds the factors of each Triangular Series Number 
        count=count+1;   
      }
      if(count==50)      //Break at the first Triangular Series Number having 50 factors
      {
       printf("%d:%d\t",a[i],count);
       break;
      }
    }
}

性能问题: 此代码在执行时间方面会产生性能问题。执行和显示输出需要“一分钟”时间。

【讨论】:

  • 谢谢阿努。感谢c版的回答。将对 C 程序员有所帮助。
  • 但是需要做一些修改。 num2 未定义,int main 不返回任何值。请修改程序,以便对 C 程序员有所帮助。
  • a[i] 值在 main 函数中打印,tri_series 函数只是在 main 函数中定义和调用。感谢陈丹的指正。
猜你喜欢
  • 1970-01-01
  • 1970-01-01
  • 2020-12-10
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 2018-03-17
  • 1970-01-01
相关资源
最近更新 更多