// Slower C++ solution:
#include <bits/stdc++.h>

using namespace std;



using VS = vector<string>;    using LL = long long;

using VI = vector<int>;       using VVI = vector<VI>;

using PII = pair<int, int>;   using PLL = pair<LL, LL>;

using VL = vector<LL>;        using VVL = vector<VL>;



#define ALL(a)  begin((a)),end((a))

#define RALL(a) (a).rbegin(), (a).rend()

#define PB push_back

#define EB emplace_back

#define MP make_pair

#define SZ(a) int((a).size())

#define SORT(c) sort(ALL((c)))

#define RSORT(c) sort(RALL((c)))

#define UNIQ(c) (c).erase(unique(ALL((c))), end((c)))

#define FOR(i, s, e) for (int(i) = (s); (i) < (e); (i)++)

#define FORR(i, s, e) for (int(i) = (s); (i) > (e); (i)--)

#define debug(x) cerr << #x << ": " << x << endl

const int INF = 1e9;                          const LL LINF = 1e16;

const LL MOD = 1000000007;                    const double PI = acos(-1.0);

int DX[8] = { 0, 0, 1, -1, 1, 1, -1, -1 };    int DY[8] = { 1, -1, 0, 0, 1, -1, 1, -1 };



/* -----  2018/03/04  Problem: 064_arc_c / Link: https://arc064.contest.atcoder.jp/tasks/arc064_c?lang=en  ----- */

/* ------問題------



xy 平面があります。 すぬけ君は座標 (xs,ys) から座標 (xt,yt) まで移動しようとしています。 すぬけ君は好きな向きへ速さ 1 で動くことができます。 なお、すぬけ君は大きさのない点と見なすことにします。

平面上には N 個の円形のバリアが張ってあります。 i 番目のバリアは中心が (xi,yi) で半径が ri です。 バリアは互いに重なっていたり、互いを含んでいたりすることがあります。

平面上の各座標について、その座標がどのバリアの内部にも含まれない場合、その座標には宇宙線が降り注いでいます。

すぬけ君は移動中にできるだけ宇宙線を浴びたくないので、宇宙線を浴びる時間が最小になるように移動します。 すぬけ君が移動中に宇宙線を浴びる時間の最小値を求めてください。



-----問題ここまで----- */

/* -----解説等-----



移動する人の気持ちになると円から円の移動は中心から中心へ移動するのが最適なのでこれに工夫なく辺を張れば良い

なぜ600点?



----解説ここまで---- */







using Type = pair<double, int>;



vector<double> Dijkstra(vector<vector<pair<int,double>>>& G, int s) {

	vector<double> dist(SZ(G), LINF);

	dist[s] = 0;

	priority_queue<Type, vector<Type>, greater<Type>> que;

	que.push(Type(0LL, s));



	while (!que.empty()) {

		Type p = que.top(); que.pop();

		int v; /* 頂点*/ double d; /* 頂点vまでの距離 */

		tie(d, v) = p;



		if (d > dist[v]) continue;



		FOR(i, 0, (int)G[v].size()) {

			int nv = G[v][i].first;

			if (dist[nv] > dist[v] + G[v][i].second) {

				dist[nv] = dist[v] + G[v][i].second;

				que.push(Type(dist[nv], nv));

			}

		}

	}

	return dist;

}



LL N;



double ans = 0LL;

struct Circle {

	double x, y, r;

	Circle() { x = y = r = 0; }

};



double distance(Circle &a, Circle &b) {

	double d =hypot(a.x - b.x, a.y - b.y)- (a.r + b.r );

	return (d > 0 ? d : 0);

}

int main() {

	cin.tie(0);

	ios_base::sync_with_stdio(false);



	double sx, sy, tx, ty;

	cin >> sx >> sy >> tx >> ty >> N;

	vector < Circle > Cir(N + 2);

	FOR(i, 0, N) {

		cin >> Cir[i].x >> Cir[i].y >> Cir[i].r;

	}

	const int S = N;

	const int T = N + 1;

	Cir[S].x = sx; Cir[S].y = sy;

	Cir[T].x = tx; Cir[T].y = ty;

	N += 2;

	

	vector<vector<pair<int,double>>> G(N);

	FOR(i, 0, N) {

		FOR(j, i + 1, N) {

			double d = distance(Cir[i], Cir[j]);

			G[i].push_back(pair<int, double>(j,d));

			G[j].push_back(pair<int, double>(i,d));

		}

	}

	vector<double> dist = Dijkstra(G, S);

	ans = dist[T];

	cout << fixed<<setprecision(10)<<ans << "\n";



	return 0;

}


// Faster C++ solution for the same problem (measured speedup: 4.2132099742498905x):
#include<cmath>

#include<cstdio>

#include<algorithm>

using namespace std;

const int MAXN = 1000;

const double EPS = 1E-10;

struct Point{

	double x, y;

	Point(double _x=0, double _y=0):x(_x), y(_y){}

};

struct Circle{

	Point p; double r;

}c[MAXN + 5];

typedef Point Vector;

Vector operator + (Vector a, Vector b) {return Vector(a.x+b.x, a.y+b.y);}

Vector operator - (Vector a, Vector b) {return Vector(a.x-b.x, a.y-b.y);}

Vector operator * (Vector a, double k) {return Vector(a.x*k, a.y*k);}

Vector operator / (Vector a, double k) {return Vector(a.x/k, a.y/k);}

double cross(Vector a, Vector b) {return a.x*b.y - a.y*b.x;}

double dot(Vector a, Vector b) {return a.x*b.x + a.y*b.y;}

double length(Vector a) {return sqrt(dot(a, a));}

double dcmp(double x) {if( fabs(x) <= EPS ) return 0; else return x > 0 ? 1 : -1;}

double pow2(double x) {return x*x;}

Vector Unit(Vector a) {double l = length(a); return Vector(a.x/l, a.y/l);}

double dist(Point a, Point b) {return sqrt(pow2(a.x - b.x) + pow2(a.y - b.y));}

bool Check(Circle a, Circle b) {return dcmp(a.r + b.r - dist(a.p, b.p)) >= 0;}

bool Check(Circle a, Point b) {return dcmp(a.r - dist(a.p, b)) >= 0;}

double dist(Circle a, Circle b) {

	Vector vb = Unit(a.p - b.p), va = Unit(b.p - a.p);

	Point A = (a.p + va*a.r), B = (b.p + vb*b.r);

	return dist(A, B);

}

double dist(Circle a, Point b) {

	Vector v = Unit(b - a.p);

	return dist(b, a.p + v*a.r);

}

Point s, t;

double G[MAXN + 5][MAXN + 5], dis[MAXN + 5];

bool vis[MAXN + 5];

int N;

void dijkstra() {

	for(int i=1;i<=N+1;i++)

		dis[i] = G[0][i];

	vis[0] = true;

	for(int i=1;i<=N+1;i++) {

		double mn = 1E18;int p;

		for(int j=0;j<=N+1;j++) {

			if( !vis[j] && dis[j] < mn ) {

				p = j;

				mn = dis[j];

			}

		}

		vis[p] = true;

		for(int j=0;j<=N+1;j++) {

			if( !vis[j] ) dis[j] = min(dis[j], dis[p] + G[p][j]);

		}

	}

}

int main() {

	scanf("%lf%lf%lf%lf", &s.x, &s.y, &t.x, &t.y);

	scanf("%d", &N);

	for(int i=1;i<=N;i++)

		scanf("%lf%lf%lf", &c[i].p.x, &c[i].p.y, &c[i].r);

	for(int i=1;i<=N;i++) {

		if( Check(c[i], s) ) G[0][i] = G[i][0] = 0;

		else G[0][i] = G[i][0] = dist(c[i], s);

		if( Check(c[i], t) ) G[i][N+1] = G[N+1][i] = 0;

		else G[i][N+1] = G[N+1][i] = dist(c[i], t);

	}

	G[0][N+1] = G[N+1][0] = dist(s, t);

	for(int i=1;i<=N;i++)

		for(int j=1;j<=N;j++) {

			if( Check(c[i], c[j]) ) G[i][j] = 0;

			else G[i][j] = dist(c[i], c[j]);

		}

	dijkstra();

	printf("%.10lf\n", dis[N+1]);

}