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

struct Point {
    int x, y;
};

double pi() {
    return atan(1) * 4;
}

double cal_angle(Point p) {
    return atan2(p.y, p.x) + pi();
}

void solve() {
    int n, k;
    cin >> n >> k;

    vector<Point> islands(n);
    for (int i = 0; i < n; i++) {
        cin >> islands[i].x >> islands[i].y;
    }

    sort(islands.begin(), islands.end(), [](const auto& a, const auto &b) {
        return cal_angle(a) < cal_angle(b);
    });

    vector<pair<double, int>> angles;
    for (int i = 0; i < islands.size(); i++) {        
        double cur_angle = cal_angle(islands[i]);
        if (angles.size() == 0 || abs(cur_angle - angles[angles.size()-1].first) > 1e-7) {
            angles.push_back(make_pair(cur_angle, 1));
        }
        else {
            angles[angles.size()-1].second++;
        }
    }
    // add 360o
    int sz = angles.size();
    for (int i = 0; i < sz; i++) {
        angles.push_back(make_pair(angles[i].first + 2 * pi(), angles[i].second));
    }
    double ans = 0.0;
    int l = 0, r = 0, cur_cnt = angles[0].second;
    while (r < angles.size()-1) {
        
        while (r < angles.size() && (cur_cnt <= k || l == r)) {
            r++;
            cur_cnt += angles[r].second;
        }

        // cout << "test " << l << " " << r << " " << angles[r].first - angles[l].first << "\n";
        if (l >= 0 && l < angles.size() && r >= l && r < angles.size()) {
            ans = max(ans, angles[r].first - angles[l].first);
        }

        while (l+1 <= r && cur_cnt > k) {
            cur_cnt -= angles[l].second;
            l++;
        }
    }
    cout << fixed << setprecision(10) << ans << "\n";
}

int main() {
    ios::sync_with_stdio(false);

    int T;
    cin >> T;
    while (T--) {
        solve();
    }

    return 0;
}