#include <stdio.h>
#include <iostream>
#include <iomanip>
 
#define NLen 32
#define MLen 512
 
using namespace std;
 
int conf[NLen][NLen];
 
double E[MLen][MLen];
long double tm[MLen];
 
int main()
{
    freopen("minesweeper.in","r",stdin);
    freopen("minesweeper.out","w",stdout);
     
    int N,M;
     
    cin>>N>>M;
     
    N*=M;
    M=0;
     
    for(int i=0;i<=N;++i)
        for(int j=0;i+j<=N;++j)
            conf[i][j]=++M;
             
    ++M;
     
    for(int b=0;b<=N;++b)
        for(int c=0;b+c<=N;++c)
        {           
            int a=N-b-c;
             
            E[conf[b][c]][conf[b][c]]=1;
             
            if(c==N) continue;
             
            E[conf[b][c]][M]=1;
             
            if(a) E[conf[b][c]][conf[b+1][c]]-=(double)a/N;
             
            if(b) E[conf[b][c]][conf[b-1][c+1]]-=(double)b/N;
             
            if(c) E[conf[b][c]][conf[b][c-1]]-=(double)c/N;
        }
         
    for(int i=1;i<M;++i)
    {
        if(E[i][i]!=0)
        {
            for(int j=i+1;j<=M;++j) E[i][j]/=E[i][i];
            E[i][i]=1;
        }
         
        for(int k=i+1;k<M;++k)
            if(E[k][i]!=0)
            {
                for(int j=i+1;j<=M;++j) E[k][j]-=E[k][i]*E[i][j];
                E[k][i]=0;
            }
    }
 
    for(int i=M;i>0;--i)
    {
        tm[i]=E[i][M];
         
        for(int j=i+1;j<M;++j)
            tm[i]-=E[i][j]*tm[j];
    }
     
    cout<<fixed<<setprecision(6)<<tm[1]<<'\n';
         
    return 0;
}