MMultiplication

#include <fstream>
#define N 1000
using namespace std;
int a[N][N],b[N][N],c[N][N];

int main()
{
    int n,m;
    ifstream f("mmultiplication.in");
    f>>n>>m;
    for(int i=0;i<n;++i)
        for(int j=0;j<m;++j)
            f>>a[i][j];
    int p,q;
    f>>p>>q;
    ofstream g("mmultiplication.out");
    if(m!=p)
    {
        g<<-1;
        return 0;
    }
    p=q;
    for(int i=0;i<m;++i)
        for(int j=0;j<p;++j)
            f>>b[i][j];
    for(int i=0;i<n;++i)
        for(int j=0;j<p;++j)
            for(int k=0;k<m;++k)
                c[i][j]+=a[i][k]*b[k][j];
    for(int i=0;i<n;++i)
    {
        for(int j=0;j<p;++j)
            g<<c[i][j]<<' ';
        g<<endl;
    }
    return 0;
}