矩阵计数

腾讯游戏开发了一款全新的编程类益智小游戏，榜首的题目是一道关于矩阵的计算，你用多久能计算出来呢？

游戏中给出一个N*M 的矩阵，若其中填入的内容是数字 1~N*M 的排列，求问有多少种不等价的矩阵？

等价矩阵：若一个矩阵 A 可以通过交换其中两行或者两列变成另一个矩阵 B，则称 A 和 B 等价。且若 A 和 B 等价，B和 C等价，则 A和 C 也等价。

答案对 998244353 取模。

说明：

当你计算一个答案需要对某大质数取模的问题时，加减乘都是可以中途取模的，例如 (A+B+C)%mod可以改为 ((A+B)%mod+C)%mod，这样可以防止运算溢出，而结果不变，注意，当你需要计算除法时，譬如计算 (A/B)%mod，也许 A 和 B 本身很大很大，但是经过取模后变成一个相对较小的数，这里再这么算是不对的，比如 mod=7 时，30/10 的结果本来是 3，但是 A 和 B 对 7 取模后变成了 2/3，直接计算得到 0，就产生了错误，你可以使用下面的代码中 inv 函数，当你需要计算 A/B%mod 时可以改写成 (A%mod)*inv(B % mod)%mod，前提是 B 不为 0（在模 mod 后不为 0），注意数据溢出问题，你可能需要使用 long long 类型。inv 函数的复杂度为O(log mod)

typedef long long ll;
ll inv(ll x){
	ll b=MOD-2,ans=1;
	while(b){
		if(b&1){
			ans=ans*x%MOD;
		}
		x=x*x%MOD;
		b>>=1;
	}
	return ans;
}

输入格式

一行两个正整数N和M，空格隔开

输出格式

一个正整数，表示答案对998244353取模的结果

数据规模

N，M<=2000

 

样例输入1：

1 2

样例输出1：

1

样例输入1：

2 2

样例输出1：

6

样例输入1：

3 3

样例输出1：

10080

 

结果=(m*n)! / m! / n!

#include <iostream>
#include <algorithm>
#include <cstring>
#include <vector>
#include <map>
#define rep1(i,s,e,c) for(int i=s;i<e;i=i+c)
#define rep2(i,s,e,c) for(int i=s;i<=e;i=i+c)
#define rep3(i,e,s,c) for(int i=e;i>=s;i=i-c)
using namespace std;
typedef long long ll;
const int MAX=1e5+1;
const int INF=0x3f3f3f3f;
const int MOD=998244353;
ll inv(ll x){
	ll b=MOD-2,ans=1;
	while(b){
		if(b&1){
			ans=ans*x%MOD;
		}
		x=x*x%MOD;
		b>>=1;
	}
	return ans;
}
ll A(ll m,ll n){
	ll sum=1;
	for(ll i=(m-n+1);i<=m;i++){
		sum*=i;
		sum%=MOD;
	}
	return sum;
} 
int main()
{
	ll n,m,ans=1;
	cin>>n>>m;
	//ans=A(m*n,m*n-n)/A(m,m);
	ans=(A(m*n,m*n-n)%MOD)*inv(A(m,m)%MOD)%MOD;
	ans%=MOD; 
	cout<<ans<<"\n";
    return 0;
}