Как реализовать вычисления по формуле без ошибки выхода за пределы массива?

int s;
                            s = p[i];
                            int summa;
                            int v = 0;
 
                            for (int n=0; n<s-1; n++)
                            {
                                int k = p[n];
                                summa= Vi[p[k]] * matrix[p[k]][j];
                                v = c[j] - summa;
                                a.Add(v);
 
                            }
                            min = a.Min();
                            Vi[i] = min;

                for (int j = 0; j < collumncount; j++)
                {
                    //первое условие
                    if (p[i] == 0)
                    {
                        if (matrix[i] [j] == 1)
                        {
                            a.Add(c[j]);
                            
    
                        }
                            min = a.Min();
                            Vi[i] = min;
                       
                    }
                    
                  //второе условие
                    if (p[i] == 1) 
                    {
                        if (matrix[i] [j] == 1)
                        {
                            int ind;
                            ind = GetNullIndex(p, 0);
 
                            a.Add(c[j] - Vi[ind]*matrix[ind][j]);
                            
                            //a.Clear();
 
                        }
                            min = a.Min();
                            Vi[i] = min;
 
                    }}

//третье условие
else
                    {
                        if (matrix[i][j] == 1)
                        {   
                            int s;
                            s = p[i];
                            int summa;
                            int v = 0;
 
                            for (int n=0; n<s-1; n++)
                            {
                                int k = p[n];
                                    
                                summa= Vi[p[k]] * matrix[p[k]][j];
                                v = c[j] - summa;
                                a.Add(v);
 
                            }
                            min = a.Min();
                            Vi[i] = min;
 
 
 
                        }