题目背景


又到了某任*堂开关中学一年一度的自主招生考试的时间了,在考试完后许多家长决定带着自己的孩子参观一下这所距千年名校还有890周年的百年学校;

题目描述


这所学校的布局非常奇怪,是一个由N 个点M 条边构成的无向图,既然来了肯定要把学校逛完.

家长们的思路也非常清奇,带孩子走过的边刚好有n-1 条,并且这n-1 条边构成了原图的一棵生成树。 每条边上所能看到的风景有两个固定的价值A 和B ,走了这条边就只能得到这条边的A价值,而没走的边就只能得到那条边的B 价值,每个家庭最终看到的总价值的具体算法为所有走过的边的A 值的乘积*所有没选的边的B 值的乘积;

zackzh当年并没有游览学校,所以他想知道所有可能走出的所有生成树的总价值和MOD1e9+7;

输入输出格式


输入格式:

第11 行两个个整数N 与M ,接下来M 每行四个整数U ,V ,A ,B 描述一条边;

输出格式:


仅一个整数,即可能走出的所有生成树的总价值和MOD1e9+7 ;

输入输出样例


输入样例#1:


输出样例#1: 


说明


N<=500,M<=200000,A,B<=10;

%10数据满足 N = 3 M = 12

另有%10数据满足N = 12 M = 50

另有%10数据满足N = 500 M = 499

另有%10数据满足N = 500 M = 1000

另有%10数据满足 N = 500 M = 100000

剩下%50 数据满足N = 500 M = 200000

对于%100数据满足1<=A,B<=10

分析:


就是考察了一下对矩阵树定理的理解。

在矩阵树定理树一棵生成树的权值求出的其实是它所有边的乘积,

矩阵树定理就求出了所有生成树的权值和。

将边权设为1就可以求出生成树个数。

这道题我们可以把生成树权值设为A / B,设所有边B乘积为ret

那么最终答案就是矩阵树定理求出来的值*ret了。

对于逆元用费马小定理就好了。

AC代码:


# include <cstdio>
# include <cstring>
# include <vector>
# include <iostream>
using namespace std;
typedef long long LL;
const LL mod = 1e9 + ;
void read(int &x)
{
x = ;char i = getchar();
while(!isdigit(i))i = getchar();
while(isdigit(i))x = x * + i - '',i = getchar();
}
LL Inc(LL x,LL y){return (x + y + mod) % mod;}
LL Dec(LL x,LL y){return (x - y + mod) % mod;}
LL mul(LL x,LL y){return (x * y % mod + mod) % mod;}
LL a[][],ans;int n,m;
LL inv(LL x)
{ LL v = ;
for(int i = mod - ;i;i >>= ,x = mul(x,x))if(i & )v = mul(v,x);
return v;
}
LL Gauss()
{
LL ret = ;
for(int i = ;i < n;i++)
{
int p = i;
for(int j = i + ;j < n;j++)if(a[j][i] > a[p][i])p = j;
if(p != i)swap(a[p],a[i]),ret = -ret;
if(!a[i][i])return ;
for(int j = i + ;j < n;j++)
{
LL t = a[j][i] * inv(a[i][i]) % mod;
for(int k = ;k < n;k++)
a[j][k] = Dec(a[j][k],mul(t,a[i][k]));
}
}
for(int i = ;i < n;i++)ret = mul(ret,a[i][i]);
return ret;
}
int main()
{
read(n);read(m);int x,y,A,B;LL ret = ;
for(int i = ;i <= m;i++)
{
read(x);read(y);read(A);read(B);
ret = mul(ret,B);
a[x][y] = Inc(a[x][y],mul(A,inv(B)));
a[y][x] = Inc(a[y][x],mul(A,inv(B)));
}
for(int i = ;i <= n;i++)
for(int j = ;j <= n;j++)
if(i != j)a[i][i] = Inc(a[i][i],a[i][j]),a[i][j] = mul(-,a[i][j]);
printf("%lld\n",mul(ret,Gauss()));
}

[洛谷U22156]未曾届到游览(矩阵树定理)的更多相关文章

  1. 洛谷4208 JSOI2008最小生成树计数(矩阵树定理+高斯消元)

    qwq 这个题目真的是很好的一个题啊 qwq 其实一开始想这个题,肯定是无从下手. 首先,我们会发现,对于无向图的一个最小生成树来说,只有当存在一些边与内部的某些边权值相同的时候且能等效替代的时候,才 ...

  2. 洛谷 P4336 黑暗前的幻想乡 —— 容斥+矩阵树定理

    题目:https://www.luogu.org/problemnew/show/P4336 当作考试题了,然而没想出来,呵呵. 其实不是二分图完美匹配方案数,而是矩阵树定理+容斥... 就是先放上所 ...

  3. 洛谷 P5406 - [THUPC2019]找树(FWT+矩阵树定理)

    题面传送门 首先看到这道题你必须要有一个很清楚的认识:这题新定义的 \(\oplus\) 符号非常奇怪,也没有什么性质而言,因此无法通过解决最优化问题的思路来解决这个问题,只好按照计数题的思路来解决, ...

  4. 洛谷 P3377 【模板】左偏树(可并堆)

    洛谷 P3377 [模板]左偏树(可并堆) 题目描述 如题,一开始有N个小根堆,每个堆包含且仅包含一个数.接下来需要支持两种操作: 操作1: 1 x y 将第x个数和第y个数所在的小根堆合并(若第x或 ...

  5. 洛谷P3377 【模板】左偏树(可并堆) 题解

    作者:zifeiy 标签:左偏树 这篇随笔需要你在之前掌握 堆 和 二叉树 的相关知识点. 堆支持在 \(O(\log n)\) 的时间内进行插入元素.查询最值和删除最值的操作.在这里,如果最值是最小 ...

  6. 【BZOJ3534】重建(矩阵树定理)

    [BZOJ3534]重建(矩阵树定理) 题面 BZOJ 洛谷 题解 这.... 矩阵树定理神仙用法???? #include<iostream> #include<cmath> ...

  7. 【BZOJ5297】【CQOI2018】社交网络(矩阵树定理)

    [BZOJ5297][CQOI2018]社交网络(矩阵树定理) 题面 BZOJ 洛谷 Description 当今社会,在社交网络上看朋友的消息已经成为许多人生活的一部分.通常,一个用户在社交网络上发 ...

  8. 【BZOJ4031】小Z的房间(矩阵树定理)

    [BZOJ4031]小Z的房间(矩阵树定理) 题面 BZOJ 洛谷 Description 你突然有了一个大房子,房子里面有一些房间.事实上,你的房子可以看做是一个包含n*m个格子的格状矩形,每个格子 ...

  9. 【SPOJ】Highways(矩阵树定理)

    [SPOJ]Highways(矩阵树定理) 题面 Vjudge 洛谷 题解 矩阵树定理模板题 无向图的矩阵树定理: 对于一条边\((u,v)\),给邻接矩阵上\(G[u][v],G[v][u]\)加一 ...

随机推荐

  1. dataTables去掉搜索框和每页多少条框体,解决Cannot reinitialise DataTable问题

    $('#example').DataTable({ searching:false, //去掉搜索框 bLengthChange:false,//去掉每页多少条框体 "language&qu ...

  2. JDBC性能优化篇

    系统性能. 少用Metadata方法     与其它的JDBC方法相比, 由ResultSet对象生成的metadata对象的相对来说是很慢的. 应用程序应该缓存从ResultSet返回的metada ...

  3. hdu 3232 Crossing Rivers 过河(数学期望)

    题意:你在点A,目的地是点B,A和B的距离为D.中间隔了好多条河(所有河不会重叠),每条河有3个参数(P,L,V),其中P表示距离A点的长度,L表示河的长度,V表示河里的船的速度.假设每条河中仅有1条 ...

  4. vijos 1772 巧妙填数

    描述 将1,2,\cdots,91,2,⋯,9共99个数分成三组,分别组成三个三位数,且使这三个三位数构成1:2:31:2:3的比例. 试求出所有满足条件的三个三位数.例如:三个三位数192,384, ...

  5. edquota - 编辑用户配额

    SYNOPSIS(总览) edquota [ -p proto-username ] [ -u | -g ] username... edquota [ -u | -g ] -t DESCRIPTIO ...

  6. 一条update语句优化小记

    遇到性能问题的sql如下: sql1: UPDATE amlclientlevel a SET    a.client_value = (SELECT l.client_value           ...

  7. ROI(Region of Interesting)

    ROI(region of interest),感兴趣区域.机器视觉.图像处理中,从被处理的图像以方框.圆.椭圆.不规则多边形等方式勾勒出需要处理的区域,称为感兴趣区域,ROI.在Halcon.Ope ...

  8. 【软件构造】第八章第三节 代码调优的设计模式和I/O

    第八章第三节 代码调优的设计模式和I/O 本节学习如何通过对代码的修改,消除性能瓶颈,提高系统性能?——代码调优.面向性 能的设计模式 Outline Java调优 代码调优的概念 单例模式(Sing ...

  9. tmp note

    cat file.txt > /dev/null 2>&1 丢弃错误和标准输出 systemctl isolate multi-user.target 切换回多用户命令行模式 sy ...

  10. 02.28 day03

    print(1 or 3 > 2 and 4 < 5 or 6 and 2 < 7)## while True:# print(11)# print(22)# # break# # ...