欧拉计划Problem 18

By starting at the top of the triangle below and moving to adjacent numbers on the row below, the maximum total from top to bottom is 23.

3
7 4
4 6
8 5 9 3

That is, 3 + 7 + 4 + 9 = 23.

Find the maximum total from top to bottom of the triangle below:

75
95 64
17 47 82
18 35 87 10
20 04 82 47 65
19 01 23 75 03 34
88 02 77 73 07 63 67
99 65 04 28 06 16 70 92
41 41 26 56 83 40 80 70 33
41 48 72 33 47 32 37 16 94 29
53 71 44 65 25 43 91 52 97 51 14
70 11 33 28 77 73 17 78 39 68 17 57
91 71 52 38 17 14 91 43 58 50 27 29 48
63 66 04 68 89 53 67 30 73 16 69 87 40 31
04 62 98 27 23 09 70 98 73 93 38 53 60 04 23

NOTE: As there are only 16384 routes, it is possible to solve this problem by trying every route. However, Problem 67, is the same challenge with a triangle containing one-hundred rows; it cannot be solved by brute force, and requires a clever method! ;o)

#define N 15

int main()
{
 clock_t te,ts;
 ts=clock();
 char t[5];
 int s[N][N]={0};
 FILE *f;
 int i,j;
 f = fopen("18.txt","r");
 for (i = 0; i < N; i++)
 {
  for (j = 0; j <= i; j++)
  {
   fgets(t,4,f);
   s[i][j] =atoi(t);
  }
 }
 fclose(f);
 for ( i = N-2; i >=0; i--)
 {
  for ( j = 0; j <= i; j++)
  {
   if (s[i+1][j] > s[i+1][j+1])
   {
    s[i][j]+=s[i+1][j];
   }
   else
   {
    s[i][j]+=s[i+1][j+1];
   }
  }
 }
 printf("\nanswer: %d",s[0][0]);
 te=clock();
 printf("\ntime difference: %ds\n",(te-ts)/CLOCKS_PER_SEC);
 return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,以下是采用Python编写的一个简单的示例程序,用于演示如何使用向前欧拉、向后欧拉和梯形欧拉方法求解初值问题,并计算相应的截断误差和收敛精度: ```python import numpy as np import matplotlib.pyplot as plt # 定义一个函数,表示初值问题的导数函数 def f(t, y): return y - t**2 + 1 # 定义一个函数,表示精确解 def exact(t): return (t+1)**2 - 0.5*np.exp(t) # 定义向前欧拉方法 def euler_forward(y0, t0, h, n): y = [y0] t = [t0] for i in range(n): y.append(y[i] + h*f(t[i], y[i])) t.append(t[i] + h) return y, t # 定义向后欧拉方法 def euler_backward(y0, t0, h, n): y = [y0] t = [t0] for i in range(n): y.append(y[i] + h*f(t[i+1], y[i+1])) t.append(t[i] + h) return y, t # 定义梯形欧拉方法 def trapezoidal(y0, t0, h, n): y = [y0] t = [t0] for i in range(n): y.append(y[i] + 0.5*h*(f(t[i], y[i]) + f(t[i+1], y[i] + h*f(t[i], y[i])))) t.append(t[i] + h) return y, t # 计算截断误差和收敛精度 def calculate_error(y, t, h): n = len(y) - 1 error = [0] for i in range(1, n+1): exact_value = exact(t[i]) error.append(abs(y[i] - exact_value)) convergence_order = np.log10(error[-1]/error[-2])/np.log10(h) return error, convergence_order # 设置初始条件和参数 y0 = 0.5 t0 = 0 h = 0.1 n = 10 # 使用向前欧拉方法求解初值问题 y_forward, t_forward = euler_forward(y0, t0, h, n) error_forward, convergence_order_forward = calculate_error(y_forward, t_forward, h) # 使用向后欧拉方法求解初值问题 y_backward, t_backward = euler_backward(y0, t0, h, n) error_backward, convergence_order_backward = calculate_error(y_backward, t_backward, h) # 使用梯形欧拉方法求解初值问题 y_trapezoidal, t_trapezoidal = trapezoidal(y0, t0, h, n) error_trapezoidal, convergence_order_trapezoidal = calculate_error(y_trapezoidal, t_trapezoidal, h) # 绘制精确解和数值解的图像 t_exact = np.linspace(t0, t0+n*h, 100) y_exact = exact(t_exact) plt.plot(t_exact, y_exact, label='Exact') plt.plot(t_forward, y_forward, 'o-', label='Euler Forward') plt.plot(t_backward, y_backward, 's-', label='Euler Backward') plt.plot(t_trapezoidal, y_trapezoidal, '^-', label='Trapezoidal') plt.legend() plt.xlabel('t') plt.ylabel('y') plt.title('Numerical Solution of Initial Value Problem') # 输出截断误差和收敛精度 print('Euler Forward: Error =', error_forward, ', Convergence Order =', convergence_order_forward) print('Euler Backward: Error =', error_backward, ', Convergence Order =', convergence_order_backward) print('Trapezoidal: Error =', error_trapezoidal, ', Convergence Order =', convergence_order_trapezoidal) plt.show() ``` 这段代码中,我们首先定义了一个初值问题的导数函数 `f` 和精确解函数 `exact`,然后分别实现了向前欧拉、向后欧拉和梯形欧拉方法的函数 `euler_forward`、`euler_backward` 和 `trapezoidal`。接着,我们定义了一个计算截断误差和收敛精度的函数 `calculate_error`,并使用这个函数计算了每种方法的截断误差和收敛精度。最后,我们绘制了精确解和数值解的图像,并输出了截断误差和收敛精度的结果。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值