算法-分数矩阵

问题概述

---

题目：算法-分数矩阵

问题描述：
我们定义如下矩阵：
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

---

1. 问题分析

---

1.1 算法核心

本题关键就是理解对称矩阵的特点，还需要注意的是上下三角矩阵实际上只用算一边就可以了，不然很容易导致超时错误！
11 12 13 14 15
12 11 12 13 14
13 12 11 12 13
14 13 12 11 12
15 14 13 12 11

---

1.2 算法分析

因为对角线都是1/1，从正对角线开始往两侧递增，依次是1/2，1/3…，且这些矩阵元素的个数也是递增的，正对角线的值可以直接作为结果初始值，因为是固定的，然后根据对称三角矩阵的特点，斜对角线元素的个数与其当前循环次数有关！所以本题解题关键就是理解对称矩阵的特点以及找出矩阵元素值的特点；
因为是对称矩阵，所以我们专注于下三角矩阵的值，最后加上每一条斜对角线值的时候乘以2就行；
当n等于5时，发现1/2的个数有4个，1/3的个数有3个，1/4的个数有2个，1/5的个数有1个，所以下三角矩阵的值也就等于8个1/2，6个1/3，4个1/4，2个1/5，以及正对角线固定值5相加；
总结规律就是:
矩阵的值等于
n + 2 * (1/n + 1/(n - 1) * 2 + 1/(n - 2) * 3 + … + 1/2 * (n - 1));
最后别忘记了保留两位小数嗷！

---

2. 解决方案

#include <iostream>
#include <vector>
using namespace std;

#define maxLength 1001

int n, t;
vector<double> result;

int main()
{
	while (scanf("%d", &n) != EOF && n != 0)
	{
		double sum = n;
		for (int i = n - 1, j = 2; i > 0; i--,j++)
			sum += i / j * 2.0;
		result.push_back(sum);
	}
	for (int i = 0; i < result.size(); i++)
		printf("%.2lf\n", result[i]);

	return 0;
}

3. 资源分享

旗木白哉のGitHub：
https://github.com/YangMengHeng
C++ / C源代码下载地址：
https://github.com/YangMengHeng/myCode/tree/master/VSCode
Java源代码下载地址：
https://github.com/YangMengHeng/myCode/tree/master/Java
旗木白哉のblog源代码下载地址：
https://github.com/YangMengHeng/myCode/tree/master/%E6%97%97%E6%9C%A8%E7%99%BD%E5%93%89%E3%81%AEblog

---