Double Patience UVA - 1637 dp+概率

问题

分析

带有记忆化搜索的状态转移，dp[n]是从现在state->空 累计成功的概率

#include <cstring>
#include <cstdio>
#include <vector>
#include <iostream>
#include <cmath>
#include <vector>
#include <utility>
using namespace std;
typedef long long ll;
const int maxn=2000000;   //状态总数有5^9=1953125种
double dp[maxn],tot;
int vis[maxn],state[15],kase;
char piles[9][5],temp[5];

inline int ID(){
    int t=state[0];
    for(int i=1;i<9;++i) t=5*t+state[i];
    return t;
}
//n是还在场的牌的数量，dp(n)是累积的概率，从n->0
double dfs(int n){
    if(n==0) return 1;
    int id=ID();
    double &ans=dp[id];
    if(vis[id]==kase) return ans;
    vis[id]=kase;
    ans=0;

    int s=0;
    for(int i=0;i<9;++i){
        if(state[i]==0) continue;
        for(int j=i+1;j<9;++j){
            if(state[j]==0) continue;
            if(piles[i][state[i]]==piles[j][state[j]]){
                state[i]--;
                state[j]--;
                ans+=dfs(n-2);
                ++s;
                state[i]++;
                state[j]++;
            }
        }
    }
    if(s>0) ans=ans/s;  //一定要s>0才能除
    return ans;
}


int main(void){
    while(true){
        ++kase;
        for(int i=0;i<9;++i){
            for(int j=1;j<5;++j){
                if(scanf("%s",temp)==EOF) return 0;
                piles[i][j]=temp[0];
            }
        }
        fill(state,state+9,4);
        printf("%.6lf\n",dfs(36));
    }
}