BZOJ4766 文艺计算姬

本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作。

本文作者:ljh2000
作者博客:http://www.cnblogs.com/ljh2000-jump/
转载请注明出处,侵权必究,保留最终解释权!

 

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

 

 

正解:$Matrix-Tree$定理

解题报告:

  学了$Matrix-Tree$之后,这道题推推式子就可以辣。

  先用小数据建图,再建出基尔霍夫矩阵,高斯消元,对答案打个表(不要取模!不然规律就GG了...)看看规律,

  发现答案就是$n^{m-1} * m^{n-1}$,用快速幂$+$快速乘法就可以了。

//It is made by ljh2000
//有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <vector>
#include <queue>
#include <map>
#include <set>
#include <string>
#include <complex>
#include <bitset>
using namespace std;
typedef long long LL;
typedef long double LB;
typedef complex<double> C;
const double pi = acos(-1);
LL n,m,p,ans;

inline LL getint(){
    LL w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
    if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
}

inline LL cheng(LL x,LL y){
	LL r=0;
	while(y>0) {
		if(y&1) r+=x,r%=p;
		x+=x; x%=p;
		y>>=1;
	}
	return r;
}

inline LL fast_pow(LL x,LL y){
	LL r=1;
	while(y>0) {
		if(y&1) r=cheng(r,x);
		x=cheng(x,x);
		y>>=1;
	}
	return r;
}

inline void work(){
	n=getint(); m=getint(); p=getint();
	ans=cheng(fast_pow(n,m-1),fast_pow(m,n-1));
	cout<<ans;
}

int main()
{
    work();
    return 0;
}
//有志者,事竟成,破釜沉舟,百二秦关终属楚;苦心人,天不负,卧薪尝胆,三千越甲可吞吴。

  

原文地址:https://www.cnblogs.com/ljh2000-jump/p/6516320.html