实现N*N矩阵的乘法，矩阵有一维数组表示

下面是一个实现3*3维数组的乘法！
#include<iostream>
#include<vector>
using namespace std;
int k=0;
vector<int> c;
void array(int a[],int n,int b[])
{

for(int i=0;i<n*n;)
{
  for(int j=0;j<n;j++)
  {
   k=a[i]*b[j]+a[i+1]*b[j+3]+a[i+2]*b[j+6];
   c.push_back(k);

  }
    i=i+n;
}
for(int m=0;m<n*n;m++)
  cout<<c[m]<<endl;

}
int main()
{
int a[]={1,2,3,4,5,6,7,8,9};
int b[]={1,2,3,4,5,6,7,8,9};
array(a,3,b);