uva11181 - Probability|Given 条件概率

Problem G
Probability|Given
Input:
Standard Input

Output: Standard Output

 

N friends go to the local super market together. The probability of their buying something from the market is respectively. After their marketing is finished you are given the information that exactly r of them has bought something and others have bought nothing. Given this information you will have to find their individual buying probability.

 

Input

The input file contains at most 50 sets of inputs. The description of each set is given below:

 

First line of each set contains two integers N (1 ≤ N ≤ 20) and r(0 ≤ r ≤ N). Meaning of N and r are given in the problem statement. Each of the next N lines contains one floating-point number  (0.1< <1) which actually denotes the buying probability of the i-th friend. All probability values should have at most two digits after the decimal point. 

 

Input is terminated by a case where the value of N and r is zero. This case should not be processes.  

 

Output
For each line of input produce N+1 lines of output. First line contains the serial of output. Each of the next N lines contains a floating-point number which denotes the buying probability of the i-th friend given that exactly r has bought something. These values should have six digits after the decimal point. Follow the exact format shown in output for sample input. Small precision errors will be allowed. For reasonable precision level use double precision floating-point numbers.

 

Sample Input                             Output for Sample Input

3 2
0.10
0.20
0.30
5 1
0.10
0.10
0.10
0.10
0.10
0 0

Case 1:

0.413043

0.739130

0.847826

Case 2:

0.200000

0.200000

0.200000

0.200000

0.200000


  给你每个人做一件事的概率,告诉你有几个人做了这件事,求每个人在这个条件下做这件事的概率。

  哎~~概率论~P(A|B) = P(AB)/P(B)

  B就是N个人里面有r人做这件事。P(B)怎么算?因为N不大,用DFS来。DFS(int cur,int r,double pi),cur是当前该判断第cur个人的状态,r是当前已经选了r个人,pi是当前已经确定的选择的概率。注意不能到r==R就返回,因为前面选到了R后面就不能再选,要到cur>N时判断r是否等于R。

#include<cstring>
#include<cstdio>
#include<iostream>
#include<climits>
#include<cmath>
#include<algorithm>
#include<queue>
#define INF 0x3f3f3f3f
using namespace std;
int N,R,cas=0;
double a[25],p[25];
double DFS(int cur,int r,double pi){
    if(cur>N) return r==R?pi:0;
    double s;
    s=DFS(cur+1,r+1,pi*p[cur]);
    a[cur]+=s;
    s+=DFS(cur+1,r,pi*(1-p[cur]));
    return s;
}
int main(){
    freopen("in.txt","r",stdin);
    while(scanf("%d%d",&N,&R),N||R){
        int i;
        for(i=1;i<=N;i++) scanf("%lf",&p[i]);
        memset(a,0,sizeof(a));
        double s=DFS(1,0,1);
        printf("Case %d:\n",++cas);
        for(i=1;i<=N;i++) printf("%.6lf\n",a[i]/s);
    }
    return 0;
}




要使用遗传算法求解函数的最小值,可以按照以下步骤进行编写代码: 1. 定义适应度函数:根据给定的函数表达式,计算适应度函数的值。适应度函数的值越小,表示个体的适应度越高。 ```python def fitness_function(x): x1, x2, x3, x4, x5, x6, x7 = x y = 100 / (abs(x1 + 1) + abs(x2) + abs(x3 - 1) + abs(x4 + 2) + abs(x5 + 3) + abs(x6 - 2) + abs(x7 - 3) + 1) return y ``` 2. 初始化种群:随机生成一组个体作为初始种群。 ```python import random def initialize_population(population_size, chromosome_length): population = [] for _ in range(population_size): individual = [random.uniform(-10, 10) for _ in range(chromosome_length)] population.append(individual) return population ``` 3. 选择操作:使用轮盘赌选择算法,根据个体的适应度选择一部分个体作为下一代的父代。 ```python def selection(population, fitness_values): total_fitness = sum(fitness_values) probabilities = [fitness / total_fitness for fitness in fitness_values] selected_population = [] for _ in range(len(population)): r = random.random() cumulative_probability = 0 for i, probability in enumerate(probabilities): cumulative_probability += probability if r <= cumulative_probability: selected_population.append(population[i]) break return selected_population ``` 4. 交叉操作:使用单点交叉算法,对父代个体进行交叉操作生成子代个体。 ```python def crossover(parent1, parent2): crossover_point = random.randint(1, len(parent1) - 1) child1 = parent1[:crossover_point] + parent2[crossover_point:] child2 = parent2[:crossover_point] + parent1[crossover_point:] return child1, child2 ``` 5. 变异操作:对子代个体进行变异操作,增加种群的多样性。 ```python def mutation(individual, mutation_rate): mutated_individual = individual.copy() for i in range(len(mutated_individual)): if random.random() < mutation_rate: mutated_individual[i] = random.uniform(-10, 10) return mutated_individual ``` 6. 遗传算法主程序:根据上述操作,进行遗传算法的迭代过程。 ```python def genetic_algorithm(population_size, chromosome_length, generations, mutation_rate): population = initialize_population(population_size, chromosome_length) for _ in range(generations): fitness_values = [fitness_function(individual) for individual in population] selected_population = selection(population, fitness_values) new_population = [] while len(new_population) < population_size: parent1 = random.choice(selected_population) parent2 = random.choice(selected_population) child1, child2 = crossover(parent1, parent2) mutated_child1 = mutation(child1, mutation_rate) mutated_child2 = mutation(child2, mutation_rate) new_population.extend([mutated_child1, mutated_child2]) population = new_population best_individual = max(population, key=fitness_function) best_fitness = fitness_function(best_individual) return best_individual, best_fitness ``` 7. 调用遗传算法函数并输出结果。 ```python population_size = 100 chromosome_length = 7 generations = 100 mutation_rate = 0.01 best_individual, best_fitness = genetic_algorithm(population_size, chromosome_length, generations, mutation_rate) print("最优个体:", best_individual) print("最优适应度:", best_fitness) ``` 这样就可以使用遗传算法求解给定函数的最小值了。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值