случайное ветвление

[VitalD]

Здравствуйте все!
Возникла вот проблема, к решению которой не знаю как подступиться. Может быть, кто сталкивался, подскажет, что и как. Так вот, представьте себе, что есть некая программа, в начале выполнения которой пользователь может указать желаемые вероятности исполнения трех различных ветвей внутри содержащегося в теле программы цикла. Например, что-то типа этого:

объявление массива М[0..99] каких-то констант;
объявление какой-то процедуры № 1;
объявление какой-то процедуры № 2;
объявление какой-то процедуры № 3;

ввод параметра Р1 (матожидание для ветви № 1);
ввод параметра Р2 (то же для ветви № 2);
ввод параметра Р3 (то же для № 3);

for і:=0 то 99 do
begin
if <какое-то условие>
then
вызов процедуры № 1 (М[i]) //первая ветвь
еlse
if <еще какое-то условие>
then
вызов процедуры № 2 (М[i]) //вторая ветвь
еlse
вызов процедуры № 3 (М[i]); //третья ветвь
end;

Как видим, отдельные константы из массива М обрабатываются не только по-разному в зависимости от условий, но еще и с вероятностью, не обязательно равной 1 (то есть некоторые могут и вообще остаться нетронутыми, если пользователь установит значение параметра , например, Р2, равным нулю, а другие могут иногда подвергнуться какой-то обработке, а могут и не подвергнуться - случай определяет). Надеюсь, что понятно объяснил задачу.

Заранее огромное спасибо тому или тем, кто хотя бы направление укажет, КАК это можно реализовать.

[BDA]

Не уверен, но можно попробовать:

randomize;
...
a := random();
if a < p1 / (p1 + p2 + p3) then
...
else if a < (p1 + p2) / (p1 + p2 + p3) then
...
else

Мог напутать

[VitalD]

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

Не уверен, но можно попробовать:

randomize;
...
a := random();
if a < p1 / (p1 + p2 + p3) then
...
else if a < (p1 + p2) / (p1 + p2 + p3) then
...
else

Мог напутать

При выполнении а получит значение либо 0, либо 1. Параметры Р1, Р2, Р3 - матожидания выполнения 1-й..3-й ветвей соответственно, могут лежать в ДИАПАЗОНЕ от 0 до 1, например, Р1=0,09, Р2=0,56, Р3=0,88 (не зависят друг от друга совершенно). Если последовать вашему коду, то в величинах <р1/(р1+р2+р3)> и <(p1 + p2)/(p1 + p2 + p3)> не вижу никакого смысла, т.к. условиями выбора той или иной ветви могут быть совершенно другие потребности, например, величина константы из массива М. Да дело и не в этих условиях, вообще-то. Просто выбрали при каком-то условии (каком-то), что для данной константы МОЖЕТ (!) быть вызвана процедура № 2 (например). И вероятность ее вызова = параметру Р2 (например). Пытаюсь вот нарыть что-то в генетических алгоритмах, но пока - глухо.

[VitalD]

Попробую чуток конкретизировать. Предположим, что ВСЕ константы, содержащиеся в массиве М, вполне соответствуют условию, при котором происходит выбор вызова ВТОРОЙ процедуры. Но вызов, не взирая на это, происходит (должен происходить, если точнее) НЕ ВСЕГДА, а согласно матожиданию, равному параметру Р2. Вот, может, так будет более ясным описание проблемы.))

[Vanta11a]

Итак, Вы перебираете все константы

for (i:=1;i<n,i++)

Далее, для каждой константы вы перебираете все возможные процедуры

for (j:=1;j<n,j++)

В каждой из итераций выполняете что-то вроде следующего:

if (rand()%100 > *заданная вероятность выполнения действия j*) then 
...
break;

Т.е. последовательно прогоняете все функции, генерируете случайный % от 100% и сравниваете с заданным % выполнения для этой функции. Если полученный % превышает заданный - функция выполняется, цикл прерывается.

Примеры приведены на С++, но, думаю, алгоритм из этого вывести можно.

[VitalD]

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

if (rand()%100 > *заданная вероятность выполнения действия j*) then 
...
break;

Т.е. последовательно прогоняете все функции, генерируете случайный % от 100% и сравниваете с заданным % выполнения для этой функции. Если полученный % превышает заданный - функция выполняется, цикл прерывается.

Кажется, что-то в этом есть... спасибо, подумаю-попробую и отпишусь сегодня или завтра.

[BDA]

1) random() вернет вещественное число 0<=a<1
2) Приведенный вариант вызывает одну из трех процедур с заданными вероятностями (обязательно будет вызвана одна), но я не учел, что Вы хотите вызвать каждую процедуру с какой-либо вероятностью.
Тогда, можно:

randomize;
...
a := random();
if a < p1 then
...
if a < p2 then
...
if a < p3 then

Что уже и предложил Vanta11a.

[VitalD]

for i2 := -WSY to WSY  do
  for i1 := -WSX to WSX  do
     case C[i1, i2] of
       1: if (Random(100) > 100-StrToFloat(ShanceM.Text)*100)
            then
              ….
       2: if (Random(100) > 100-StrToFloat(ShanceF.Text)*100)
            then
              ….
       3: if (Random(100) > 100-StrToFloat(ShanceC.Text)*100)
            then
              ….
      end

Опробовал - работает. Спасибо за помощь.