此题建图稍显麻烦,我先固定一个点,然后对斜率排序,斜率则写成了分数形式,然后嘛就是上模板
#include <utility>
#include <algorithm>
#include <string>
#include <cstring>
#include <cstdio>
#include <iostream>
#include <iomanip>
#include <set>
#include <vector>
#include <cmath>
#include <queue>
#include <bitset>
#include <map>
using namespace std;
#define clr(a,v) memset(a,v,sizeof(a))
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
const int INF = 0x7f7f7f7f;
const int maxn = 311;
const int mod = 10007;
const double pi = acos(-1.0);
const double eps = 1e-10;
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> VI;
typedef vector<VI> VVI;
typedef vector<VVI> VVVI;
typedef pair<int, int> pii;
struct node {
int x, y, dis, id;
node() {
}
node(int x, int y, int dis, int id) :
x(x), y(y), dis(dis), id(id) {
}
} p[maxn << 2];
inline bool same(long x, long y) {
return (x >> 31) == (y >> 31);
}
inline int cross(const node& a, const node& b) {
return a.x * b.y - a.y * b.x;
}
bool cmp(const node& a, const node& b) {
int t = cross(a, b);
if (!t)
return a.dis < b.dis;
return same(t, a.y * b.y);
}
int f[maxn], n, m;
int r[maxn], x[maxn], y[maxn];
int det[maxn][maxn];
int Det(int n) {
int ans = 1;
int i, j, k, t;
for (i = 0; i < n; ++i) {
for (j = i + 1; j < n; ++j) {
while (det[j][i]) {
t = det[i][i] / det[j][i];
for (k = i; k < n; ++k) {
det[i][k] -= det[j][k] * t;
det[i][k] %= mod;
swap(det[i][k], det[j][k]);
}
ans *= -1;
}
}
if (!det[i][i])
return 0;
}
for (i = 0; i < n; ++i)
ans = ans * det[i][i] % mod;
ans = (ans % mod + mod) % mod;
return ans;
}
int find(int n) {
return f[n] == n ? n : f[n] = find(f[n]);
}
void Union(int x, int y) {
int a = find(x);
int b = find(y);
if (a != b)
f[a] = b, r[b] += r[a];
}
int main() {
int R, i, j, cnt, T, dis, dx, dy;
scanf("%d", &T);
while (T--) {
clr(det,0);
scanf("%d%d", &n, &R);
R *= R;
for (i = 0; i < n; ++i)
f[i] = i, r[i] = 1;
cnt = 0;
for (i = 0; i < n; ++i)
scanf("%d%d", x + i, y + i);
for (i = 0; i < n; ++i) {
cnt = 0;
for (j = i + 1; j < n; ++j) {
dx = x[i] - x[j];
dy = y[i] - y[j];
dis = dx * dx + dy * dy;
if (dis <= R)
p[cnt++] = node(dx, dy, dis, j);
}
sort(p, p + cnt, cmp);
for (j = 0; j < cnt; ++j) {
if (!j || cross(p[j], p[j - 1])) {
det[i][i]++;
det[p[j].id][p[j].id]++;
det[i][p[j].id] = -1;
det[p[j].id][i] = -1;
Union(i, p[j].id);
}
}
}
int u = find(0);
if (r[u] < n)
puts("-1");
else
printf("%d\n", Det(n - 1));
}
return 0;
}