Множество точек и окружностей

#include "stdafx.h"
#include <iostream>
#include <cstdlib>
#include <clocale>
#include <ctime>
#include <cmath>
#include <iostream>
#include <string>
 
using namespace std;
 
 
struct point {
    double x,y;
};
struct okr {
    double x,y,r;
};
 
const int nmax=30;
 
/* функция определяет пересекается прямая с окржуностью или нет */
bool Peres(point a, point b, okr c) 
{
    double h; /*строим треугольник с концами в координатах прямой и вершиной в центре окружности. Ищем высоту от координаты окружности к прямой*/
    h = (abs((b.y-a.y)*c.y + (a.x-b.x)*c.x + (a.x*b.y - a.y*b.x)))/(sqrt((b.y-a.y)*(b.y-a.y)+(a.x-b.x)*(a.x-b.x)));
    return (h<c.r); // если высота меньше радиуса, пересекаются
}
 
 
int main() {
    point a[nmax];
    okr b[nmax];
    int n,m,i,j,k,p,mx,imx,jmx;
 
    do {
        cout << "Количество точек до " << nmax << " n=";
        cin >> n;
    }while(n<2 || n>30);
 
    do {
        cout << "Количество окружностей до " << nmax << " m=";
        cin >> m;
    }while(m<1 || m>30);
 
    for(i=0; i<n; ++i) {
       cout << "Введите координаты точки("<<i+1<<"):\nX=";
    cin >> a[i].x;
    cout << "Y=";
    cin >> a[i].y;
    }
 
    for(i=0; i<m; ++i) {
        cout << "Введите координаты окружности("<<i+1<<"):\nX=";
    cin >> b[i].x;
    cout << "Y=";
    cin >> b[i].y;
    cout << "R=";
    cin >> b[i].r;
    
    }
 
    cout << "Координаты точек:\nX:";
    for(i=0;i<n;++i)
        cout << a[i].x << " ";
    cout << "\nY:";
    for(i=0;i<n;++i)
        cout << a[i].y << " ";
    cout << "\n\nПараметры окружностей:\nX:";
    for(i=0;i<m;++i)
        cout << b[i].x << " ";
    cout << "\nY:";
    for(i=0;i<m;++i)
        cout << b[i].y << " ";
    cout << "\nR:";
    for(i=0;i<m;++i)
        cout << b[i].r << " ";
    cout << "\n\n";
    /* ищем прямую с максимальными пересечениями */
    mx=0;
    imx=0;
    jmx=0;
 
    for(i=0;i<n-1;++i)
        for(j=i+1;j<n;++j) 
        {
            k=0;
            for(p=0;p<m;++p)
                if(Peres(a[i],a[j],b[p]))
                    ++k;
            if(k>mx) 
            {
                mx=k;
                imx=i;
                jmx=j;
            }
        }
 
    if(mx==0)
        cout << "Нет пересекающихся прямых с окружностями\n";
    else
        cout << "Максимальное число пересечений прямой с окружностями="<< mx
            << "\nЭта прямая проходит через точки ("<<a[imx].x<<';'<<a[imx].y<<") и ("<<a[jmx].x<<';'<<a[jmx].y<<")\n";
    cout << "\n";
    system("pause");
}