Content-Type : Application/json

		static double[] Slow(double[] pb, double[] pc)
		{
			int n = pb.Length;
			double[] ev = new double[n];

			for (int a = 0; a < n; a++)
			{
				for (int b = 0; b < n; b++)
				{
					for (int c = 0; c < n; c++)
					{
						if (a == b || a == c || b == c) continue;
						if (a > b && a > c) ev[a] += pb[b] * pc[c];
						if (a < b || a < c) ev[a] -= pb[b] * pc[c];
					}
				}
			}

			return ev;
		}

  static double[] cumulative(double[] d){
    double[] r = d.clone();
    for(int i=1;i<r.length;i++)
      r[i]+=r[i-1];
    return r;
  }
  
  static double[] fast(double[] pb, double[] pc){
    int n = pb.length;
    double[] wmax = new double[n];
    double[] wmin = new double[n];
    
    double[] pcc = cumulative(pc);
    double[] pbc = cumulative(pb);
    
    for (int i = 1; i < n; i++){  
      wmax[i]+=pb[i]*pcc[i-1];
      wmax[i]+=pc[i]*pbc[i-1];
    }
    for (int i = 0; i < n-1; i++){  
      wmin[i]+=pb[i]*(pcc[n-1]-pcc[i]);
      wmin[i]+=pc[i]*(pbc[n-1]-pbc[i]);
    }
    
    double[] wmaxc = cumulative(wmax);
    double wtotal = wmaxc[n-1];
    
    double[] ev = new double[n];
    for(int a=0;a<n;a++){
      if(a>0)
        ev[a]+=(wmaxc[a-1]);
      if(a<n-1)
        ev[a]-=(wtotal-wmaxc[a])*(1.0-wmin[a]/(wtotal-wmaxc[a]));
    }
    
    return ev;
  }