HDU - 4305

题意:

    比较裸的一道生成树计数问题,构造Krichhoof矩阵,求解行列式即可。但是这道题还有一个限制,就是给定的坐标中,两点连线中不能有其他的点,否则这两点就不能连接。枚举点,用叉积计算是否共线即可。

#include <algorithm>
#include  <iterator>
#include  <iostream>
#include   <cstring>
#include   <cstdlib>
#include   <iomanip>
#include    <bitset>
#include    <cctype>
#include    <cstdio>
#include    <string>
#include    <vector>
#include     <stack>
#include     <cmath>
#include     <queue>
#include      <list>
#include       <map>
#include       <set>
#include   <cassert>using namespace std;
//#pragma GCC optimize(3)
//#pragma comment(linker, "/STACK:102400000,102400000")  //c++
// #pragma GCC diagnostic error "-std=c++11"
// #pragma comment(linker, "/stack:200000000")
// #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
// #pragma GCC optimize("-fdelete-null-pointer-checks,inline-functions-called-once,-funsafe-loop-optimizations,-fexpensive-optimizations,-foptimize-sibling-calls,-ftree-switch-conversion,-finline-small-functions,inline-small-functions,-frerun-cse-after-loop,-fhoist-adjacent-loads,-findirect-inlining,-freorder-functions,no-stack-protector,-fpartial-inlining,-fsched-interblock,-fcse-follow-jumps,-fcse-skip-blocks,-falign-functions,-fstrict-overflow,-fstrict-aliasing,-fschedule-insns2,-ftree-tail-merge,inline-functions,-fschedule-insns,-freorder-blocks,-fwhole-program,-funroll-loops,-fthread-jumps,-fcrossjumping,-fcaller-saves,-fdevirtualize,-falign-labels,-falign-loops,-falign-jumps,unroll-loops,-fsched-spec,-ffast-math,Ofast,inline,-fgcse,-fgcse-lm,-fipa-sra,-ftree-pre,-ftree-vrp,-fpeephole2",3)#define lson (l , mid , rt << 1)
#define rson (mid + 1 , r , rt << 1 | 1)
#define debug(x) cerr << #x << " = " << x << "\n";
#define pb push_back
#define pq priority_queuetypedef long long ll;
typedef unsigned long long ull;typedef pair<ll ,ll > pll;
typedef pair<int ,int > pii;
typedef pair<int,pii> p3;//priority_queue<int> q;//这是一个大根堆q
//priority_queue<int,vector<int>,greater<int> >q;//这是一个小根堆q
#define fi first
#define se second
//#define endl '\n'#define OKC ios::sync_with_stdio(false);cin.tie(0)
#define FT(A,B,C) for(int A=B;A <= C;++A)  //用来压行
#define REP(i , j , k)  for(int i = j ; i <  k ; ++i)
#define max3(a,b,c) max(max(a,b), c);
//priority_queue<int ,vector<int>, greater<int> >que;const ll mos = 0x7FFFFFFF;  //2147483647
const ll nmos = 0x80000000;  //-2147483648
const int inf = 0x3f3f3f3f;
const ll inff = 0x3f3f3f3f3f3f3f3f; //18
const int mod = 10007;
const double esp = 1e-8;
const double PI=acos(-1.0);
const double PHI=0.61803399;    //黄金分割点
const double tPHI=0.38196601;template<typename T>
inline T read(T&x){x=0;int f=0;char ch=getchar();while (ch<'0'||ch>'9') f|=(ch=='-'),ch=getchar();while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();return x=f?-x:x;
}/*-----------------------showtime----------------------*/const int maxn = 500;int n,r;int a[maxn][maxn],fa[maxn];struct node{int x,y;                }p[maxn];int find(int x){if(fa[x] == x)return x;return fa[x] = find(fa[x]);}void uni(int x,int y){int px = find(x);int py = find(y);fa[px] = py;}void cal(){ll ans = 1;int sign = 0;for(int i=1; i<=n; i++){        //当前行for(int j=i+1; j<=n; j++){int x = i, y = j;while(a[y][i]){ //利用gcd的方法,不停地进行辗转相除,达到消去其他行对应列元素的目的ll t = a[x][i] / a[y][i];for(int k=i; k<=n; k++)a[x][k] = (a[x][k] - a[y][k]*t)%mod;swap(x,y);}if(x != i){     //奇数次交换,则D=-D'整行交换for(int k = 1; k<=n; k++){swap(a[i][k], a[x][k]);}sign ^= 1;}}if(a[i][i] == 0){   //斜对角中有一个0,则结果为0puts("0");return;}else ans = ans * a[i][i] %mod;}if(sign) ans *= -1;if(ans < 0) ans += mod;printf("%lld\n", ans);}   double dis(int i,int j){return sqrt(1.0*(p[i].x - p[j].x)*(p[i].x - p[j].x) + 1.0*(p[i].y - p[j].y)*(p[i].y - p[j].y));}bool check(int i,int k,int j){return ((p[j].x - p[k].x)*(p[j].y - p[i].y) == (p[j].x - p[i].x)*(p[j].y - p[k].y)) \&&(max(p[i].x,p[j].x) >= p[k].x) &&(min(p[i].x,p[j].x) <= p[k].x)\&&(max(p[i].y,p[j].y) >= p[k].y) &&(min(p[i].y,p[j].y) <= p[k].y);    }
int main(){int t;  scanf("%d", &t);while(t--){memset(a,0,sizeof(a));scanf("%d%d", &n, &r);for(int i=1; i<=n; i++)fa[i] = i;for(int i=1; i<=n; i++){scanf("%d%d", &p[i].x, &p[i].y);}for(int i=1; i<=n; i++){for(int j=1; j<i; j++){int ok = 1;for(int k=1; k<=n; k++){if(k==i||k==j)continue;if(check(i,k,j)){ok=0;break;}}if(ok && dis(i,j) <= r){a[i][j] = a[j][i] = -1;a[i][i]++,a[j][j]++;uni(i,j);}}}int c = 0;for(int i=1; i<=n; i++){if(fa[i] == i)c++;}if(c!=1)puts("-1");else {n--;cal();}}return 0;
}

HDU - 4305

转载于:https://www.cnblogs.com/ckxkexing/p/9704052.html

HDU - 4305 - Lightning 生成树计数 + 叉积判断三点共线相关推荐

  1. 【CCCC】L3-018 森森美图 (30分),计算几何+判断三点共线+bfs最短路

    problem L3-018 森森美图 (30分) 森森最近想让自己的朋友圈熠熠生辉,所以他决定自己写个美化照片的软件,并起名为森森美图.众所周知,在合照中美化自己的面部而不美化合照者的面部是让自己占 ...

  2. HDU 4305 Lightning (高斯消元解kirchhoff矩阵+逆元)

    题意是:给一些坐标点,如果两点之间的距离小于R,并且两点之间没有其他点,则这两个点保持连通,这样构成了一个图.问这个图中生成树的个数. 因为数据量并不大,O(N^3)的建图没有问题. 建好图以后就可以 ...

  3. 笔试面试题目:三点共线的判断

    看腻了认知鸡汤,你一无所获,味同嚼蜡?看累了万字长文,你收藏了事,心事重重? 那么,我们一起,上点算法小菜,斟点程序小酒,品味一下短小的趣味算法和程序,别有一番风味. 面试题目 在本文中,我们要讨论的 ...

  4. HDU 5304 Eastest Magical Day Seep Group's Summer(状压DP+环缩点+生成树计数)

    题意:给出一个图m条边和n个顶点,要求删除m-n条边使得原图联通,问有多少种方法. 思路:显然删除后的图是由一棵树和额外的一条边组成的,如果只有树那么可以用MATRIX-TREE定理求出生成树的数量, ...

  5. Matrix-Tree (生成树计数)

    生成数计数 对于一个无向简单图,nnn个点和n−1n-1n−1条边构成一个生成树,生成树计数就是求这个无向图中一共有几种不同的生成树(任意两边不同) 度矩阵 n×nn × nn×n 的矩阵,统计每个点 ...

  6. 牛客多校3 - Operation Love(几何+叉积确定三点顺逆)

    题目链接:点击查看 题目大意:在平面直角坐标系中给出一只手的轮廓,需要判断是左手还是右手 题目分析:因为这些点在给出时不能确定时顺时针还是逆时针给出,所以我们需要自己来判断,这里存一下给出三点判断顺时 ...

  7. 最短路径生成树计数+最短路径生成树

    最短路径生成树计数. 我们应该先明白什么是最短路径生成树,不会戳这里. 计数方法明显是要使用乘法原理计数,也就是说我们可以得出每一步的方案数再乘进答案中. 接下来考虑如何的出每一步的方案数,所谓方案数 ...

  8. HDU 4380 Farmer Greedy(叉积和三角形知识的综合应用)

    题目链接:http://acm.hust.edu.cn/vjudge/contest/view.action?cid=115760#problem/A 题目意思大致为由n个点(n小于100)和m个金矿 ...

  9. 蓝桥杯java生成树计数_第七届蓝桥杯决赛杂记

    反正最近不想干什么正经事,回忆一下这些东西也好. 总的来说,这是我第二次去北京参加蓝桥杯决赛.去年获得了优秀奖,换句话来说就是重在参与.不过,今年成绩比去年好一点. (订的宾馆一如既往的温馨就不说了) ...

最新文章

  1. python爬去学校_python爬取学校教务系统
  2. 《Effective STL》学习笔记(第四部分)
  3. go+cookie+angular踩过的坑
  4. 这或才是产品人「年终奖」的真实情况!
  5. ssm mysql增删改查_SSM配置并实现简单的数据库增删改查操作
  6. virtualbox 网络配置
  7. swoft php怎么样,[原创]Swoole和Swoft的那些事 (Http/Rpc服务篇)
  8. 计算机专业基础857考试大纲,2018年哈尔滨工业大学854计算机基础考研大纲
  9. Broadcast receiver之我见
  10. linux运维和python运维_Linux运维对比:Shell运维 VS Python运维
  11. 基于linux的触摸屏组态软件,基于LINUX的控制系统组态软件研究与开发
  12. 运营进阶:打造好文案的万能公式
  13. 开关型直流电源比线性直流电源效率高的原因?
  14. weui popup的使用
  15. 机器学习任务的一般流程及必要步骤
  16. Java狐仙ol,狐仙八大职业介绍 简要分为四类
  17. 英语中 distinct 与 distinctive 的区别
  18. 如何用大数据做行为预测的?
  19. javaweb商城实现在线支付
  20. InfiniBand RDMA

热门文章

  1. 使用代理进行界面之间的反向传值
  2. Elasticsearch template(待续...)
  3. 虚拟机安装 xp步骤(参照百度文库)
  4. 《Spring参考手册》中定义了以下几个AOP的重要概念
  5. IT从业者几大职位薪酬对比
  6. ​多大分辨率图像做分类更适合?浙大华为国科大等提出Dynamic Resolution Network,降低计算量还提性能!...
  7. CVPR2021 Oral《Seeing Out of the Box》北科大中山大学微软提出端到端视觉语言表征预训练方法...
  8. 天池OCR大赛前排方案总结!
  9. EfficientPS:目前排名第一的高效高精度全景分割算法
  10. 华为云发起美食图片分类大赛!奖品丰厚还可免费使用云资源