uva 1637 - Double Patience

uva 1637 Double Patience (全概率)

根据全概率公式： 每次取得的概率是一样的，所以进行操作的时候，可以将所有的概率加起来，然后除以他们发生的种类数，递归中所有的情况都是全概率的情况，所以每个概率都需要计算出来，然后慢慢往上推出来，最后获取一开始的情况，乘以条件发生的概率1/k，k为一开始可以取得方法数。

WRONG ANSWER:

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
int d[10][5];
double d_status[6][6][6][6][6][6][6][6][6];
int vis[6][6][6][6][6][6][6][6][6];

double solve(int d1,int d2,int d3,int d4,int d5,int d6,int d7,int d8,int d9){
	
	if(vis[d1][d2][d3][d4][d5][d6][d7][d8][d9])return d_status[d1][d2][d3][d4][d5][d6][d7][d8][d9];
	
	int d_fal[10]={
		0,d1,d2,d3,d4,d5,d6,d7,d8,d9
	};
	
	int sum = 0;
	
	for(int i = 1;i<=9;i++){
		sum+=d_fal[i];
	}
	
	if(sum >= 5*9)return 1;
	int ans = 0;
	double s = 0;
	for(int i = 1;i<=9;i++)if(d_fal[i]<5)
	    for(int j = i+1;j<=9;j++){
    		if(d_fal[j]<5&&d[i][d_fal[i]] == d[j][d_fal[j]]){
		    	ans++;
		    	d_fal[i]++;
		    	d_fal[j]++;
		    	s+=solve(d_fal[1],d_fal[2],d_fal[3],d_fal[4],d_fal[5],d_fal[6],d_fal[7],d_fal[8],d_fal[9]);
		    	d_fal[i]--;
		    	d_fal[j]--;
		    }
    	}
	vis[d1][d2][d3][d4][d5][d6][d7][d8][d9] = 1;
	return ans == 0?0:d_status[d1][d2][d3][d4][d5][d6][d7][d8][d9] = s/ans;

}

void fs(char c,int i,int j){
	
	if(c == 'A')d[i][j] = 1;
	else if(c == 'T')d[i][j] = 10;
	else if(c =='J')d[i][j] = 11;
	else if(c == 'Q')d[i][j] = 12;
	else if(c == 'K')d[i][j] = 13;
	else d[i][j] = c-'0';
	
}

int main(){
	char a[4];
	while(~scanf("%s",a)){
	    memset(vis,0,sizeof(vis));
	    memset(d,0,sizeof(d));
		fs(a[0],1,4);
		
		for(int i = 3;i>=1;i--){
			scanf("%s",a);
			fs(a[0],1,i);
		}
		for(int i = 2;i<=9;i++)
		    for(int j = 4;j>=1;j--){
    			scanf("%s",a);
    			fs(a[0],i,j);
    		}
        printf("%.6f\n",solve(1,1,1,1,1,1,1,1,1));
	}
	
} 
/*
1 9 6 13
11 12 1 13
7 12 11 13
12 10 11 9
6 10 1 8
12 10 13 8
8 9 10 7
9 7 11 7
8 6 1 6
*/
/*
13 6 9 1
13 1 12 11
13 11 12 7
9 11 10 12
8 1 10 6
8 13 10 12
7 10 9 8
7 11 7 9
6 1 6 8
*/

ACCEPTED:

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
int d[10][5];
double d_status[6][6][6][6][6][6][6][6][6];
int vis[6][6][6][6][6][6][6][6][6];

double solve(int d1,int d2,int d3,int d4,int d5,int d6,int d7,int d8,int d9){
	
	if(vis[d1][d2][d3][d4][d5][d6][d7][d8][d9])return d_status[d1][d2][d3][d4][d5][d6][d7][d8][d9];
	
	int d_fal[10]={
		0,d1,d2,d3,d4,d5,d6,d7,d8,d9
	};
	
	int sum = 0;
	
	for(int i = 1;i<=9;i++){
		sum+=d_fal[i];
	}
	
	if(sum >= 5*9)return 1;
	int ans = 0;
	double s = 0;
	for(int i = 1;i<=9;i++)if(d_fal[i]<5)
	    for(int j = i+1;j<=9;j++){
    		if(d_fal[j]<5&&d[i][d_fal[i]] == d[j][d_fal[j]]){
		    	ans++;
		    	d_fal[i]++;
		    	d_fal[j]++;
		    	d_status[d_fal[1]][d_fal[2]][d_fal[3]][d_fal[4]][d_fal[5]][d_fal[6]][d_fal[7]][d_fal[8]][d_fal[9]] = solve(d_fal[1],d_fal[2],d_fal[3],d_fal[4],d_fal[5],d_fal[6],d_fal[7],d_fal[8],d_fal[9]);
		    	s+=d_status[d_fal[1]][d_fal[2]][d_fal[3]][d_fal[4]][d_fal[5]][d_fal[6]][d_fal[7]][d_fal[8]][d_fal[9]];
		    	vis[d_fal[1]][d_fal[2]][d_fal[3]][d_fal[4]][d_fal[5]][d_fal[6]][d_fal[7]][d_fal[8]][d_fal[9]] = 1;
		    	d_fal[i]--;
		    	d_fal[j]--;
		    }
    	}
	vis[d1][d2][d3][d4][d5][d6][d7][d8][d9] = 1;
	return ans == 0?0:d_status[d1][d2][d3][d4][d5][d6][d7][d8][d9] = s/ans;

}

void fs(char c,int i,int j){
	
	if(c == 'A')d[i][j] = 1;
	else if(c == 'T')d[i][j] = 10;
	else if(c =='J')d[i][j] = 11;
	else if(c == 'Q')d[i][j] = 12;
	else if(c == 'K')d[i][j] = 13;
	else d[i][j] = c-'0';
	
}

int main(){
	char a[4];
	while(~scanf("%s",a)){
	    memset(vis,0,sizeof(vis));
	    memset(d,0,sizeof(d));
		fs(a[0],1,4);
		
		for(int i = 3;i>=1;i--){
			scanf("%s",a);
			fs(a[0],1,i);
		}
		for(int i = 2;i<=9;i++)
		    for(int j = 4;j>=1;j--){
    			scanf("%s",a);
    			fs(a[0],i,j);
    		}
        printf("%.6f\n",solve(1,1,1,1,1,1,1,1,1));
	}
	
} 
/*
1 9 6 13
11 12 1 13
7 12 11 13
12 10 11 9
6 10 1 8
12 10 13 8
8 9 10 7
9 7 11 7
8 6 1 6
*/
/*
13 6 9 1
13 1 12 11
13 11 12 7
9 11 10 12
8 1 10 6
8 13 10 12
7 10 9 8
7 11 7 9
6 1 6 8
*/