Dev C++ или VS построение фракталов Драконовы ключи

#define _CRT_SECURE_NO_WARNINGS  // чтоб не жаловался на всякие небезопасные функции
 
#include<windows.h>
#include<string.h>
 
#include <vector>
#include <math.h>
 
#define M_PI        3.14159265358979323846   // стянул из GCC С++
 
using namespace std;
 
LRESULT CALLBACK WndProc(HWND, UINT, WPARAM, LPARAM);
char szProgName[] = "Имя программы";
 
 
/*
Построение ломанной удобно разбить на две функции.
Первая возвращает массив точек, соответсвующих углам ломаной:
*/
vector<POINT> GetDragonPolyline(int order, int cxClient, int cyClient);
 
/*
Вторая (главная) запрашивает ломанную в виде массива точек и выводит её на экран:
*/
void Dragon(int order, HDC hdc, int cxClient, int cyClient);
 
vector<POINT> GetDragonPolyline(int order, int cxClient, int cyClient)
{
    if (order == 0)
    {
        // ломанная нулевого порядка состоит из одного сегмента
        vector<POINT> line;
        POINT p;
        
        p.x = cxClient / 4;
        p.y = cyClient / 2;
        
        // добавить первую точку
        line.push_back(p);
 
        p.x = 3 * cxClient / 4;
        p.y = cyClient / 2;
 
        // добавить последнюю точку
        line.push_back(p);
 
        return line;
    }
 
    // получить ломаную предыдущего порядка
    vector<POINT> prevLine = GetDragonPolyline(order - 1, cxClient, cyClient);
    vector<POINT> line;
 
    // знак направления угла (1 - влево, -1 - вправо)
    int DirSign = 1;
 
    // начальная точка ломаной не изменяется
    line.push_back(prevLine[0]);
    for (size_t i = 0; i < prevLine.size() - 1; ++i)
    {
        POINT p1 = prevLine[i];
        POINT p2 = prevLine[i + 1];
 
        double alpha = atan2(p2.y - p1.y, p2.x - p1.x) - DirSign * M_PI / 4;
        double R = sqrt((p1.x - p2.x) * (p1.x - p2.x) + (p1.y - p2.y) * (p1.y - p2.y)) / sqrt(2);
 
        // найти новую точку ломаной (pc)
        POINT pc;
        pc.x = p1.x + R * cos(alpha);
        pc.y = p1.y + R * sin(alpha);
 
        // добавить pc и p2 в список точек текущей ломаной
        line.push_back(pc);
        line.push_back(p2);
 
        // изменить направление на противоположное
        DirSign = -DirSign;
    }
    return line;
}
 
void Dragon(int order, HDC hdc, int cxClient, int cyClient)
{
    vector<POINT> line = GetDragonPolyline(order, cxClient, cyClient);
    for (size_t i = 0; i < line.size() - 1; ++i)
    {
        MoveToEx(hdc, line[i].x, line[i].y, NULL);
        LineTo(hdc, line[i + 1].x, line[i + 1].y);      
    }
}
 
 
int WINAPI WinMain(HINSTANCE hInstance, HINSTANCE hPrevInstance, LPSTR lpszCmdLine, int nCmdShow)
{
    HWND hWnd;
    MSG lpMsg;
    WNDCLASS w;
 
    w.lpszClassName = szProgName; //имя программы - объявлено выше
    w.hInstance = hInstance; //идентификатор текущего приложения
    w.lpfnWndProc = WndProc; //указатель на функцию окна
    w.hCursor = LoadCursor(NULL, IDC_ARROW); //загружаем курсор
    w.hIcon = 0;
    w.lpszMenuName = 0;
    w.hbrBackground = (HBRUSH)GetStockObject(WHITE_BRUSH); //цвет фона окна
    w.style = CS_HREDRAW | CS_VREDRAW;
    w.cbClsExtra = 0;
    w.cbWndExtra = 0;
 
    //Если не удалось зарегистрировать класс окна - выходим
    if (!RegisterClass(&w))
        return 0;
 
    //Создадим окно в памяти, заполнив аргументы CreateWindow
    hWnd = CreateWindow(szProgName, //Имя программы
        "Грфические возможности Win32 API", //Заголовок окна
        WS_OVERLAPPEDWINDOW, //Стиль окна - перекрывающееся
        100, //положение окна на экране по х
        100, //положение по у
        1024, //ширина
        768, //высота
        (HWND)NULL, //идентификатор родительского окна
        (HMENU)NULL, //идентификатор меню
        (HINSTANCE)hInstance, //идентификатор экземпляра программы
        (HINSTANCE)NULL); //отсутствие дополнительных параметров
 
                          //Выводим окно из памяти на экран
    ShowWindow(hWnd, nCmdShow);
    //Обновим содержимое окна
    UpdateWindow(hWnd);
 
    //Цикл обработки сообщений
 
    while (GetMessage(&lpMsg, NULL, 0, 0)) {
        TranslateMessage(&lpMsg);
        DispatchMessage(&lpMsg);
    }
    return(lpMsg.wParam);
}
 
 
//Функция окна
LRESULT CALLBACK WndProc(HWND hWnd, UINT messg,
    WPARAM wParam, LPARAM lParam)
{
    // будут за координаты x и у
    static int cxClient, cyClient;
 
    HDC hdc; //создаём контекст устройства
    PAINTSTRUCT ps; //создаём экземпляр структуры графического вывода
    LOGFONT lf;
    HFONT hFont;
    RECT r;
    HBRUSH hBrush = NULL;   // инициализировал нулом так как студия ругается
    HPEN hPen = NULL;       // типа - использована неинициализированная локальная переменная 
                            // можно вынести в глобальную область - там она сама проинициализируется 
 
    //Цикл обработки сообщений
    switch (messg)
    {
    case WM_SIZE:
        cxClient = LOWORD(lParam);
        cyClient = HIWORD(lParam);
        break;
 
        //сообщение рисования
    case WM_PAINT:
        hdc = BeginPaint(hWnd, &ps);
 
        hBrush = CreateSolidBrush(RGB(10, 200, 100));
        SelectObject(hdc, hBrush);
 
        Dragon(10, hdc, cxClient, cyClient);
 
        
        EndPaint(hWnd, &ps);
        break;
 
        //сообщение выхода - разрушение окна
    case WM_DESTROY:
        PostQuitMessage(0); //Посылаем сообщение выхода с кодом 0 - нормальное завершение
        DeleteObject(hPen);
        DeleteObject(hBrush);
        break;
 
    default:
        return(DefWindowProc(hWnd, messg, wParam, lParam)); //освобождаем очередь приложения от нераспознаных
    }
    return 0;
}

...
    void build(carr<int>& a)
    {   int w = a.column_();
        int h = a.row_   ();
 
        int R = w / 2 - 2 ;
 
        int r  = carr<int>::rrand(  10,  48);
        int dx = carr<int>::rrand(-200, 200);
        int dy = carr<int>::rrand(-200, 200);
 
        ucolor u;
               u.I.i = get_rand_color();
 
        FORi(r)
        {   float X = i - r;
            float Y = sqrt(r * r - i * i);
            a[int(Y) + R + dy][i + R + dx] = u.I.i;
            a[R - int(Y) + dy][i + R + dx] = u.I.i;
            a[int(Y) + R + dy][R - i + dx] = u.I.i;
            a[R - int(Y) + dy][R - i + dx] = u.I.i;
 
            if(u.R.r > 3) u.R.r -= 3;
            if(u.R.g > 3) u.R.g -= 3;
            if(u.R.b > 3) u.R.b -= 3;
        }
    }
...

...
    void duplicator()
    {   const float A12 = PI_2/12;
        const float A6  = PI_2/6;
 
        float step = 0.f;
        for(int i = 0.f; i < 6; i++, step += A6)
        {
            for(float A = 0.f; A < A12; A += 0.002f)
            {
                const float D = a.row_()/2;
                for(float r = 0.f; r < D; r += 0.5f)
                {
                    float aa = A + step;
 
                    int x = r * sinf(A) + D;
                    int y = r * cosf(A) + D;
                    int v = a[y][x];
 
 
                    int fx1  = r * sinf(aa) + D;
                    int fy1  = r * cosf(aa) + D;
                    f[fy1][fx1] = v;
 
                    int fx2  = r * sinf(2 * A12 - aa) + D;
                    int fy2  = r * cosf(2 * A12 - aa) + D;
                    f[fy2][fx2] = v;
                }
            }
        }
    }
...