分数矩阵（找规律/DP）

题目描述
我们定义如下矩阵：

1/1 1/2 1/3
1/2 1/1 1/2
1/3 1/2 1/1

矩阵对角线上的元素始终是1/1，对角线两边分数的分母逐个递增。
请求出这个矩阵的总和。
输入

输入包含多组测试数据。每行给定整数N（N<50000），表示矩阵为N*N。当N=0时，输入结束。 输出

输出答案，结果保留2位小数。

样例输入

1
2
3
4
0

样例输出

1.00
3.00
5.67
8.83

思路一：动态规划（DP大法不超时）
当n == 1时，矩阵为：

1/1

dp[ 1 ] = 1.00

当n == 2时，矩阵为：

1/1 1/2
1/2 1/1

dp[ 2 ] = 3.00

当n == 3时，矩阵为：

1/1 1/2 1/3
1/2 1/1 1/2
1/3 1/2 1/1

容易发现红色部份为dp[ 2 ]且不难发现
dp[ 3 ] = 2 * dp[ 2 ] - dp[ 1 ] + 2 * (1 / 3)

当n == 4时，矩阵为：

1/1 1/2 1/3 1/4
1/2 1/1 1/2 1/3
1/3 1/2 1/1 1/2
1/4 1/3 1/2 1/1

容易发现红色部份为dp[ 3 ]且红色部份 = 蓝色部分 + dp[ 2 ]
dp[ 4 ] = 2 * dp[ 3 ] - dp[ 2 ] + 2 * (1 / 4)

… …

推导出dp[ n ] = 2 * dp[n - 1] - dp[n - 2] + 2 / n

AC代码

#include <cstdio>
#include <iostream>

using namespace std;

double dp[50010];

int main()
{
    int n;

    dp[1] = 1.00;
    dp[2] = 3.00;
    
    for (int i = 3; i < 50010; i++)
    {
        //打表
        dp[i] = 2 * dp[i - 1] - dp[i - 2] + (double)2 / i; 
    }
    while(scanf("%d",&n) && n) printf("%.2lf\n",dp[n]);
    return 0;
}