【BZOJ1016】【Luogu P4208】 [JSOI2008]最小生成树计数 最小生成树,矩阵树定理
蛮不错的一道题,遗憾就遗憾在数据范围会导致暴力轻松跑过。
最小生成树的两个性质:
不同的最小生成树,相同权值使用的边数一定相同。
不同的最小生成树,将其都去掉同一个权值的所有边,其连通性一致。
这样我们随便跑一个\(MST\),就可以知道所有\(MST\)边的构造情况。由于性质二,我们可以考虑枚举每一种权值的所有边,保留所有非此权值的树边,看可以连出来多少种不同的最小生成树。也就是按照权值构造最小生成树,这个过程满足乘法原理。
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int N = 100 + 5;
const int M = 1000 + 5;
const int Mod = 31011;
struct Len {
int u, v, w;
bool operator < (Len rhs) const {
return w < rhs.w;
}
}l[M];
vector <Len> v;
int n, m, S[N];
int find (int x) {
return S[x] == x ? x : S[x] = find (S[x]);
}
vector <int> val, use, tot;
vector <int> :: iterator it;
void kruskal () {
sort (l + 1, l + 1 + m);
for (int i = 0; i <= n; ++i) S[i] = i;
for (int i = 1; i <= m; ++i) {
int fu = find (l[i].u);
int fv = find (l[i].v);
it = lower_bound (val.begin (), val.end (), l[i].w);
if (it == val.end ()) {
val.push_back (l[i].w);
use.push_back (0);
tot.push_back (1);
} else {
tot[it - val.begin ()]++;
}
if (fu != fv) {
S[fu] = fv;
it = lower_bound (val.begin (), val.end (), l[i].w);
use[it - val.begin ()]++;
v.push_back (l[i]);
}
}
}
int mat[N][N];
int gauss (int n) {
int ret = 1;
for (int i = 1; i <= n; ++i) {
for (int k = i + 1; k <= n; ++k) {
while (mat[k][i]) {
int d = mat[i][i] / mat[k][i];
for (int j = i; j <= n; ++j) {
(((mat[i][j] -= d * mat[k][j]) %= Mod) += Mod) %= Mod;
}
swap (mat[i], mat[k]); ret = -ret;
}
}
(((ret *= mat[i][i]) %= Mod) += Mod) %= Mod;
}
return abs (ret);
}
void add_edge (int u, int v) {
mat[u][u]++;
mat[v][v]++;
mat[u][v]--;
mat[v][u]--;
}
int sep[N];
int solve () {
kruskal ();
if (v.size () < n - 1) return 0;
int ans = 1;
for (int i = 0; i < val.size (); ++i) {
memset (mat, 0, sizeof (mat));
if (use[i] == 0 || tot[i] == use[i]) continue;
for (int j = 0; j <= n; ++j) S[j] = j;
for (int j = 0; j < v.size (); ++j) {
if (v[j].w != val[i]) {
S[find (v[j].u)] = find (v[j].v);
}
}
int cnt = 0;
for (int i = 1; i <= n; ++i) {
sep[++cnt] = find (i);
}
sort (sep + 1, sep + 1 + cnt);
cnt = unique (sep + 1, sep + 1 + cnt) - sep - 1;
for (int j = 1; j <= m; ++j) {
if (l[j].w == val[i]) {
int fu = find (l[j].u);
int fv = find (l[j].v);
fu = lower_bound (sep + 1, sep + 1 + cnt, fu) - sep;
fv = lower_bound (sep + 1, sep + 1 + cnt, fv) - sep;
add_edge (fu, fv);
}
}
(ans *= gauss (use[i])) %= Mod;
}
return ans;
}
signed main () {
// freopen ("data.in", "r", stdin);
cin >> n >> m;
for (int i = 1; i <= m; ++i) {
static int u, v, w;
cin >> u >> v >> w;
l[i] = (Len) {u, v, w};
}
cout << solve () << endl;
}
【BZOJ1016】【Luogu P4208】 [JSOI2008]最小生成树计数 最小生成树,矩阵树定理的更多相关文章
- bzoj1016 [JSOI2008]最小生成树计数——Kruskal+矩阵树定理
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1016 从 Kruskal 算法的过程来考虑产生多种方案的原因,就是边权相同的边有一样的功能, ...
- BZOJ 1016 最小生成树计数(矩阵树定理)
我们把边从小到大排序,然后依次插入一种权值的边,然后把每一个联通块合并. 然后当一次插入的边不止一条时做矩阵树定理就行了.算出有多少种生成树就行了. 剩下的交给乘法原理. 实现一不小心就会让程序变得很 ...
- 洛谷4208 JSOI2008最小生成树计数(矩阵树定理+高斯消元)
qwq 这个题目真的是很好的一个题啊 qwq 其实一开始想这个题,肯定是无从下手. 首先,我们会发现,对于无向图的一个最小生成树来说,只有当存在一些边与内部的某些边权值相同的时候且能等效替代的时候,才 ...
- 【BZOJ 1016】 1016: [JSOI2008]最小生成树计数 (DFS|矩阵树定理)
1016: [JSOI2008]最小生成树计数 Description 现在给出了一个简单无向加权图.你不满足于求出这个图的最小生成树,而希望知道这个图中有多少个不同的最小生成树.(如果两颗最小生成树 ...
- [spoj104][Highways] (生成树计数+矩阵树定理+高斯消元)
In some countries building highways takes a lot of time... Maybe that's because there are many possi ...
- spoj104 highways 生成树计数(矩阵树定理)
https://blog.csdn.net/zhaoruixiang1111/article/details/79185927 为了学一个矩阵树定理 从行列式开始学(就当提前学线代了.. 论文生成树的 ...
- @总结 - 7@ 生成树计数 —— matrix - tree 定理(矩阵树定理)与 prüfer 序列
目录 @0 - 参考资料@ @0.5 - 你所需要了解的线性代数知识@ @1 - 矩阵树定理主体@ @证明 part - 1@ @证明 part - 2@ @证明 part - 3@ @证明 part ...
- 【Luogu】P3317重建(高斯消元+矩阵树定理)
题目链接 因为这个专门跑去学了矩阵树定理和高斯消元qwq 不过不是很懂.所以这里只放题解 玫葵之蝶的题解 某未知dalao的矩阵树定理 代码 #include<cstdio> #inclu ...
- Luogu P4336 [SHOI2016]黑暗前的幻想乡 矩阵树定理+容斥原理
真是菜到爆炸....容斥写反(反正第一次写qwq) 题意:$n-1$个公司,每个公司可以连一些边,求每个边让不同公司连的生成树方案数. 矩阵树定理+容斥原理(注意到$n$不是很大) 枚举公司参与与否的 ...
- 【BZOJ4596】【Luogu P4336】 [SHOI2016]黑暗前的幻想乡 矩阵树定理,容斥
同样是矩阵树定理的裸题.但是要解决它需要能够想到容斥才可以. \(20\)以内的数据范围一定要试试容斥的想法. #include <bits/stdc++.h> using namespa ...
随机推荐
- MSSQL字符串分割
CREATE FUNCTION dbo.f_splitstr( @str varchar(8000) )RETURNS @r TABLE(id int IDENTITY(1, 1), value va ...
- 【POJ - 3087】Shuffle'm Up(模拟)
Shuffle'm Up 直接写中文了 Descriptions: 给定两个长度为len的字符串s1和s2, 接着给出一个长度为len*2的字符串s12. 将字符串s1和s2通过一定的变换变成s12, ...
- mariadb(第二章)增删改 MariaDB 数据类型
MariaDB 数据类型 MariaDB数据类型可以分为数字,日期和时间以及字符串值. 使用数据类型的原则:够用就行, 尽量使用范围小的,而不用大的 常用的数据类型 整数:int, bit 小数:de ...
- 【DSP开发】TI第二代KeyStone SoC诠释德仪的“云”态度
11月14日,期盼已久的德州仪器基于ARM Cortex-A15的产品终于新鲜出炉.伴随着TIKeyStone II多核 SoC系列产品的发布,结合了ARM Cortex-A15 处理器.C66x D ...
- MySQL Explain命令详解--表的读取顺序,数据读取操作的类型等
表示索引中使用的字节数,可通过该列计算查询中使用的索引的长度(key_len显示的值为索引字段的最大可能长度,并非实际使用长度,即key_len是根据表定义计算而得,不是通过表内检索出的) 不损失精确 ...
- Redis缓存击穿
缓存击穿 缓存击穿,是指一个key非常热点,在不停的扛着大并发,大并发集中对这一个点进行访问,当这个key在失效的瞬间,持续的大并发就穿破缓存,直接请求数据库,就像在一个屏障上凿开了一个洞. 比如在做 ...
- 掌握MyBatis的核心对象
一.获取SqlSessionFactoryBuilder对象 1.SqlSessionFactoryBuilder的作用 所有的MyBatis应用都是以SqlSessionFactory实例为中心.S ...
- C++ 的关键字(保留字)完整介绍
转载至:https://www.runoob.com/w3cnote/cpp-keyword-intro.html 1. asm asm (指令字符串):允许在 C++ 程序中嵌入汇编代码. 2. a ...
- 2019中山纪念中学夏令营-Day2[JZOJ]
博客的开始,先聊聊代码实现: 每次比赛以后,要有归纳错误的习惯. 错误小结: 1.读入:scanf(“%c”)会读入回车和空格,但cin不会. 2.对于二维数组的输入,不能把m,n搞混了,会引起严重的 ...
- (转)当margin-top、padding-top的值为百分比时是如何计算的?
本文链接:https://blog.csdn.net/qq_27437967/article/details/72625900问题:当margin-top.padding-top的值是百分比时,分别是 ...