Генератор псевдослучайных чисел

@melanu · Регистрация: 10.03.2012

Данный генератор работает по конгруэнтному методу, формула X_{n} = ( a X_{n-1}^3 + b X_{n-1}^2 + c X_{n-1} + d )\mod 2^e.
При компиляции выдаётся следующая ошибка:
[C++ Error] Unit1.cpp(34): E2060 Illegal use of floating point
Поискав в интернете номер ошибки, выяснил что проблема в %m. Но как её исправить так и не понял

C++

 float x[1000];
        int a=0, b=0, c=0, d=0, m=0, i=0 ,y[1000];
        a = StrToInt(Edit1->Text);
        b = StrToInt(Edit2->Text);
        c = StrToInt(Edit3->Text);
        d = StrToInt(Edit4->Text);
        m = StrToInt(Edit5->Text);
        y[0] = StrToInt(Edit6->Text);
        for(i = 1; i<=1000; i++) {
                y[i] = ( (a*pow(y[i-1],3) + b*pow(y[i-1],2) + c*(y[i-1])) + d ) % m;
                x[i] = (y[i]) / double(m);
                Memo1->Lines->Add(FloatToStrF(x[i],ffFixed,8,2));
        }

@~~IrineK~~ · 22.02.2015, 15:36

А если так:

C++

1	y[i] = (float)((int)( (apow(y[i-1],3) + bpow(y[i-1],2) + c*(y[i-1])) + d ) % m);

@melanu · 22.02.2015, 15:46 **[ТС]**

Так программа запускается, но при вводе тестовых значений выдаёт не корректные результаты!
Пример: а=8,б=0,с=0,d=65,m=100,y[0]=2; Результат x[1] = 0,29. Хотя посчитав на калькуляторе должен быть x[1] = 0,81.
y[i] = (int)((a * y[i-1]) + d) % m; - а вот эта похожая функция выдаёт корректные результаты! Не могу понять в чём проблема

@~~IrineK~~ · 22.02.2015, 16:16

Сообщение от melanu

Хотя посчитав на калькуляторе

На калькуляторе есть деление нацело?

@melanu · 22.02.2015, 16:19 **[ТС]**

Ещё раз всё перепроверил, вроде бы правильно считает. Можете объяснить, (float)((int)выражение) вот эту часть? Не совсем понимаю приведение типов

@~~IrineK~~ · 22.02.2015, 16:21

У вас y[i] - интовые, поэтому

C++

1	y[i] = (int)((a * y[i-1]) + d) % m

это правильно.
Перевод по float не нужен.

Добавлено через 1 минуту
А вот для x[i] - нужен.

@melanu · 22.02.2015, 16:23 **[ТС]**

Т.е. y[i] целочисленные, а для вычисления x[i] только float т.к. деление, поэтому ошибка я правильно понял?

@~~IrineK~~ · 22.02.2015, 16:31

pow не нужен. Тогда и приведение типов не нужно

C++

1	y[i] = (ay[i-1]y[i-1]y[i-1] + by[i-1]y[i-1] + cy[i-1] + d ) % m;

По x[i] - вы объявили массив как float, а в строке 11 пытаетесь получить double. Вы уж как-то определитесь.
Скажем, так:

C++

1
2
3

double x[1000];
//...
x[i] =  1.0* y[i] / m;

@melanu · 23.02.2015, 00:38 **[ТС]**

Исправил на 1 тип, спасибо!

Добавлено через 8 часов 2 минуты
Возникла 1 проблема, при больших значениях констант а b c d m y[0] выдаются некорректные результаты(отрицательные числа), подсчитав числа на калькуляторе получились значения в триллионах, а максимальное число типа int может быть 2 147 483 647. Как можно решить данную проблему?

@~~IrineK~~ · 23.02.2015, 09:54

long long int

@melanu 1 / 1 / 0 Регистрация: 10.03.2012 Сообщений: 120
	22.02.2015, 15:46 [ТС]	3
	Так программа запускается, но при вводе тестовых значений выдаёт не корректные результаты! Пример: а=8,б=0,с=0,d=65,m=100,y[0]=2; Результат x[1] = 0,29. Хотя посчитав на калькуляторе должен быть x[1] = 0,81. y[i] = (int)((a * y[i-1]) + d) % m; - а вот эта похожая функция выдаёт корректные результаты! Не могу понять в чём проблема 0

@~~IrineK~~ Заблокирован
	22.02.2015, 16:16	4
	Сообщение от melanu Хотя посчитав на калькуляторе На калькуляторе есть деление нацело? 1

@melanu 1 / 1 / 0 Регистрация: 10.03.2012 Сообщений: 120
	22.02.2015, 16:19 [ТС]	5
	Ещё раз всё перепроверил, вроде бы правильно считает. Можете объяснить, (float)((int)выражение) вот эту часть? Не совсем понимаю приведение типов 0

@melanu 1 / 1 / 0 Регистрация: 10.03.2012 Сообщений: 120
	22.02.2015, 16:23 [ТС]	7
	Т.е. y[i] целочисленные, а для вычисления x[i] только float т.к. деление, поэтому ошибка я правильно понял? 0

@melanu 1 / 1 / 0 Регистрация: 10.03.2012 Сообщений: 120
	23.02.2015, 00:38 [ТС]	9
	Исправил на 1 тип, спасибо! Добавлено через 8 часов 2 минуты Возникла 1 проблема, при больших значениях констант а b c d m y[0] выдаются некорректные результаты(отрицательные числа), подсчитав числа на калькуляторе получились значения в триллионах, а максимальное число типа int может быть 2 147 483 647. Как можно решить данную проблему? 0

@~~IrineK~~ Заблокирован
	23.02.2015, 09:54	10
	long long int 0

Опции темы