博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ4766:文艺计算姬(矩阵树定理)
阅读量:7012 次
发布时间:2019-06-28

本文共 1441 字,大约阅读时间需要 4 分钟。

Description

"奋战三星期,造台计算机"。小W响应号召,花了三星期造了台文艺计算姬。文艺计算姬比普通计算机有更多的艺术细胞。
普通计算机能计算一个带标号完全图的生成树个数,而文艺计算姬能计算一个带标号完全二分图的生成树个数。
更具体地,给定一个一边点数为n,另一边点数为m,共有n*m条边的带标号完全二分图K_{n,m},计算姬能快速算出其生成树个数。
小W不知道计算姬算的对不对,你能帮助他吗?

Input

仅一行三个整数n,m,p,表示给出的完全二分图K_{n,m}
1 <= n,m,p <= 10^18

Output

仅一行一个整数,表示完全二分图K_{n,m}的生成树个数,答案需要模p。

Sample Input

2 3 7

Sample Output

5

Solution

首先先把(度数矩阵-邻接矩阵)搞出来,这里以样例为例。

$
\left\{
\begin{matrix}
3&0&-1&-1&-1\\
0&3&-1&-1&-1\\
-1&-1&2&0&0\\
-1&-1&0&2&0\\
-1&-1&0&0&2\\
\end{matrix}
\right\}
$
按照求矩阵树的方法随便删掉一行一列,这里删掉了最后一行和最后一列。
$
\left\{
\begin{matrix}
3&0&-1&-1\\
0&3&-1&-1\\
-1&-1&2&0\\
-1&-1&0&2\\
\end{matrix}
\right\}
$
把前$n-1$行和后$m-1$行都加到第$n$行
$
\left\{
\begin{matrix}
3&0&-1&-1\\
1&1&0&0\\
-1&-1&2&0\\
-1&-1&0&2\\
\end{matrix}
\right\}
$
用第$n$行的去加到后面$m-1$行上,把$-1$给消掉。
$
\left\{
\begin{matrix}
3&0&-1&-1\\
1&1&0&0\\
0&0&2&0\\
0&0&0&2\\
\end{matrix}
\right\}
$
这样的话这个矩阵的行列式显然就是$m^{n-1}n^{m-1}$了。
记得快速乘。

Code

1 #include
2 #define LL long long 3 using namespace std; 4 5 LL n,m,p; 6 7 LL Mul(LL a,LL b) 8 { 9 LL tmp=a*b-(LL)((long double)a*b/p+0.1)*p;10 return tmp<0?tmp+p:tmp;11 }12 13 LL Qpow(LL a,LL b)14 {15 LL ans=1;16 while (b)17 {18 if (b&1) ans=Mul(ans,a);19 a=Mul(a,a); b>>=1;20 }21 return ans;22 }23 24 int main()25 {26 scanf("%lld%lld%lld",&n,&m,&p);27 printf("%lld",Mul(Qpow(n,m-1),Qpow(m,n-1)));28 }

转载于:https://www.cnblogs.com/refun/p/10175635.html

你可能感兴趣的文章
构建千万级别用户的后台数据库架构
查看>>
学习SLF4J,LOG4J
查看>>
模拟实现C语言strncpy.strncat.strncmp函数
查看>>
右侧扩展栏效果实现
查看>>
//输出电脑系统的年月日星期以及时间
查看>>
Mysql主从复制及读写分离的实现
查看>>
哈希表的静态,动态,以及key/value形式
查看>>
基于memcached-session-manager的tomcat session共享集群
查看>>
关于Apache的一些配置
查看>>
清除Windows Server Backup 备份副本
查看>>
MySQL数据库性能优化之一(缓存参数优化)
查看>>
常见Sqlite管理工具
查看>>
Web工程加入日志组件:slf4j+logback
查看>>
MySQL中创建、删除表和库,添加、删除记录
查看>>
cacti snmp error
查看>>
Cloudflare能检测HTTPS并拦截MITMEngine
查看>>
独家解密:阿里大规模数据中心性能分析
查看>>
一线实践 | 借助混沌工程工具 ChaosBlade 构建高可用的分布式系统
查看>>
智能合约就是程序
查看>>
python项目实战:实现游戏敏感词汇过滤器
查看>>