Question

我试图解决项目Euler问题＃2 ：查找斐波纳契系列偶数项的总和达4,000,000。 Ex：高达1000，总和是2 + 8 + 34 + 144 + 610 = 798。

现在，我的算法是从2开始添加每第三个斐波那契数 - 因为每三个斐波纳契数都必然遵循逻辑：

偶：E，奇数：O，

O+E=O    E+O=O    O+O=E
1+2=3    2+3=5    5+3=8    (example)

所以，我写了下面的代码来找出答案..

#include<stdio.h>
#define LT 4000000
int main()
{
    double i0,i1,sum=0,cycle,eSum=2,status=1;
    i0 = 1;
    i1 = 2;

    while(i1<LT && status == 1)
    {
        for(cycle=3;cycle>0;cycle--)
        {
            sum=i0+i1;
            i0=i1;
            i1=sum;

            if((i1+i0)>LT)
            {
                status = 0;
                break;
            }
        }
        eSum+=(status == 1)?sum:0;
    }
    printf("\nThe required Answer: %8.0f\n",eSum);
    return 0;
}

现在，正常用于 LT = 1000 ，但对于需要LT = 4,000,000的问题，程序显示错误的值 1089154 < / strong>，而不是 4613732 的正确值。

我无法弄清楚这段代码有什么问题。另外，我不明白它如何能够正常工作LT = 1000，但不是更多的数字..我错过了一些令人尴尬的明显的东西吗？请帮助..

Answer 1

OP过早退出循环。将测试更改为sum而不是i1+i0 next sum。

    for(cycle=3;cycle>0;cycle--) {
        sum=i0+i1;
        i0=i1;
        i1=sum;
        // if((i1+i0)>LT)
        if(sum>LT)
        {
            status = 0;
            break;
        }
    }

由于您的号码在[1 ... 2 * 4,000,000]范围内，因此您可以使用long，unsigned long，float，double获得可接受的结果。由于这是一个整数问题，建议使用整数。

Answer 2

内心看起来很奇怪，你不应该使用double来表示自然数。尝试这样的事情：

#include<stdio.h>
#define LT 4000000
int main()
{
  long i0,i1,sum=0,eSum=2;
  int counter = 0;
  i0 = 1;
  i1 = 1;

  while(sum<LT)
  {
    sum=i0+i1;
    i0=i1;
    i1=sum;

    eSum+=((counter++ % 3 == 0) && ( sum <= LT))?sum:0;
  }
  printf("\nThe required Answer: %8.0f\n",eSum);
  return 0;
}

Answer 3

为什么不简单呢？

#include <stdio.h>
#include <inttypes.h>
#define LT 4000000
typedef uint64_t Numtype;
int main( void )
{
    Numtype n0 = 1;
    Numtype n1 = 1;
    Numtype n2 = 0;
    Numtype sum = 0;
    for(;;)
    {
        n2 = n0 + n1;
        if ( n2 > LT ) break;
        if ( n2 % 2 == 0 ) sum += n2;
        n0 = n1;
        n1 = n2;
    }

    printf("\nThe required Answer: %ld\n", sum);
    return 0;
}

这个算法有什么问题？

3 个答案: