题目传送门
这题一看就知道是期望的dp。
我们只关心目前还有多少牌没翻过来,其中多少已经知道了。
所以我们设
d
p
i
,
j
dp_{i,j}
dpi,j表示有i张牌未反过来,j张已经知道了。
然后就分情况讨论,就行了。
边界判断要小心!!
/*
{By GWj
*/
#pragma GCC optimize(2)
#include<bits/stdc++.h>
#define rb(a,b,c) for(int a=b;a<=c;++a)
#define rl(a,b,c) for(int a=b;a>=c;--a)
#define LL long long
#define IT iterator
#define PB push_back
#define II(a,b) make_pair(a,b)
#define FIR first
#define SEC second
#define FREO freopen("check.out","w",stdout)
#define rep(a,b) for(int a=0;a<b;++a)
#define SRAND mt19937 rng(chrono::steady_clock::now().time_since_epoch().count())
#define random(a) rng()%a
#define ALL(a) a.begin(),a.end()
#define POB pop_back
#define ff fflush(stdout)
#define fastio ios::sync_with_stdio(false)
#define R(a) cin>>a
#define R2(a,b) cin>>a>>b
#define check_min(a,b) a=min(a,b)
#define check_max(a,b) a=max(a,b)
using namespace std;
const int INF=0x3f3f3f3f;
typedef pair<int,int> mp;
/*}
*/
double dp[2600][2600];//剩余i张不知道的牌,有j张知道的牌的期望
class PerfectMemory{
public:
double getExpectation(int n, int m){
n*=m;
rb(i,1,n){
rl(j,i/2,0){
if((i)&1) continue;
double pp=0;
if(i>=2&&j*2<=i-2){
dp[i][j]+=(dp[i-2][j]+1.0)*double(i-j-j)/double(i-j)/double(i-j-1);
pp+=double(i-j-j)/double(i-j)/double(i-j-1);
// cout<<"1."<<pp<<endl;
}
if(i!=j&&j&&(i-2)>=(j-1)*2){
dp[i][j]+=(dp[i-2][j-1]+1.0)*double(j)/(double)(i-j);
pp+=double(j)/(double)(i-j);
// cout<<"2."<<pp<<endl;
}
if(i>=2&&j&&i-2>=j*2){
dp[i][j]+=(dp[i-2][j]+2.0)*(double)(i-j-j)/(double)(i-j)*double(j)/(double)(i-j-1);
pp+=(double)(i-j-j)/(double)(i-j)*double(j)/(double)(i-j-1);
// cout<<"3."<<pp<<endl;
}
if(i-j>=2&&i>=2*(j+2)){
dp[i][j]+=(dp[i][j+2]+1.0)*(double)(i-j-j)/double(i-j)*(double)(i-2-j-j)/(double)(i-j-1);
pp+=(double)(i-j-j)/double(i-j)*(double)(i-2-j-j)/(double)(i-j-1);
// cout<<"4."<<pp<<endl;
}
}
}
return dp[n][0];
}
};