// fannhkn
#include<bits/stdc++.h>
#define forinc(i,a,b) for(int i=a;i<=b;i++)
#define int long long
#define fi first
#define se second
#define pii pair<int,int>
using namespace std;
const int N=1e5+10;
int n,k;
pii a[N];
double xuli(pii x,pii y)
{
    return x.fi * y.se - x.se * y.fi;
}
void sol()
{
    cin>>n>>k;
    forinc(i,0,n-1) cin>>a[i].fi>>a[i].se;
    double res=0;
    int b=0;
    int c=k;
    int x=-1;
    double tmp=0;
    for(int i=1;i<=c;i++) tmp+=xuli(a[i-1],a[i]);
    for(int i=n-1 ; i != c ; i--)
    {
        double S = tmp + xuli(a[c],a[i]) + xuli(a[i] , a[b]);
        if(S>res)
        {
            res=S;
            x=i;
        }
    }
    forinc(i,1,n-1)
    {
        c=(c+1)%n;
        tmp+=xuli(a[(c - 1 + n) % n] , a[c]) - xuli(a[b] , a[(b+1)%n]);
        b = (b + 1) % n;
        if(x==c) x=(c+1)%n;
        double pre = tmp + xuli(a[c],a[x]) + xuli(a[x],a[b]);
        while(x!=b)
        {
            x=(x+1)%n;
            if(x==b)
            {
                x=(b-1+n)%n;
                break;
            }
            if(pre<tmp + xuli(a[c],a[x]) + xuli(a[x],a[b])) pre=tmp + xuli(a[c],a[x]) + xuli(a[x],a[b]);
            else
            {
                x=(x-1+n)%n;
                break;
            }
        }
        res=max(res,pre);
    }
    cout << setprecision(12) << res/2 << '\n';
    //forinc(i,(l+1)%n , r)
}
int32_t main()
{
    #define task "task"
  if(fopen(task".inp","r"))
  {
    freopen(task".inp","r",stdin);
    freopen(task".out","w",stdout);
  }
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    int t;cin>>t;
    while(t--)
    {
        sol();
    }
}