国科大——高性能计算系统习题1

环境：visual stdio 2022
时间：2022.3.13
——2021届沈计所 Sen.W

1. 矩阵加法（C++一维数组实现）

   描述：在主程序中定义了矩阵A和B (比如C语言, A[50][67],| B[73][49]),请写一个通用子程序完成矩阵相加，使得对于任何m x n的子矩阵A与子矩阵B相加都能调用该子程序，用C或Fortran语言(比如,matradd(m, n,… A, … B, …C)。

   • 进行了运算符的重载
   • 使用一维数组模拟二维数组（这个方式行列需要自己计算可能出现失误）
     可以使用二维数组，c语言也支持，详情见https://blog.csdn.net/qq_43840665/article/details/128003164

#include<iostream>
using namespace std;


class matirx {

private:
	int rows;// 矩阵行数
	int cols;// 矩阵列数
	int* mValue;// 矩阵的值

public: 
	// 标准构造函数格式
	
	matirx() {}
	// 获取输入
	void initMatirx() {
		// 获取矩阵的行列和值
		cout << "请输入行数：";
		cin >> this->rows;
		cout << endl;
		cout << "请输入列数：";
		cin >> this->cols;
		cout << endl;
		this->mValue = new int[rows * cols];
		for (int i = 0; i < rows; i++)
			for (int j = 0; j < cols; j++) {
				cout << "请输入m[" << i + 1 << "][" << j + 1 << "]：";
				cin >> mValue[i * cols + j];
				cout << endl;
			}
	}
	void deleteMatirx() {
		delete[] this->mValue;
	}
	// 打印矩阵
	void printMatirx() {
		for (int i = 0; i < rows; i++){
			for (int j = 0; j < cols; j++) {
				cout << mValue[i * cols+ j]<<"\t";
			}
			cout << endl;	
		}
	}
	// 重载加法运算符进行矩阵相乘
	matirx operator+(const matirx &m2) {
		// 行和列
		if (this->cols != m2.rows) {
			cout << "input error!!!";
			return matirx();// 返回空对象
		}
		
		matirx m3;
		m3.mValue = new int[rows * m2.cols];
		m3.rows = rows;
		m3.cols = m2.cols;
		// 矩阵乘法
		int sum = 0;
		for (int i = 0; i < rows; i++) {
			for (int j = 0; j < m2.cols; j++) {
				for (int k = 0; k < cols; k++) {
					sum += mValue[i * cols + k] * m2.mValue[j + k * m2.cols];
				}
				m3.mValue[i*m2.cols + j ] = sum;// 离大普，卡在这一下午
				sum = 0;
			}
			
		}	
		return m3;
	}

	int getRows() {
		return this->rows;
	}
	int getCols() {
		return this->cols;
	}

};

int main()
{
	
	matirx m1,m2,m3;
	m1.initMatirx();
	m2.initMatirx();

	cout << "m1:" << endl;
	m1.printMatirx();
	cout << "m2:" << endl;
	m2.printMatirx();

	m3 = m1 + m2;
	m3.printMatirx();
}

2. 设多项式$P_n(x)=a_n{x}^n+a_{n-1}{x}^{n-1}+\cdot \cdot \cdot +a_{1}x+a_0$，请给出计算$P_n(x_0)$的方法并写出计算它的子程序

#include<iostream>
#include<math.h>
using namespace std;

// 定义链表节点存储每一项的值。
//不使用数组存储，因为避免高次项数而低项数对于存储的浪费
typedef struct L{
	int n; // 第n项
	double a_n; // 常数项a_n的值
	struct L* next; // 下一项的链表指针
}linkList;

double Add(linkList* p, int x) {
	double result = 0;
	while (p != NULL) {
		result += p->a_n * pow(x, p->n);
		p = p->next;
	};
	return result;
}
 // 打印公式
void show(linkList* pHead) {
	linkList* p = pHead->next;
	cout << "多项式为：" << endl;
	while (p != NULL) {
		cout << p->a_n << "*a_" << p->n << "*X^" << p->n << " + "  ;
		p = p -> next;
	};
	printf("\b\b ");// 删除最后一个加号和X^0
}


int main() 
{

	linkList* pHead = new linkList;// 头指针
	pHead->n = -1;
	pHead->a_n = 0;
	pHead->next = NULL;
	double x = 0;
	linkList* p;// 工作指针
	int maxTerm = 0;// 最大项次数
	double temp = 0; // 暂存常数项a_n的值
	cout << "请输入最大项幂指数：";
	cin >> maxTerm;
	cout << "请输入X的值：";
	cin >> x;
	// 获取输入，使用头插法形成链表
	for (int i = 0; i < maxTerm+1; i++) {
		cout << "请输入常数项a[" << i << "]：";
		cin >> temp;
		if (temp > 0) {// 常数项为0的节点不计入链表中
			p = new linkList;
			p->n = i;
			p->a_n = temp;
			p->next = pHead->next;
			pHead->next = p;
		}
	}
	show(pHead);
	cout <<endl<< "结果为：" << Add(pHead, x) << endl;
	return 0;
}

3. 对任意给定的x和n，用算法1写一个程序计算$x^n$(使用矩阵计算多项式)

少年，等待你的答案！！！！