Random с неоднородной вероятностью выпадения событий

[Vago]

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

Давайте посмотрим в задание автора:

Давайте.

Вот Вам ТЗ:

Написать датчик случайных чисел, генерирующий случайную величину, являющуюся натуральным числом из диапазона [1,6], вероятности появления которой определяются так:
1 - 0,08
2 - 0,08
3 - 0,32
4 - 0,32
5 - 0,10
6 - 0,10

Я, как видите, сократил количество интервалов до трёх, оставив в каждом интервале по 2 равновероятных значения, но суть задачи - сохранил. С громадным интересом посмотрю на Ваш код! Естественно, ничего против использования внутри него ф-ии rand() не имею.

[EC.cpp]

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

Давайте.

Вот Вам ТЗ:

Написать датчик случайных чисел, генерирующий случайную величину, являющуюся натуральным числом из диапазона [1,6], вероятности появления которой определяются так:
1 - 0,08
2 - 0,08
3 - 0,32
4 - 0,32
5 - 0,10
6 - 0,10

Я, как видите, сократил количество интервалов до трёх, оставив в каждом интервале по 2 равновероятных значения, но суть задачи - сохранил. С громадным интересом посмотрю на Ваш код! Естественно, ничего против использования внутри него ф-ии rand() не имею.

Ну собственно на первый взгляд ничего трудного нет.
Попробую кратко описать: из исходных вероятностей нам надо сделать 100 бросков чтобы получить соответствующие распределения значений, т.е. известно заранее что после 100 бросков: в 8 случаях выпадет 1, в 32 случаях выпадет 3 и т.д.

//Надо описать структуру для хранения данных:

struct Card
{
	unsigned int value;		//Ценность элемента, в данном случае принимает значение от 1-6
	Card* next;			//Указатель на след элемент, односвязныый список
}

//Описать методы реализующие стандартную работу с односвязным списком:

discInfo* add_elem (...);	//Функция добавления элемента списка
void show_elem ();			//Процедура вывода одного элемента списка на консоль
void shake_me_up(...);		//Процедура случайного перемешивания списка

Последняя процедура самая интересная, в ней то и будем смешивать нашу "Колоду" чисел от 1-6. Алгоритм может выглядеть приблизительно так:

	unsigned int Unsort_cards = 100;	//Текущее количество неперемешанных карт
	unsigned int i;						//Просто счетчик
	unsigned int rnd;					//Случайная величина, указывающая на номер смешиваемой в текущий момент карты
	Card* Base;							//Заголовок списка

	for (i=99; i > -1; i--)
	{
		rnd = rand() %i;				//Получаем номер перемещаемого в конец списка элемента структуры
		{...}							//Обработка списка относительно смещенной карты(величины)
	}

	ShowList();							//Показать уже смешанную колоду чисел с заданными вероятностями выпадения

[Vago]

Уважаемый (? ) EC.cpp! Повторяю, я буду дико благодарен, если увижу здесь код. Не "на первый взгляд", не "приблизительно", не обрывки, а полнофункциональный КОД, обеспечивающий то задание, которое я в #11 написал! Чтобы я вызвал Ваш датчик тыщу раз, построил гистограммку, увидел три ступеньки, удовлетворённо хмыкнул и сказал бы себе "О! С этим парнем (? ) можно поговорить дальше!.."

И ещё. Одним из требований, с необходимостью предъявляемых к Д(П)СЧ, является быстродействие. Ибо в реальных задачах моделирования этот датчик вызывается миллиарды и миллиарды раз. Поэтому использовать в ПО этого типа связанные списки - ну, как бы это помягче выразиться... непрофессионально. Впрочем, ладно. Я готов посмотреть на код, использующий внутри любые структуры данных, лишь бы он соответствовал моему ТЗ.

Заранее спасибо.

[EC.cpp]

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

... а полнофункциональный КОД, обеспечивающий то задание, которое я в #11 написал! ...

Хм, знаете вот что было плохо в Советском Союзе? :-)
Говорят - отсутствие стимула.
Так вот собственно чтобы качественно реализовать и оттестировать это в коде мне понадобится около дня - зачем мне этим заниматься?
Я Вам предоставил алгоритм - ответ на вопрос КАК.
А реализовывать это в коде любого языка программирования - работа программиста

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

И ещё. Одним из требований, с необходимостью предъявляемых к Д(П)СЧ, является быстродействие. Ибо в реальных задачах моделирования этот датчик вызывается миллиарды и миллиарды раз. Поэтому использовать в ПО этого типа связанные списки - ну, как бы это помягче выразиться... непрофессионально. Впрочем, ладно. Я готов посмотреть на код, использующий внутри любые структуры данных, лишь бы он соответствовал моему ТЗ.

Заранее спасибо.

Списки - один и самых эффективных и быстрых методов для реализации задач подобного типа. Если Вы знаете ( знаете? :-) ) какие-либо другие, поделитесь пожалуйста с окружающими.

[Сtrl]

Вот Вам ТЗ:

Написать датчик случайных чисел, генерирующий случайную величину, являющуюся натуральным числом из диапазона [1,6], вероятности появления которой определяются так:
1 - 0,08
2 - 0,08
3 - 0,32
4 - 0,32
5 - 0,10
6 - 0,10

int getrand(void)
{
	if (!(rand()%5)) //события выполняются с вероятностью 1:5, то есть 0,20
	{
		return rand()%2+5; //	либо 0+5, либо 1+5, то есть возвращает 5 или 6
	}					   //	вероятность выпадения 0,20 / 2 = 0,10, как и заказывали
	else //в этой ветке шансы уже равны 0,80
	{
		if (!(rand()%5)) //события выполняются с вероятностью 0,80 * (1:5) = 0,16
		{
			return rand()%2+1; // либо 0+1, либо 1+1, то есть возвращает 1 или 2
		}					   // вероятность выпадения 0,16 / 2 = 0,08, как и заказывали
		else //в этой ветке шансы равны 0,80-0,16 = 0,64
		{
			return rand()%2+3; // возвращает либо 0+3, либо 1+3, то есть 3 или 4
		}					   // вероятность выпадения 0,64 / 2 = 0,32, как и заказывали
	}
}

int main()
{
	srand(time(0));
	printf("rand nums\n");
	for (int i = 0; i < 10000; i++)
		printf("%d\t", getrand());
	getchar();
}

Здесь правда я сильно обобщил - лучше для получения вероятностей и самих чисел использовать разные генераторы сл.чис., потому что иначе нормальный разброс не получится.

[EC.cpp]

Вероятность выпадания значений фукции rand() не поддается параметризации или систематизации, т.е. нельзя сказать с какой вероятностью какое значение выпадет - это функция генерации воистину случайного числа.
Чуть подробнее о методах порождения случайных величин написано тут (ток на английском)

Вышесказанное означает, что Вы обычным применением метода rand() никогда не получите заданное распределение вероятностей.

[Сtrl]

это функция генерации воистину случайного числа.

На это мой алгоритм и рассчитан, проблема одна - на самом-то деле оно, как ни крути, псевдослучайное.

[EC.cpp]

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

На это мой алгоритм и рассчитан, проблема одна - на самом-то деле оно, как ни крути, псевдослучайное.

Вот именно поэтому, я предлагаю использовать этот самый rand() только для выборки 1-го из 100 чисел с последующий его перестановкой в списке - в результате мы получим перемешанную последовательность, но при выводе их на экран вероятности будут соответствовать строго заданным параметрам.

[Vago]

Сообщение от EC.cpp

знаете вот что было плохо в Советском Союзе? :-)

Знаю Но не об этом сейчас речь...

Сообщение от EC.cpp

зачем мне этим заниматься?

Незачем. Если это - не является предметом Вашей профессиональной специализации (так что подходящие исходники - всегда под рукой), то, кроме удовлетворённого тщеславия, - незачем. Но эти Ваши слова только подтверждают мою реакцию на слово "просто", прозвучавшее (не из Ваших, правда, уст) в этом обсуждении.

Сообщение от EC.cpp

Я Вам предоставил алгоритм - ответ на вопрос КАК.

Да?

Сообщение от EC.cpp

Если Вы знаете ( знаете? :-) ) какие-либо другие, поделитесь пожалуйста с окружающими.

Хорошо. Я вечером в завалах посмотрю. Если окажется, что они остались на трёхдюймовых дискетах - тогда завтра.

[Vago]

Уважаемый Сtrl, я использовал Вашу программку из #15, вызывая её вот таким образом.

int main() {
	
   double f[6] ;
   const int nSeeds = 100000 ;
   int   k ;

   srand(time(0)) ;
   printf("rand nums\n");

   for ( k = 0; k < 6; k++ )
      f[k] = 0. ;

   for (int i = 0; i < nSeeds; i++) {
      k = getrand() ;
      f[k-1] += 1. / double( nSeeds ) ;
   }

   for ( k = 0; k < 6; k++ ) 
      printf( "%1d - %.3f\n", k+1, f[k] ) ;

   getchar();
   return 0 ;

}

И получил вот что:
110410.jpg

Ну ладно. До вечера!