uva 1637 Double Patience

#include<cstdio>
using namespace std;
char s[10][5];
char ss[5];
int left[10];
int v[2000001];
double dp[2000001];
int t;
int state()
{
    int tmp=0;
    for(int i=1;i<=9;i++) tmp=tmp*5+left[i];
    return tmp;  
}
double solve()
{
    int x=state();
    if(!x) return 1.0;
    if(v[x]==t) return dp[x];
    double ans=0; int cnt=0;
    for(int i=1;i<9;i++)
     for(int j=i+1;j<=9;j++) 
      if(left[i]&&left[j]&&s[i][left[i]]==s[j][left[j]])
      {
           left[i]--; left[j]--;
          ans+=solve();
          cnt++;
          left[i]++; left[j]++;
      }
    if(cnt) dp[x]=1.0*ans/cnt;
    else dp[x]=0;
    v[x]=t;
    return dp[x];
}
int main()
{
    while(scanf("%s",ss)!=EOF)
    {
        s[1][1]=ss[0];
        for(int j=2;j<=4;j++) 
        {
            scanf("%s",ss);
            s[1][j]=ss[0];
        }
        for(int i=2;i<=9;i++)
         for(int j=1;j<=4;j++)
          {
              scanf("%s",ss);
              s[i][j]=ss[0];
          }
        t++;
        for(int i=1;i<=9;i++) left[i]=4;
        printf("%.6lf\n",solve());
    }
}