题意

给定一个 \(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]最小生成树计数的更多相关文章

  1. P4208 [JSOI2008]最小生成树计数

    现在给出了一个简单无向加权图.你不满足于求出这个图的最小生成树,而希望知道这个图中有多少个不同的最小生成树.(如果两颗最小生成树中至少有一条边不同,则这两个最小生成树就是不同的)输出方案数对31011 ...

  2. [洛谷P4208][JSOI2008]最小生成树计数

    题目大意:有$n$个点和$m$条边(最多有$10$条边边权相同),求最小生成树个数 题解:对于所有最小生成树,每种边权的边数是一样的.于是就可以求出每种边权在最小生成树中的个数,枚举这种边的边集,求出 ...

  3. 洛谷P4208 [JSOI2008]最小生成树计数——题解

    题目传送 前置知识:对于同一个图的所有最小生成树,权值相等的边的数量相同. 可以简单证明一下: 我们可以从kruskal的过程考虑.这个算法把所有边按权值大小从小到大排序,然后按顺序看每条边,只要加上 ...

  4. bzoj1016 [JSOI2008]最小生成树计数

    1016: [JSOI2008]最小生成树计数 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 3517  Solved: 1396[Submit][St ...

  5. BZOJ 1016: [JSOI2008]最小生成树计数( kruskal + dfs )

    不同最小生成树中权值相同的边数量是一定的, 而且他们对连通性的贡献是一样的.对权值相同的边放在一起(至多10), 暴搜他们有多少种方案, 然后乘法原理. ----------------------- ...

  6. 1016: [JSOI2008]最小生成树计数

    1016: [JSOI2008]最小生成树计数 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 6200  Solved: 2518[Submit][St ...

  7. 【BZOJ 1016】 1016: [JSOI2008]最小生成树计数 (DFS|矩阵树定理)

    1016: [JSOI2008]最小生成树计数 Description 现在给出了一个简单无向加权图.你不满足于求出这个图的最小生成树,而希望知道这个图中有多少个不同的最小生成树.(如果两颗最小生成树 ...

  8. 【bzoj1016】[JSOI2008]最小生成树计数

    1016: [JSOI2008]最小生成树计数 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 4863  Solved: 1973[Submit][St ...

  9. bzoj1016: [JSOI2008]最小生成树计数(kruskal+dfs)

    1016: [JSOI2008]最小生成树计数 题目:传送门 题解: 神题神题%%% 据说最小生成树有两个神奇的定理: 1.权值相等的边在不同方案数中边数相等  就是说如果一种方案中权值为1的边有n条 ...

随机推荐

  1. spring ioc 源码分析(三)--finishBeanFactoryInitialization(beanFactory)

    之前的博客专门分析了一个类到BeanDefinition的创建过程,现在分析BeanDefinition 到一个bean的创建过程:从refresh() 方法的---->finishBeanFa ...

  2. Apache shiro权限基本使用

    l shiro框架的核心功能:认证.授权.会话管理.加密 Application Code:应用程序代码,由开发人员负责开发的 Subject:框架提供的接口,代表当前用户对象 SecurityMan ...

  3. 嵌入式arm-linux mips-linux 交叉编译GDB,结合vscode图形化调试使用,coredump定位段错误

    第一部分:使用GDB GDB源码下载路径:http://ftp.gnu.org/gnu/gdb/ 遇到的主要难点: 选择合适的GDB源码版本 我的mips-linux交叉编译器不支持C++11特性,所 ...

  4. Java知识日常收集整理001Java获取变量的数据类型的实现方法

    一.具体情况区分 对于简单类型变量,是无法直接获得变量类型的:要想获取,必须自定义函数进行返回. 对于包装类型变量,是可以直接获得的,变量名称.getClass().getName(); 二.代码实现 ...

  5. sysfs是什么??

    来源:https://blog.csdn.net/qq_36412526/article/details/83751520 第一次接触:sysfs, 这里记录过程: 原文:Documenttation ...

  6. ASP。netcore,Angular2 CRUD动画使用模板包,WEB API和EF 1.0.1

    下载Angular2ASPCORE.zip - 1 MB 介绍 在本文中,让我们看看如何创建一个ASP.NET Core CRUD web应用程序与Angular2动画使用模板包,web API和EF ...

  7. Oracle 存储过程解锁及表解锁和停止执行

    查看进程: select * from v$process; 根据存储过程名称查找是否被锁: select * FROM dba_ddl_locks where name =upper('sp_1') ...

  8. 在Windows7系统中设置虚拟内存大小

    当我们的电脑物理内存空间不够用时,操作系统就会自动从硬盘空间上分出一块空间来当内存使用,这就是虚拟内存.可以说虚拟内存是物理内存的补充,是备用的物理内存.一般来说,如果电脑里的程序不多,占用内存资源不 ...

  9. golang Gin framework with websocket

    概述 golang websocket 库 示例 后端 前端 结论 概述 对于 golang 的 web 开发, 之前写过 2 篇 blog, 分别介绍了: 在 Gin 框架下, 各类 http AP ...

  10. Node.js安装及环境配置 for winer

    Node.js安装及环境for Windows 一.安装环境 1.本机系统:Windows 10 Pro(64位) 2.Node.js:v6.9.2LTS(64位) 二.安装Node.js步骤 1.下 ...