F - Make Pai(区间dp+计数)

F - Make Pai(区间dp+计数)

传送门

思路:

状态划分: d p [ i ] [ j ] dp[i][j] dp[i][j]​为i到j配对完的的方案数

我们让i去找个k配对,那么其中间的一定也要配对完才能选[i+1,k-1],即可 将其划分[i+1,k-1]和[k+1,j]两个集合,然后进行状态转移。

状态转移: d p [ i ] [ j ] = d p [ i + 1 ] [ k − 1 ] ∗ d p [ k + 1 ] [ j ] ∗ C x + y + 1 x dp[i][j]=dp[i+1][k-1]*dp[k+1][j]*C_{x+y+1}^{x} dp[i][j]=dp[i+1][k1]dp[k+1][j]Cx+y+1x​​其中 x = k − i − 1 > > 1 , y = j − k > > 1 x=k-i-1>>1,y=j-k>>1 x=ki1>>1,y=jk>>1​。​​​​​​​

C j − i + 1 j − k + 1 C_{j-i+1}^{j-k+1} Cji+1jk+1​​​​​​​​,的意义是两个集合中有x和y+1个配对放x+y+1个位置的方案。

代码:
#include <iostream>
#include <cstring>
using namespace std;
typedef long long ll;
const ll mod=998244353;
const int N=410;
int n,m,e[N][N];
ll f[N][N],fac[2000010];
ll ksm(ll x,ll p){
    ll res=1;
    while(p){
        if(p%2==1) res=res*x%mod;
        p/=2;
        x=x*x%mod;
    }
    return res;
}
ll inv(ll a) {
	return ksm(a,mod-2)%mod;
}
ll C(ll n,ll k){
	if(k>n)return 0;
	if(k==1)return n;
	return (fac[n]*inv(fac[k])%mod*inv(fac[n-k])%mod);
}
ll dfs(int l,int r){
    if(l>r) return 1;
    if(f[l][r]!=-1) return f[l][r];
    ll ans=0;
    for(int i=l+1;i<=r;i+=2){
        int x=i-l-1>>1,y=r-i>>1;
        if(!e[l][i]) continue;
        ans=(ans+1ll*dfs(l+1,i-1)*dfs(i+1,r)%mod*1ll*C(x+y+1,y)%mod)%mod;
    }
    return f[l][r]=ans;
}
int main(){
    cin>>n>>m;
    memset(f,-1,sizeof f);
    fac[0] = 1;
	for(int i = 1;i < 2000006; i++) {
		fac[i] = (fac[i-1]*i)%mod;
	}
    for(int i=1;i<=m;i++){
        int u,v;
        cin>>u>>v;
        e[u][v]=e[v][u]=1;
    }
    cout<<dfs(1,2*n)<<"\n";
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值