оригинальный алгоритм сложения элементов матрицы

[Ledom]

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

[Ledom]

если долго мучатся что-нибудь получиться)
написал алгоритм считает правильно, возможно не самый простой вариант, но увы до другого не дошел...
написал через копирование элементов одного из треугольника, затем транспонирование основной матрицы и сложение с матрицей с копией треугольником, далее нахожу сумму элементов матрицы лежащих на и ниже главной диагонали....
вот так... может кто-нибудь предложит как-то упростить?

#include <stdio.h>
#include <conio.h>
#include <math.h>
void main()
{
unsigned char i,j,N,K;
float S=0,A[10][10],B[10][10],At[10][10],C[10][10];
clrscr();
printf("Vvedite razmernost' kvadratnoy matricu: ");
scanf ("%d",&N);
if(N%2==0)K=N/2;
else K=(N+1)/2;
for(i=1;i<=N;i++)
for(j=1;j<=N;j++)
	{
	printf("vvedite element matricu A[%d,%d]: ",i,j);
	scanf("%f",&A[i][j]);
	}
for(i=1;i<=N;i++)
for(j=1;j<=N;j++)
if(i>=j) B[i][j]=A[i][j];
else B[i][j]=0;

for(i=1;i<=N;i++)
for(j=1;j<=N;j++)
if(i+j>N+1)B[i][j]=0;

for(i=1;i<=N;i++)
for(j=1;j<=N;j++)
At[i][j]=A[j][i];

for(i=1;i<=N;i++)
for(j=1;j<=N;j++)
{
if(i+j<N+1) At[i][j]=0;
if((i==K)&&(i==j)) At[i][j]=0;
}

for(i=1;i<=N;i++)
{C[i][j]=0;
for(j=1;j<=N;j++)
C[i][j]=B[i][j]+At[i][j];
}

for(i=1;i<=N;i++)
for(j=1;j<=N;j++)
if(i>=j)S=S+C[i][j];
printf("\nS=%-5.2f",S);
getch();
}

[Hacker19_90]

А мож так попробовать

sum := 0;
//это для левого
for (j = 0; j < n/2; j ++)
{
     for(i = j; i < (n - j); i++)
     sum += matrix[i][j];
}
//это для правого
for (j = n-1; j > (n/2+1); j--)
{
     for(i = (n-j -1); i < (j + 1); i++)
     sum += matrix[i][j];
}

Ну примерно глянь мож доработаешь

[Ledom]

спасибо, так сразу не пришло... но еще погляжу... может ченить придумаю)

[Ledom]

нашел более рациональный вариант может кому поможет:

if(N1%2!=0) { //N1 размерность квадратной матрицы
	     K=(N1-1)/2;//проверка на четность для нахождения индекса элемента на пересечении диагоналей 
	      for(j=0;j<N1;j++)
		for (i=0;i<N1;i++)
		    if(j<=K)
			{
			if((i>=j)&&((i+j)<=(N1-1)))
			S+=X[i][j];
			}
		    else if((i<=j)&&((i+j)>=(N1-1))) S+=X[i][j];
	    }

  else {
	K=N1/2; // 
	for(j=0;j<N1;j++)
	 for (i=0;i<N1;i++)
	      if(j<K)
		    {
		     if((i>=j)&&(i+j<=N1-1))
		     S+=X[i][j];
		    }
	      else if((i<=j)&&(i+j>=N1-1)) S+=X[i][j];
       }