HDU --- 4305 Lighting 【生成树计数 + 向量 】

传送门

题意: 有一道闪电, 如果两个点之间的距离小于R ， 并且他们直接没有点, 则这两个点之间就有一条边, 最后问这个图中生成树有几颗.

思路 : 矩阵树 + 判断向量是否共线既可以解决!!!
注意 : 判断完了共线后, 还要判断那个点是不是在他们中间, 因为有可能三个点共线, 但是与之共线的那个点不一定在判断的那两个点之间. 这道题难在判断两个点之间是否还有其他点.

调BUG调到吐血

AC 代码 :

#include<bits/stdc++.h>
#define mod 10007
#define ll long long
using namespace std;
const int maxn=3e2+5;
int du[maxn];
int a[maxn][maxn];
int n,r;
struct node
{
    int x,y;
}s[maxn];

ll det(int n)  //矩阵树定理.
{
    ll res=1;
    for(int i=2; i<=n; i++)
    {
        for(int j=i+1; j<=n; j++){
            while(a[j][i])
            {
                int t=a[i][i] / a[j][i] % mod;
                for(int k=i; k<=n; k++)
                    a[i][k]=(a[i][k]-a[j][k]*t) % mod;
                for(int k=i; k<=n; k++)
                    swap(a[i][k],a[j][k]);
                res=-res;
            }
        }
        if(a[i][i] == 0)  //图无法联通.
            return 0;
        res *= a[i][i];
        res %= mod;
    }
    res = ((res % mod ) + mod ) % mod;
    return res;
}

bool check(int i,int j,int k)
{   //向量共线并且坐标的对应关系.
    if(((s[i].x-s[k].x)*(s[j].y-s[k].y) == (s[i].y-s[k].y)*(s[j].x-s[k].x) && 
        (s[k].x <= max(s[i].x,s[j].x) && s[k].x >= min(s[i].x,s[j].x) 
        && s[k].y <= max(s[i].y,s[j].y) && s[k].y >= min(s[i].y,s[j].y) ) )  || 
       ((double)sqrt(((s[i].x-s[j].x) * (s[i].x-s[j].x) + 
                    (s[i].y-s[j].y) * (s[i].y-s[j].y))*1.0)) > r)
        return true;
    return false;
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--){
        memset(a,0,sizeof(a));
        memset(du,0,sizeof(du));
        scanf("%d%d",&n,&r);
        for(int i=1;i<=n;i++){
            scanf("%d%d",&s[i].x,&s[i].y);
        }
        int flag = 0;
        for(int i=1;i<=n;i++){
            for(int j=i+1;j<=n;j++){
                flag = 1;
                for(int k=1;k<=n;k++){
                    if(k == i || k == j) continue;
                    if(check(i,j,k))
                        flag = 0;
                }
                if(flag){  //这两个点符合要求.
                    a[i][j] = a[j][i] = -1;
                    du[i] ++;
                    du[j] ++;
                }
            }
        }
        for(int i=1;i<=n;i++)
            a[i][i] = du[i];
        ll res = det(n);
        if(res)
            printf("%lld\n",res);
        else
            printf("-1\n");
    }
}