Minimum Spanning Tree

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)

【Problem Description】
XXX is very interested in algorithm. After learning the Prim algorithm and Kruskal algorithm of minimum spanning tree, XXX finds that there might be multiple solutions. Given an undirected weighted graph with n (1<=n<=100) vertexes and m (0<=m<=1000) edges, he wants to know the number of minimum spanning trees in the graph.
【Input】
There are no more than 15 cases. The input ends by 0 0 0. For each case, the first line begins with three integers --- the above mentioned n, m, and p. The meaning of p will be explained later. Each the following m lines contains three integers u, v, w (1<=w<=10), which describes that there is an edge weighted w between vertex u and vertex v( all vertex are numbered for 1 to n) . It is guaranteed that there are no multiple edges and no loops in the graph.
【Output】
For each test case, output a single integer in one line representing the number of different minimum spanning trees in the graph. The answer may be quite large. You just need to calculate the remainder of the answer when divided by p (1<=p<=1000000000). p is above mentioned, appears in the first line of each test case.
【Sample Input】

【Sample Output】


【题意】

一张无向图,要求求出其中最小生成树的棵树。
【分析】

生成树计数可以使用Matrix-Tree定理解决,本题最主要的区别是有了一个最小生成树的额外条件。

首先考虑一下如何得到最小生成树。

Kruskal算法的基本思想是,按照边长排序,然后不断将短边加入集合,最终一步如果能成功把n-1条边都加入同一个集合,则找到了最小生成树。在维护集合时,可以使用并查集来快速处理。

如果把Kruskal的过程按照边长划分成多个阶段,实际上是处理了所有短边的连通性之后继续处理下一个长度的边的连通性,并依次继续处理剩下边的连通性。然后我们可以发现,不同长度的边之间的连通性互不影响!!!

假设存在n1条长度为c1的边,n2条长度为c2的边...则Kruskal首先处理c1边的连通性,然后处理c2边的连通性,对于c1边的连通性的处理可能有多种方案,即从n1条边中取出一定数量的边构成最大连通图,但是最终处理完之后的结果对于c2来说是完全一样的。因此算法就出来了,在Kruskal的基础上,使用Matrix-Tree定理处理每个阶段生成树的种数,最后将所有阶段的结果相乘即可。

具体实现为:

在Kruskal的基础上,每完成一个阶段(检查完一个长度),就将所有遍历过的点缩成一个点,然后用Matrix-Tree定理计算该点与下一组点组成的连通图中生成树的个数。最终把每一个阶段的结果相乘即可。

 /* ***********************************************
MYID : Chen Fan
LANG : G++
PROG : Counting_MST_HDU4408
************************************************ */ #include <vector>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <bitset>
#define N 405
#define M 4005 using namespace std; typedef struct nod
{
int a,b,c;
} node;
node edge[M]; bool op(node a,node b)
{
return a.c<b.c;
} int n,m,o,fa[N],ka[N];
long long ans,mod,gk[N][N],kir[N][N]; bitset<N> flag;
vector<int> gra[N]; int getfather(int x,int father[])
{
if (father[x]!=x) father[x]=getfather(father[x],father);
return father[x];
} long long det(long long a[][N],int n) //Matrix-Tree定理求Kirchhoff矩阵
{
for (int i=;i<n;i++)
for (int j=;j<n;j++) a[i][j]%=mod;
long long ret=;
for (int i=;i<n;i++)
{
for (int j=i+;j<n;j++)
while (a[j][i])
{
long long t=a[i][i]/a[j][i];
for (int k=i;k<n;k++) a[i][k]=(a[i][k]-a[j][k]*t)%mod;
for (int k=i;k<n;k++) swap(a[i][k],a[j][k]);
ret=-ret;
}
if (a[i][i]==) return ;
ret=ret*a[i][i]%mod;
//ret%=mod;
}
return (ret+mod)%mod;
} void matrix_tree()
{
for (int i=;i<=n;i++) //根据访问标记找出连通分量
if (flag[i])
{
gra[getfather(i,ka)].push_back(i);
flag[i]=;
}
for (int i=;i<=n;i++)
if (gra[i].size()>) //枚举连通分量
{
memset(kir,,sizeof(kir)); int len=gra[i].size();
for (int a=;a<len;a++)
for (int b=a+;b<len;b++)
{
int la=gra[i][a],lb=gra[i][b];
kir[b][a]-=gk[la][lb];
kir[a][b]=kir[b][a];
kir[a][a]+=gk[la][lb];
kir[b][b]+=gk[la][lb];
} //构造矩阵 long long ret=det(kir,len);
ret%=mod;
ans=(ans*ret%mod)%mod; for (int a=;a<len;a++) fa[gra[i][a]]=i;
}
for (int i=;i<=n;i++) //连通图缩点+初始化
{
fa[i]=getfather(i,fa);
ka[i]=fa[i];
gra[i].clear();
}
} int main()
{
freopen("4408.txt","r",stdin); while (scanf("%d%d%lld",&n,&m,&mod)==)
{ if (n==&&m==&&mod==) break;
for (int i=;i<=m;i++) scanf("%d%d%d",&edge[i].a,&edge[i].b,&edge[i].c);
sort(&edge[],&edge[m+],op); for (int i=;i<=n;i++) gra[i].clear();
for (int i=;i<=n;i++)
{
fa[i]=i;
ka[i]=i;
}
flag.reset();
memset(gk,,sizeof(gk));
ans=;
o=edge[].c;
for (int i=;i<=m;i++)
{
int pa=getfather(edge[i].a,fa),pb=getfather(edge[i].b,fa);
if (pa!=pb)
{
flag[pa]=;
flag[pb]=; //访问标记
ka[getfather(pa,ka)]=getfather(pb,ka);
gk[pa][pb]++;
gk[pb][pa]++; //邻接矩阵
}
if (i==m||edge[i+].c!=o) //所有相同的边并成一组
{
matrix_tree();
o=edge[i+].c;
}
} bool done=true;
for (int i=;i<=n;i++)
if(ka[i]!=ka[i-])
{
done=false;
break;
}
if (!done) printf("0\n");
else
{
ans%=mod;
printf("%lld\n",ans);
}
} return ;
}

HDU 4408 Minimum Spanning Tree 最小生成树计数的更多相关文章

  1. hdu 4408 Minimum Spanning Tree

    Problem Description XXX is very interested in algorithm. After learning the Prim algorithm and Krusk ...

  2. 数据结构与算法分析–Minimum Spanning Tree(最小生成树)

    给定一个无向图,如果他的某个子图中,任意两个顶点都能互相连通并且是一棵树,那么这棵树就叫做生成树(spanning tree). 如果边上有权值,那么使得边权和最小的生成树叫做最小生成树(MST,Mi ...

  3. 【算法】关于图论中的最小生成树(Minimum Spanning Tree)详解

    本节纲要 什么是图(network) 什么是最小生成树 (minimum spanning tree) 最小生成树的算法 什么是图(network)? 这里的图当然不是我们日常说的图片或者地图.通常情 ...

  4. 【HDU 4408】Minimum Spanning Tree(最小生成树计数)

    Problem Description XXX is very interested in algorithm. After learning the Prim algorithm and Krusk ...

  5. Educational Codeforces Round 3 E. Minimum spanning tree for each edge 最小生成树+树链剖分+线段树

    E. Minimum spanning tree for each edge time limit per test 2 seconds memory limit per test 256 megab ...

  6. 说说最小生成树(Minimum Spanning Tree)

    minimum spanning tree(MST) 最小生成树是连通无向带权图的一个子图,要求 能够连接图中的所有顶点.无环.路径的权重和为所有路径中最小的. graph-cut 对图的一个切割或者 ...

  7. 最小生成树(Minimum Spanning Tree)——Prim算法与Kruskal算法+并查集

    最小生成树——Minimum Spanning Tree,是图论中比较重要的模型,通常用于解决实际生活中的路径代价最小一类的问题.我们首先用通俗的语言解释它的定义: 对于有n个节点的有权无向连通图,寻 ...

  8. 多校 HDU - 6614 AND Minimum Spanning Tree (二进制)

    传送门 AND Minimum Spanning Tree Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 131072/131072 ...

  9. CF# Educational Codeforces Round 3 E. Minimum spanning tree for each edge

    E. Minimum spanning tree for each edge time limit per test 2 seconds memory limit per test 256 megab ...

随机推荐

  1. 安卓无法生成R文件原因

    原因个人总结出来: 清单文件报错,则无法生成R文件 gen和bin目录可以删除

  2. linux 配置tomcat服务器

    1. 找到tomcat安装包    find / -name apache-tomcat* 2. 解压包 tar zxvf apache-tomcat-7.0.67.tar.gz rpm -ivh j ...

  3. Unity3d 物体沿着正七边形轨迹移动

    不对之处,敬请谅解. 1.圆内接正七边形半径 public static float r = 10; 2.存储七个顶点的值 Vector3[] ar = new Vector3[7]; 3.圆心角 s ...

  4. Inno Setup入门(二十)——Inno Setup类参考(6)

    存储框 存储框也是典型的窗口可视化组件,同编辑框类似,可以输入.显示文本,但是和编辑框不同的是,编辑框只能编辑.显示单行文本,而存储框则可以对多行文本进行操作.存储框的类定义如下:< xmlna ...

  5. byte数组如何转为short数组 (转)

    最近在搞毕业设计,做的是有关语音识别的手机应用.在处理音频的过程中,发现需要用short数组处理音频,可能光用byte会越界.但是java读文件没有一次性读到short数组中的api,要么是一个一个读 ...

  6. WPF中ListBox控件选择多个数据项

    XAML: <Window x:Class="WpfApplication1.MainWindow" xmlns="http://schemas.microsoft ...

  7. jsp提交表单问题

    以form形式提交的话 String usernameInForm = hreq.getParameter("username");String passwordInForm = ...

  8. oracle恢复一个数据表的方法

    今天提交给客户方一个sql脚本去跟新历史数据,结果客户那边的部署人员犯了一个错误,直接拿系统账号去部署,结果第一段代码没有执行成功,结果第二段代码却执行成功了,并且已经提交了的,....由于事前没有备 ...

  9. weak和assign区别

    weak比assign多了一个功能,当对象消失后自动把指针变成nil haofanazenmeban[4002:406590] controller:<SecondViewController: ...

  10. 文本格式ANSI,Unicode等有什么区别

    首先DBCS是亚洲的字符集,包含了ANSI,ANSI也就是ASCII值为0-255之间的字符,当字符为ANSI时,存放于文件中占用的是一个字节.如果是非ANSI的呢,则占用两字节.用VB的ASC函数可 ...