Особый цикл

#include "stdafx.h"
#include <iostream>
 
using namespace std;
 
#define dim 100 
 
using namespace std; 
 
int N; //количество точек первой группы
int K; //количество точек второй группы
 
//динамический массив первой группы точек
double points1_x[dim]; 
double points1_y[dim]; 
 
 
int main()
{
    setlocale(LC_ALL,"rus");
 
    cout << "\nВведите количество точек в группе, N = 3, 6, i + 3: \n";
    cin >> N;
 
    cout << "\nВведите координаты точек в группе: \n";
 
    for(int i = 0; i < N; i++ )
    {
        
            cout << "\nКоординаты " << i + 1 <<"-й точки: \n" ;
            cout << " x[" << i + 1 << "] = ";
            scanf_s("%lf",&points1_x[i]);
            cout << " y[" << i + 1 << "] = ";
            scanf_s("%lf",&points1_y[i]);
        
    }
 
    cout << "\nКоординаты имеют точки: \n\n";
 
        for(int i = 0; i < N; i++)
        {
            
            cout << " x[" << i + 1 << "] = ";
            printf_s("%3.2lf",points1_x[i]);
            cout << " y[" << i + 1 << "] = ";
            printf_s("%3.2lf",points1_y[i]);
            printf("\n");
 
        }
 
 
            for(int i = 0; i < N; i++)
            {
                for(int j = 0; j < N; j++)
                {
                    if( i == j) //Проход цикла, исключая n = 1, ... , j член 
                    {
                        continue;
                    }
 
                                         //Скорее всего здесь что-то нужно добавить
                    
                }
            }
 
 
    return 0;
}

bool do_step_for_i( Smth& result, const YourType& array, vector& placesArray, int place = 0, int i = K )
{
if ( i == 0 )
{
//Делаешь, что хочешь - выбранные элементы маркируются значениями в placesArray
return false;
}
 
for ( int k = place; k < N - i; ++k )
{
placesArray.push_back( k );
do_step_for_i(  result, array, placesArray, k + 1, i - 1 );
placesArray.pop_back( k );
}
}

#include <cmath>
#include <ctime>
#include <iostream>
using namespace std;
 
struct my_point
{
    double x;
    double y;
};
 
my_point getCentoid(my_point p1, my_point p2, my_point p3);
my_point * getCentr(my_point * pCoord, int n);
 
int main()
{
    int N = 15;
    int i = 0;
    srand(time(0));//Устанавливаем нач значение ген-ра случ чисел
    my_point * pCoords = new my_point[N];
    cout<<"Input points : "<<endl;
    for(i = 0; i < N; i++)
    {
        pCoords[i].x = (rand()%20 - 10) / 10.0;
        pCoords[i].y = (rand()%20 - 10) / 10.0;
        cout<<"("<<pCoords[i].x<<";"<<pCoords[i].y<<")"<<endl;
    }
    do
    {
        pCoords = getCentr(pCoords, N);
    }
    while(1 < (N /= 3));
    cout<<"Center of figure : "<<"("<<pCoords[0].x<<";"<<pCoords[0].y<<")"<<endl;
    return 0;
}
 
my_point getCentoid(my_point p1, my_point p2, my_point p3)
{
    my_point pCentr = {0};
    pCentr.x = (p1.x + p2.x + p3.x) / 3;
    pCentr.y = (p1.y + p2.y + p3.y) / 3;
    return pCentr;
}
 
my_point * getCentr(my_point * pCoord, int n)
{
    my_point * pCentr = new my_point[n / 3 + 1];
    my_point p1 = {0};
    my_point p2 = {0};
    my_point p3 = {0};
    int i, j = 0;
    for(i = 2; i < n; i += 2, j++)
    {
        p1 = pCoord[i - 2];
        p2 = pCoord[i - 1];
        p3 = pCoord[i - 0];
        pCentr[j] = getCentoid(p1, p2, p3);
    }
    return pCentr;
}

Point ConvexCenter( const Point * src, int num )
{
 Point cntr(0, 0);
 float sumW = 0.0f;
 for (int i = 0; i < num - 2; ++i) {
   float w = cross(src[0], src[i + 1], src[i + 2]);
   cntr += (src[0] + src[i + 1] + src[i + 2]) * w / 3;
   sumW += w;
 }
 cntr /= sumW;
 return cntr;
}

for(i = 2; i < n; i += 2, j++)
    {
        p1 = pCoord[i - 2];
        p2 = pCoord[i - 1];
        p3 = pCoord[i - 0];
        pCentr[j] = getCentoid(p1, p2, p3);
    }

 float sumW = 0.0f; //0.0f - что это?)

Point CenterMonteCarlo( const Point bounds[2], const Point * src, int num, int numTry )
{
 srand(0);
 Point cntr(0, 0);
 int numInside = 0;
 for (int i = 0; i < numTry; ++i) {
  float x = float(rand()) / RAND_MAX;
  float y = float(rand()) / RAND_MAX;
  Point test(bounds[0].x + x * (bounds[1].x - bounds[1].x),
                  bounds[0].y + y * (bounds[1].y - bounds[1].y));
  if (!PtInside(test, src, num)) continue;
  cntr += test;
  ++numInside;
 }
 if (!numInside) return (bounds[0] + bounds[1]) / 2;  // too small
 return cntr / numInside;
}