csp-s模擬測試101的T3代碼+註釋

由於題目過於大神因此單獨拿出來講。並且既然下發std了頹代碼貌似也不算可恥233

很難講啊，因此仍是寫在代碼註釋裏面吧

由於比較認真的寫了很多註釋，因此建議縮放到80%觀看，或者拿到gedit上

 1 //不要在乎我寫的是「氣」，輸入法找不到那個字。。。  2 #include<cstdio>  3 #include<algorithm>  4 long double C[205][205];int n,g,t,a[205],N;  5 struct Ex{//用於存儲指望值  6 long double cnt,tot;//cnt表示到達這個狀態的總方案數，tot表示這些方案所得到的神的總和  7 void add(long double e,long double c){tot+=e*c;cnt+=c;}//這裏是正推的，因此機率加和的結果不爲1,因此指望不能直接相加  8 //指望其實相似於加權平均數。如今你有cnt個方案，和爲tot。要加入c個方案，這些方案的平均數是e，因此這c個方案的和爲e*c  9 //那麼總的方案數就是cnt+c。總的和就是tot+e*c 10 long double E(){return cnt?tot/cnt:0;}//指望值就是總量/方案數啊，能夠重載成+=，可是由於參數有2個還得寫pair什麼的因此就沒有重載 11 }w[205][205],f[205][205];//全局變量初值爲0,因此剛開始的方案數和指望值都是0 12 long double cal(int l,int r,int ord){return ord>n?0:(l+r+1)/2.0;}//計算平均數，表示用了一份l的氣所獲得的神的指望。 13 //若是用的是原來不存在的氣，就是你加的那些編號爲[n+1,n+t]的氣，這些氣是不存在的，因此其實你放棄了這一份精，因此並不會獲得神，返回0 14 int main(){ 15 freopen("surmount.in","r",stdin);freopen("surmount.out","w",stdout); 16 scanf("%d%d%d",&n,&g,&t);//方便查閱：n是氣的份數，g是單次最大的精量，t是輪數 17 for(int i=1;i<=n;++i)scanf("%d",&a[i]); 18 for(int i=1;i<=n;++i)if(a[i]>g)a[i]=g;//大於g的氣是沒有用的，直接重置爲g 19 std::sort(a+1,a+1+n);N=n+t;//排序方便處理連續的區間 20 for(int i=n+1;i<=N;++i)a[i]=g;//添加一些足夠大的氣，這樣在下面處理的時候就能保證每一份精都能被壓制 21 for(int i=0;i<=t;++i)C[i][0]=1; 22 for(int i=1;i<=t;++i)for(int j=1;j<=i;++j)C[i][j]=C[i-1][j-1]+C[i-1][j];//處理組合數 23 //w[i][j]表示對於[i,j)這段區間的全部氣，若是使用了它們所得到的神的指望值 24 for(int i=0;i<=N;++i)w[i][i].add(0,1);//賦初值：空區間的方案數是1,對神的貢獻爲0 25 for(int i=N;~i;--i)for(int j=i+1;j<=N&&j-i<=t;++j)for(int k=i;k<j;++k) 26 w[i][j].add(w[i][k].E()+w[k+1][j].E()+cal(a[k+1],a[i],k+1),w[i][k].cnt*w[k+1][j].cnt*(a[k+1]-a[i])*C[j-i-1][k-i]); 27 //如今的操做是要把[i,k)與[k+1,j)合併成[i,j)。假設最後一個使用的精是k。注意區間開閉。 28 //add的兩個參數：指望值就是[i,k)和[k+1,j)這兩段的指望值的和，然而還有使用k所獲得的神，你使用它是用來壓制[a[k+1]，a[i]]的精，獲得對應的神，詳見cal函數 29 //總方案數就是兩邊的方案相乘（乘法計數原理）。還要乘上組合數，由於要考慮這些氣的順序能夠顛倒。由於左右兩邊內部已經有序了，因此只須要把兩個混合起來。 30 //就是相似於歸併排序兩個內部有序的數組，總方案數就是C[j-i-1][k-i]，表示在全部j-i-1個位置裏選出左邊所用的k-i個元素的位置，注意第k個被你強制定爲最後一個了因此是j-i還要-1 31 for(int i=0;i<=t;++i)f[i][i]=w[0][i]; 32 //初值，表示前i-1份氣沒有用。超過t沒有意義。 33 for(int i=0;i<=N;++i)for(int k=0;k<=t;++k)if(f[i][k].cnt)for(int j=i+1;j<=N&&k+j-i-1<=t;++j) 34 f[j][k+j-i-1].add(f[i][k].E()+w[i+1][j].E(),f[i][k].cnt*w[i+1][j].cnt*C[k+j-i-1][k]); 35 //i表示轉移以前已經用了幾份氣（含），j表示轉移以後用的氣，k表示轉移了幾輪，因此k+j-i-1就是轉移以後已經通過了的輪數 36 //f[a][b]表示用了前a份氣，已經通過了b輪以後獲得的神的指望，因此最後答案是f[N][t] 37 //轉移就是，先看指望，就是已經作完的f[i][k]再加上新的一輪你用了[i+1,j)這一段的氣，這樣所得到的神 38 //總方案數其實也同樣，和w的轉移比較相似，也是排序了兩個內部有序的過程 39 printf("%.10Lf\n",f[N][t].E()); 40 }//By DeepinC究級壓行+垃圾註釋（硬核25行）

累。。。