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
n和m较小时可以用矩阵树定理
现在打表找规律
发现 Kn,m=nm−1∗mn−1
记得使用快速乘
#include<bits/stdc++.h> using namespace std; #define For(i,n) for(int i=1;i<=n;i++) #define Fork(i,k,n) for(int i=k;i<=n;i++) #define ForkD(i,k,n) for(int i=n;i>=k;i--) #define Rep(i,n) for(int i=0;i<n;i++) #define ForD(i,n) for(int i=n;i;i--) #define RepD(i,n) for(int i=n;i>=0;i--) #define Forp(x) for(int p=pre[x];p;p=next[p]) #define Forpiter(x) for(int &p=iter[x];p;p=next[p]) #define Lson (o<<1) #define Rson ((o<<1)+1) #define MEM(a) memset(a,0,sizeof(a)); #define MEMI(a) memset(a,0x3f,sizeof(a)); #define MEMi(a) memset(a,128,sizeof(a)); #define MEMx(a,b) memset(a,b,sizeof(a)); #define INF (0x3f3f3f3f) #define pb push_back #define mp make_pair #define fi first #define se second #define vi vector<int> #define pi pair<int,int> #define SI(a) ((a).size()) #define Pr(kcase,ans) printf("Case #%d: %lld\n",kcase,ans); #define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl; #define PRi2D(a,n,m) For(i,n) { \ For(j,m-1) cout<<a[i][j]<<' ';\ cout<<a[i][m]<<endl; \ } #pragma comment(linker, "/STACK:102400000,102400000") #define ALL(x) (x).begin(),(x).end() typedef long long ll; typedef long double ld; typedef unsigned long long ull; ll n,m,F; ll add(ll a,ll b){return (a+b)%F;} ll sub(ll a,ll b){return ((a-b)%F+F)%F;} void upd(ll &a,ll b){a=(a%F+b%F)%F;} inline int read() { int x=0,f=1; char ch=getchar(); while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();} while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();} return x*f; } ll mul(ll x,ll y) { ll p=0; while(y) { if (y&1) p=(p+x)%F; x=(x*2)%F; y>>=1; } return p; } ll pow2(ll a,ll b) { if (!b) return 1%F; ll p=pow2(a,b/2); p=mul(p,p); if (b&1) p=mul(p,a); return p; } int main() { // freopen("bzoj4766.in","r",stdin); // freopen(".out","w",stdout); cin>>n>>m>>F; cout<<mul(pow2(n,m-1),pow2(m,n-1))<<endl; return 0; }