matprod

#include <fstream>
using namespace std;
int a[1001][1001],s[1001][1001];

int main()
{
    int n;
    ifstream fin("matprod.in");
    fin>>n;
    for(int i=1;i<=n;++i)
        for(int j=1;j<=n;++j)
            fin>>a[i][j];
    for(int i=1;i<=n;++i)
        s[i][1]+=s[i-1][1]+a[i][1],a[i][1]=a[i-1][1]+a[i][1]*s[i-1][1];
    for(int j=2;j<=n;++j)
        for(int i=1;i<=n;++i)
            s[i][j]=s[i][j-1]+a[i][j],a[i][j]=a[i][j-1]+a[i][j]*s[i][j-1];
    ofstream fout("matprod.out");
    for(int i=1;i<=n;++i,fout<<'\n')
        for(int j=1;j<=n;++j)
            fout<<a[i][j]<<' ';
    return 0;
}