Wannafly冬のキャンプ2020日6J K重排列 - DP

探している\（Kは\）数である\（N \）サイクルの順列。\（T \当量100、N \当量50、K \当量10 ^ {18} \）

解決

置換が提供し、DPを考慮し、グラフを有向いくつかのリングを試みることができる\（F [i]が\）を表す（\ N = I）\次に、応答時間を
\ [F [I] = \ sum_ { J = 1} ^ N [J | K] \ CDOT C_ {I-1} ^ {J-1} \ CDOT（J-1）\ CDOT F [IJ] \]！

#include <bits/stdc++.h>
using namespace std;

#define int long long

const int mod = 998244353;
const int N = 105;
int t,n,k,f[N],c[N][N],frac[N];

signed main() {
    frac[0]=1;
    for(int i=1;i<=100;i++) frac[i]=(frac[i-1]*i)%mod;
    for(int i=0;i<=100;i++) c[i][0]=1;
    for(int i=1;i<=100;i++) {
        for(int j=1;j<=i;j++) c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
    }
    cin>>t;
    while(t--) {
        cin>>n>>k;
        memset(f,0,sizeof f);
        f[0]=1;
        for(int i=1;i<=n;i++) {
            for(int j=1;j<=i;j++) {
                if(k%j==0) {
                    f[i]+=c[i-1][j-1]*frac[j-1]%mod*f[i-j]%mod;
                    f[i]%=mod;
                }
            }
        }
        cout<<f[n]<<endl;
    }
}