#include<bits/stdc++.h>
using namespace std;
const int N=2e5;
struct ed{int x,y,d;}e[2000*2001/2+5];
int pr[N+5],p[N+5],cnt=0;
int x[2005],y[2005];
int d[2005];
int T,n,m,tot;
bitset<2005> gr[2005],le[2005];
int main(){
    // freopen("a.in","r",stdin);
    pr[1]=1;
    for (int i=2;i<=N;++i)
        if (!pr[i]){
            p[++cnt]=i;
            for (int j=1;j<=cnt && 1ll*i*p[j]<=N;++j){
                pr[i*p[j]]=1;
                if (i%p[j]==0) break;
            }
        }
    for (scanf("%d",&T);T;--T){
        scanf("%d%d",&n,&m);tot=0;
        for (int i=1;i<=n;++i) scanf("%d%d",&x[i],&y[i]),gr[i].reset(),le[i].reset();
        for (int i=1;i<=n;++i)
            for (int j=1;j<=n;++j)
                if (j!=i)
                    gr[i].set(j);
        for (int i=1;i<=n;++i)
            for (int j=i+1;j<=n;++j)
                e[++tot]=(ed){i,j,abs(x[i]-x[j])+abs(y[i]-y[j])};
        sort(e+1,e+1+tot,[](ed a,ed b){return a.d<b.d;});
        int ans=0;
        for (int i=1,it1=1,it2=1;i<=tot;++i){
            while (it1<=tot && e[it1].d<e[i].d){
                le[e[it1].x].set(e[it1].y);
                le[e[it1].y].set(e[it1].x);
                ++it1;
            }
            while (it2<=tot && e[it2].d<=e[i].d){
                gr[e[it2].x].flip(e[it2].y);
                gr[e[it2].y].flip(e[it2].x);
                ++it2;
            }
            if (!pr[e[i].d]){
                bitset<2005> t=(le[e[i].x]&gr[e[i].y])|(gr[e[i].x]&le[e[i].y]);
                ans+=t.count();
            }
        }
        for (int i=1,j=1;i<=tot;){
            int tmp=0;
            while (j<=tot && e[j].d==e[i].d){
                tmp+=d[e[j].x]+d[e[j].y];
                ++d[e[j].x],++d[e[j].y];
                ++j;
            }
            if (!pr[e[i].d]) ans+=tmp;
            while (i<=tot){
                --d[e[i].x],--d[e[i].y];
                ++i;
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}