腾讯游戏开发了一款全新的编程类益智小游戏,榜首的题目是一道关于矩阵的计算,你用多久能计算出来呢?
游戏中给出一个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;
}