правило Саррюса (правило треугольников) C++ или Delphi

[carrie bradshaw]

привет, всем активным участникам форума. помогите, пожалуйста, написать программу на языке C++ или Delphi .

задание:
нужно написать программу, чтобы она решала матрицу по правилу Саррюса (правило треугольников)
вот картинка этого правила:
k3_image001.jpg k3_image002.jpg
вот так нужно, чтобы считалась матрица:
=(a11*a22*a33)+(a12*a23*a31)+(a13*a 21*a32)-(a13*a22*a31)-(a11*a32*a23)-(a12*a21*a33)

на форме компонент StrinGrid в который нужно вводить матрицу и кнопка для её решения (чтобы выводилось решение полностью, а не только ответ).

всем огромное спасибо за помощь и внимание!!!

[carrie bradshaw]

Если кто то может помочь, то напишите, пожалуйста!

[carrie bradshaw]

есть программа на Pascal, возможно кто-то знает как это можно реализовать на Delphi или C++

program zd;
Const n=3; m=3;

Type
matr  = array[1..n,1..m] of integer;

Const
A: matr = ( ( 2, 10, 1 ),
            ( 3, 4, -1 ),
            ( 1, 6, 2 )
           );
Var
S:Real;
Begin
S:=(   ( ( A[1,1] )*( A[2,2] )*( A[3,3] ) )+
       ( ( A[1,2] )*( A[2,3] )*( A[3,1] ) )+
       ( ( A[1,3] )*( A[2,1] )*( A[3,2] ) )-
       ( ( A[1,3] )*( A[2,2] )*( A[3,1] ) )-
       ( ( A[1,1] )*( A[3,2] )*( A[2,3] ) )-
       ( ( A[1,2] )*( A[2,1] )*( A[3,3] ) )  );

Writeln('( A[1,1] )*( A[2,2] )*( A[3,3] ) + ( A[1,2] )*( A[2,3] )*( A[3,1] ) + ( A[1,3] )*( A[2,1] )*( A[3,2] ) - ( A[1,3] )*( A[2,2] )*( A[3,1] ) - ( A[1,1] )*( A[3,2] )*( A[2,3] ) - ( A[1,2] )*( A[2,1] )*( A[3,3] ) = ',
 A[1,1],'*',A[2,2],'*',A[3,3],'+',A[1,2],'*',A[2,3],'*',A[3,1],'+',A[1,3],'*',A[2,1],'*',A[3,2],'-',A[1,3],'*',A[2,2],'*',A[3,1],'-',A[1,1],'*',A[3,2],'*',A[2,3],'-',A[1,2],'*',A[2,1],'*',A[3,3],'=', S:0,4 );
END.

надеюсь на помощь!

[carrie bradshaw]

вот интерфейс программы:
123.JPG

всем огромное спасибо за внимание и надеюсь за помощь!

[carrie bradshaw]

вот программа на C++

#pragma hdrstop
#pragma argsused

#include <tchar.h>
#include <stdio.h>
#include <iostream>
#include<cstdio>
using std::cout;
using std::cin;
using std::endl;

void qwadMatrixPrint(int* *pM, int n)
{
    for (int str=0;str<n;str++)
	{
        for (int col=0;col<n;col++)
        {
			std::cout<<pM[col][str]<<"\t";
        }
		std::cout<<"\n";
	}
}
//функция вывода произвольной матрицы nxn в файл
void MatrixMovF(int n)
{
	for (int str=0;str<n;str++)
	{
        for (int col=0;col<n;col++)
		{
			std::cout<<((rand() % 100)*((rand() % 3)-1))<<("\t");
		}
		std::cout<<"\n";
	}
}
//функция нахождения определителя матрицы по правилу треугольника
int determMatrix(int* *pM, int n,int deep=0)
{
    if(n==2)
		return pM[0][0+deep]*pM[1][1+deep]-pM[0][1+deep]*pM[1][0+deep];

	int znak= 1;//переменная определяющая знак алгебраического дополнения Aij
//  ( Aij=(-1)^(i+j)*Mij где Mij-минор;минор произвольного элемента квадратной
//  матрицы Pij  равен определителю матрицы, полученной из исходной
//вычеркиванием i-ой строки и j-го столбца.)
	int sum_det= 0;

    for(int i= 0; i< n; i++)
    {
		int* *pI= new int* [n-1];//создание нового масива указатей
        //указывающих только на столбцы минора(верхние ненужные строки(ка)
        //отсекаються за щет переменной deep(глубина рекурсии))

		//создание очередного минора
            int p=0;
        for(int ii= 0; ii< n; ii++ )
        {
			if(ii==i)continue;

            pI[p]=pM[ii];
            p++;
		}

        sum_det+= znak*pM[i][deep]*determMatrix(pI,(n-1),(deep+1));//переменная deep(глубина рекурсии)
        //указывает первую строку с которой начинаетья данный минор

        delete [] pI;//освобождение памяти
        znak*= -1;

    }
	return sum_det;
}


int main(int argc, char * argv[])
{
    setlocale(LC_ALL, "rus");

    int n=0;
    cout<<"Введите порядок матрицы : "<<endl;
	cin>>n;

    freopen("test.in", "w", stdout);//открыть файл test.in для записи
                                    //связать с потоком stdout(вывода)
    MatrixMovF(n);
	freopen("test.in", "r", stdin);//открыть файл test.in для чтения
                                    //связать с потоком stdin (ввода)

    int* *pMat=new int* [n];//pMat присваеваеться указатель на масив указателей на столбцы матрицы
    //Инициализация элементов матрицы
	for(int col=0; col<n; col++) pMat[col]=new int[n];


        for (int str=0;str<n;str++)
        {
			for (int col=0;col<n;col++)
            {
                cin>>pMat[col][str];//инициализация элемента столбца матрицы
                                    //прочитаным из файла test.in значением
			}
        }
    freopen("test.out", "w", stdout);//открыть файл test.in для записи
									//связать с потоком stdout(вывода)

    qwadMatrixPrint(pMat, n);
	int determ =determMatrix(pMat, n);
    cout<<"Определитель матрицы "<<n<<"x"<<n<<endl;
    cout<<"равен "<<determ<<endl;



//освобождение памяти
for (int i=0;i<n;i++)
    {
        delete [] pMat[i];
	}
    delete [] pMat;

    return 0;
}

ошибка в строчке:

		std::cout<<((rand() % 100)*((rand() % 3)-1))<<("\t");

[reflexx]

На Delphi 7

[carrie bradshaw]

reflexx, огромное вам спасибо!