#include<bits/stdc++.h>
#define int long long
#define double long double
#define endl '\n'
using namespace std;
const double pi=acos(-1);
void solve()
{
    int n,m;
    cin>>n>>m;
    vector<double> sum(n+1),tot(n+1);
    auto check=[&](int x,int now,int pos)->double
    {
        if(x>now) return 1e18;
        int tmp;
        if(pos>m) tmp=2*m-pos+1;
        else tmp=pos-1;
        double len=(now-x)*2+tmp*(pi*x)/m;
        return len;
    };
    for(int i=1;i<=n;i++)
    {
        vector<double> tmp;
        for(int j=2;j<=2*m;j++)
        {
            int l=0,r=i+1;
            while (l < r) {
                int mid=(l+r)/2;
                if (check(mid,i,j) <= check(mid + 1,i,j)) r = mid;
                else l = mid + 1;
            }
            tmp.push_back(check(l,i,j));
            tot[i]+=check(l,i,j);
            // cout<<i<<' '<<j<<endl;
            // cout<<check(l,i,j)<<endl;
        }
        
        // for(int j=0;j<tmp.size();j++) cout<<tmp[j]<<' ';
        // cout<<endl;
        for(int j=tmp.size()-2;j>=0;j--)
            tmp[j]+=tmp[j+1];
        // for(int j=0;j<tmp.size();j++) cout<<tmp[j]<<' ';
        // cout<<endl;
        for(int j=0;j<tmp.size();j++) sum[i]+=tmp[j];
    }
    // for(int i=1;i<=n;i++) cout<<sum[i]<<' ';
    // cout<<endl;
    // for(int i=1;i<=n;i++) cout<<tot[i]<<' ';
    // cout<<endl;
    double ans=0;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<i;j++)
        {
            // cout<<i<<' '<<j<<endl;
            ans+=((i-j)*2*m+tot[j])*2*m;
        }
        // cout<<ans<<endl;
        ans+=sum[i];
    }
    if(m!=1)
    {
        for(int i=1;i<=n;i++)
        ans+=i*2*m;
    }
    cout<<fixed<<setprecision(10)<<ans<<endl;
}
signed main()
{
    cin.tie(0)->sync_with_stdio(0);
    int t=1;
    while(t--) solve();
    return 0;
}