#include<bits/stdc++.h>
#define ll long long
#define ld double
using namespace std;
struct node{
	ld x,y;
}t[1000010];
ld cal(node a,node b,node c,node d){
	ld dx1=a.x-b.x,dy1=a.y-b.y;
	ld dx2=c.x-d.x,dy2=c.y-d.y;
	return dx1*dy2-dx2*dy1;
}
bool cmp1(node a,node b){
	node c;
	c.x=0;
	c.y=0;
	ld num=cal(a,c,b,c);
	if(num==0){
		if(a.x>0) return a.x<b.x;
		else return a.x>b.x; 
	}
	else if(num>0) return 1;
	else return 0;
}
vector<node> st1,st2;
vector<node> st;
ld ccal(node a,node b){
//		printf(" x1=%lld y1=%lld x2=%lld y2=%lld \n",a.x,a.y,b.x,b.y);
	node c;
	c.x=0;
	c.y=0;
	ld len=cal(a,c,b,c);
	if(len==0){
		if(a.x!=0){
			if(a.x*b.x<0) return 3.1415926536;
			else if(abs(a.x)<abs(b.x)){
				return 0;
			}
			else return 6.2831853072;
		}
		else{
			if(a.y*b.y<0) return 3.1415926536;
			else if(abs(a.y)<abs(b.y)){
				return 0;
			}
			else return 6.2831853072;
		}
	}
	ld t1,t2;
	t1=a.x*b.x+a.y*b.y;
	t2=sqrt(a.x*a.x+a.y*a.y)*sqrt(b.x*b.x+b.y*b.y);
//	printf("t1/t2=%Lf acos=%Lf\n",t1/t2,acos(t1/t2));
	if(len>0) return acos(t1/t2);
	else return 6.2831853072-acos(t1/t2);
}
void slove(){
	st1.clear();
	st2.clear();
	st.clear();
	ll n,k;
	scanf("%lld%lld",&n,&k);
	for(int i=1;i<=n;i++){
		scanf("%lf%lf",&t[i].x,&t[i].y);
		if(t[i].y==0){
			if(t[i].x>0) st1.push_back(t[i]);
			else st2.push_back(t[i]);
		}
		else if(t[i].y>0){
			st1.push_back(t[i]);
		}
		else{
			st2.push_back(t[i]);
		}
	}
	sort(st1.begin(),st1.end(),cmp1);
	sort(st2.begin(),st2.end(),cmp1);
	for(int i=0;i<st1.size();i++){
		st.push_back(st1[i]);
	}
	for(int i=0;i<st2.size();i++){
		st.push_back(st2[i]);
	}
	for(int i=0;i<k;i++){
		st.push_back(st[i]);
	}
//	printf("\n\n\n");
//	for(int i=0;i<st.size();i++){
//		printf("%d %lf %lf\n",i,st[i].x,st[i].y);
//	}
//	printf("\n\n\n");
	ld ans=0;
	for(int i=0;i<n;i++){
		ans=max(ans,ccal(st[i],st[i+k]));
	}
	printf("%.11lf\n",ans);
	return ;
}
int main(){
	int T;
	scanf("%d",&T);
	while(T--) slove();
	return 0;
}