Как привести матрицу к унитреульной матрице или верхней треугольной матрице?

Как привести матрицу к унитреульной матрице или верхней треугольной матрице?
#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <iostream>
#define N 5 
void main (){
int mt, nt, i,j,k,c;
int A[N+3][N]={
1,1,0,0,0,
1,1,0,1,0,
1,0,1,0,0,
1,0,1,1,1,
1,0,0,0,0,
0,0,1,0,0,
0,1,0,0,1,
0,1,0,1,0}, B[N+3]={0,1,2,3,4,5,6,7};
mt=8; nt=5;
for(int u=0; u<mt; u++){
    for(int y=0; y<nt; y++){
        printf("%d\t",A[u][y]);
    if(y==(nt-1))printf("\n");}
}
printf("\n");
 
 
k=0;
while(k<mt){
j=k;
if(A[B[j]][k]==0 && j<mt){j++;}
if(j>(mt-1)){k++;}
if(j>(k-1)){c=B[j];B[j]=B[k];B[k]=c;}
i=k;
while(i<(mt-1)){
    i++;
if(A[B[i]][k]==1){
    for(j=0;j<nt;j++){
        A[B[i]][j]=(A[B[i]][j]+A[B[k]][j])%2;}
}
}
for(int u=0; u<N+3; u++){
    for(int y=0; y<N; y++){
        printf("%d\t",A[u][y]);
    if(y==(N-1))printf("\n");}
}
printf("\n");
k++;
}
}
написал код, но работает не так. может кто посмотрит
7 ответов

а вот лови прямо здесь на форуме лежит я правдо немного отредактировал что бы сдать преподу мог вот функция
#define FOR_k  for (int k = 0; k < dim; ++k)
#define FOR_j  for (int j = 0; j < dim; ++j)
#define FOR_i  for (int i = k + 1; i < dim; ++i)
void swap(****** &a, ****** &b)
{
    ****** temp = a;
    a = b;
    b = temp;                         
 
}
 
void Gauss(****** **A, ****** *B, int dim)
{
    FOR_k
    {
        if (A[k][k] == 0.0)
        {
           FOR_i
            {
                if (A[i][k] != 0.0)
                {
                    FOR_j
                        swap(A[i][j], A[k][j]);
 
                    swap(B[i], B[k]);
 
                    break;
                }
            }
        }
 
        ****** div = A[k][k];
 
        FOR_j
            A[k][j] /= div;
 
        B[k] /= div;
 
       FOR_i
        {
            ****** multi = A[i][k];
 
            for (int j = 0; j < dim; ++j)
                A[i][j] -= multi * A[k][j];
 
            B[i] -= multi * B[k];
        }
    }
 
    for (int k = dim - 1; k > 0; --k)
    {
        for (int i = k - 1; i + 1 > 0; --i)
        {
            ****** multi = A[i][k];
 
            for (int j = 0; j < dim; ++j)
                A[i][j] -= multi * A[k][j];
 
            B[i] -= multi * B[k];
        }
    }                                                    
}
решение матриц методом Гауса


dim что за переменная? размерность?


dim что за переменная? размерность?
да размерность


у меня есть матрица из нулей и единиц, нужно получить унитреугольную матрицу. что тогда я должен в качестве матрицы B[] подавать?


int A[N+3][N]={
1,1,0,0,0,
1,1,0,1,0,
1,0,1,0,0,
1,0,1,1,1,
1,0,0,0,0,
0,0,1,0,0,
0,1,0,0,1,
0,1,0,1,0}, B[N+3]={0,1,2,3,4,5,6,7};
я так понимаю это матрица и столбец свободных членов ?? так если так то программа сделает еденичную матрицу треугольную


нет мне нужно только совершить эквивалентные преобразования над матрицей А, матрица В тут не при чём


нет мне нужно только совершить эквивалентные преобразования над матрицей А, матрица В тут не при чём
матрица) Б у тебя там массив когда вижу такое сразу понимаю что это матрица со столбоцом свободных членов