欧拉计划Problem_18

C#

#include <stdio.h>
#define num 105
int arr[num][num] ={{75},
                    {95,64},
                    {17,47,82},
                    {18,35,87,10},
                    {20,4,82,47,65},
                    {19,1,23,75,3,34},
                    {88,2,77,73,7,63,67},
                    {99,65,4,28,6,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,4,68,89,53,67,30,73,16,69,87,40,31},
                    {4,62,98,27,23,9,70,98,73,93,38,53,60,4,23}};
int arr1[num][num] = {0};
//从上往下arr1={{75},{170,139},{187,217,221},....}
int main(){
    arr1[0][0] = arr[0][0];
    for (int i = 1; i <= 14; i++){
        for(int j = 0; j <= i;j++){
            if(j != 0 && arr1[i-1][j-1] > arr1[i-1][j])
                arr1[i][j] = arr1[i-1][j-1];
            if(j != i && arr1[i-1][j] > arr1[i-1][j-1])
                arr1[i][j] = arr1[i-1][j];
            arr1[i][j] += arr[i][j];
            printf("%d\n",arr1[i][j]);
        }
    }
    int max = 0;
    for(int i = 0;i < 16;i++){
        if(arr1[14][i] > max) 
            max = arr1[14][i];
    }
    printf("%d\n",max);
    return 0;
}

JAVA

package exercise;

public class _最大路径和1 {

	public static void main(String[] args) {
		int arr[][] = {{75},
                {95,64},
                {17,47,82},
                {18,35,87,10},
                {20,4,82,47,65},
                {19,1,23,75,3,34},
                {88,2,77,73,7,63,67},
                {99,65,4,28,6,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,4,68,89,53,67,30,73,16,69,87,40,31},
                {4,62,98,27,23,9,70,98,73,93,38,53,60,4,23}};
		int arr1[][] = new int[105][105];
		//从下往上存放和
		for(int i = 0;i < 15; i++) {
			arr1[14][i] = arr[14][i];
		}
		for(int i = 13; i >=0; i--) {
	       for (int j = 0; j <= i; j++) {
	    	   arr1[i][j] = arr[i][j];
	          if(arr1[i+1][j] > arr1[i+1][j+1]) {
	                arr1[i][j]+=arr1[i+1][j];
	            } else {
	                arr1[i][j]+=arr1[i+1][j+1];
	            }
	          System.out.println(arr1[i][j]);
	        }
	    }
		System.out.println(arr1[0][0]);
	}

}

结果:1074

  • 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、付费专栏及课程。

余额充值