#include <bits/stdc++.h>
using namespace std;

#define endl '\n'
using ll = long long;
const double PI = acos(-1);
const double eps = 1e-9;

struct point {
	double x, y;
};

// 计算叉积
double cross(const point &a, const point &b) {
	return a.x * b.y - a.y * b.x;
}

// 计算点积
double dot(const point &a, const point &b) {
	return a.x * b.x + a.y * b.y;
}

// 判断点的象限
int quadrant(const point &p) {
	if (p.x >= 0 && p.y >= 0) return 1; // 第一象限
	if (p.x < 0 && p.y >= 0) return 2;  // 第二象限
	if (p.x < 0 && p.y < 0) return 3;   // 第三象限
	return 4;                           // 第四象限
}

// 极角排序：先按象限排，再用叉积
// 时间快，精度低
bool cmp1(point a,point b)
{
	if(atan2(a.y,a.x)!=atan2(b.y,b.x))
		return atan2(a.y,a.x)<atan2(b.y,b.x);
	else return a.x<b.x;
}


// 计算向量长度
double len(const point &a) {
	return sqrt(a.x * a.x + a.y * a.y);
}

void solve() {
	int n, k;
	cin >> n >> k;
	vector<point> a(n);
	for (int i = 0; i < n; i++) {
		cin >> a[i].x >> a[i].y;
	}
	sort(a.begin(), a.end(), cmp1);
	double ans = 0;
	for (int i = 0, j = k; i < n; i++, j++) {
		double res = cross(a[i], a[j % n]);
		double res1 = dot(a[i], a[j % n]) / (len(a[i]) * len(a[j % n]));
		double angle = acos(min(1.0, max(-1.0, res1))); // 防止精度误差导致的问题
		if (res < 0 || (fabs(res) < eps && j % n <= i)) {
			angle = 2 * PI - angle;
		}
		ans = max(ans, angle);
	}
	cout << fixed << setprecision(10) << ans << endl;
}

int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr), cout.tie(nullptr);
	
	int T = 1;
	cin >> T;
	while (T--) {
		solve();
	}
	
	return 0;
}