Умножение матриц по Винограду

Пытаюсь сделать умножение матриц по Винограду, но как-то не очень получается. Прошу помощи у вас. То что есть :

#include <iostream>
#include <vector>
int main()
{
    const int n = 5;
    const int max = 10;
    const int min =1;
    int arr1[n][n];
    int arr2[n][n];
    int arr3[n][n];
    std::vector <int> rowFactor;
    std::vector <int> columnFactor;
    int d = n/2;
    for(int i = 0; i < n; ++i)
    {
      for(int j = 0; j < n; ++j)
      {
        arr1[i][j] = rand()%(max-min+1) + min;
        arr2[i][j] = rand()%(max-min+1) + min;
      }
    }

    for(int i = 0 ;i<n;++i)
    {
      rowFactor.insert(rowFactor.end(), arr1[i][1]*arr1[i][2]);
      for(int j =0;j<d;++j)
      {
        rowFactor.insert(rowFactor.end(),rowFactor[i] + arr1[i][2*j - 1] * arr1[i][2*j]);
      }
    } 
    for(int i = 0 ;i<n;++i)
    {
      columnFactor.insert(columnFactor.end(), arr2[i][1]*arr2[i][2]);
      for(int j =0;j<d;++j)
      {
        columnFactor.insert(columnFactor.end(), columnFactor[i] +  arr2[2*j - 1][ i] * arr2[2*j][ i]);
      }
    }
    for(int i =0; i<n;++i)
    {
      for(int j =0;j<n;++j)
      {
        arr3[i][j] = -rowFactor[i] - columnFactor[j];
        for(int k = 0; k< d;k++)
        {
            arr3[i][j]=arr3[i][j]+(arr1[i][ 2*k-1]+arr2[2*k][j])*(arr1[i][ 2*k] + arr2[2*k-1][j]);
        }
      }
    }
    std::cout<<"\nresult\n"<<std::endl;
    for (int i=0; i<n; i++) 
    {
      for (int j=0; j<n; j++)
         std::cout << " " << arr3[i][j];
      std::cout << std::endl; 
    }
}

Даже что-то считает, но не правильно. Думаю в этой строке ошибка

arr3[i][j]=arr3[i][j]+(arr1[i][ 2*k-1]+arr2[2*k][j])*(arr1[i][ 2*k] + arr2[2*k-1][j]);