HDU - 4305 - Lightning 生成树计数 + 叉积判断三点共线
HDU - 4305
题意:
比较裸的一道生成树计数问题,构造Krichhoof矩阵,求解行列式即可。但是这道题还有一个限制,就是给定的坐标中,两点连线中不能有其他的点,否则这两点就不能连接。枚举点,用叉积计算是否共线即可。
![](/assets/blank.gif)
![](/assets/blank.gif)
#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 生成树计数 + 叉积判断三点共线相关推荐
- 【CCCC】L3-018 森森美图 (30分),计算几何+判断三点共线+bfs最短路
problem L3-018 森森美图 (30分) 森森最近想让自己的朋友圈熠熠生辉,所以他决定自己写个美化照片的软件,并起名为森森美图.众所周知,在合照中美化自己的面部而不美化合照者的面部是让自己占 ...
- HDU 4305 Lightning (高斯消元解kirchhoff矩阵+逆元)
题意是:给一些坐标点,如果两点之间的距离小于R,并且两点之间没有其他点,则这两个点保持连通,这样构成了一个图.问这个图中生成树的个数. 因为数据量并不大,O(N^3)的建图没有问题. 建好图以后就可以 ...
- 笔试面试题目:三点共线的判断
看腻了认知鸡汤,你一无所获,味同嚼蜡?看累了万字长文,你收藏了事,心事重重? 那么,我们一起,上点算法小菜,斟点程序小酒,品味一下短小的趣味算法和程序,别有一番风味. 面试题目 在本文中,我们要讨论的 ...
- HDU 5304 Eastest Magical Day Seep Group's Summer(状压DP+环缩点+生成树计数)
题意:给出一个图m条边和n个顶点,要求删除m-n条边使得原图联通,问有多少种方法. 思路:显然删除后的图是由一棵树和额外的一条边组成的,如果只有树那么可以用MATRIX-TREE定理求出生成树的数量, ...
- Matrix-Tree (生成树计数)
生成数计数 对于一个无向简单图,nnn个点和n−1n-1n−1条边构成一个生成树,生成树计数就是求这个无向图中一共有几种不同的生成树(任意两边不同) 度矩阵 n×nn × nn×n 的矩阵,统计每个点 ...
- 牛客多校3 - Operation Love(几何+叉积确定三点顺逆)
题目链接:点击查看 题目大意:在平面直角坐标系中给出一只手的轮廓,需要判断是左手还是右手 题目分析:因为这些点在给出时不能确定时顺时针还是逆时针给出,所以我们需要自己来判断,这里存一下给出三点判断顺时 ...
- 最短路径生成树计数+最短路径生成树
最短路径生成树计数. 我们应该先明白什么是最短路径生成树,不会戳这里. 计数方法明显是要使用乘法原理计数,也就是说我们可以得出每一步的方案数再乘进答案中. 接下来考虑如何的出每一步的方案数,所谓方案数 ...
- HDU 4380 Farmer Greedy(叉积和三角形知识的综合应用)
题目链接:http://acm.hust.edu.cn/vjudge/contest/view.action?cid=115760#problem/A 题目意思大致为由n个点(n小于100)和m个金矿 ...
- 蓝桥杯java生成树计数_第七届蓝桥杯决赛杂记
反正最近不想干什么正经事,回忆一下这些东西也好. 总的来说,这是我第二次去北京参加蓝桥杯决赛.去年获得了优秀奖,换句话来说就是重在参与.不过,今年成绩比去年好一点. (订的宾馆一如既往的温馨就不说了) ...
最新文章
- python爬去学校_python爬取学校教务系统
- 《Effective STL》学习笔记(第四部分)
- go+cookie+angular踩过的坑
- 这或才是产品人「年终奖」的真实情况!
- ssm mysql增删改查_SSM配置并实现简单的数据库增删改查操作
- virtualbox 网络配置
- swoft php怎么样,[原创]Swoole和Swoft的那些事 (Http/Rpc服务篇)
- 计算机专业基础857考试大纲,2018年哈尔滨工业大学854计算机基础考研大纲
- Broadcast receiver之我见
- linux运维和python运维_Linux运维对比:Shell运维 VS Python运维
- 基于linux的触摸屏组态软件,基于LINUX的控制系统组态软件研究与开发
- 运营进阶:打造好文案的万能公式
- 开关型直流电源比线性直流电源效率高的原因?
- weui popup的使用
- 机器学习任务的一般流程及必要步骤
- Java狐仙ol,狐仙八大职业介绍 简要分为四类
- 英语中 distinct 与 distinctive 的区别
- 如何用大数据做行为预测的?
- javaweb商城实现在线支付
- InfiniBand RDMA
热门文章
- 使用代理进行界面之间的反向传值
- Elasticsearch template(待续...)
- 虚拟机安装 xp步骤(参照百度文库)
- 《Spring参考手册》中定义了以下几个AOP的重要概念
- IT从业者几大职位薪酬对比
- ​多大分辨率图像做分类更适合?浙大华为国科大等提出Dynamic Resolution Network,降低计算量还提性能!...
- CVPR2021 Oral《Seeing Out of the Box》北科大中山大学微软提出端到端视觉语言表征预训练方法...
- 天池OCR大赛前排方案总结!
- EfficientPS:目前排名第一的高效高精度全景分割算法
- 华为云发起美食图片分类大赛!奖品丰厚还可免费使用云资源