题意:
几何处理后+求生成树个数+高斯消元+求逆元
知识点分析:
1.这题最主要的就是矩阵树了吧。。。学习新的图论知识好开心啊~~~
Matrix-Tree定理(Kirchhoff矩阵-树定理)。
论文参见:周冬《生成树的计数及其应用》
Matrix-Tree定理是解决生成树计数问题最有力的武器之一。
1、G的度数矩阵D[G]是一个n*n的矩阵,并且满足:当i≠j时,dij=0;当i=j时,dij等于vi的度数。
2、G的邻接矩阵A[G]也是一个n*n的矩阵, 并且满足:如果vi、vj之间有边直接相连,则aij=1,否则为0。
我们定义G的Kirchhoff矩阵(也称为拉普拉斯算子)C[G]为C[G]=D[G]-A[G],则Matrix-Tree定理可以描述为:G的所有不同的生成树的个数等于其Kirchhoff矩阵C[G]任何一个n-1阶主子式的行列式的绝对值。所谓n-1阶主子式,就是对于r(1≤r≤n),将C[G]的第r行、第r列同时去掉后得到的新矩阵,用Cr[G]表示。
2.逆元:知识点参考
一、模运算
模运算的性质总结如下:
(1) (a + b) % n = (a % n + b % n) % n
(2) (a - b) % n = (a % n - b % n) % n
(3) (a * b) % n = (a % n * b % n) % n
(4) ab % n = ((a % n)b) % n
定义 如果ab≡1(mod m), 则称b是a的模m逆,记作a的模m逆是方程ax≡1(mod m)的解. 两个数互质一定有逆元。
求逆元可以使用辗转相除法,但是只有两个数都是质数的时候才有逆元,举例如下:
例:求5的模7逆
做辗转相除法, 求得整数b,k使得 5b+7k=1, 则b是5的模7逆.
计算如下:
7=5+2, 5=2×2+1.
回代 1=5-2×2=5-2×(7-5)= 3×5-2×7,
得 5 -1≡3(mod7).
三、此处使用的模版(大神模板)long long inv(long long a,long long m)
{
if(a == 1)return 1;
return inv(m%a,m)*(m-m/a)%m;
}
反思:(坑到自己的地方)
一、几何处理:
应该要判断的是是否在线段上,这里学到了个姿势。(叉积为0+不是端点+判断x,y范围+可以有等号)
求距离不要开方,精度保证
二、建图:
用O(n^3)处理,可以有常数优化。
但是我优化挫了。。。多建边了。。。
就是在判断是否有中间阻挡点的时候之需要删除边而不需要加边,这样会把已经删掉的边再加回去。
三、高斯消元:
这题还是在套模版。。。Orz。。。但要注意传进去的参数n应该是(n-1)
AC代码
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
#define maxn 305
#define pi acos(-1)
#define eps 1e-6
using namespace std;
int T,n,r;
struct Point
{
int x,y;
}P[maxn];
int map[maxn][maxn],ma[maxn][maxn];
int cross(Point a,Point b,Point s)
{
int x1=a.x-s.x,y1=a.y-s.y;
int x2=b.x-s.x,y2=b.y-s.y;
return x1*y2-x2*y1;
}
int dis(Point a,Point b)
{
//cout<<"a.x = "<<a.x<<" b.x = "<<b.x<<" a.y = "<<a.y<<" b.y = "<<b.y<<endl;
return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
bool judge(int a,int b,int i)
{
if(P[i].x<=max(P[a].x,P[b].x) && P[i].y<=max(P[a].y,P[b].y) && P[i].x>=min(P[a].x,P[b].x) && P[i].y>=min(P[a].y,P[b].y))
return true;
else return false;
}
void check(int a,int b)
{
int i;
for(i=0;i<n;i++)
{
if(cross(P[a],P[i],P[b])==0 && a!=i &&b!=i)
{
if(judge(a,b,i))
{
map[a][b]=map[b][a]=0;
//map[a][i]=map[i][a]=1;
}
}
}
}
const int MOD = 10007;
int INV[MOD];
//求ax = 1( mod m) 的x值,就是逆元(0<a<m)
long long inv(long long a,long long m)
{
if(a == 1)return 1;
return inv(m%a,m)*(m-m/a)%m;
}
void init()
{
int i,j;
memset(map,0,sizeof(map));
memset(ma,0,sizeof(ma));
for(i=0;i<n;i++)
for(j=i+1;j<n;j++)
{
if(j==i) continue;
else
{
if(dis(P[i],P[j])<=r*r )
{
map[i][j]=map[j][i]=1;
check(i,j);
}
}
}
for(i=0;i<n;i++)
for(j=0;j<n;j++)
{
if(i!=j && map[i][j])
ma[i][i]++,ma[i][j]=-1;
}
}
int guss(int n)
{
int i,j,k;
for(i=0;i<n;i++)
for(j=0;j<n;j++)
ma[i][j] = ( ma[i][j]%MOD+MOD ) % MOD;
int res=1;
for(i=0;i<n;i++)
{
for(j=i;j<n;j++)
if(ma[j][i]!=0)
{
for(k=i;k<n;k++)
swap(ma[i][k],ma[j][k]);
if(i!=j)
res=(-res+MOD)%MOD;
break;
}
if(ma[i][i]==0)
{
res=-1;
break;
}
for(j=i+1;j<n;j++)
{
int tmp= (ma[j][i]*inv(ma[i][i],MOD))%MOD;
for(k=i;k<n;k++)
ma[j][k] = (ma[j][k]-(ma[i][k]*tmp)%MOD+MOD)%MOD;
}
res = (res * ma[i][i]) %MOD;
}
return res;
}
int main()
{
freopen("in","r",stdin);
//freopen("out","w",stdout);
scanf("%d",&T);
while(T--)
{
scanf("%d %d",&n,&r);
int i,j;
for(i=0;i<n;i++)
{
scanf("%d %d",&P[i].x,&P[i].y);
}
init();
// for(i=0;i<n;i++)
// {
// for(j=0;j<n;j++)
// printf("%d ",map[i][j]);
// putchar('\n');
// }
int ans=guss(n-1);
printf("%d\n",ans);
}
return 0;
}