草稿1

Posted 2021-01-22 junk-yao-blog

#include<cstdio>
#include<iostream>
#include<Eigen/Dense>
using namespace std;
using namespace Eigen;
int n, m;
MatrixXd input;
MatrixXd A;
MatrixXd B;
VectorXd b;
VectorXd c;
void presolve(){    
    //***读取输入大小
    
    cin >> m >> n;
    
    //***输入系数矩阵并建立初始A, B矩阵
    
    input.resize(m, n);
    for(int i = 0; i < m; ++i)
        for(int j = 0; j < n; ++j)
            cin >> input(i, j);
    n += m;
    A.resize(m, n);
    B.resize(m, m);
    B = MatrixXd::Identity(m, m);
    A << input, B;
//    cout << "Now A is: " << A << endl;
    
    //***输入b, c矩阵
    
}
int main(void){
    presolve();    
}