题目背景


又到了某任*堂开关中学一年一度的自主招生考试的时间了,在考试完后许多家长决定带着自己的孩子参观一下这所距千年名校还有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. InvocationTargetException异常的深入研究-servlet的setAttribute与getAttribute

    在某项目中,前端jsp传入的用户id数据通过session域传入后台servlet进行处理的过程中,无意间出现了InvocationTargetException异常 前端部分代码如下:测试代码,非原 ...

  2. 一个PHP开发APP接口的视频教程

    感觉php做接口方面的教程很少,无意中搜到了这个视频教程,希望能给一些人带来帮助http://www.imooc.com/learn/163

  3. Animate.css_css3动画库介绍

    插件描述:Animate.css内置了很多典型的css3动画,兼容性好使用方便. Animate.css是一个有趣的,跨浏览器的css3动画库.很值得我们在项目中引用. 用法 1.首先引入animat ...

  4. 联玛客(W 笔试)

    纸质算法题 1. 输入数据:1.3.2.4.8... 输出数据:3.1.4.2.8... 找出规律,写出一个程序求解,并附上时间复杂度和空间复杂度 我的答案: 规律一:奇偶位互换 假设输入数据长度为5 ...

  5. How to Slove MB SD C4 Cannot Access OBD2 Port

    When using china clone mb sd connect compact 4 Multiplexer, it could not link to the car computer, M ...

  6. Open Cascade:使用鼠标画线

    Open Cascade:使用鼠标画线 在View类文件中创建以下代码: 1.创建鼠标消息: afx_msg void OnLButtonDown(UINT nFlags, CPoint point) ...

  7. 20道必须掌握的C++面试题

    20道必须掌握的C++面试题 在面试C++方面的工作时,经常会遇到各种面试题,这对应聘人员的知识掌握能力要求较高.本文将为大家带来的就是20道必须掌握的C++面试题,不要错过哦! 问1:请用简单的语言 ...

  8. QT +菜单栏和工具栏

    #include "mainwindow.h" #include <QMenuBar>//菜单栏需要的头文件 #include <QMenu>//菜单 #i ...

  9. python之bool (布尔值)

    用途: ​ 判断真假 识记: 空的字符串是False,非空的就是True 0 是False,非0的都是True 格式: True False 布尔值转换成字符串: print(type(str(Tru ...

  10. django 127.0.0.1 将您重定向的次数过多

    "GET /?next=/%3Fnext%3D/%253Fnext%253D/ HTTP/1.1" 302 0 solution reference from django.con ...