【NOIP2016】换教室（DP，期望）

题意：

对于刚上大学的牛牛来说, 他面临的第一个问题是如何根据实际情况中情合适的课程。

在可以选择的课程中,有2n节课程安排在n个时间段上。在第 i ( 1≤ i≤n)个时同段上, 两节内容相同的课程同时在不同的地点进行, 其中, 牛牛预先被安排在教室 ci上课, 而另一节课程在教室 di进行。

在不提交任何申请的情况下,学生们需要按时间段的顺序依次完成所有的n节安排好的课程。如果学生想更换第i节课程的教室,则需要提出中情。若申请通过,学生就可以在第 i个时间段去教室 di上课, 否则仍然在教室 ci上课。

由于更换教室的需求太多, 申请不一定能获得通过。 通过计算, 牛牛发现申请更换第 i节课程的教室时, 中情被通过的概率是一个已知的实数 ki, 并且对于不同课程的申请, 被通过的概率是互相独立的。

学校规定, 所有的申请只能在学期开始前一次性提交, 并且每个人只能选择至多m节课程进行申请。 这意味着牛牛必须一次性决定是否申请更换每节课的教室, 而不能根据某些课程的申请结果来决定其他课程是否申请; 牛牛可以申请白己最希望更换教室的 m门课程,也可以不用完这m个中情的机会,甚至可以一门课程都不申请。

因为不同的课程可能会被安排在不同的教室进行, 所以牛牛需要利用课问时间从一间教室赶到另一间教室。

牛牛所在的大学有 v个教室,有 e条道路。每条道路连接两间教室, 并且是可以双向通行的。 由于道路的长度和拥;i者程度不同, 通过不同的道路耗费的体力可能会有所不同。当第i ( 1≤i≤n-1 )节课结束后,牛牛就会从这节课的教室出发,选择一条耗费体力最少的路径前往下一节课的教室。

现在牛牛想知道,申请哪几门课程可以使他因在教室问移动耗费的体力值的总和的期望值最小,请你帮他求出这个最小值

保证1≤n≤2000, 0≤m≤2000, 1≤v≤300, 0≤ e≤90000。

思路：萎靡的期望DP

        dp[i,j,0/1]表示前i门申请j门（不一定通过），最后一门是否申请的最小期望

        分四种情况讨论

        V打成N毁一生

var dp:array[0..2100,0..2100,0..1]of double;
    f:array[0..2100,0..2100]of double;
    c,d:array[0..2100]of longint;
    g:array[0..2100]of double;
    n,m,v1,e1,i,j,k,x,y,z:longint;
    ans:double;

function min(x,y:double):double;
begin
 if x<y then exit(x);
 exit(y);
end;

begin
 
 readln(n,m,v1,e1);
 for i:=1 to n do read(c[i]);
 for i:=1 to n do read(d[i]);
 for i:=1 to n do read(g[i]);
 for i:=0 to v1 do
  for j:=0 to v1 do f[i,j]:=1<<62;
 for i:=0 to v1 do f[i,i]:=0;
 for i:=1 to v1 do begin f[0,i]:=0; f[i,0]:=0; end;
  for i:=1 to e1 do
 begin
  read(x,y,z);
  f[x,y]:=min(f[x,y],z);
  f[y,x]:=min(f[y,x],z);
 end;
 for i:=1 to v1 do
  for j:=1 to v1 do
   for k:=1 to v1 do
    if(i<>j)and(i<>k)and(j<>k) then f[j,k]:=min(f[j,k],f[j,i]+f[i,k]);
 for i:=0 to n do
  for j:=0 to m do
  begin
   dp[i,j,0]:=1<<62; dp[i,j,1]:=1<<62;
  end;
 dp[0,0,0]:=0; //dp[0,0,1]:=0;
 g[0]:=0;
 for i:=0 to n-1 do
  for j:=0 to m do
  begin
   dp[i+1,j,0]:=min(dp[i+1,j,0],dp[i,j,0]+f[c[i],c[i+1]]);

   if j+1<=m then dp[i+1,j+1,1]:=min(dp[i+1,j+1,1],dp[i,j,0]+g[i+1]*f[c[i],d[i+1]]
                                   +(1-g[i+1])*f[c[i],c[i+1]]);

   dp[i+1,j,0]:=min(dp[i+1,j,0],dp[i,j,1]+g[i]*f[d[i],c[i+1]]
                               +(1-g[i])*f[c[i],c[i+1]]);

   if j+1<=m then dp[i+1,j+1,1]:=min(dp[i+1,j+1,1],dp[i,j,1]+
                      g[i]*g[i+1]*f[d[i],d[i+1]]+
                      g[i]*(1-g[i+1])*f[d[i],c[i+1]]+
                      (1-g[i])*g[i+1]*f[c[i],d[i+1]]+
                      (1-g[i])*(1-g[i+1])*f[c[i],c[i+1]]);


  end;
 ans:=1<<62;
 for i:=0 to m do ans:=min(ans,min(dp[n,i,0],dp[n,i,1]));
 writeln(ans:0:2);

end.