#include<bits/stdc++.h>
using namespace std;
#define int long long
const long double pai=acos(-1);
long double n,m,p,q,sum1;
int down,up;
void init(){
	down=2.0*m/pai,up=ceil(2.0*m/pai);
}
long double calc(int x){
	return (x+1)*x/2;
}
pair<long double,long double>sum(int l,int r){
	if(l>r){
		return {0,0};
	}
	long double A=0,B=0,limit=down;
	if(limit<l){
		B=2*(r-l+1);
	}
	if(limit>r){
		A=calc(r)-calc(l-1);
	}
	if(l<=limit&&limit<=r){
		A=calc(limit)-calc(l-1),B=(r-limit)*2;
	}
	return {A,B};
}
pair<long double,long double>sum2(int l,int r){
	if(l>r){
		return {0,0};
	}
	long double A=0,B=0,limit=up;
	if(limit<l){
		A=r-l+1;
	}
	if(limit>r){
		B=calc(r)-calc(l-1);
	}
	if(l<=limit&&limit<=r){
		A=r-limit+1,B=calc(limit-1)-calc(l-1);
	}
	return {A,B};
}
signed main(){
	cin>>n>>m,init();
	if(m==1){
		int ans=0;
		for(int i=1;i<=2*n;i++){
			ans+=i*(2*n-i+1);
			if(i<=n){
				ans-=2*i;
			}
		}
		cout<<ans;
		return 0;
	}
	for(int i=1;i<=n;i++){
		sum1=(sum1+i*(n-i)),q=(q+4*m*calc(n-i)*m);
	}
	for(int i=1;i<=2*m;i++){
		q=(q+2*(sum2(i,-(i+m+1-i-2*m)).first+sum2(1,min(2*m,i+m)-i).first)*calc(n)),p=(p+(sum2(i,-(i+m+1-i-2*m)).second+sum2(1,min(2*m,i+m)-i).second)*calc(n)),p=(p+(sum(2*m+1-i,m-1).first+2*sum(1,min(i+m,2*m)-i).first+sum(i,m-1).first)*sum1),q=(q+(sum(2*m+1-i,m-1).second+2*sum(1,min(i+m,2*m)-i).second+sum(i,m-1).second)*sum1);
	}
	for(int i=1;i<=n;i++){
		q=(q+2*m*i);
	}
	p/=m;
	cout<<setprecision(100)<<p*pai+q;
}