Răspuns :
#include <iostream>
#include<fstream>
using namespace std;
int a[100][100],n;
void citire()
{
ifstream f("matrice.in");
f>>n;
int i,j;
for(i=1; i<=n; i++)
for(j=1; j<=n; j++)
f>>a[i][j];
}
int cmmdc(int a,int b)
{
int r;
while(b!=0)
{
r=a%b;
a=b;
b=r;
return a;
}
}
int suma(int x)
{
int i,s=0;
for(i=1; i<=n; i++)
s=s+a[x][i];
return s;
}
int produs(int y)
{
int i,p=1;
for(i=1; i<=n; i++)
p=p*a[i][y];
return p;
}
int main()
{
citire();
ofstream g("matrice.out");
for(int i=1; i<=n; i++)
for(int j=1; j<=n; j++)
if(i==j&&a[i][j]==cmmdc(suma(i),produs(j)))
g<<a[i][j]<<" ";
}