UVALive 6346 Wormtongue's Mind

题意：x 是０~１的随机变量，问最后表达式的期望

变量 x 的概率分布为　F(x)=x,

H(x)=max(F(x),G(x))

则 H(x)=F(x)*G(x)

 

H(x)=min(F(x),G(x))

则 H(x)=1-(1-F(x))*(1-G(x))

用多项式表示

我代码是 java 的

用 double 开始没过，就用BigDecimal
 

import java.math.BigDecimal;
import java.util.Scanner;


public class Main{
	static String a;
	static int num;
	static int maxn=105;
	static long[] work()
	{
			num++;
		if(a.charAt(num)=='x')
		{
			long[] z=new long[maxn];
			z[1]=1;
			return z;
		}
		if(a.charAt(num)=='M')
		{
			long[] z=work();
			long[] x=work();
			long[] c=new long[maxn];
			for(int q=0;q<maxn;q++)
				for(int w=0;w+q<maxn;w++)
					c[q+w]+=z[q]*x[w];
			return c;
		}
		else
		{
			long[] z=work();
			long[] x=work();
			long[] c=new long[maxn];
			for(int q=0;q<maxn;q++)
				for(int w=0;w+q<maxn;w++)
					c[q+w]-=z[q]*x[w];
			for(int q=0;q<maxn;q++)
				c[q]+=z[q]+x[q];
			return c;
		}
	}
	public static void main(String[] as)
	{
		Scanner scan=new Scanner(System.in);
		int T=scan.nextInt();
		for(int idx=0;idx<T;idx++)
		{
			a=scan.next();
			num=-1;
			long[] z=work();
			double ans=0;
			BigDecimal b=BigDecimal.ZERO,v,c;
			for(int q=1;q<maxn;q++)
			{
				ans+=((double)z[q])*q/(q+1);
				v=new BigDecimal(z[q]*q+"");
				c=v.divide(new BigDecimal(q+1+""),30,BigDecimal.ROUND_HALF_UP);
				b=b.add(c);
				
			}
			System.out.println(b);
			
			
			
			
			
			
			
			
			
			
			
			
		}
	}
}