題意：比較裸的生成樹計數問題。

  如何處理生成樹計數問題？

基爾霍夫矩陣：

if i==j  Kir[i][j] = i的度數

if i!=j   Kir[i][j] = i到j的平行邊的個數的負數

即，基爾霍夫矩陣 = 度數矩陣 - 鄰接矩陣

將基爾霍夫矩陣刪去第i行和第i列，餘下i-1階的行列式的值即為生成樹個數。（證明略）

求行列式的值可以將行列式轉為上三角陣，求對角線上的積即為行列式的值。

可以用高斯消元的方法轉換上三角陣，過程中要mod p，除法改為逆元。

程式碼：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int N = 400+10;
typedef long long ll;
const ll mod = 10007;

struct point{
    int x,y;
    int id;
}p[N];
bool d[N][N],f[N][N];
int T,n,r;
int flag = 0;
ll dis(point a,point b){
    return (a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y);
}
int online(point a,point b,point c){
    return (ll)(a.x-c.x)*(b.y-c.y) == (ll)(b.x-c.x)*(a.y-c.y);
}
void build(){
    memset(d,0,sizeof(d));

    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++) if(i!=j){
        if(dis(p[i],p[j]) <= (ll)r*r){
            int ok = 1;
            for(int k=0;k<n;k++)if(k!=i&&k!=j){
                if( (p[i].x<=p[k].x && p[k].x<=p[j].x)||(p[j].x<=p[k].x && p[k].x<=p[i].x) ){
                    if(online(p[i],p[j],p[k])){
                        ok = 0;
                        break;
                    }
                }
            }
            if(ok) d[i][j]=1;
        }
    }

    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++) f[i][j]=d[i][j];

    for(int k=0;k<n;k++)
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
        f[i][j] = f[i][j] | (f[i][k]&&f[k][j]);

    flag = 1;
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    if(!f[i][j]) flag = 0;
}

ll exgcd(ll a,ll b,ll &x,ll &y)//乘法逆元返回的d是a,b的公約數，x是a mod b的逆元
{
    if(b==0)
    {
        x=1ll;y=0;
        return a;
    }
    ll d=exgcd(b,a%b,x,y);
    ll t=x;
    x=y;
    y=t-a/b*y;
    return d;
}
ll Gauss(int C[][N],int n)//計算n階行列式的絕對值 % mod
{
    ll ans=1ll;
    int flag=1;//行列交換的次數
    int i,j,k;
    for(i=0;i<n;i++)
    {
        if(C[i][i]==0)
        {
            for(j=i+1;j<n;j++)
                if(C[j][i])break;
            if(j==n)return 0;//某列的值全是0的ans=0；
            flag=!flag;
            for(int k=i;k<n;k++)
                swap(C[i][k],C[j][k]);//i和j行交換
        }
      ans=ans*C[i][i]%mod;//對角線相乘
      ll x,y;
      int tp=exgcd(C[i][i],mod,x,y);//x為逆元


      for(k=i+1;k<n;k++)
        C[i][k]=C[i][k]*x%mod;


      for(int j=i+1;j<n;j++)
      for(int k=i+1;k<n;k++)
      {
          C[j][k]=(C[j][k]-(ll)C[j][i]*C[i][k])%mod;
          if(j==k)
          C[j][k]=(C[j][k]+mod)%mod;
      }
      for(k=i+1;k<n;k++)
        C[i][k]=(ll)C[i][k]*C[i][i]%mod;

    }


    ans=(ans%mod+mod)%mod;
    if(flag) return ans;
    else  return mod-ans;
}

int Kir[N][N];
int solve(){
    memset(Kir,0,sizeof(Kir));
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
        Kir[i][i] += d[i][j];
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
        Kir[i][j] -= d[i][j];

    return Gauss(Kir,n-1);
}
int main(){
    cin >> T;
    while(T--){
//        cin >> n >> r;
        scanf("%d%d",&n,&r);
        for(int i=0;i<n;i++){
            scanf("%d%d",&p[i].x,&p[i].y);
            p[i].id = i;
        }
        build();
        if(flag==0){
            puts("-1");
            continue;
        }
        ll ans = solve();
//        cout<<ans<<endl;
        printf("%d\n",(int)ans);
    }
    return 0;
}