Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
 
Рейтинг: Рейтинг темы: голосов - 14, средняя оценка - 5.00
Woody-krsk
1 / 1 / 0
Регистрация: 20.12.2010
Сообщений: 62
#1

Какая-то ошибка с памятью, где накосячил? - C++

18.08.2011, 18:21. Просмотров 1841. Ответов 24
Метки нет (Все метки)

Всех приветствую! Собственно программа - зачаток простого генетического алгоритма. Что делает: создается популяция с заданным количеством индивидов, которые есть бинарные строки, нужно найти строчку с максимальным количеством переходов 1-0.
1) создаем популяцию и заполняем случайно;
2) вычисляем пригодность;
3) проверяем нет ли в популяции нужного решения (если есть, то все);
4) селекция равномерная, методом "рулетки". Заполняем рулетку с зависимости от пригодности индивида (чем больше пригодность, тем больше ячеек на поле рулетки он получает);
5) Выполняется скрещивание, выбирая индивидов на которых указало "колесо рулетки"
6) применяем оператор мутации
7) повторяем шаги 2-6 заданное число итераций.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
#include <iostream>
#include <cstdlib>
#include <stdlib.h>
#include <conio.h>
#include <time.h>
using namespace std;
 
void showmass(int**, int, int);
void showmass(int*, int);
void showmass(double*, int);
int fitness (int*, int);
int fitness_1 (int*, int);
void cross(int*&, int*&, int);
void cross(int*&, int*&, int*&, int);
void mutation(int*& gen, int len, int ver);
int round (double);
int pow (int, unsigned int);
int main()
{
    srand(time(NULL));
    setlocale(0,"");
    int dm1 = 0, dm2 = 10, itr, mut;
    cout << "Введите размер популяции" << endl;
    cin >> dm1;
    cout << endl;
    cout << "Введите число прогонов" << endl;
    cin >> itr;
    cout << endl;
    cout << "Введите вероятность мутаци" << endl;
    cin >> mut;
    cout << endl;
    int** popul = new int*[dm1];//создается популяция
    for (int i = 0; i < dm1; i++){
        popul[i] = new int[dm2];
    }
    for (int i = 0; i < dm1; i++){
        for (int j = 0; j < dm2; j++){
            popul[i][j] = rand()%2;
        }   
    }
    int* fitn = new int[dm1];//массив для пригодности каждого индивида
for (int gener = 0; gener < itr; gener++){
    int av_fit = 0, best = 0, best_pos = 0;
           //------------------------------------
           //Вычисляем пригодность каждого в популяции
           //---------------------------------------
    for (int i = 0; i < dm1; i++){
        fitn[i] = fitness_1(popul[i], dm2);
        av_fit= av_fit + fitn[i];
        if(best<fitn[i]){
            best=fitn[i];
            best_pos = i;
        }
    }
    cout << gener+1 <<" поколение" << endl;
    showmass(popul, dm1, dm2);
    cout << endl;
    cout << "Его пригодность" << endl;
    showmass(fitn, dm1);
    cout << endl;   
    //---------------
    //TEST (Проверка не нашлось ли решение)
    //---------------
    if(best==dm2/2){
        cout << "DONE!" << endl << "Find in " << gener+1 << " iteration" << endl;
        cout << "Индивид номер " << best_pos+1 << endl;
        break;
    }
    //------------------------------
    //Селекция
    //------------------------------
    int* possible = new int[100];//массив "рулетка"
    for (int i = 0; i < 100; i++){
        possible[i]=0;
    }
    int pos = 0;    
    for (int i = 0; i < dm1; i++){
        double p = 0;
        p = (double)fitn[i]/av_fit;
        p = round(p);
                      //----------------------------
                      //Крайнему отдадим остатки от поля рулетки
                      //-----------------------------
        if(i == (dm1-1)){
            for(int j = pos; j < 100; j++){
            possible[j] = i;
            }
        } else {//здесь заполняем поля рулетки в зависимости пригодности
                for(int j = pos; j < (pos + p); j++){
                possible[j] = i;
            }
        }
        pos += p;
    }   
    cout << "Рулетка" << endl;
    showmass(possible,100); 
    //----------------------------------
    //Промежуточная популяция (из нее будем брать индивидов на скрещивание)
    //----------------------------------
    int** inter_popul = new int*[dm1];
    for (int i = 0; i < dm1; i++){
        inter_popul[i] = new int[dm2];
    }
    for (int i = 0; i < dm1; i++){
        for (int j = 0; j < dm2; j++){
            inter_popul[i][j] = popul[i][j];
        }   
    }
    cout << "Промежуточное поколение" << endl;
    showmass(inter_popul, dm1, dm2);
    cout << endl;
    //-----------------------------------
    //Скрещивание
    //-----------------------------------
    for(int i = 0; i < dm1; i++){
        cross(inter_popul[possible[rand()%100]],inter_popul[possible[rand()%100]],popul[i], dm2);
    }
    //-----------------------
    //Мутация
    //-----------------------
    for(int i = 0; i < dm1; i++){
        mutation(popul[i],dm2,mut);
    }
    for (int i = 0; i < dm1; i++){
        delete[] inter_popul[i];    
    }
    delete[] inter_popul;
}
    cout << "Последняя популяция" << endl;
    showmass(popul, dm1, dm2);
    _getch();
 
    for (int i = 0; i < dm1; i++){
        delete[] popul[i];  
    }
    delete[] popul;
    delete[] fitn;
    return 0;
}
int fitness (int* ind, int lng){
    int fitness = 0;
    for (int i = 0; i < lng; i++){
        fitness = fitness + ind[i];
    }
    return fitness;
}
int fitness_1 (int* ind, int lng){
    int fitness = 0;
    for(int i = 0; i < lng; i++){
        if(ind[i]==1){
            if(ind[i+1]==0){
            fitness++;
            }
        }
    }
    return fitness;
}
void showmass(int** mass, int dm1, int dm2){
    for (int i = 0; i < dm1; i++){
        for (int j = 0; j < dm2; j++){
            cout << mass[i][j] << " ";
        }
    cout << endl;
    }
}
void showmass(int* mass, int dm1){
    for (int i = 0; i < dm1; i++){
        cout << mass[i]<< " ";
        
    }
    cout << endl;
}
void showmass(double* mass, int dm1){
    for (int i = 0; i < dm1; i++){
        cout << mass[i]<< " ";
        
    }
    cout << endl;
}
void cross(int*& gen1, int*& gen2, int len){
    int pos = rand()%len;
    cout << pos << endl;
    int buff = 0;
    for(int i = 0; i < (len-pos); i++){
       buff = gen1[i];
       gen1[i] = gen2[i];
       gen2[i] = buff;
   }
}
void cross(int*& gen1, int*& gen2, int*& gen3, int len){
    int pos = rand()%len;
    for(int i = 0; i < (len-pos); i++){
       gen3[i] = gen1[i];
   }
    for(int i = pos; i < len; i++){
       gen3[i] = gen2[i];
   }
}
//---------------------------------------------------------------------------
//Мутация
//---------------------------------------------------------------------------
void mutation(int*& gen, int len, int ver){
    int i = rand()%len;
    int p = rand()%100;
    if(p < ver){
        switch (gen[i]){
            case 0: {
                gen[i] = 1;
                break;
            }
            case 1: {
                gen[i] = 0;
            }
        }
    }
 
}
int round (double digit){
    int integer = 0;
    double diff = 0;
    int ex = 100000;
    for (int i = 0; i < 3; i++){
        digit *= ex;
        integer = digit;
        diff = digit - integer;
        if(diff > .5){
            integer++;
        }
        digit = (double)integer/ex;
        ex /= pow(10,i+1);
    }
    return integer;
}
int pow (int x, unsigned int y){
    int result = x;
    switch (y){
        case 0 : {
            x = 1;
            break;
        }
        case 1: {
            break;
        }
        default: {
            for(int i = 1; i < y; i++){
                result*=x;
            }
        }
    }
    return result;
}
Программа вроде как работает, но если взять большую популяцию (от 25 и выше), вылетает какая-то совсем непонятная мне ошибка:
Необработанное исключение в "0x7575b9bc" в "listen.exe": Исключение Microsoft C++: std::bad_alloc по адресу 0x0015f7b0..
которая отсылает меня на домашнем компьютере (win7) на файл mlock.c, на рабочем ноуте (winxp) постоянно было new.cpp. Пишу из дома, указывает мне вот на эту функцию в файле (зеленая стрелочка показывает по крайнюю фигурную скобку):
C++
1
2
3
4
5
6
7
8
9
void __cdecl _unlock (
        int locknum
        )
{
        /*
         * leave the critical section.
         */
        LeaveCriticalSection( _locktable[locknum].lock );
}
Спросил у товарища, который мается той же фигней, сказал ищи косяки мелкие, где-то вылазишь за границу массива, все излазил - не нашел.
Путем установки точек останова выяснил, что косяк, возможно, возникает где то в блоке с селекцией, но что там не так понять не могу.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
18.08.2011, 18:21
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Какая-то ошибка с памятью, где накосячил? (C++):

Задача про ГАИ и файлы (где накосячил в коде?) - C++
Добрый день! Задача звучит так: Вдоль шоссе в точках X1,X2,...,XN расположены посты ГАИ. В точке X произошло мелкое ДТП...

Трай , где неправильная работа с памятью? - C++
class Trie { struct Node { int count ; Node **r; Node() { count = 0; ...

Ошибка с памятью - C++
Добрый день, Уважаемые форумчане! Следующий код заточен под параллельный алгоритм с использованием библиотеки mpi. Сам...

Ошибка с памятью - C++
Здравствуйте! Пытаюсь преобразовать картинки одного из векторов в серые. (имеются 2 вектора с одинаковыми изображениями) ...

Ошибка в работе с памятью - C++
После использования программы и последующего выхода из нее возникает ошибка HEAP CORRUPTION DETECTED. Хотелось бы узнать, в чем причина и...

Ошибка с динамической памятью - C++
Очень прошу помочь исправить ошибку в области дин.памяти. Проблема: 1. Шаг Создаю массив для записи 2. Шаг удаляю элемент под...

24
Deviaphan
Делаю внезапно и красиво
Эксперт С++
1306 / 1221 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
19.08.2011, 16:06 #16
Я конкретно про функцию round спрашиваю. До ближайшего целого округлять должна или что это вообще? Не понимаю.(

Добавлено через 45 секунд
Цитата Сообщение от Woody-krsk Посмотреть сообщение
округляя это число 0,012345 -> 12
Ну не округление это.)))
0
Woody-krsk
1 / 1 / 0
Регистрация: 20.12.2010
Сообщений: 62
19.08.2011, 16:07  [ТС] #17
Deviaphan, ааа, округляет до ближайшего двузначного целого грубо говоря. В общем вытаскивает первые 3 старших разряда из дробной части в целую.

Добавлено через 16 секунд
Deviaphan, а как это назвать? Пусть будет увеличение порядка с округлением, например
0
Deviaphan
Делаю внезапно и красиво
Эксперт С++
1306 / 1221 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
19.08.2011, 16:09 #18
То же самое делает.
C++
1
2
3
4
5
6
int round (double digit)
{
      digit *= 1e3;
      digit += 0.5;
      return (int)floor(digit);
}
Добавлено через 58 секунд
Цитата Сообщение от Woody-krsk Посмотреть сообщение
ааа, округляет до ближайшего двузначного целого грубо говоря. В общем вытаскивает первые 3 старших разряда из дробной части в целую.
Ближайшее целое для этого будет ноль.
В общем, название не правильное. Функционал чуть выше.
0
Woody-krsk
1 / 1 / 0
Регистрация: 20.12.2010
Сообщений: 62
19.08.2011, 16:10  [ТС] #19
Deviaphan, ой можно по косточкам пожалуйста? Что делает floor, в какой оно библиотеке, и что значет 1е3?
0
Deviaphan
Делаю внезапно и красиво
Эксперт С++
1306 / 1221 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
19.08.2011, 16:13 #20
Цитата Сообщение от Woody-krsk Посмотреть сообщение
1е3
Это 1000

Цитата Сообщение от Woody-krsk Посмотреть сообщение
floor
Округляет до меньшего целого. Есть ещё ceil. MSDN

Т.к. округляет до меньшего, а не до ближайшего, то добавляем 0,5

Можно и без floor, сразу в int приводить. Я сперва думал, что double должен возвращаться.
0
Woody-krsk
1 / 1 / 0
Регистрация: 20.12.2010
Сообщений: 62
19.08.2011, 16:21  [ТС] #21
Цитата Сообщение от Deviaphan Посмотреть сообщение
MSDN
Это библиотека такая что ли? С правке MSDN вбил
Returns the largest integer less than or equal to the given numeric expression.

Текст
Syntax
FLOOR ( numeric_expression )

Arguments
numeric_expression

Is an expression of the exact numeric or approximate numeric data type category, except for the bit data type.

Return Types
Returns the same type as numeric_expression.

Examples
This example shows positive numeric, negative numeric, and currency values with the FLOOR function.

SELECT FLOOR(123.45), FLOOR(-123.45), FLOOR($123.45)
The result is the integer portion of the calculated value in the same data type as numeric_expression.

--------- --------- -----------
123 -124 123.0000

Что делает понятно, что надо подключить ни слова, а с ней код не компилируется
В любом случае проблему с погрешностью такой вариант не решает.
0
Deviaphan
Делаю внезапно и красиво
Эксперт С++
1306 / 1221 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
19.08.2011, 16:27 #22
Ну так ты справку по С++, а не по SQL смотри.)))
Значит, проблема не в погрешности, а в неправильном алгоритме.
0
Woody-krsk
1 / 1 / 0
Регистрация: 20.12.2010
Сообщений: 62
19.08.2011, 16:38  [ТС] #23
Deviaphan, вроде делал поиск по с++
Да я уже тоже понял, что это тупой способ получение количества ячеек, другого пока в голову не приходит, буду думать. Если есть желание можно какую нибудь идейку подкинуть
0
Deviaphan
Делаю внезапно и красиво
Эксперт С++
1306 / 1221 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
19.08.2011, 16:48 #24
Измени алгоритм заполнения вектора пригодности.
Сперва вычисли количество ечеек для каждого индивида (пригодность одного вида в процентах от всех пригодностей). А потом заполни массив possible в соответствии с известным распределением.

Добавлено через 1 минуту
По С++ вот так хэлп выглядит.)
http://msdn.microsoft.com/en-us/library/x39715t6%28v=vs.71%29.aspx
0
Woody-krsk
1 / 1 / 0
Регистрация: 20.12.2010
Сообщений: 62
19.08.2011, 17:20  [ТС] #25
Deviaphan, хех, что-то подобное в голову и пришло, буду пробовать так, спасибо огромное за подсказку. Не прощаюсь, вопросов будет еще миллион, я уверен
0
19.08.2011, 17:20
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
19.08.2011, 17:20
Привет! Вот еще темы с ответами:

Ошибка при работе с памятью - C++
Есть класс, в котором я предпринимаю безрезультатные попытки работать с динамической памятью из-за того, что экземпляров класса может быть...

Ошибка при работе с памятью - C++
Написан следующий код: #include &quot;stdafx.h&quot; #include &lt;iostream&gt; #include &lt;stdio.h&gt; #include &lt;math.h&gt; #include &lt;cstdlib&gt; #include...

Ошибка при работе с памятью - C++
Если ввести например 6, 0, чтобы выбросить исключение. То по завершение программы выскакивает ошибка. В чём может быть дело? Так что то про...

При перемножении матриц ошибка с памятью - C++
Написал класс матрица с алгоритмом умножения, но при выходе из проги выбивает ошибку. Не пойму как написать деструктор. Подскажите как...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
25
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru