memset(c,0,sizeof(c));for(i=1;i<=n;i++)for(j=1;j<=n;j++){if(a[i][j]==0)continue;for(k=1;k<=n;k++)c[i][k]=c[i][k]+a[i][j]*b[j][k];}