Luogu P4208 [JSOI2008]最小生成树计数
题意
给定一个 \(n\) 个点 \(m\) 条边的图,求最小生成树的个数。
\(\texttt{Data Range:}1\leq n\leq 100,1\leq m\leq 10^4\)
题解
一道好题。
根据本题后面提供的与那题正解没什么关联的方法可知,这个操作过程是这样的:
首先求出原图的某一个最小生成树,接下来考虑从小到大枚举最小生成树上边的边权 \(w\)。
将最小生成树上边权不为 \(w\) 的边保留下来进行缩点,接下来再连上不在最小生成树中边权为 \(w\) 的边。
这个时候会建出一个无向图,对每一个可能的 \(w\) 建出的图求一下生成树的个数乘起来即可。
证明的话可以利用 Kruskal 的性质,求生成树的时候使用 Matrix-Tree 定理即可。
模合数求行列式的方法是辗转相消,每一次需要交换两行,行列式要乘上 \(-1\),实在不理解可以看我代码
容易看出这个东西的复杂度是 \(O(n^3\log n)\) 的。
代码
#include<bits/stdc++.h>
using namespace std;
typedef int ll;
typedef long long int li;
const ll MAXN=251,MOD=31011;
struct EdgeForKruskal{
ll from,to,dist;
inline bool operator <(const EdgeForKruskal &rhs)const
{
return this->dist<rhs.dist;
}
};
EdgeForKruskal ed[2051],tree[MAXN];
ll n,m,x,y,z,kk,res=1,totw;
ll ffa[MAXN],bel[MAXN],mat[MAXN][MAXN],wt[MAXN];
inline ll read()
{
register ll num=0,neg=1;
register char ch=getchar();
while(!isdigit(ch)&&ch!='-')
{
ch=getchar();
}
if(ch=='-')
{
neg=-1;
ch=getchar();
}
while(isdigit(ch))
{
num=(num<<3)+(num<<1)+(ch-'0');
ch=getchar();
}
return num*neg;
}
inline void add(ll x,ll y)
{
mat[x][y]--,mat[y][x]--,mat[x][x]++,mat[y][y]++;
}
inline ll find(ll x)
{
return x==ffa[x]?x:ffa[x]=find(ffa[x]);
}
inline void setup(ll n)
{
for(register int i=1;i<=n;i++)
{
ffa[i]=i;
}
}
inline void merge(ll x,ll y)
{
ll fx=find(x),fy=find(y);
fx!=fy?ffa[fy]=fx:1;
}
inline ll Kruskal()
{
ll tott=0;
for(register int i=1;i<=m;i++)
{
if(find(ed[i].from)!=find(ed[i].to))
{
merge(ed[i].from,ed[i].to),tree[++tott]=ed[i];
if(wt[totw]!=ed[i].dist)
{
wt[++totw]=ed[i].dist;
}
if(tott==n-1)
{
break;
}
}
}
return tott==n-1;
}
inline void mergePoint(ll wt)
{
for(register int i=1;i<n;i++)
{
tree[i].dist!=wt?merge(tree[i].from,tree[i].to):(void)1;
}
}
inline ll det(ll n)
{
ll res=1,sgn=1,cof;
for(register int i=1;i<=n;i++)
{
for(register int j=1;j<=n;j++)
{
mat[i][j]=(mat[i][j]+MOD)%MOD;
}
}
for(register int i=1;i<=n;i++)
{
for(register int j=i+1;j<=n;j++)
{
while(mat[j][i])
{
cof=mat[i][i]/mat[j][i];
for(register int k=i;k<=n;k++)
{
mat[i][k]=(mat[i][k]-(li)cof*mat[j][k]%MOD+MOD)%MOD;
}
swap(mat[i],mat[j]),sgn*=-1;
}
}
}
for(register int i=1;i<=n;i++)
{
res=(li)res*mat[i][i]%MOD;
}
return sgn==1?res:MOD-res;
}
inline ll calc(ll wt)
{
ll blk=0;
memset(mat,0,sizeof(mat)),setup(n),mergePoint(wt);
for(register int i=1;i<=n;i++)
{
find(i)==i?bel[i]=++blk:1;
}
for(register int i=1;i<=n;i++)
{
bel[i]=bel[find(i)];
}
for(register int i=1;i<=m;i++)
{
ed[i].dist==wt?add(bel[ed[i].from],bel[ed[i].to]):(void)1;
}
return det(blk-1);
}
int main()
{
n=read(),m=read();
for(register int i=1;i<=m;i++)
{
x=read(),y=read(),z=read(),ed[i]=(EdgeForKruskal){x,y,z};
}
sort(ed+1,ed+m),setup(n);
if(!Kruskal())
{
return puts("0"),0;
}
for(register int i=1;i<=totw;i++)
{
res=res*calc(wt[i])%MOD;
}
printf("%d\n",res);
}
Luogu P4208 [JSOI2008]最小生成树计数的更多相关文章
- P4208 [JSOI2008]最小生成树计数
现在给出了一个简单无向加权图.你不满足于求出这个图的最小生成树,而希望知道这个图中有多少个不同的最小生成树.(如果两颗最小生成树中至少有一条边不同,则这两个最小生成树就是不同的)输出方案数对31011 ...
- [洛谷P4208][JSOI2008]最小生成树计数
题目大意:有$n$个点和$m$条边(最多有$10$条边边权相同),求最小生成树个数 题解:对于所有最小生成树,每种边权的边数是一样的.于是就可以求出每种边权在最小生成树中的个数,枚举这种边的边集,求出 ...
- 洛谷P4208 [JSOI2008]最小生成树计数——题解
题目传送 前置知识:对于同一个图的所有最小生成树,权值相等的边的数量相同. 可以简单证明一下: 我们可以从kruskal的过程考虑.这个算法把所有边按权值大小从小到大排序,然后按顺序看每条边,只要加上 ...
- bzoj1016 [JSOI2008]最小生成树计数
1016: [JSOI2008]最小生成树计数 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 3517 Solved: 1396[Submit][St ...
- BZOJ 1016: [JSOI2008]最小生成树计数( kruskal + dfs )
不同最小生成树中权值相同的边数量是一定的, 而且他们对连通性的贡献是一样的.对权值相同的边放在一起(至多10), 暴搜他们有多少种方案, 然后乘法原理. ----------------------- ...
- 1016: [JSOI2008]最小生成树计数
1016: [JSOI2008]最小生成树计数 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 6200 Solved: 2518[Submit][St ...
- 【BZOJ 1016】 1016: [JSOI2008]最小生成树计数 (DFS|矩阵树定理)
1016: [JSOI2008]最小生成树计数 Description 现在给出了一个简单无向加权图.你不满足于求出这个图的最小生成树,而希望知道这个图中有多少个不同的最小生成树.(如果两颗最小生成树 ...
- 【bzoj1016】[JSOI2008]最小生成树计数
1016: [JSOI2008]最小生成树计数 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 4863 Solved: 1973[Submit][St ...
- bzoj1016: [JSOI2008]最小生成树计数(kruskal+dfs)
1016: [JSOI2008]最小生成树计数 题目:传送门 题解: 神题神题%%% 据说最小生成树有两个神奇的定理: 1.权值相等的边在不同方案数中边数相等 就是说如果一种方案中权值为1的边有n条 ...
随机推荐
- tomcat在linux下安装
1.下载地址: https://tomcat.apache.org/download-90.cgi 2.上传linux 3.查看是否上传成功 4.解压: 5.进入后,查看README.md文件,可以查 ...
- 普转提Day2
T1 给定一个区间,求这个区间中只有一个数字与其他数组不相同的数的个数. 给出的区间范围较大,但是要求的数比较少.所以我的想法是这样的:因为这些数只有一个数字和每个数字都相同的数不同,所以考虑将所有数 ...
- spring初始(介绍、核心架构)
1.spring介绍 Spring是个java企业级应用的开源开发框架.主要用来开发Java应用,但是有些扩展是针对构建J2EE平台的web应用.Spring框架目标是简化Java企业级应用开发,并通 ...
- RT-Thread小知识点收集
- 把python文件打包成可执行文件(win10实验成功)
总是有人来找我帮看下工单状态,又懒得写页面展示出来,干脆打包成exe文件好啦 打包很简单,难点在于安装pyinstaller这个依赖包,主要是网络问题~ 我也是参考别人的博文,别人的文章写得很详细,我 ...
- Python 中 pip 工具的安装与使用
pip 是 Python 包管理工具,该工具提供了对Python 包的查找.下载.安装.卸载的功能. 目前如果你在 python.org 下载最新版本的安装包,则是已经自带了该工具. Python 2 ...
- day14 Pyhton学习
一.迭代器-概念 可迭代协议:内部含有__iter__方法的值/变量都是可迭代的 如何得到一个迭代器:可迭代变量.__iter__()返回一个迭代器 迭代器协议:内部含有__iter__方法和__ne ...
- C语言实现和 *.ini文件。
本文之前由于技术不到位,写的比较挫,最近花了大半天时间写了一个高级点的版本. 这里是我写的代码,已经上传到github了.跳转到github 主要是使用了链表保存ini文件的内容,在程序运行最初会初始 ...
- Oracle函数总结
<Trunc()> 描 述(实际应用):截取小数或者日期整数 简 介:https://baike.baidu.com/item/trunc/9657216?fr=al ...
- 飞翔---------双重线性dp
题目: 鹰最骄傲的就是翱翔,但是鹰们互相都很嫉妒别的鹰比自己飞的快,更嫉妒其他的鹰比自己飞行的有技巧.于是,他们决定举办一场比赛,比赛的地方将在一个迷宫之中. 这些鹰的起始点被设在一个N*M矩阵的左下 ...