Sicily 3712. Matrix multiplication

#include<iostream>using namespace std;int max1[1001][1001];int max2[1002][1001];int sum[1001][1001];int main(){    int n;    while(cin >> n){    for(int i=0; i < n; i++)        for(int j=0; j < n; j++)            cin >> max1[i][j];    for(int i=0; i < n; i++)        for(int j=0; j < n; j++)            cin >> max2[i][j];    for(int i=0; i < n; i++)        for(int j=0; j < n; j++)            for(int f=0; f < n; f++){                sum[i][j]+=max1[i][f]*max2[f][j];                }    for(int i=0; i < n; i++)        for(int j=0; j < n; j++)            if(j<n-1)cout << sum[i][j] <<" ";            else cout << sum[i][j] <<endl;    for(int i=0; i < n; i++)        for(int j=0; j < n; j++)            sum[i][j]=0;     }    return 0;}