#include <bits/stdc++.h>

using namespace std;
const double PI=acos(-1);
typedef long long ll;
int n,m;

double  c[514];


double val[514];

vector<pair<int,double>> e[500010];

double d[500010];
bool tf[500010];


void dijkstra(int x){
    priority_queue<pair<double,int>,vector<pair<double,int>>,greater<pair<double,int>>> q;
    q.push({0,x});
    d[x]=0;
    while(q.size()){
        auto t=q.top();
        q.pop();
        if(tf[t.second]) continue;
        tf[t.second]=1;
        for(auto [j,w]:e[t.second]){
            if(d[j]>d[t.second]+w)
            {
                d[j]=d[t.second]+w;;
                q.push({d[j],j});
            }
        }
    }
}

signed main() {
    cin>>n>>m;

    double ans=0;
    m*=2;
    int j=1;
    // double ans=0;
    for(int i=1;i<=n;i++){
        
        j+=m;
        for(int k=0;k<=j;k++) d[k]=1e9,tf[k]=0;
        
        d[j-1]=0;
        int base=j-1;
        base--;
        int R=base/m+1,C=base%m;
        for(int k=j-2;k>0;k--){
            int r=(k-1)/m+1,c=(k-1)%m;
            double d1=R-r+PI*2*r/m*min(m-abs(C-c),abs(C-c));
            double d2=R+r;
            // printf("%.15lf %.15lf\n",d1,d2);
            d[k]=min(d1,d2);
        }

        for(int k=1;k<j-m;k++) ans+=m*d[k];
        if(m!=2)
        ans+=m*i;
        double res=0;
        // cout<<j-m+1<<endl;
        for(int k=j-m;k<j;k++) res+=d[k];
        // for(int k=0;k<j;k++) printf("%.15f ",d[k]);
        // cout<<endl;
        res/=2;
        ans+=m*res;
    }
    printf("%.15f",ans);
}