首页 > 其他 > 详细

AGC030F - Permutation and Minimum

时间:2019-06-11 11:07:22      阅读:144      评论:0      收藏:0      [点我收藏+]

https://atcoder.jp/contests/agc030/tasks/agc030_f

题解

我们先把这个排列从\(1 \sim 2n\)表达出来,然后题面中的每一对数我们可以用一条线把他们连起来,那么在序列中表达出的值是这条线的左端点。

如果一开始每个数都没有限制的话,我们则需要求有哪些数会成为左端点,这个其实就是卡特兰数,求答案的话还需要求一个阶乘。

现在有一些位置有了限制,我们就把有限制的位置成为特殊位置。

还是从\(1\sim 2n\)这个序列上考虑,因为我们的贡献是在左端点产生,所以我们设\(dp[i][j][k]\)表示做到\(i\),有\(j\)个普通点没有被匹配,有\(k\)个特殊点没有被匹配。

普通点可以任意匹配,特殊点只能匹配普通点。

如果普通点匹配了特殊点,那么它的位置已经确定了,所以最后我们要乘上普通配普通的对数的阶乘。

代码

#include<bits/stdc++.h>
#define N 602
using namespace std;
typedef long long ll;
int n,now,a[N],sm;
int tag[N];
ll jie[N],dp[2][302][302];
const int mod=1e9+7;
inline void MOD(ll &x){x=x>=mod?x-mod:x;}
inline ll rd(){
    ll x=0;char c=getchar();bool f=0;
    while(!isdigit(c)){if(c=='-')f=1;c=getchar();}
    while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
    return f?-x:x;
}
int main(){
    n=rd();
    sm=n*2;
    int cnt=0;
    for(int i=1;i<=sm;++i){
        a[i]=rd();
        if(a[i]!=-1)tag[a[i]]=1;
        if(i%2==0&&a[i]!=-1&&a[i-1]!=-1){
            tag[a[i]]=-1;
            tag[a[i-1]]=-1;
        }
        if(i%2==0&&a[i]==-1&&a[i-1]==-1)cnt++;
    }
    jie[0]=1;
    for(int i=1;i<=sm;++i)jie[i]=jie[i-1]*i%mod;
    dp[1][0][0]=1;
    int now=1,pre=0;
    for(int i=sm;i>=1;--i)if(tag[i]!=-1){
        swap(now,pre);
        memset(dp[now],0,sizeof(dp[now]));
        for(int j=0;j<=n;++j)
            for(int k=0;j+k<=n;++k)if(dp[pre][j][k]){
                if(tag[i]){
                  MOD(dp[now][j][k+1]+=dp[pre][j][k]);
                  if(j)MOD(dp[now][j-1][k]+=dp[pre][j][k]);
                }
                else{
                  MOD(dp[now][j+1][k]+=dp[pre][j][k]);
                  if(j)MOD(dp[now][j-1][k]+=dp[pre][j][k]);
                  if(k)MOD(dp[now][j][k-1]+=dp[pre][j][k]*k%mod);
                }
            }
    }
    printf("%lld",dp[now][0][0]*jie[cnt]%mod);
    return 0;
}

AGC030F - Permutation and Minimum

原文:https://www.cnblogs.com/ZH-comld/p/11002302.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!